91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 5KB
    文件類(lèi)型: .rar
    金幣: 2
    下載: 0 次
    發(fā)布日期: 2021-05-18
  • 語(yǔ)言: 其他
  • 標(biāo)簽: 預(yù)測(cè)??GPC??

資源簡(jiǎn)介

GPC有一定幫助 想學(xué)習(xí)的 可以下載

資源截圖

代碼片段和文件信息

%?%?clear;
%?%?clc;
%?%?c=[2;3;-5];
%?%?a=[-25-1;131];
%?%?b=[-10;12];
%?%?aeq=[111];
%?%?beq=7;
%?%?x=linprog(-cabaeqbeqzeros(31));
%?%?value=c‘*x
%?
%?clear?;
%?%?na=2;%多項(xiàng)式A的階
%?%?nb=1;%多項(xiàng)式B的階
%?%?nc=0;
%?%?N=4;%預(yù)測(cè)時(shí)間長(zhǎng)度
%?%?Nu=1;%控制時(shí)間長(zhǎng)度
%?%?L=zeros(NNu);%依據(jù)P82公式初始化矩陣L,H,G,E,Q,R
%?%?H=zeros(Nnb);
%?%?G=zeros(Nna+1);
%?%?E=zeros(NN);
%?%?Q=0.4*eye(N);
%?%?%?R=1*zeros(NuNu);
%?%?R=5*eye(NuNu);
%?[EGLH]=ComputationEGLH(41210[-0.9;?0;?-0.1][0.001;0.0015]);
%?%?[EGLH]=ComputationEGLH(NNunanbncAGanBVec);




%?function?[EGLH]=ComputationEGLH1(NNunanbncAGanBVecCVec)
???clear;
???clc;
???N=6;
???Nu=1;
???na=2;
???nb=1;
???nc=0;
???n=2;
%????AGan=[0.9?-0.2?0?0;0?0?0.8?-0.1]‘;
???AGan=[-1.9?1.1?-0.2?0?0?0;0?0?0?-1.8?0.9?-0.1]‘;
???BVec=[0.6?0.3?0.5?0.3;0.4?0.3?0.5?0.3]‘;
???%注:AGan第一項(xiàng)是z(z-1)的系數(shù),AGan首項(xiàng)1省略
???L=zeros(n*Nn*Nu);%依據(jù)P82公式初始化矩陣L,H,G,E,Q,R
???H=zeros(n*Nn*nb);
???G=zeros(n*Nn*(na+1));
???E=zeros(n*Nn*N);
???
???E(1:n1:n)=eye(n);?%依據(jù)P79公式計(jì)算E1
???G(1:n:)=-AGan‘;?%依據(jù)P79公式計(jì)算G1
???L(1:n1:n)=BVec(1:n1:n);%依據(jù)P80公式計(jì)算L1
???H(1:n1:n)=BVec(n+1:end1:n);%依據(jù)P80公式計(jì)算H1
???for?j=2:N
???????EVecNow=E(n*(j-1)-1:n*(j-1)1:n*(N-1));
???????GVecNow=G(n*(j-1)-1:n*(j-1):);%分離G矩陣的第(j-1)行,用以計(jì)算新的E矩陣公式
???????E(n*j-1:n*j1:n)=GVecNow(1:n1:n);%公式4-48確定矩陣第j行第一個(gè)值
???????E(n*j-1:n*j2*n-1:n*N)=EVecNow;%確定E矩陣第j行的第2-N的值,見(jiàn)P80公式
???????%確定G矩陣的第j行?
???????for?i=1:na
%??????????G(ji)=G(j-1i+1)-G(j-11)*AGan(i1);?%公式(4-49)注意AGan(i1)的標(biāo)號(hào)??
???????????G(n*j-1:n*jn*i-1:n*i)=G(n*(j-1)-1:n*(j-1)n*(i+1)-1:n*(i+1))-G(n*(j-1)-1:n*(j-1)1:n)*AGan(n*i-1:n*i1:n);?%公式(4-49)注意AGan(i1)的標(biāo)號(hào)
???????end
%????????G(jna+1)=-G(j-11)*AGan(na+11);%公式(4-49)
???????G(n*j-1:n*jn*(na+1)-1:n*(na+1))=-G(n*(j-1)-1:n*(j-1)1:n)*AGan(n*(na+1)-1:n*(na+1)1:n);%公式(4-49)
???????L(n*j-1:n*j1:n)=E(n*j-1:n*j1:n)*BVec(1:n1:n)+H(n*(j-1)-1:n*(j-1)1:n);%公式(4-51)計(jì)算L矩陣第j行
???????if?Nu>1
???????????L(j2:Nu)=L(j-11:Nu-1);
???????end
???????nh=max(nb-1nc-1);
???????if?nh>0
???????????for?i=1:nh
???????????????H(n*j-1:n*jn*i-1:n*i)=E(n*j-1:n*j1:n)*BVec(n*i-1:n*i1:n)+H(n*(j-1)-1:n*(j-1)n*(i+1)-1:n*(i+1))-L(n*(j+1)-1:n*(j+1)1:n)*CVec(n*i-1:n*i1:n);?%公式(4-52)注意CVec(i1)不存在
???????????end
???????????H(n*j-1:n*jn*(nh+1)-1:n*(nh+1))=E(n*j-1:n*j1:n)*BVec(n*(nb+1)-1:n*(nb+1)1:n);%公式(4-52),另一部分?-L(j1)*CVec(nb1)置0
???????else
???????????H(n*j-1:n*jn*(nh+1)-1:n*(nh+1))=E(n*j-1:n*j1:n)*BVec(n*(nb+1)-1:n*(nb+1)1:n);%對(duì)于c=1的情形,由公式(4-52),另一部分?-L(j1)*CVec(nb1)置0
???????end
???end

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????2805??2011-05-25?19:14??ComputationEGLH.m

?????文件???????4947??2008-06-04?13:46??GPC.m

?????文件???????4171??2011-05-19?21:07??mimo.m

?????文件????????433??2011-01-05?21:15??mimopoleplacement.m

-----------?---------??----------?-----??----

????????????????12356????????????????????4


評(píng)論

共有 條評(píng)論

相關(guān)資源