stepresponse.m

来自「预测控制中动态矩阵控制（DMC）的Matlab代码」· M 代码 · 共 19 行

%----------------------------subfunction----------------------------------- 
function a=stepresponse(P,AA,BB) %求阶越响应a 
a=zeros(P,1); %初始化阶越响应 
nu=length(BB)-1; 
ny=length(AA)-1; 
Y=0; 
YY(1:ny+1)=0; 
UU(1:nu+1)=0; 
for t=1:P 
UU(2:end)=UU(1:end-1); 
UU(1)=1; 
a(t)=Y; 
Y=BB*UU(1:nu+1)'-AA(2:ny+1)*YY(1:ny)'; 
if ny>=1 
YY(2:end)=YY(1:end-1); 
end 
YY(1)=Y; 
end