数字信号处理实验

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) 按实验内容要求，上机实验，并写出实验报告。 >> %