数字信号处理实验五数字信号理实验五处理实验五.doc

下载文档 降价啦

21
0
约2.32千字
约 9页
2017-03-18 发布于贵州
举报
版权申诉
保障服务

数字信号处理实验五数字信号理实验五处理实验五.doc

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

数字信号处理实验五数字信号理实验五数字信号处理实验五数字信号处理实验五

实验clear all N=100; n=0:N-1; xn=(0.9).^n; XK=fft(xn,N); magXK=abs(XK); phaXK=angle(XK); subplot(1,2,1); plot(n,xn) xlabel(n); ylabel(x(n)); title(x(n)N=100); subplot(1,2,2); k=0:length(magXK)-1; stem(k,magXK,.); xlabel(k);ylabel(|X(k)|); title(X(k)N=100); 2) 周期为N 的余弦序列：v2(n) = cos(2*pi*n/N)，且0≤n≤N-1 程序 clear all N=100; n=0:N-1; xn=cos(2*pi*n/N); XK=fft(xn,N); magXK=abs(XK); phaXK=angle(XK); subplot(1,2,1); plot(n,xn) xlabel(n); ylabel(x(n)); title(x(n)N=100); subplot(1,2,2); k=0:length(magXK)-1; stem(k,magXK,.); xlabel(k);ylabel(|X(k)|); title(X(k)N=100); 3)复合函数序列：v3(n) = 0.9sin(2*pi*n/N) + 0.6sin(2*pi*n*3/N) 程序 clear all N=100; n=0:N-1; xn=0.9*sin(2*pi*n/N)+0.6*sin(2*pi*n/(N/3)); XK=fft(xn,N); magXK=abs(XK); phaXK=angle(XK); subplot(1,2,1); plot(n,xn) xlabel(n); ylabel(x(n)); title(x(n)N=100); subplot(1,2,2); k=0:length(magXK)-1; stem(k,magXK,.); xlabel(k);ylabel(|X(k)|); title(X(k)N=100); （3）N=32 clear all N=32; n=0:N-1; xn=(0.9).^n; XK=fft(xn,N); magXK=abs(XK); phaXK=angle(XK); subplot(1,2,1); plot(n,xn) xlabel(n); ylabel(x(n)); title(x(n)N=32); subplot(1,2,2); k=0:length(magXK)-1; stem(k,magXK,.); xlabel(k);ylabel(|X(k)|); title(X(k)N=32); N=64: N=128: (4)取N=32; clear all N=32; n=0:N-1; xn=cos(2*pi*n/N); XK=fft(xn,N); magXK=abs(XK); phaXK=angle(XK); subplot(1,2,1); plot(n,xn) xlabel(n); ylabel(x(n)); title(x(n)N=32); subplot(1,2,2); k=0:length(magXK)-1; stem(k,magXK,.); xlabel(k);ylabel(|X(k)|); title(X(k)N=32); 2N=64; N+2=34; (5) ① 32 点FFT: clear all N=32; n=0:N-1; xn=0.15*sin(2*pi*n/32)+sin(4*pi*n/32)-0.1*sin(6*pi*n/32); XK=fft(xn,N); magXK=abs(XK); phaXK=angle(XK); subplot(1,2,1); plot(n,xn) xlabel(n); ylabel(x(n)); title(x(n)N=32); subplot(1,2,2); k=0:length(magXK)-1; stem(k,magXK,.); xlabel(k);ylabel(|X(k)|); title(X(k)N=32); ②64 点FFT; clear all N=64; n=0:N-1; xn=0.15*sin(2*pi*n/32)+sin(4*pi*n/32)-0.1*sin(6*pi*n/32); XK=fft(xn,N); magXK=abs(XK); phaXK=angle(XK); subplot(1,2,1); plot(n,xn) xlabel(n); ylabel(x(n)); title(x(n)N=64