利用DFT计算模拟信号与离散信号.docVIP

实验2利用DFT计算模拟信号和离散信号 实验目的： 通过应用DFT分析模拟信号的频谱，加深对DFT的理解； 通过应用DFT进一步加深对离散信号DTFT和DFT的概念及其相互关系的理解。 一、利用DFT计算连续周期信号的频谱 1. 实验原理 连续周期信号在满足一定条件下，可以展开为傅里叶级数： 傅里叶级数的系数按下式求得： 式中是周期信号的周期，是基本角频率（通称基频），为次谐波频率。是离散频率的复函数，称为频谱函数，可表示为： 反映了组成周期信号的不同频率谐波分量的幅度随频率变化的特性，即幅度谱；反映了不同谐波分量的初相角随频率变化的特性，即相位谱。 连续周期信号傅里叶变换的示意图如下图所示。可以看出，时域的连续函数造成频域是非周期的频谱函数，而频域的离散频谱就与时域的周期时间函数相对应。 连续周期信号的频谱求解步骤： ① 根据抽样定理，确定时域抽样间隔T; 抽样频率: 抽样间隔: 式中 为信号的最高频率 ② 计算一个周期内的抽样点数N; 式中 为频率分辨率 ③ 使用fft命令作N点FFT计算，求得; ④ 最后求得连续周期信号的频谱为 2. 实验内容 （1） 已知一连续信号为 其中。试利用DFT分析其频谱。 解：信号的最高频率，抽样频率，取抽样频率；最低的频率分辨率为，最少的信号样点数为。 的MATLAB程序如下： % program exa_1_1.m,利用矩形窗计算有限长余弦信号频谱 N=30; %数据的长度 L=512; %DFT的点数 f1=100; f2=120; fs=600; %抽样频率 T=1/fs; %抽样间隔 t=(0:N-1)*T; x=cos(2*pi*f1*t)+cos(2*pi*f2*t); y=fft(x,L); mag=abs(y); f=(0:length(y)-1)*fs/length(y); plot(f(1:L/2),mag(1:L/2)); xlabel(频率(Hz)) ylabel(幅度谱) 程序运行结果如下图所示。 由图可见，频谱图显示出两个较为明显的峰值（对应）。结论：当截取信号样点时，频率分辨率，刚好能够分辨出和两个频谱分量，但频谱泄漏较严重。 若取（达不到最低的频率分辨率）。 的MATLAB程序如下： % program exa_1_2.m,利用矩形窗计算有限长余弦信号频谱 N=20; %数据的长度 L=512; %DFT的点数 f1=100; f2=120; fs=600; %抽样频率 T=1/fs; %抽样间隔 t=(0:N-1)*T; x=cos(2*pi*f1*t)+cos(2*pi*f2*t); y=fft(x,L); mag=abs(y); f=(0:length(y)-1)*fs/length(y); plot(f(1:L/2),mag(1:L/2)); xlabel(频率(Hz)) ylabel(幅度谱) 程序运行结果如下图所示。 由图可见，频谱图显示不出对应的两个明显峰值。 结论：当截取信号样点时，频率分辨率，达不到最低的频率分辨率，故分辨不出和两个频谱分量，且频谱泄漏更为严重。 若取频率分辨率，则对应的信号样点数为。 的MATLAB程序如下： % program exa_1_3.m,利用矩形窗计算有限长余弦信号频谱 N=600; %数据的长度 f1=100; f2=120; fs=600; %抽样频率 T=1/fs; %抽样间隔 t=(0:N-1)*T; x=cos(2*pi*f1*t)+cos(2*pi*f2*t); y=fft(x); mag=abs(y); f=(0:length(y)-1)*fs/length(y); plot(f(1:length(y)/2),mag(1:length(y)/2)); xlabel(频率(Hz)) ylabel(幅度谱) 程序运行结果如下图所示。 由图可见，频谱图显示出两个特别明显的峰值（对应）。结论：当截取信号样点时，频率分辨率，高分辨率的频谱图具有较高的质量，频谱分析时必须保证获取足够的信号数据长度。 （2） 改变有关参数，进一步观察结果的变化，并加以分析说明。 二、离散信号的DFT和DTFT 1．实验原理： 序列x[n] 的DTFT定义： N点序列x[n] 的DFT定义： 在MAT