如果您无法下载资料,请参考说明:
1、部分资料下载需要金币,请确保您的账户上有足够的金币
2、已购买过的文档,再次下载不重复扣费
3、资料包下载后请先用软件解压,在使用对应软件打开
Homework王明皓11001006251、利用Matlab设计Butterworth模拟滤波器,(一)通带留余量程序Wp=2*pi*5000;Rp=2;Ws=2*pi*12000;Rs=30;[N,Wc]=buttord(Wp,Ws,Rp,Rs,'s');%通带留余量[z,p,k]=buttap(N);[Bp,Ap]=zp2tf(z,p,k);%归一化函数的系数[Bs,As]=lp2lp(Bp,Ap,Wc);%解归一化f=0:100:16000;HK=freqs(Bs,As,f*2*pi);subplot(2,1,1);plot(f/1000,abs(HK));axis([01601]);xlabel('f/KHZ');ylabel('幅度');title('归一化图');subplot(2,1,2);plot(f/1000,20*log10(abs(HK)));axis([016-500]);gridonxlabel('f/KHZ');ylabel('幅度/dB');title('dB图');gtext('f=5Khzy=-0.6354dB');gtext('f=12Khzy=-30dB');结果:Bs=7.7094e+22As=[11.2230e+057.4785e+092.8263e+146.6014e+187.7094e+22];(二)阻带留余量Wp=2*pi*5000;Rp=2;Ws=2*pi*12000;Rs=30;[N,Wc]=buttord(Wp,Ws,Rp,Rs,'s');Wc=Wp/((10^(.1*abs(Rp))-1)^(1/(2*N)));%阻带留余量[z,p,k]=buttap(N);[Bs,As]=lp2lp(Bp,Ap,Wc);%解归一化f=0:100:16000;HK=freqs(Bs,As,f*2*pi);subplot(2,1,1);plot(f/1000,abs(HK));%axis([01601]);xlabel('f/KHZ');ylabel('归一化幅度');title('归一化图');subplot(2,1,2);plot(f/1000,20*log10(abs(HK)));set(gca,'XTick',[0:1:16]);axis([016-500]);gridonxlabel('f/KHZ');ylabel('幅度/dB');title('dB图');gtext('f=5Khzy=-2dB');gtext('f=12Khzy=-35.69dB');:Bs=4.4414e22As=[11.0727e055.7529e091.9069e143.9065e184.0014e22];作业2画出下图低通滤波器的幅度响应,频率范围取0-1kHz,其中RC=0.001。传递函数:程序:>>f=0:10:1000;>>w=2*pi*f;>>hw=1./(1+0.001*j*w);>>xlabel('f/HZ');>>ylabel('|h(jw)|');>>title('一阶RC网络幅频特性曲线');结果: