原创力文档- 1、本文档共6页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
实验3-离散信号的离散频域分析
数字信号处理A
实验项目名称: 3离散信号的离散频域分析
学 院:______计算机与通信工程____
专 业:______ _通信工程 _________
学 号:______201xxxxxxxxxx_______
班 级:______ 通信1x-x_________
报 告 人:_____ ___xx__________
指 导 老 师:___ 胡双红 _ _______
实 验 时 间:_______201x-1x-2x________
实验目的：
1、掌握离散时间系统的DFT的MATLAB实现；
2、熟悉DTFT和DFT之间的关系。
3、了解信号不同变形的DFT与原信号DFT之间的关系
实验内容：
：
实验：
2、在信号每两个相邻样本之间插入一个零值，扩充为16点序列，作DFT，画出幅度谱和相位谱，并与原序列的DFT进行比较。
3、将信号以8为周期扩展，得到长为16的两个周期，作DFT，画出幅度谱和相位谱，并与原序列的DFT进行比较。
实验：
n=0:7;
x=(0.9*exp(j*pi/3)).^n;
x1=[x zeros(1,8)];x2=[x1 zeros(1,16)];
Xk=dft(x,8);Xk1=dft(x1,16);Xk2=dft(x2,32);
w=0:pi/200:2*pi;X=dtft(x,n,w);
magX=abs(X);angX=angle(X);
magXk=abs(Xk);angXk=angle(Xk);magXk1=abs(Xk1);angXk1=angle(Xk1);
magXk2=abs(Xk2);angXk2=angle(Xk2);
subplot(4,2,1);plot(w/pi,magX);xlabel(w/pi);ylabel(幅度|X|);grid on;
subplot(4,2,2);plot(w/pi,angX);xlabel(w/pi);ylabel(相位(rad/π));
subplot(4,2,3);stem(0:7,magXk);xlabel(k);ylabel(幅度|X(k)|);axis([0,8,0,6]);
subplot(4,2,4);stem(0:7,angXk);xlabel(k);ylabel(相位(rad/π));axis([0,8,-2,2]);
subplot(4,2,5);stem(0:15,magXk1);xlabel(k);ylabel(幅度|X1(k)|);axis([0,16,0,6]);
subplot(4,2,6);stem(0:15,angXk1);xlabel(k);ylabel(相位(rad/π));axis([0,16,-2,2]);
subplot(4,2,7);stem(0:31,magXk2);xlabel(k);ylabel(幅度|X2(k)|);axis([0,32,0,6]);
subplot(4,2,8);stem(0:31,angXk2);xlabel(k);ylabel(相位(rad/π));axis([0,32,-2,2]);
n=0:7;
x=(0.9*exp(j*pi/3)).^n;
for i=1:8 %将原序列每两个相邻样本之间插入一个零值，扩展为长为16的序列
y1(2*i-1:2*i)=[x(i) 0];
end
Yk1=dft(y1,16);
Warning: Function call dft invokes inexact match D:\MATLAB7\work\DFT.M.
magYk1=abs(Yk1);angYk1=angle(Yk1);
y2=[x x]; %将信号以8为周期扩展，得到长为16的两个周期
Yk2=dft(y2,16);
magYk2=abs(Yk2);angYk2=angle(Yk2);
Xk=dft(x,8);
magXk=abs(Xk);angXk=angle(Xk);
subplot(3,2,1);stem(0:7,magXk);xlabel(k);ylabel(幅度|X(k)|);
subplot(3,2,2);stem(0:7,angXk);xlabel(k);ylabel(相位(rad/π));
subplot(3,2,3);stem(0:15,magYk1);xlabel(k);ylabel(幅度|Y1(k)|);
subplot(3,2,4);stem(0:15,angYk1);xlabel(k);ylabel(相位(rad/π));
subplot(3,2,5);stem(
文档评论（0）