. 3
, . , , L2, L2. , L1. > L2
, , - ( L2) - ( > L2).
. 1
System #3 | |
;
;
, ;
, .
, .
- .
5 . , , :
, , i - .
Matlab.
gisterN(lam, m_1, m_2, L_1, L_2) , , :
. 4 GisterN
gisterQ(lam, m_1, m_2, L_1, L_2) , , :
. 5 GisterQ
function Q = gisterQ(lam, m_1, m_2, L_1, L_2)
%GISTERP returns average Q value for hysteresys control system
% k state number = queue length
% lam input intensity
% m_1 output intensity for 1st mode
% m_2 output intensity for 2nd mode
% L_1 hysteresis control 1st parameter
% L_2 hysteresis control 2nd parameter
r=lam./m_1;
r_1=lam./m_2;
q=L_2-L_1-1;
P_0=(1./(1-r)-((q+1).*r.^(L_1+q).*(r-r_1))./((1-r.^(q+1)).*(1-r_1))).^(-1);
|
|
Q=P_0.*(((r-(L_1+1).*r.^(L_1+1)+L_1.*r.^(L_1+2))./(1-r).^2 -(r-r.^(L_1+1))./(1-r))+(r.^(L_1).*(1-r))./(1-r.^(q+1)).*((L_1-1).*(((r-r.^(q+1))./(1-r)-q.*r.^(q+1))./(1-r)+(r.^q.*r_1.*(q-(r_1-r_1.^(q+1))./(1-r_1)))./(1-r_1))+(((r-(q+1).*r.^(q+1)+q.*r.^(q+2))./(1-r).^2 -(r.^(q+1).*q.*(q+1))./2)./(1-r)+(r.^q.*r_1.*((q.*(q+1))./2-(r_1-(q+1).*r_1.^(q+1)+q.*r_1.^(q+2))./(1-r_1).^2))./(1-r_1)))+((1-r_1.^(q+1)).*(1-r).*r.^(L_1+q))./((1-r.^(q+1)).*(1-r_1)).*((r_1.*(L_1+q-1))./(1-r_1)+r_1./(1-r_1).^2));
end
. 6 MatLab
. 5.
. 7
. 7 i - . , . , 1 , .
L , . ; , . , .
:
.
. 8 GisterW
function W = gisterW(lam, m_1, m_2, L_1, L_2)
% k state number = queue length
% lam input intensity
% m_1 output intensity for 1st mode
% m_2 output intensity for 2nd mode
% L_1 hysteresis control 1st parameter
% L_2 hysteresis control 2nd parameter
r=lam./m_1;
r_1=lam./m_2;
q=L_2-L_1-1;
P_0=(1./(1-r)-((q+1).*r.^(L_1+q).*(r-r_1))./((1-r.^(q+1)).*(1-r_1))).^(-1);
W=P_0./lam.*(((r-(L_1+1).*r.^(L_1+1)+L_1.*r.^(L_1+2))./(1-r).^2 -(r-r.^(L_1+1))./(1-r))+(r.^(L_1).*(1-r))./(1-r.^(q+1)).*((L_1-1).*(((r-r.^(q+1))./(1-r)-q.*r.^(q+1))./(1-r)+(r.^q.*r_1.*(q-(r_1-r_1.^(q+1))./(1-r_1)))./(1-r_1))+(((r-(q+1).*r.^(q+1)+q.*r.^(q+2))./(1-r).^2 -(r.^(q+1).*q.*(q+1))./2)./(1-r)+(r.^q.*r_1.*((q.*(q+1))./2-(r_1-(q+1).*r_1.^(q+1)+q.*r_1.^(q+2))./(1-r_1).^2))./(1-r_1)))+((1-r_1.^(q+1)).*(1-r).*r.^(L_1+q))./((1-r.^(q+1)).*(1-r_1)).*((r_1.*(L_1+q-1))./(1-r_1)+r_1./(1-r_1).^2));
end
. 9 GisterW
. 10 .
. 10
. 10 i - . , ; . , . . , . ,
|
|
.
. 11 GisterS
function S = gisterS(lam, m_1, m_2, L_1, L_2)
% k state number = queue length
% lam input intensity
% m_1 output intensity for 1st mode
% m_2 output intensity for 2nd mode
% L_1 hysteresis control 1st parameter
% L_2 hysteresis control 2nd parameter
r=lam./m_1;
r_1=lam./m_2;
q=L_2-L_1-1;
P_0=(1./(1-r)-((q+1).*r.^(L_1+q).*(r-r_1))./((1-r.^(q+1)).*(1-r_1))).^(-1);
S=P_0./lam.*(r./(1-r).^2 -((q+1).*r.^(L_1+q).*(r-r_1))./((1-r.^(q+1)).*(1-r_1)).*((2.*L_1+q)./2+(1-r.*r_1)./((1-r).*(1-r_1))))-P_0./lam.*(((r-(L_1+1).*r.^(L_1+1)+L_1.*r.^(L_1+2))./(1-r).^2 -(r-r.^(L_1+1))./(1-r))+(r.^(L_1).*(1-r))./(1-r.^(q+1)).*((L_1-1).*(((r-r.^(q+1))./(1-r)-q.*r.^(q+1))./(1-r)+(r.^q.*r_1.*(q-(r_1-r_1.^(q+1))./(1-r_1)))./(1-r_1))+(((r-(q+1).*r.^(q+1)+q.*r.^(q+2))./(1-r).^2 -(r.^(q+1).*q.*(q+1))./2)./(1-r)+(r.^q.*r_1.*((q.*(q+1))./2-(r_1-(q+1).*r_1.^(q+1)+q.*r_1.^(q+2))./(1-r_1).^2))./(1-r_1)))+((1-r_1.^(q+1)).*(1-r).*r.^(L_1+q))./((1-r.^(q+1)).*(1-r_1)).*((r_1.*(L_1+q-1))./(1-r_1)+r_1./(1-r_1).^2));
end
. 12 GisterS
. 13 .
. 13
. 13 i - . , , . , . ,
.
, 1, . 2.
. 2