Fs300Hz和Fs200Hz的程序与上面Fs1000Hz完全相同。
实验结果运行如下：
（2）频域采样理论的验证。给定信号如下：

10
13
x
27
14
26
0
其它
编写程序分别对频谱函数XejFTx
在区间02上等间隔采样32
《数字信号处理》实验报告
3
f计算机科学与工程学院
和16点，得到X32k和X16k：
X32kXej2kk0123132
X16kXej2kk0121516
再分别对X32k和X16k进行32点和16点IFFT，得到x32
和x16
：x32
IFFTX32k32
01231x16
IFFTX16k16
01215
分别画出Xej、X32k和X16k的幅度谱，并绘图显示x
、x32
和x16
的波形，
进行对比和分析，验证总结频域采样理论。
频域采样理论验证程序exp2bmM27N32
0M产生M长三角波序列x
xa0floorM2xbceilM2110x
xaxbXkfftx
10241024点FFTx
用于近似序列x
的TFX32kfftx
3232点FFTx
x32
ifftX32k32点IFFTX32k得到x32
X16kX32k12N隔点抽取X32k得到X16Kx16
ifftX16kN216点IFFTX16k得到x16
subplot322stem
x
boxo
titleb三角波序列x
xlabel
ylabelx
axis032020k01023wk2k1024subplot321plotwkabsXktitleaFTx
xlabelomegapiylabelXejomegaaxis010200k0N21subplot323stemkabsX16kboxo
titlec16点频域采样xlabelkylabelX_1_6kaxis080200
10N21subplot324stem
1x16
boxo
titled16点IDFTX_1_6kxlabel
ylabelx_1_6
axis032020k0N1subplot325stemkabsX32kboxo
titlee32点频域采样xlabelkylabelX_3_2kaxis0160200
10N1subplot326stem
1x32
boxo

titlef32点IDFTX_3_2kxlabel
ylabelx_3_2
axis032020
《数字信号处理》实验报告
4
f试验运行结果如下：
计算机科学与工程学院
三、思考题
如果序列x
的长度为M，希望得到其频谱Xej在02上的N点等间隔采样，当
NM时，如何用一次最少点数的DFT得到该频谱采样？先对原序列x
以N为周期进行周期延拓后取主值区序列，
xN
x
iNRN
i
《数字信号处理》实验报告
5
f计算机科学与工程学院
再计算N点DFT则得到N点频域采样：
XNkDFTxN
NXej2kk012N实验总结
N1
通过这次实验，学习了线性卷积的计算机编程方法，利用卷积的方法观察、分析时域采样与频域采样特性。对matlab有了更多的了解。
《数字信号处理》实验报告
6
fr