mimo_noise.m

Software for design and tuninig of SISO and MIMO contol systems

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;