%calcolo dell'errore percentule della formula approssimata del limite
%superiore di banda trovata dal Millman, nel caso di un sistema con 3 poli
% f1, f2 ed f3
clear all
% f1 ed f2 variano, mentre f3 e' fissato
f3=10000
%f1 ed f2 variano tra 1 ed 10000 Hz
f1min=1;
f2min=1;
f1max=10000;
f2max=10000;
f1=linspace(f1min,f1max,20);
f2=linspace(f2min,f2max,20);
[F1,F2]=meshgrid(f1,f2);
nx=length(F1)
F3=f3*ones(nx,nx);
%calcolo della formula approssimata
RHS=sqrt(1./F1.^2.+1./F2.^2.+1./(F3.^2));
FAPPR=1/1.1./RHS;
for I = 1:nx,
    for J = 1:nx,
        b=tf([1],[1 f1(I)+f2(J)+f3 (f1(I)*f2(J)+f2(J)*f3+f1(I)*f3) f1(I)*f2(J)*f3]);
        fb2=bandwidth(b);
        FB(I,J)=fb2;
    end
end
PERCENTO=abs(FAPPR-FB)./FB*100;
mesh(F1,F2,PERCENTO)