资源简介
GPC有一定帮助 想学习的 可以下载

代码片段和文件信息
% % 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;%多项式A的阶
% % nb=1;%多项式B的阶
% % nc=0;
% % N=4;%预测时间长度
% % Nu=1;%控制时间长度
% % L=zeros(NNu);%依据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第一项是z(z-1)的系数,AGan首项1省略
L=zeros(n*Nn*Nu);%依据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); %依据P79公式计算E1
G(1:n:)=-AGan‘; %依据P79公式计算G1
L(1:n1:n)=BVec(1:n1:n);%依据P80公式计算L1
H(1:n1:n)=BVec(n+1:end1:n);%依据P80公式计算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)行,用以计算新的E矩阵公式
E(n*j-1:n*j1:n)=GVecNow(1:n1:n);%公式4-48确定矩阵第j行第一个值
E(n*j-1:n*j2*n-1:n*N)=EVecNow;%确定E矩阵第j行的第2-N的值,见P80公式
%确定G矩阵的第j行
for i=1:na
% G(ji)=G(j-1i+1)-G(j-11)*AGan(i1); %公式(4-49)注意AGan(i1)的标号
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)的标号
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)计算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);%对于c=1的情形,由公式(4-52),另一部分 -L(j1)*CVec(nb1)置0
end
end
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 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
- 上一篇:苏昱CSS2.0中文手册
- 下一篇:JPG图片文件捆绑工具
相关资源
- LPC线性预测分析及编码
- 基于经济时间序列预测的露天矿开采
- 基于数据驱动的故障预测模型
- 提高风出力预测精度的储能系统模糊
- 电网负荷预测论文集9篇
- 基于蚁群算法优化SVM的瓦斯涌出量预
- 基于模糊聚类和SVM的瓦斯涌出量预测
- 最小预测跷跷板模型中的脂肪生成
- 中微子混合矩阵和无中微子双β衰
- 基于时序SAR技术的采空区上方高速公
- 约化分子电距矢量用于茎用莴苣花挥
- 基于CAN总线与ZigBee的瓦斯实时监测及
- 基于熵权法的PHC管桩承载力组合预测
- 基于动态分类器集成选择和GM(21)
- 缓冲算子理论在矿井瓦斯涌出量预测
- 瓦斯涌出量预测的GM(11)模型改进
- 矿井瓦斯涌出量建模预测
- 基于灰色系统理论的矿井瓦斯涌出量
- 矿井回采工作面瓦斯涌出量预测新途
- 矿井瓦斯涌出量的GM(11)新陈代谢模
- GM(11)模型与线性回归组合方法在矿
- 灰色系统理论在矿井瓦斯涌出量预测
- 等维新息模型在矿井瓦斯涌出量预测
- 基于改进灰色模型的矿井瓦斯涌出量
- GM(11)与ARIMA模型在中国一次能源消
- 基于时间加权-新陈代谢GM(11)模型在
- 基于新陈代谢GM(11)-ARMA(pq)模型的
- 改进的GM(11)模型在城市规划货运量
- 大型强子对撞机中单个矢量样夸克生
- 超外围碰撞中的排他和/ψ光产生:
评论
共有 条评论