完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
扫一扫,分享给好友
附上代码
wp=0.4*pi;ws=0.6*pi;deltap=0.01;deltas=0.001; deltaH=max(deltap,deltas);deltaL=min(deltap,deltas); weights=[deltas/deltap 1] deltaf=(ws-wp)/(2*pi) M=ceil((-20*log10(sqrt(deltap*deltas))-13)/(14.6*deltaf)+1) f=[0 wp/pi ws/pi 1] m=[1 1 0 0]; h=remez(M-1,f,m,weights); [db,mag,pha,grd,w]=freqz_m(h,[1]); delta_w=2*pi/1000;wsi=ws/delta_w+1;wpi=wp/delta_w; Asd=-max(db(wsi:1:501)) M=M+1; h=remez(M-1,f,m,weights); [db,mag,pha,grd,w]=freqz_m(h,[1]); Asd=-max(db(wsi:1:501)) M=M+1; h=remez(M-1,f,m,weights); [db,mag,pha,grd,w]=freqz_m(h,[1]); Asd=-max(db(wsi:1:501)) M=M+1; h=remez(M-1,f,m,weights); [db,mag,pha,grd,w]=freqz_m(h,[1]); Asd=-max(db(wsi:1:501)) M=M+1; h=remez(M-1,f,m,weights); [db,mag,pha,grd,w]=freqz_m(h,[1]); Asd=-max(db(wsi:1:501)) M [Hr,omega,P,L]=ampl_res(h); %画图 figure(1);subplot(1,1,1) subplot(2,2,1);stem([0:1:M-1,h]);title('脉冲响应') axis([0 M-1-0.1 0.3]);ylabel('h(n)');text(M,-0.1,'n') set(gca,'XTickMode','manual','XTick',[0,M-1]) set(gca,'YTickMode','manual','YTick',[-0.1:0.1:0.3]) subplot(2,2,2);plot(w/pi,db);title('幅度响应(单位:dB)'); axis([0,1,-80,10]);xlabel('');ylabel('分贝数') set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]) set(gca,'YTickMode','manual','YTick',[-50,0]); set(gca,'YTickLabelMode','mamual','YTickLabels',['50';'0']);grid subplot(2,2,3);plot(omega/pi,Hr);title('振幅响应'); axis([0 1 -0.1 1.1]);xlabel('频率(单位:pi)');ylabel('Hr(w)') set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]) set(gca,'YTickMode','manual','YTick',[0,1]);grid subplot(2,2,4); pbw=omega(1:1:wpi+1)/pi;pbe=Hr(1:1:wpi+1)-1; ***w=omega(wsi+1:501)/pi;***e=Hr(wsi+1:501); plot(pbw,pbe,***w,***e); axis([0,1,-deltaH,deltaH]);title('误差响应'); xlabel('频率(单位:pi)');ylabel('Hr(w)') set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]) set(gca,'YTickMode','manual','YTick',[-deltaH,-deltaL,0,deltaL,deltaH]); set(gca,'XGrid','on') |
|
|
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问,请使用“评论”功能。
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2024-4-20 06:32 , Processed in 0.408932 second(s), Total 40, Slave 34 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191 工商网监 湘ICP备2023018690号