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

如题所述

第1个回答推荐于2016-02-01

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('Butterworth BS DF');
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('ChebyshevI BS DF');
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('ChebyshevII BS DF');
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('Tuoyuan BS DF');
ylabel('Gain,dB');本回答被网友采纳

相似回答

求一个matlab程序设计,双线性变换法的低通滤波器的设计答：a]=lp2lp(Bp,Ap,Wn); %模拟低通滤波器到模拟低通滤波器的转换 %用双线性变换法得到数字低通滤波器系数[bz,az]=bilinear(b,a,fs); figure(1);

利用MATLAB仿真软件系统结合双线性变换法设计一个数字切比雪夫带阻II...答：利用双线性变换法设计IIR带阻滤波器 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]...

联系matlab用双线性变换法设计Butterworth低通滤波器m答：1.用Butterworth低通模拟滤波器原型%%%%%%%%%%%%%%%%%%%%%%%%%%%%数字滤波器指标 Wp=0.3*pi; %通带边界频率(归一化)：wp=fp*2*pi/fs Ws=0.5*pi; %阻带边界频率(归一化)：ws=fr*2*pi/fs Rp=0.5; %通带波纹 As=40; %阻带衰减转换为模拟滤波器指标(预畸)；Ts=0....

关于matlab的双线型变换答：1、你手算的有错误吧？按照z=(2-s)/(2+s)替换，不应该是你写的结果：>> syms z s>> g = (3*z-4)/(z^2-2*z+1);>> g1=simple(subs(g,z,(2+s)/(2-s))) g1 = 1/4*(-4+16*s-7*s^2)/s^2 2、直接用d2c(g)并不是使用双线性变换法，而是使用默认的零阶保持法...