- 1、本文档共3页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
程序清单及波形显示:实验(1)
Clear; clc
x1=[ones(1,4)];
x2=[1 2 3 4 4 3 2 1];
x3=[4 3 2 1 1 2 3 4];
X1K8=fft(x1,8); X1K16=fft(x1,16);
X2K8=fft(x2,8); X2K16=fft(x2,16);
X3K8=fft(x3,8); X3K16=fft(x3,16);
subplot(3,2,1);stem(0:length(X1K8)-1,abs(X1K8),'.');title('(1a) 8点DFT[x_1(n)]');
subplot(3,2,2);stem(0:length(X1K16)-1,abs(X1K16),'.');title('(1b) 16点DFT[x_1(n)]');
subplot(3,2,3);stem(0:length(X2K8)-1,abs(X2K8),'.');title('(2a) 8点DFT[x_2(n)]');
subplot(3,2,4);stem(0:length(X2K16)-1,abs(X2K16),'.');title('(2b) 16点DFT[x_2(n)]');
subplot(3,2,5);stem(0:length(X3K8)-1,abs(X3K8),'.');title('(3a) 8点DFT[x_3(n)]');
subplot(3,2,6);stem(0:length(X3K16)-1,abs(X3K16),'.');title('(3b)16DFT[x_3(n)]');
%实验内容(2) 周期序列谱分析
clear; clc
N=8;n=0:N-1;
x4n=cos(pi*n/4);
x5n=cos(pi*n/4)+cos(pi*n/8);
X4k8=fft(x4n); X5k8=fft(x5n);
N=16;n=0:N-1;
x4n=cos(pi*n/4);
x5n=cos(pi*n/4)+cos(pi*n/8);
X4k16=fft(x4n); X5k16=fft(x5n);
subplot(2,2,1);stem(0:length(X4k8)-1,abs(X4k8),'.');title('(4a) 8点DFT[x_4(n)]');
subplot(2,2,2);stem(0:length(X4k16)-1,abs(X4k16),'.');title('(4b)16点DFT[x_4(n)]');
subplot(2,2,3);stem(0:length(X5k8)-1,abs(X5k8),'.');title('(5a) 8点DFT[x_5(n)]');
subplot(2,2,4);stem(0:length(X5k16)-1,abs(X5k16),'.');title('(5b)16点DFT[x_5(n)]');
%实验内容(3) 模拟周期信号谱分析
Fs=64;T=1/Fs; N=16;n=0:N-1;
x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);
X6k16=fft(x6nT);
Tp=N*T;F=1/Tp; k=0:N-1;fk=k*F;
subplot(3,1,1);stem(fk,abs(X6k16),'.');title('(6a) 16点|DFT[x_6(nT)]|');
axis([0,N*F-1,0,1.2*max(abs(X6k16))])
N=32;n=0:N-1;
x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);
X6k32=fft(x6nT);
Tp=N*T;F=1/Tp; k=0:N-1;fk=k*F;
subplot(3,1,2);stem(fk,abs(X6k32),'.');title('(6b) 32点|DFT[x_6(nT)]|');
axis([0,N*F-1,0,1.2*max(abs(X6k32))])
N=64;n=0:N-1;
x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);
X6k64=fft(x6nT);
Tp=N*T;F=1/Tp; k=0:N-1;fk=k*F;
subplot(3,1,3);stem(fk,abs(X6k64),'.');title('(6a) 64点|DFT[x_6(nT)]|');
axis([0,N*F-1,0,1.2*max(abs(X6k64))])
文档评论(0)