mimo_noise.m

来自「内模控制器（IMC）工具箱。包括参数整定、PID控制器参数转换等」· M 代码 · 共 22 行

function f=mimo_noise(e,qs,order,Na,y)
n=length(order);
s=1e-5;
fil=zeros(1,n);
for i=1:n
   fil(i)=1/(e(i)*s+1)^order(i);
end
fil=diag(fil);
q0=inv(eval(qs))*fil;
f=ones(1,n);
while s < 1e5
s=s*10;
fil=zeros(1,n);
for i=1:n
   fil(i)=1/(e(i)*s+1)^order(i);
end
fil=diag(fil);
qinf=inv(eval(qs))*fil;
f=max([f;qinf./q0]);
end
f=f-Na;