| [UorXmin,UorXmax]=xu2minmax(POLright,POLgur,prostoPOL,NofS1,TofS1)
|
function [UorXmin,UorXmax]=xu2minmax(POLright,POLgur,prostoPOL,NofS1,TofS1)
%<u^2 min> , <x^2 max>
alpha=roots(POLright);
POLstar=zeros(1,length(alpha));
for k=1:(length(alpha))
aj=alpha(k);
POLstar=POLstar+deconv(POLright,[-1,aj])*polyval(POLright,-1*aj)*polyval(NofS1,aj)/polyval(TofS1,aj)/polyval(prostoPOL,aj)/polyval(polyder(POLright),aj);
end;
drPOLright=POLright(1)*((-1)^(length(POLright)+1))*poly(-1*roots(POLright));
UorXmin=(norm(tf(POLstar,drPOLright),2))^2*pi;
vspomPOL=tf( polsum(conv( conv(prostoPOL,TofS1),POLstar),conv(NofS1,drPOLright)),conv(conv(TofS1,drPOLright),conv(POLgur,drPOLright)) );
UorXmax=(norm(vspomPOL,2))^2*pi;
|
|