chl.m

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

%----------------------------subfunction----------------------------------- 
function HL=CHL(na,nb,z,u) 
L=length(z); 
HL=zeros(L,na+nb+1); 
for i=1:L 
for j=1:na+nb+1 
if j<=na 
if i-j<0 
HL(i,j)=0; 
else 
HL(i,j)=-z(i-j+1); 
end 
else 
if i-j+na<0 
HL(i,j)=0; 
else 
HL(i,j)=u(i-j+na+1); 
end 
end 
end 
end 
HL=HL(na:L-1,:);