Code covered by the BSD License  

Highlights from
Logiciels temps réel

Logiciels temps réel

by

 

Real-time software concerning the domains of the signal processing, feedback and regulation

test_passe_bande.m
n=[1]
d=[1 1 1]
 [nt,dt]=lp2bp(n,d,0.7,0.2)
 nt(2:end)=0

[z,p,k]=tf2zp(nt,dt)
 w=logspace(-2,2,1000);
semilogx(w,20*log10(abs(polyval(nt,1i*w)./polyval(dt,1i*w))))
 abs(p)

 kpl=sqrt(k)*abs(p(1))

 [nl,dl]=zp2tf([],p(1:2),kpl)

 hold on
semilogx(w,20*log10(abs(polyval(nl,1i*w)./polyval(dl,1i*w))),'r')
 kph=sqrt(k)/abs(p(1))

 [nh,dh]=zp2tf([0;0],p(3:4),kph)
semilogx(w,20*log10(abs(polyval(nh,1i*w)./polyval(dh,1i*w))),'m')

Contact us