matlab求双线性变换法,如何用MATLAB编写双线性变换法

最新推荐文章于 2024-01-09 08:52:41 发布

海燕老师

最新推荐文章于 2024-01-09 08:52:41 发布

阅读量2.2k

点赞数

文章标签： matlab求双线性变换法

wp1=0。3*pi;wp2=0。6*pi;ws1=0。4*pi;ws2=0。5*pi;Ap=3;As=10;Fs=1;B=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2));ws=1;[N,wc]=buttord(wp,ws,Ap,As,'s');[num,den]=butter(N,wc,'s');[numt,dent]=lp2bs(num,den,w0,B);[numd,dend]=bilinear(numt,dent,Fs);w=linspace(0,pi,512);h=freqz(numd,dend,w);norm=max(abs(h));numd=numd/norm;subplot(2,2,1);plot(w/pi,20*log(abs(h)/norm));grid;xlabel('ButterworthBSDF');ylabel('Gain,dB');wp1=0。

3*pi;wp2=0。6*pi;ws1=0。4*pi;ws2=0。5*pi;Ap=3;As=10;Fs=1;B=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2));ws=1;[N,wc]=cheb1ord(wp,ws,Ap,As,'s');[num,den]=cheby1(N,Ap,wc,'s');[numt,dent]=lp2bs(num,den,w0,B);[numd,dend]=bilinear(numt,dent,Fs);w=linspace(0,pi,512);h=freqz(numd,dend,w);norm=max(abs(h));numd=numd/norm;subplot(2,2,2);plot(w/pi,20*log(abs(h)/norm));grid;xlabel('ChebyshevIBSDF');ylabel('Gain,dB');wp1=0。

3*pi;wp2=0。6*pi;ws1=0。4*pi;ws2=0。5*pi;Ap=3;As=10;Fs=1;B=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2));ws=1;[N,wc]=cheb2ord(wp,ws,Ap,As,'s');[num,den]=cheby2(N,As,wc,'s');[numt,dent]=lp2bs(num,den,w0,B);[numd,dend]=bilinear(numt,dent,Fs);w=linspace(0,pi,512);h=freqz(numd,dend,w);norm=max(abs(h));numd=numd/norm;subplot(2,2,3);plot(w/pi,20*log(abs(h)/norm));grid;xlabel('ChebyshevIIBSDF');ylabel('Gain,dB');wp1=0。

3*pi;wp2=0。6*pi;ws1=0。4*pi;ws2=0。

5*pi;Ap=3;As=10;Fs=1;B=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2));ws=1;[N,wc]=ellipord(wp,ws,Ap,As,'s');[num,den]=ellip(N,Ap,As,wc,'s');[numt,dent]=lp2bs(num,den,w0,B);[numd,dend]=bilinear(numt,dent,Fs);w=linspace(0,pi,512);h=freqz(numd,dend,w);norm=max(abs(h));numd=numd/norm;subplot(2,2,4);plot(w/pi,20*log(abs(h)/norm));grid;xlabel('TuoyuanBSDF');ylabel('Gain,dB');。

全部