数字信号处理MATLAB上机实验2014年题稿.ppt

下载文档 降价啦

5
0
约2.68万字
约 127页
2017-02-14 发布于湖北
举报
版权申诉
保障服务

数字信号处理MATLAB上机实验2014年题稿.ppt

1、本文档共127页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

* kwp=fix(wp*M /(2*pi))+1 kws=fix(ws*M /(2*pi))+1 k=[0:1:M-1]; wk = (2*pi/M)*k; Hrs=[ones(1,kwp),zeros(1,M-2*kwp+1),ones(1,kwp-1)] % Hrs = [1,1,1,1,0.5,zeros(1,11),0.5,1,1,1] wd = [0,0.3,0.6,1]; Hrd = [1,1,0,0]; k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1; angH = [-alpha*(2*pi)/M*k1, alpha*(2*pi)/M*(M-k2)]; * H = Hrs.*exp(j*angH); h = real(ifft(H,M)); [Hm,wm] = freqz(h,1,512); mag=abs(Hm); pha=angle(Hm); db=20*log10((mag+eps)/max(mag)) n=[0:255]; fs=2000; x1=sin(2*pi/fs*100*n); x2=sin(2*pi/fs*800*n); x=x1+x2; * y=filter(h,1,x); figure(4); subplot(2,2,1); plot(wk(1:11)/pi,Hrs(1:11),o,wd,Hrd);grid on axis([0,1,-0.1,1.1]); title(Frequency Samples: M=20) xlabel(frequency in pi units); ylabel(Hr(k)); set(gca,XTickMode,manual,XTick,[0,0.3,0.4,0.6,1]) set(gca,YTickMode,manual,YTick,[0,0.5,1]); * subplot(2,2,2); stem(k,h);grid on title(Impulse Response); xlabel(n); ylabel(h(n)); subplot(2,2,3); plot(wk(1:11)/pi,Hrs(1:11),o,wm/pi,mag); grid on axis([0,1,-0.2,1.2]); title(Amplitude Response) xlabel(frequency in pi units); ylabel(Hr(w)) set(gca,XTickMode,manual,XTick,[0,0.3,0.4,0.6,1]) set(gca,YTickMode,manual,YTick,[0,0.5,1]); * subplot(2,2,4); plot(wm/pi,db);grid on axis([0,1,-60,10]);title(Magnitude Response); xlabel(frequency in pi units);ylabel(Decibels); set(gca,XTickMode,Manual,XTick,[0;0.3;0.6;1]); set(gca,YTickMode,Manual,YTick,[-29;0]); set(gca,YTickLabelMode,manual,YTickLabels,[29; 0]) % figure(5); subplot(2,1,1); plot(wm,mag); * title(Designed Lowpass Filter Magnitude Response) xlabel(w/rad); ylabel(|H[exp(jw)]|); grid subplot(2,1,2);plot(wm,pha); title(Designed Lowpass Filter Phase Response in radians); xlabel(w/rad); ylabel(arg(H[exp(jw)]); grid % figure(6); subplot(2,1,1);plot(n,x); title(Souce signal:x=sin(2*pi*100*t)+sin(2*pi*800*t) ); xlabel(time in s); ylabel(x); grid subplot(2,1,2);plot(n,y);title(Filtered signal:y); xlabel(time in s); ylabel(y); grid * * * * % Frequency_Sample2.m % Linear Phase Lowpass Filter Design % Digital Filter Speci