1、往往是选取有限项的取值，当然这里会产生一个截尾误差，但只要这个误差在相对小一个范围里，可以忽略不计。另外，在一些实际的离散系统中，往往不是事先就能得到描述系统的差分方程的，而是通过得到系统的某些相应值，则此时系统的分析就需借助计算机的数值处理来进行，得到描述系统的某些特征，甚至进而得到描述系统的数学模型。本实验首先给出描述系统的差分方程，通过迭代的方法求得系统的单位样值响应，进而求得该离散系统的频率响应。限于试验条件，虽然给出了系统方程，但处理的方法依然具有同样的实际意义。具体的方法是：1在给定系统方程的条件下，选取激励信号为（n），系统的起始状态为零 状态，通过迭代法，求得系统的单位样值响应
2、 h（n）（n=0,N）。2利用公式 北航信号与系统实验报告 计算系统的频率响应。三、三、实验内容实验内容 1 已知系统的差分方程为：y（n）-1.3y（n-1）+0.4y（n-2）=x（n-1）利用迭代法求得系统的单位样值响应，令 N=10 2 利用公式，其中 的取值范围是 0-2,计算系统的频率响应，计算时令 的步长为 0.1。四、四、实验程序源代码、流程图实验程序源代码、流程图和实验结果和实验结果#include#include#define PI 3.14159#define N 10 double h11=0;double x11=0;double y11=0;double find
3、（int n）if（n0）return 0;else return（xn-1+1.3*find（n-1）-0.4*find（n-2）;int main（）int i;double re,im,fp,xp,w;FILE*fp2;fp2=fopen（05.xls,w）;x0=1;/激励初值 for（i=0;i=N;i+）hi=find（i）;for（i=0;i+）printf（h%d=%fn,i,hi）;printf（n）;for（w=0;w=2*PI;w+=0.1*PI）re=0;im=0;i+）re=re+hi*cos（-i*w）;im=im+hi*sin（-i*w）;fp=sqrt（re*re+im*im）;xp=atan（im/re）;if（re0）xp+=PI;if（re0&im0）xp-=PI;printf（w=%f,Re（H）=%f,Im（H）=%f,H=%f,F=%fn北航信号与系统实验报告,w,re,im,fp,xp）;fprintf（fp2,%ft%ft%fn,w,fp,xp）;fclose（fp2）;getchar（）;流程图：运行结果：北航信号与系统实验报告 幅频特性：相频特性：01234567891001234567北航信号与系统实验报告 -4-3-2-10123401234567
升级会员 