dy.m

预测控制中动态矩阵控制（DMC）的Matlab代码

%----------------------------subfunction----------------------------------- 
function [A,s,e1]=DY(P,M,a) %%赋值动态矩阵A,移位矩阵s,首1矩阵 e1 
%%%%%%%%%% 赋值动态矩阵A %%%%%%%%%%%%% 
for j=1:M 
cc=1; 
for i=j:P %i,j是A的下标,cc是a的下标 
A(i,j)=a(cc); 
cc=cc+1; 
end 
end 
%%%%%%%%%% 赋值移位矩阵s %%%%%%%%%%%%%% 
for i=1:(P-1) 
j=i+1; 
s(i,j)=1; 
s(P,P)=1; 
end 
%%%%%%%%%% 赋值首1矩阵 %%%%%%%%%%%%%% 
temp1=eye(M); 
e1=temp1(1,:); %首 1 矩阵