北航信号与系统实验五离散时间系统频响.docxVIP

实验三?离散时间系统特性分析140222班刘嘉晖一、?实验目的?1．深入理解单位样值响应，离散系统的频率响应的感念。?2．掌握借助计算机进行求离散系统单位样值响应，以及离散系统的频率响应的方法。?二、?实验原理?对于离散系统的单位样值而言，在实际处理过程中，不可能选取无穷多项的取值。往往是选取有限项的取值，当然这里会产生一个截尾误差，但只要这个误差在相对小一个范围里，可以忽略不计。?另外，在一些实际的离散系统中，往往不是事先就能得到描述系统的差分方程的，而是通过得到系统的某些相应值，则此时系统的分析就需借助计算机的数值处理来进行，得到描述系统的某些特征，甚至进而得到描述系统的数学模型。?本实验首先给出描述系统的差分方程，通过迭代的方法求得系统的单位样值响应，进而求得该离散系统的频率响应。限于试验条件，虽然给出了系统方程，但处理的方法依然具有同样的实际意义。?具体的方法是：?在给定系统方程的条件下，选取激励信号为δ(n)，系统的起始状态为零?状态，通过迭代法，求得系统的单位样值响应h(n)（n=0,…,N）。利用公式其中 Ω 的取值范围为 0- 2 π 。计算系统的频率响应。三、 实验内容1 已知系统的差分方程为利用迭代法求得系统的单位样值响应，令N=102 利用公式其中 Ω 的取值范围为 0 2 π 。计算系统的频率响应，计算时 Ω 的步长为 0.1π。用 C 语言编制相应的计算程序。四、 实验程序源代码、流程图和实验结果1、流程图2、代码#include stdio.h#include math.h#includestdlib.h#define PI 3.1415927#define N 10double h[11]={0};double x[11]={0};double y[11]={0};double find(int n)//定义求解的迭代函数{if (n0) return 0;else return (x[n-1]+1.3*find(n-1)-0.4*find(n-2));}int main(){ int i; double re,im,fp,xp,j; FILE *fp2; fp2 = fopen(trail05.xls,w); if(fp2==NULL){printf(error!\n);exit(0);} x[0]=1; //激励初值设定 for (i=0;i=N;i++) h[i]=find(i); for (i=0;i=N;i++) { printf(h[%d]=%f\n,i,h[i]); } printf(\n); for (j=0;j=2*PI;j+=0.1*PI){ re=0; im=0; for (i=0;i=N;i++) { re=re+h[i]*cos(-i*j); im=im+h[i]*sin(-i*j); } fp=sqrt(re*re+im*im); xp=atan(im/re);if (re0 im0) xp+=PI;if (re0 im0) xp-=PI;printf(j=%f\tRe(H)=%f\tIm(H)=%f\tH=%f\tF=%f\n,j,re,im,fp,xp);fprintf(fp2,%f\t%f\t%f\t%f\n,j,fp,xp,h[i]);} fclose(fp2); getchar();}五、实验数据：5.1 单位样值响应激励X[I] 响应h[I]1h[0]=0.0000000h[1]=1.0000000h[2]=1.3000000h[3]=1.2900000h[4]=1.1570000h[5]=0.9881000h[6]=0.8217300h[7]=0.6730090h[8]=0.5462200h[9]=0.4408820h[10]=0.3546595.2幅频特性和相频特性083141596.11825-1.364510.6283182.137043-2.114030.9424771.773064-2.252311.2566360.742628-2.735591.5707950.90273-2.59921.8849540.375449-3.021732.1991130.630436-2.83312.5132720.246655-3.13562.8274310.5384-3.040373.141590.2123833.141573.4557490.53843.0403873.7699080.2466513.1355694.0840670.6304372.8331164.3982260.3754413.0217124.7123850.902732.5992145.0265440.7426142.7355755.34