数字信号处理实验2013教程.docVIP

实验一 离散傅里叶变换的性质 实验目的 掌握离散傅里叶变换的性质，包括线性特性、时移特性、频移特性、对称性和循环卷积等性质； 通过编程验证傅里叶变换的性质，加强对傅里叶变换性质的认识。 实验原理和方法 线性特性 时移特性 频移特性 对称性 设由x(n) 延拓成的周期序列为 则 共轭对称序列 共轭反对称序列 将 和 截取主周期，分别得 则 x(n)序列的实部和虚部的离散立叶变换 当x(n)为实数序列 循环卷积 有限长序列线性卷积与循环卷积的关系 x1(n)和x2(n)的线性卷积： 将x1(n)和x2(n)延拓成以N为周期的周期序列 则它们的周期卷积为 x1(n)和x2(n)周期延拓后的周期卷积等于他们的线性卷积的周期延拓。 实验内容和步骤 x1[n]，x2[n] 为长度N=8的实序列，x1[n]=[1 3 5 3 6 8 3 9]，x2[n]=[2 4 3 6 7 9 0 2 ]， 采用MATLAB编程验证傅里叶变换的如下性质 线性特性 给出序列x1[n]的傅里叶变换X1[k]，并画出其幅度谱和相位谱 给出序列x2[n]的傅里叶变换X2[k] ，并画出其幅度谱和相位谱 给出序列Z=2*X1[k]+6*x2[k],并与序列2*x3[n]+6*x4[n]的傅里叶变换比较， 提示：可以用fft函数实现DFT，用abs实现幅度谱，angle获得相位谱，画图用stem 时移特性 给出序列x1[n]右移3位(循环移位)后的傅里叶变换的幅度谱和相位谱，并和原始序列的幅度谱和相位谱相比较 提示：可以用b=circshift(x1,[0,3])实现x1序列的右移3位循环移位，或者用如下语句实现循环移位a=1:10 %实现a =1 2 3 4 5 6 7 8 9 10n=3;b1=[a(n+1:end),a(1:n)] %循环移位，左移3位；实现b1 =4 5 6 7 8 9 10 1 2 3 n0=length(a)-3 b2=[a(n0+1:end),a(1:n0)] %循环移位，右移3位；实现b2 = 8 9 10 1 2 3 4 5 6 7 对称性 （1）利用x1[n]构造圆周共轭对称序列和圆周共轭反对称序列（第三章学习），讨论如下问题 画出圆周共轭对称序列的傅里叶变换的幅度谱和相位谱 画出圆周共轭对称序列的傅里叶变换的实部和虚部 利用x1[n]构造共轭对称序列和共轭反对称序列（第二章学习），讨论与圆周共轭对称圆周共轭反对称的区别 提示：应先构造出圆周共轭对称或共轭反对称序列。可以用 x1=[1 3 5 3 6 8 3 9] x2=fliplr(x1) %翻褶序列，x2=[9 3 8 6 3 5 3 1] x3=circshift(x2,[1,1]) % x2的右移一位序列，注意是循环移位，x3=[1 9 3 8 6 3 5 3] xe=(x1+conj(x3))/2 %构造出圆周共轭对称序列 xo=(x1-conj(x3))/2 %构造出圆周共轭反对称序列 函数real、imag、abs、angle分别为求实部、虚部、幅度和相位。 思考：共轭对称和共轭反对称序列如何构造？ 循环卷积与线性卷积 计算序列x1[n]和x2[n]的线性卷积y[n] 计算x1[n] 和x2[n]的傅里叶变换X1[k]和X2[k]，Y[k]=X1[k]*X2[k]，求Y[k]的反傅里叶变换y2[n]（y2[n]即为x1[n]和x2[n]的循环卷积），比较y[n]与y2[n]的异同. 补零 用MATLAB计算如下N点序列的M点DFT： ??? （1）取N=8，M=8 （2）取N=8，M=16 （3）取N=8，M=32 根据实验结果，分析同一序列不同点数的离散傅里叶变换的特点。 实验报告要求 按照实验设计完成实验内容，要求保存对应实验内容的程序和图形； 记录实验过程中遇到的主要问题，并解释问题产生原因和解决办法； 实验二 时域采样与频域采样 1. 时域采样定理的验证，分析采样序列的特性。 对xa(t)=Ae-atsin(Ω0t)u(t)进行采样，可得到采样序列 xa(n)=xa(nT)=Ae-anT sin(Ω0nT)u(n), 其中A为幅度因子，a为衰减因子，Ω0是模拟角频率，T为采样间隔。 a. 取采样频率fs=1 kHz, 即T=1 ms，观察|X(ejω)|。 b. 改变采样频率，fs=300 Hz，观察|X(ejω)|的变化； c.进一步降低采样频率，fs=200 Hz，观察频谱混叠是否明显存在， 说明原因， 并记录|X(ejω)|曲线。 2. 频域采样定理的验证。 给定信号 编写程序分别对频谱函数在区间上等间隔采样32点和16点，得到X32(k)和X16(k)，再分别对其进行32点和16点的