《数字信号处理实验报告》由会员分享，可在线阅读，更多相关《数字信号处理实验报告（15页珍藏版）》请在金锄头文库上搜索。
1、实 验 报 告实验名称 _课程名称 院 系 部: 专业班级： 学生姓名： 学 号: 同 组 人: 实验台号: 指导教师： 成 绩： 实验日期: 华北电力大学（北京）一：实验目的及要求 深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义，熟练掌握利用MATLAB分析离散系统的时域响应、频响特性和零极点的方法。掌握利用DTFT和DFT确定系统特性的原理和方法。 二：实验仪器仪器名称规格/型号数量备注计算机一台装有Matlab三：实验原理 MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。1.
2、 离散系统的时域响应离散时间LTI系统可用如下的线性常系数差分方程来描述： 已知输入信号xk以及系统初始状态y-1，y-2 .，就可以求出系统的响应。MATLAB提供了计算差分方程的数值解的函数，可以计算上述差分方程描述的离散系统的单位脉冲响应、零状态响应和完全响应等。在调用MATLAB函数时，需要利用描述该离散系统的系数函数。对差分方程进行Z变换即可得系统函数： 在MATLAB中可使用向量a和向量b分别保存分母多项式和分子多项式的系数： 这些系数均从z0按z的降幂排列。离散系统的单位冲激响应hk的计算 h,k = impz(b, a)：计算系统的单位脉冲响应hk和相应的时间向量k；也可简写为
3、：h = impz(b, a)。其中： h,k = impz(b, a, n)：计算n点单位脉冲响应hk； 也可简写为：h = impz(b, a, n)。impz(b, a)：绘制单位脉冲响应hk的图形离散系统响应yk的计算 y = filter(b, a, x)： 计算系统在输入x作用下的零状态响应yk； y = filter(b, a, x, zi)：计算系统在输入x和初始状态作用下的完全响应yk。zi是由系统的初始状态经过filtic函数转换而得到的初始条件：zi= filtic (b, a, Y0) ，Y0为系统的初始状态，Y0= y-1，y-2，y-3，.。2离散系统的系统函数零极
4、点分析离散LTI系统的系统函数H(z)可以表示为零极点形式：使用MATLAB提供的roots函数计算离散系统的零极点；使用zplane函数绘制离散系统的零极点分布图。注意：在利用这些函数时，要求H(z)的分子多项式和分母多项式的系数的个数相等，若不等则需要补零。3离散系统的频率响应 当离散因果LTI系统的系统函数H(z)的极点全部位于z平面单位圆内时，系统的频率响应可由H(z)求出，即H, w=freqz(b, a, n)： 计算系统的n点频率响应H，w为频率点向量。H=freqz(b, a, w) ：计算系统在指定频率点向量w上的频响；freqz(b,a)： 绘制频率响应曲线。其中：b和a分
5、别为系统函数H(z)的分子分母系数矩阵；4利用DTFT和DFT确定离散系统的特性 在很多情况下，需要根据LTI系统的输入和输出对系统进行辨识，即通过测量系统在已知输入xk激励下的响应yk来确定系统的特性。 若系统的脉冲响应为hk，由于存在yk=xk* hk，因而可在时域通过解卷积方法求解hk。 在实际应用中，进行信号解卷积比较困难。因此，通常从频域来分析系统，这样就可以将时域的卷积转变为频域的乘积，从而通过分析系统输入序列和输出序列的频率特性来确定系统的频率特性 ，再由 得到系统的脉冲响应hk。四：实验步骤第一：实验内容1.已知某LTI系统的差分方程为：（1）初始状态 ，输入 计算系统的完全响
6、应。（2）当以下三个信号分别通过系统时，分别计算离散系统的零状态响应：（3）该系统具有什么特性？2.已知某因果LTI系统的系统函数为：（1）计算系统的单位脉冲响应。（2）当信号 通过系统时，计算系统的零状态响应。3.已知LTI系统的输入输出序列分别为（1）利用系统辨识原理从频域分别求解系统的单位脉冲响应 。（2）利用解析方法确定 ，并求出系统的单位脉冲响应 。比较解析方法与系统辨识方法得到的系统单位脉冲响应，分析误差原因。4.已知某离散系统的输入输出序列。输入序列：2，0.8333，0.3611，0.162，0.0748，0.0354，0.017，0.0083，0.0041，0.002，0.0
7、01，0.0005，0.0002，0.0001，0.0001，后面的数值均趋于0；输出序列：0.0056，-0.0259，0.073，-0.1593，0.297， -0.4974，0.7711，-1.1267，1.5702，-2.1037，2.724，-3.4207，4.174，-4.9528，5.7117，-6.3889，6.9034，-7.1528，7.012，-6.3322，4.9416，-2.648，-0.7564，5.4872，-11.7557，19.7533，-29.6298，41.4666，-55.2433，70.7979，-87.7810（1）绘出输入输出信号的波形。（2）计算
8、该系统的频率响应 ，并绘出其幅频特性。（3）计算该系统的单位脉冲响应 ，并绘出其波形。第二：实验结果，波形：第一题（1）%xk=uk求全响应N=100; b=1,-1.143,0.412; a=0.0675,0.1349,0.0675;x=ones(1,N);zi=filtic(b,a,1,2);y=filter(b,a,x,zi)plot(y)xlabel(k);ylabel(幅值);title(N=100 xk=uk时的全响应)（2）X1 作用时%N=100时在X1k作用下的响应N=100; k=0:N-1;b=1,-1.143,0.412; a=0.0675,0.1349,0.0675;
9、x1=cos(pi*k/10);zi=filtic(b,a,1,2);y=filter(b,a,x1,zi)plot(y)xlabel(k);ylabel(幅值);title(N=100时 在X1k作用下的响应)X2 作用时%N=100时在X2k作用下的响应N=100; k=0:N-1;b=1,-1.143,0.412; a=0.0675,0.1349,0.0675;x2=cos(pi*k/5);zi=filtic(b,a,1,2);y=filter(b,a,x2,zi)plot(y)xlabel(k);ylabel(幅值);title(N=100时 在X2k作用下的响应)X3 作用时%N=1
10、00时在X3k作用下的响应N=100; k=0:N-1;b=1,-1.143,0.412; a=0.0675,0.1349,0.0675;x3=cos(7*pi*k/10);zi=filtic(b,a,1,2);y=filter(b,a,x3,zi)plot(y)xlabel(k);ylabel(幅值);title(N=100时 在X3k作用下的波形)第二题(1)%单位脉冲响应N=40;k=0:N-1;a=1,-1.035,0.8246,-0.2605,0.04033;b=0.03571,0.1428,0.2143,0.1428,0.3571; y=impz(b,a,N); stem(y)xl
11、abel(k);ylabel(幅值);title(N=40 单位脉冲响应hk)波形图为：(2)%零状态响应N=40;k=0:N-1;a=1,-1.035,0.8246,-0.2605,0.04033;b=0.03571,0.1428,0.2143,0.1428,0.3571; x=ones(N); y=filter(b,a,x) plot(y);xlabel(k);ylabel(幅值);grid;title(N=40 零状态响应y0k曲线)波形图为：第三题%3.1(a)求解单位脉冲响应N=50; k=0:N-1;x=(0.5).k;y=0.25*(0.5).k+(0.5).k;X=fft(x,
12、N);Y=fft(y,N);H=Y./X;h=ifft(H,N);plot(h,k);xlabel(k);ylabel(幅值Y);title(单位脉冲响应3.1(a)（2）%3.1(b)求解单位脉冲响应N=50; k=0:N-1;x=(0.25).k;y=(0.25).k-(0.25).(k-1);X=fft(x,N);Y=fft(y,N);H=Y./X;h=ifft(H,N);plot(h);xlabel(k);ylabel(幅值Y);title(单位脉冲响应3.1(b)第四题(1)%第四题 绘制输入信号波形clearN=30;k=0:N;x=2.0000,0.8333,0.3611,0.1
13、620,0.0748,0.0354,0.0170,0.0083,0.0041,0.0020,0.0010,0.0005,0.0002,0.0001,0.0001,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000;plot(x);xlabel(k);ylabel(幅值);title(输入信号波形);波形图为：(2) %第四题 系统频域响应的波形%系统频域响应的波形clearN=30;k=0:N;x= 2.0000,0.8333,0.3611,0.1620,0.0748,0.0354,0.0170,0.0083, 0.0
