1、百度文库下载几篇本科的现成论文 你就知道了2、仿真就秒杀吧低通巴特沃斯模拟滤波器设计。通带截至频率3400 Hz,通带最大衰减3dB阻带截至频率4000 Hz,阻带最小衰减40dBIir2:模拟低通滤波器转换为数字低通滤波器,脉冲响应不变法和双线性变换法。Iir3:切比雪夫二型低通数字滤波器设计通带边界频率π,通带最大衰减1dB阻带截至频率π,阻带最小衰减80dBIir4:椭圆带通数字滤波器设计Iir5:高通和带通巴特沃思数字滤波器设计双线性变换%低通巴特沃斯模拟滤波器设计clear; close allfp=3400; fs=4000; Rp=3; As=40;[N,fc]=buttord(fp,fs,Rp,As,'s')[B,A]=butter(N,fc,'s');[hf,f]=freqs(B,A,1024);plot(f,20*log10(abs(hf)/abs(hf(1))))grid, xlabel('f/Hz'); ylabel('幅度(dB)')axis([0,4000,-40,5]);line([0,4000],[-3,-3]);line([3400,3400],[-90,5])%用脉冲响应不变法和双线性变换法将模拟滤波器离散化clear; close allb=1000;a=[1,1000];w=[0:1000*2*pi];[hf,w]=freqs(b,a,w);subplot(2,3,1); plot(w/2/pi,abs(hf)); grid;xlabel('f(Hz)'); ylabel('幅度'); title('模拟滤波器频响特性')Fs0=[1000,500];for m=1:2Fs=Fs0(m)[d,c]=impinvar(b,a,Fs)[f,e]=bilinear(b,a,Fs)wd=[0:512]*pi/512;hw1=freqz(d,c,wd);hw2=freqz(f,e,wd);subplot(2,3,2); plot(wd/pi,abs(hw1)/abs(hw1(1))); grid on; hold ontitle('脉冲响应不变法')subplot(2,3,3); plot(wd/pi,abs(hw2)/abs(hw2(1))); grid on; hold ontitle('双线性变换法')end%切比雪夫Ⅱ型低通数字滤波器设计clear; close allwp=; ws=; Rp=1; Rs=80;[N,wc]=cheb2ord(wp,ws,Rp,Rs)[B,A]=cheby2(N,Rs,wc)freqz(B,A)%直接设计带通数字椭圆滤波器clear; close allWp=[]; Ws=[];Rp=; Rs=60;[N,wc]=ellipord(Wp,Ws,Rp,Rs)[b,a]=ellip(N,Rp,Rs,wc)[hw,w]=freqz(b,a);subplot(2,1,1); plot(w/pi,20*log10(abs(hw))); gridaxis([0,1,-80,5]); xlabel('w/π'); ylabel('幅度(dB)')subplot(2,1,2); plot(w/pi,angle(hw)); gridaxis([0,1,-pi,pi]); xlabel('w/π'); ylabel('相位(rad)')%用双线性变换法设计数字高通和带通滤波器clear; close allT=1; wch=pi/2;wlc=*pi; wuc=*pi;B=1; A=[1,];[h,w]=freqz(B,A,512);subplot(2,2,1); plot(w,20*log10(abs(h))); grid%axis([0,10,-90,0]); xlabel('w/π'); title('模拟低通幅度(dB)')%高通omegach=2*tan(wch/2)/T;[Bhs,Ahs]=lp2hp(B,A,omegach);[Bhz,Ahz]=bilinear(Bhs,Ahs,1/T);[h,w]=freqz(Bhz,Ahz,512);subplot(2,2,3); plot(w/pi,20*log10(abs(h))); gridaxis([0,1,-150,0]); xlabel('w/π'); title('数字高通幅度(dB)')%带通omegalc=2*tan(wlc/2)/T;omegauc=2*tan(wuc/2)/T;wo=sqrt(omegalc*omegauc); Bw=omegauc-omegalc;[Bbs,Abs]=lp2bp(B,A,wo,Bw);[Bbz,Abz]=bilinear(Bbs,Abs,1/T);[h,w]=freqz(Bbz,Abz,512);subplot(2,2,4); plot(w/pi,20*log10(abs(h))); gridaxis([0,1,-150,0]); xlabel('w/π'); title('数字带通幅度(dB)')