数字信号处理实验教案.docx

数字信号处理实验教案 数字信号处理试验教案 信息工程学院-通信工程教研室 数字信号处理是一门理论和实际亲密结合的课程，为深化驾驭课程内容，最好在学习理论的同时，做习题和上机试验。上机试验不仅可以帮助读者深化的理解和消化基本理论，而且能熬炼同学们的独立解决问题的实力。本讲义在第三版的基础上编写了五个试验，前2个试验属基础性的验证性试验，第 3、 4、5个试验属基本应用综合性试验。 试验一 离散时间信号的MATLAB实现 试验二 线性卷积与循环卷积的原理及应用 试验三 频率采样定理 试验四 离散系统的因果性和稳定性及频率响应特性 试验五 基于MATLAB的快速傅里叶变换 依据教学进度，理论课结束后进行相关试验。 试验一 时域离散信号的产生 一 试验目的 (1)了解常用的时域离散信号及其特点。 (2)驾驭MATLAB产生常用时域离散信号的方法。 二 试验内容 (1) 编写程序，产生下列离散序列： A.f(n)=(n) (-3 Bf(n)=e(0.1+j1.6)n (0 （2） 一个连续的周期性三角波信号频率为50Hz,信号幅度在0+2V之间，在窗口上显示2个周期信号波形，对信号的一个周期进行16点采样来获得离散信号。试显示原连续信号和采样获得的离散信号波形。 (3)一个连续的周期性方波信号频率为200Hz,信号幅度在-1+1V之间，在窗口上显示2个周期信号波形，用Fs=4kHz的频率对连续信号进行采样，试显示原连续信号和采样获得的离散信号波形。 三 试验步骤 (1) 在matlab吩咐窗口中逐行输入下列语句 >> n1=-3;n2=4;n0=0; %在起点n 1、终点n2的范围内，于n0处产生冲激 >> n=n1:n2; %生成离散信号的时间序列 >> x=n=n0; %生成离散信号x(n) >> stem(n,x,'filled'); %绘制杆状图，且圆心处用实心圆表示 >> title('单位脉冲序列'); >> xlabel('时间（n）');ylabel('幅度x（n）'); 在上述语句输入完成之后，敲击回车键，弹出图形窗口，显示出如下图形，即已经满意题干所述条件，产生了 f(n)=(n)，(-3 >> n1=16;a=0.1;w=1.6*pi; >> n=0:n1; >> x=exp(a+j*w)*n); >>subplot(2,1,1),stem(n,real(x); %在指定位置描绘图像 >> title('复指数序列的实部'); >> subplot(2,1,2),stem(n,imag(x); >> title('复指数序列的虚部'); 在上述语句输入完成之后，敲击回车键，弹出图形窗口，显示出如下图形，即已经满意题干所述条件，产生了f(n)=e(0.1+j1.6)n，(0 >> f=50;Um=1;nt=2; %输入信号频率、振幅、显示周期 >> N=16;T=1/f; %N为信号一个采样周期的采样点数，T为信号周期 >> dt=T/N; %采样时间间隔 >> n=0:nt*N-1; %建立离散时间的时间序列 >> tn=n*dt; %确定时间序列样点在时间轴上的位置 >> f=Um*sawtooth(2*f*pi*tn)+1; >> subplot(2,1,1),stem(tn,f); %显示经采样的信号 >> title('离散信号'); >> subplot(2,1,2),plot(tn,f); %显示原连续信号 >> title('连续信号'); 在上述语句输入完成之后，敲击回车键，弹出图形窗口，显示出如下图形，即已经满意题干所述条件，显示了原连续信号和采样获得的离散信号波形 (4) 在matlab吩咐窗口中逐行输入下列语句 >> f=200;Um=1;nt=2; %输入信号频率、振幅、显示周期 >> Fs=4000;N=Fs/f;T=1/f; %输入采样频率、求采样点数N、T为信号周期 >> dt=T/N; %采样时间间隔 >> n=0:nt*N-1; %建立离散时间的时间序列 >> tn=n*dt; %确定时间序列样点在时间轴上的位置 >> f=Um*sin(2*f*pi*tn); >> subplot(2,1,2),plot(tn,f); %显示原连续信号 >> title('连续信号'); >> subplot(2,1,1),stem(tn,f); %显示经采样的信号 >> title('离散信号'); 在上述语句输入完成之后，敲击回车键，弹出图形窗口，显示出如下图形，即已经满意题干所述条件，显示了原连续信号和采样获得的离散信号波形 四 思索题 (1) 如何在matlab下生产f(n)=3sin(n/4)(0 （2）变更试验步骤中最终两个试验的频率参数，分别重新生成相关的信号？ 试验二 线性卷积与循环卷积的原理及应用 一、试验目的 （1）驾驭两种卷积的原理和两者的异同。 （2）驾驭MATLAB实现两种卷积的计算和比较。 二、试验内容 （1）用MATLAB设计线性卷积； （2）调试写出线性卷积和源代码； （3）用MATLAB设计循环卷积； 三 试验步骤 1 线性卷积定理 1）线性卷积的引入 在实际应用中，为了分析时域离散线性非移变系统或者对序列进行滤波处理等，须要计算两个序列的线性卷积。线性卷积既可以在时域中干脆计算，也可以通过变换在频域中计算得到。 2）线性卷积的时域计算方法 计算卷积的基本运算是翻转、移位、相乘和相加，这类卷积称为序列的线性卷积。假如两个序列的长度为N和M，那么卷积结果的长度为N+M-1。线性卷积有四步运算：卷积运算时，y(n)要先反折得到y(-n);m>0表示y(-n)序列右移，m 假设h(n)和x(n)都是有限长序列，长度分别为N和M，它们的线性卷积可以表示如下： y2 循环卷积定理 l=h(n)*x(n)=m=0åh(m)x(n-m)N-1 MATLAB信号处理工具箱供应了conv函数，该函数用于计算两个有限序列的卷积。 1）循环卷积的引入 为了提高线性卷积的速度，希望用DFT(FFT)计算线性卷积。从而引入循环卷积来运用DFT快速计算线性卷积。循环卷积运用到离散傅立叶变换的循环移位性质，即时域循环移位定理。 2）循环卷积的时域计算方法 假设h(n)和x(n)都是有限长序列，长度分别为N和M，它们的L点循环卷积可以表示如下： M, L³maxNL称为循环卷积区间长度。n和m的改变区间均是0,L-1，干脆计算该式比较麻烦。计cm=0LLy=h(n)*x(n)=åh(m)L-1x(n-m)R(n)算机中采纳矩阵相乘或快速傅里叶变换（FFT）的方法计算循环卷积。用矩阵相乘的方法计算两个序列的循环卷积，这里关键是先形成循环卷积矩阵。假如h(n)的长度N 3）线性卷积与循环卷积的关系 y上式说明，cc(n)=q=-¥åy(n+qL)Rl¥L(n) ly(n)等于y(n)以L为周期的周期延拓序列的主值序列。y(n)的长度 l为N+M-1，因此只有当循环卷积长度L³N+M-1时， cly(n)以L为周期进行周期延 l(n)y(n)y拓才无混叠现象。因此取其主值序列就满意=。即满意条件L³N+M-1。 四 思索题 (1)比较线性卷积和循环卷积在序列长度不同时两者的联系？ (2)试着写出循环卷积的源代码？ 试验三 时域采样理论与频域采样定理验证 一、试验目的 1时域采样理论与频域采样理论是数字信号处理中的重要理论。要求驾驭模拟信号采样前后频谱的改变，以及如何选择采样频率才能使采样后的信号不丢失信息；要求驾驭频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。 二、试验原理及方法 时域采样定理的要点是： (a)对模拟信号 (jW)xa(t)以间隔T进行时域等间隔志向采样，形成的采样信号的频谱X是原模拟信号频谱为： Xa(jW)以采样角频率Ws（Ws=2p/T）为周期进行周期延拓。公式 1¥=åXa(jW-jnWs) X(jW)=FTxa(t)Tn=-¥ a （b）采样频率Ws必需大于等于模拟信号最高频率的两倍以上，才能使采样信号的 频谱不产生频谱混叠。 利用计算机计算上式并不便利，下面我们导出另外一个公式，以便用计算机上进行试验。 志向采样信号a(t)xx(t)之间的关系为： 和模拟信号a¥ a(t)=xa(t)åd(t-nT)xn=-¥ 对上式进行傅立叶变换，得到： ¥Xa(jW)=òxa(t)åd(t-nT)e-jWtdt-¥n=-¥¥ å¥n=-¥ò¥-¥xa(t)d(t-nT)e-jWtdt 在上式的积分号内只有当t=nT时，才有非零值，因此： (jW)=Xa 上式中，在数值上 n=-¥åx¥a(nT)e-jWnT xa(nT)x(n)，再将w=WT代入，得到： (jW)=Xa n=-¥åx(n)e¥-jwn jwX(e)，即 上式的右边就是序列的傅立叶变换 (jW)=X(ejw)Xaw=WT 上式说明志向采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变量用WT代替即可。 频域采样定理的要点是： a) 对信号x(n)的频谱函数X(ej)在0，2上等间隔采样N点，得到 XN(k)=X(ejw) 则N点IDFT列，公式为： w=2pkN , k=0,1,2,L,N-1 XN(k)得到的序列就是原序列x(n)以N为周期进行周期延拓后的主值区序 ¥ xN(n)=IDFTXN(k)N=åx(n+iN)RN(n)i=-¥ (b)由上式可知，频域采样点数N必需大于等于时域离散信号的长度M(即NM)，才能使时域不产生混叠，则N点IDFT果N>M， XN(k)得到的序列xN(n)就是原序列x(n),即xN(n)=x(n)。如xN(n)比原序列尾部多N-M个零点；假如N 了时域混叠失真，而且 在数字信号处理的应用中，只要涉刚好域或者频域采样，都必需听从这两个采样理论的要点。 对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓”。因此放在一起进行试验。 三 试验步骤 （1）时域采样理论的验证。 -atx(t)=Aesin(W0t)u(t) a给定模拟信号， W式中A=444.128，a=502，0=502rad/s，它的幅频特性曲线如图10.2.1 图10.2.1 xa(t)的幅频特性曲线 现用DFT(FFT)求该模拟信号的幅频特性，以验证时域采样理论。 安照xa(t)的幅频特性曲线，选取三种采样频率，即Fs=1kHz，300Hz，200Hz。观测时间选Tp=50ms。为运用DFT，首先用下面公式产生时域离散信号，对三种采样频率，采样序列按依次用x1(n)，x2(n)，x3(n)表示。 -anTx(n)=x(nT)=Aesin(W0nT)u(nT) a 因为采样频率不同，得到的x1(n)，x2(n)，公式 x3(n)的长度不同， 长度（点数）用N=Tp´Fs计算。选FFT的变换点数为M=64，序列长度不够64的尾部加零。 X(k)=FFTx(n) ， k=0,1,2,3,-,M-1 式中k代表的频率为 wk=2pkM。 要求： 编写试验程序，计算x1(n)、x2(n)和析频谱混叠失真。 （2）频域采样理论的验证 给定信号如下： x3(n)的幅度特性，并绘图显示。视察分 ìn+10£n£13ïx(n)=í27-n14£n£26ï0其它î jwX(e)=FTx(n)在区间0,2p上等间隔采样32 编写程序分别对频谱函数和16点，得到X32(k)和X16(k)： X32(k)=X(ejw) w=2pk32 , k=0,1,2,L31 X16(k)=X(ejw) 再分别对 w=p216k , k=0,1,2,L15 X32(k)和X16(k)进行32点和16点IFFT，得到x32(n)和x16(n)： x32(n)=IFFTX32k(3)2 =n , L0, 1,2,x16(n)=IFFTX16k(1)6 =n , L0, 1,2,31,15jwX(e)、X32(k)和X16(k)的幅度谱，并绘图显示x(n)、x32(n)和x16(n)的波形，分别画出进行对比和分析，验证总结频域采样理论。 提示：频域采样用以下方法简单变程序实现。 jwX(k)=FFTx(n)X(e)在0,2p的3232 干脆调用MATLAB函数fft计算就得到32点频率域采样 抽取X32(k)的偶数点即可得到X(ejw)在0,2p的16点频率域采样X16(k)，即X16(k)=X32(2k) , k=0,1,2,L,15。 3 当然也可以根据频域采样理论，先将信号x(n)以16为周期进行周期延拓，取其主值区（16 jwX(e)在0,2p的16点频率域采样X16(k)。点），再对其进行16点DFT(FFT),得到的就是 四 思索题 假如序列x(n)的长度为M，希望得到其频谱()jXe在2,0上的N点等间隔采样，当N 试验四 离散系统的零极点分析 一、试验目的 1 熟识MATLAB的仿真及应用环境 2 在MATLAB的环境下探讨限制系统稳定性 二、试验内容和要求 1.了解离散系统的零极点与系统因果性和稳定性的关系。 2.视察离散系统零极点对系统冲激响应的影响。 3.熟识MATLAB中进行离散系统零极点分析的常用子函数。 三 试验步骤 一）MATLAB子函数 1.zplane 功能：显示离散系统的零极点分布图。 调用格式： zplane(z，p)；绘制由列向量z确定的零点、列向量p确定的极点构成的零极点分布图。 zplane(b，a)；绘制由行向量b和a构成的系统函数确定的零极点分布图。 hz，hp，htzplane(z，p)；执行后可得到3个句柄向量：hz为零点线句柄，hp为极点线句柄，ht为坐标轴、单位圆及文本对象的句柄。 2.roots 功能：求多项式的根。 调用格式： rroots(a)；由多项式的分子或分母系数向量求根向量。其中，多项式的分子或分母系数按降幂排列，得到的根向量为列向量。 二）试验原理 1.离散系统的因果性和稳定性 1)因果系统 由理论分析可知，一个离散系统的因果性在时域中必需满意的充分必要条件是： h(n)0 n 即系统的冲激响应必需是右序列。 在变换域，极点只能在z平面上一个有界的以原点为中心的圆内。假如系统函数是一个多项式，则分母上z的最高次数应大于分子上z的最高次数。 2)稳定系统 在时域中，离散系统稳定的充分必要条件是：它的冲激响应肯定可加，即 ¥ h(n)<¥ n=0 在变换域，则要求全部极点必需在z平面上以原点为中心的单位圆内。 3)因果稳定系统 综合系统的因果性和稳定性两方面的要求可知，一个因果稳定系统的充分必要条件是：系统函数的全部极点必需在z平面上以原点为中心的单位圆内。 2.系统极点的位置对系统响应的影响 系统极点的位置对系统响应有着特别明显的影响。下面举例说明系统的极点分别是实数 å和复数时的状况，运用MATLAB供应的zplane子函数制作零极点分布图进行分析。 3.系统的因果稳定性实例分析 在MATLAB中供应了roots子函数，用于求多项式的根。协作运用zplane子函数制作零极点分布图，可以帮助我们进行系统因果稳定性的分析。 4.试验任务 (z-0.3)H3(z)= (z+1-j)(z+1+j) 4-1.6z-1-1.6z-2+4z-3H2(z)= 1+0.4z-1+0.35z-2-0.4z-3 求该系统的零极点及零极点分布图，并推断系统的因果稳定性。 四 思索题 1结合本次试验与书本上相关原理，对书本后面的习题进行相关的matlab软件仿真？ 2因果稳定的离散系统必需满意的充分必要条件是什么？MATLAB供应了哪些进行零极点求解的子函数？如何运用？ 试验五 基于MATLAB的快速傅里叶变换 一 试验目的 学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的误差及其缘由，以便正确应用FFT。 二 试验原理 用FFT对信号作频谱分析是学习数字信号处理的重要内容。常常须要进行谱分析的信号是模拟信号和时域离散信号。对信号进行谱分析的重要问题是频谱辨别率D和分析误差。频谱辨别率干脆和FFT的变换区间N有关，因为FFT能够实现的频率辨别率是N/2，因此要求DN/2。可以依据此式选择FFT的变换区间N。误差主要来自于用FFT作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N较大时离散谱的包络才能靠近于连续谱，因此N要适当选择大一些。 周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT，得到的离散谱才能代表周期信号的频谱。假如不知道信号周期，可以尽量选择信号的视察时间长一些。 对模拟信号进行谱分析时，首先要根据采样定理将其变成时域离散信号。假如是模拟周期信号，也应当选取整数倍周期的长度，经过采样后形成周期序列，根据周期序列的谱分析进行。 三 试验步骤及内容 （1）对以下序列进行谱分析。 x1(n)=R4(4) n+1 0n3 x2(n)= 8-n 4n7 0 其它n 4-n 0n3 X3(n)= n-3 4n7 0 其它n 选择FFT的变换区间N为8和16 两种状况进行频谱分析。分别打印其幅频特性曲线。 并进行对比、分析和探讨。 （2）对以下周期序列进行谱分析。 x3(n)=cospn4 选择FFT的变换区间N为8和16 两种状况分别对以上序列进行频谱分析。分别打印其幅频特性曲线。并进行对比、分析和探讨。 （3）对模拟周期信号进行谱分析 x4(n)=cos8pt+cos16pt+cos20pt 选择 采样频率HzFs64=，变换区间N=16,32,64 三种状况进行谱分析。分别打印其幅频特性，并进行分析和探讨。 四 思索题 （1）对于周期序列，假如周期不知道，如何用FFT进行谱分析？ （2）如何选择FFT的变换区间？（包括非周期信号和周期信号） 数字信号处理试验教案 数字信号处理试验 数字信号处理教案 数字信号处理试验5 数字信号处理试验4 数字信号处理试验讲稿 数字信号处理试验二 数字信号处理试验报告 数字信号处理试验报告 数字信号处理试验报告 本文来源：网络收集与整理，如有侵权，请联系作者删除，谢谢！第21页 共21页第 21 页 共 21 页第 21 页 共 21 页第 21 页 共 21 页第 21 页 共 21 页第 21 页 共 21 页第 21 页 共 21 页第 21 页 共 21 页第 21 页 共 21 页第 21 页 共 21 页第 21 页 共 21 页