发布时间 : 星期五 文章数字信号处理实验更新完毕开始阅读c093dd45ed3a87c24028915f804d2b160a4e860e
title('|X(k)|');
xn200xn|X(k)|10005000-2000204060n,fs=1000Hzxn8000204060k,fs=1000Hz|X(k)|80200xn2001000510n,fs=200Hzxn1500510k,fs=200Hz|X(k)|150-200200xn5000-2000102030n,fs=500Hz4000102030k,fs=500Hz40
>> %频域采样定理验证
M=26;N=32;n=0:M;n1=0:13;x1=n1+1; n2=14:26;x2=27-n2; x=[x1,x2];Xk=fft(x,512); X32k=fft(x,32); k=0:511;w=(pi/512)*k;
subplot(321);stem(n,x);xlabel('n');
ylabel('xn');axis([0,31,0,15]);
subplot(322);plot(w,abs(Xk));xlabel('k'); ylabel('|X(k)|');axis([0,1,0,200]) X16k=X32k(1:2:N);
x32n=ifft(X32k);x16n=ifft(X16k,16); k1=0:31;k2=0:15;
subplot(323);stem(k1,abs(X32k));xlabel('k'); ylabel('X32k');axis([0,31,0,200]);
subplot(325);stem(k2,abs(X16k));xlabel('k'); ylabel('|X(k)|');axis([0,15,0,200]) n=0:31;
subplot(324);stem(n,abs(x32n));xlabel('n'); ylabel('|x(n)|');axis([0,31,0,15]) n1=0:15;
subplot(326);stem(n1,abs(x16n));xlabel('n');
ylabel('|x(n)|');axis([0,31,0,15])
15|X(k)|xn20010010500102030n200k32X10000102030k200|)(kX100|0051015k000.51k15|)10(n|x500102030n15|)10(n|x500102030n
实验二： 用FFT作谱分析
1.实验目的
(1) 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法， 所以FFT的运算结果必然满足DFT的基本性质)。
(2) 熟悉FFT算法原理和FFT子程序的应用。
(3) 学习用FFT对连续信号和时域离散信号进行谱分析的方法， 了解可能出现的分析误差及其原因， 以便在实际中正确应用FFT。
2.实验步骤
(1) 复习DFT的定义、 性质和用DFT作谱分析的有关内容。 (2) 复习FFT算法原理与编程思想， 并对照DIT-FFT运算流图和程序框图，读懂本实验提供的FFT子程序。
(3) 编制信号产生子程序， 产生以下典型信号供谱分析用：
(4) 编写程序。
(5) 按实验内容要求，上机实验，并写出实验报告。 >> %