2022年2022年离散时间系统特性分析 .pdf

《2022年2022年离散时间系统特性分析 .pdf》由会员分享，可在线阅读，更多相关《2022年2022年离散时间系统特性分析 .pdf（5页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1 实验五实验报告实验名称 ：离散时间系统特性分析名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 5 页 - - - - - - - - - 2 一、实验目的：1 。深入理解单位样值响应，离散系统的频率响应的概念；2。 掌握通过计算机进行求得离散系统的单位样值响应，以及离散系统的频率响应的方法。二、实验原理：对于离散系统的单位样值而言，在实际处理过程中，不可能选取无穷多项的取值。往往是选取有限项的取值， 当然这里会产生一个截尾误差，但只要这个误差在相对小一个范围里，可以忽

2、略不计。另外，在一些实际的离散系统中，往往不是事先就能得到描述系统的差分方程的，而是通过得到系统的某些相应值，则此时系统的分析就需借助计算机的数值处理来进行，得到描述系统的某些特征，甚至进而得到描述系统的数学模型。本实验首先给出描述系统的差分方程，通过迭代的方法求得系统的单位样值响应，进而求得该离散系统的频率响应。限于试验条件，虽然给出了系统方程，但处理的方法依然具有同样的实际意义。具体的方法是：1 在给定系统方程的条件下，选取激励信号为(n) ，系统的起始状态为零状态，通过迭代法，求得系统的单位样值响应h(n) （n=0, ,N ） 。2 利用公式其中 的取值范围为02 。计算系统的频率响应

3、。三、实验内容1 已知系统的差分方程为利用迭代法求得系统的单位样值响应，取N10。2 利用公式其中 的取值范围为02 。四、程序流程图和代码T n=n+1 hn=1.3hn-1-0.4hn-2 +x n-1 n0) 开始定义 h20=0,x20 输出F 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 5 页 - - - - - - - - - 3 #include #include #define N 10 #define M 20 #define pi 3.141592

4、6 struct pinlv double fu; double xiang; double hN+1,xN+1; struct pinlv PL(double w) double a=0, b=0,fu,xiang; int k; struct pinlv FX; for(k=0;k=N;k+) a=a+hk*cos(-k*w); b=b+hk*sin(-k*w); fu=sqrt(a*a+b*b); xiang=atan(b/a); if(a0) xiang=xiang+pi; if(a0)&(b0) xiang=xiang-pi; FX.fu=fu; FX.xiang=xiang; re

5、turn(FX); main() int i,j; double w0; struct pinlv FXM+1; FILE *fp1,*fp2; 单位样值响应 .txt,w); 频率特性 .txt,w); h-1=0;h-2=0; for(i=-1;i=N;i+) xi=0; x0=1; for(i=0;i=N;i+) hi=1.3*hi-1-0.4*hi-2+xi-1; printf(系统的单位样值响应为n); fprintf(fp1,系统的单位样值响应( 从 x0 开始 ) 为n); fprintf(fp1,激励 xi 响应 yin); for(i=0;i=N;i+) 名师资料总结 - -

6、 -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 5 页 - - - - - - - - - 4 printf(x%d=%6.3fth%d=%6.3fn,i,xi,i,hi); fprintf(fp1,%6.3ft%6.3fn,xi,hi); printf(n); printf(系统的幅频和相频特性为n); fprintf(fp2,系统的幅频和相频特性为n); fprintf(fp2,w0t幅度值 t 相位值 n); for(j=0;j=M;j+) w0=j*0.1*pi; FXj.fu=PL

7、(w0).fu; FXj.xiang=PL(w0).xiang; printf(w0=%6.3ft fu=%6.3ft xiang=%6.3fn,w0,FXj.fu,FXj.xiang); fprintf(fp2,%6.3ft%6.3ft%6.3fn,w0,FXj.fu,FXj.xiang); printf(n); 五、实验所得数据如下：单位样值响应激励 xi 响应 yi 1 0 0 1 0 1.3 0 1.29 0 1.157 0 0.988 0 0.822 0 0.673 0 0.546 0 0.441 0 0.355 幅频特性w0 幅度值相位值0.000 8.572 0.000 0.31

8、4 6.118 -1.365 0.628 2.137 -2.114 0.942 1.773 -2.252 1.257 0.743 -2.736 1.571 0.903 -2.599 1.885 0.375 -3.022 2.199 0.630 -2.833 2.513 0.247 -3.136 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 5 页 - - - - - - - - - 5 2.827 0.538 -3.040 3.142 0.212 3.142 3.456

9、 0.538 3.040 3.770 0.247 3.136 4.084 0.630 2.833 4.398 0.375 3.022 4.712 0.903 2.599 5.027 0.743 2.736 5.341 1.773 2.252 5.655 2.137 2.114 5.969 6.118 1.365 6.283 8.572 0.000 六、幅频特性和相频特性曲线并对系统进行分析。幅频特性曲线相频特性曲线虽然离散系统的信号是离散的，但是其频谱图是连续的，并且 0 时对应的是低频，Pi（3.14 ）时对应的为高频。根据幅频图可以看出该系统为低通滤波器。01234567891002468系列 1-4-3-2-10123402468系列 1名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 5 页 - - - - - - - - -