北航信号与系统实验报告三资料下载.pdf

北航信号与系统实验报告三资料下载.pdf

《北航信号与系统实验报告三资料下载.pdf》由会员分享，可在线阅读，更多相关《北航信号与系统实验报告三资料下载.pdf（5页珍藏版）》请在冰豆网上搜索。

往往是选取有限项的取值，当然这里会产生一个截尾误差，但只要这个误差在相对小一个范围里，可以忽略不计。

另外，在一些实际的离散系统中，往往不是事先就能得到描述系统的差分方程的，而是通过得到系统的某些相应值，则此时系统的分析就需借助计算机的数值处理来进行，得到描述系统的某些特征，甚至进而得到描述系统的数学模型。

本实验首先给出描述系统的差分方程，通过迭代的方法求得系统的单位样值响应，进而求得该离散系统的频率响应。

限于试验条件，虽然给出了系统方程，但处理的方法依然具有同样的实际意义。

具体的方法是：

1在给定系统方程的条件下，选取激励信号为（n），系统的起始状态为零状态，通过迭代法，求得系统的单位样值响应h（n）（n=0,N）。

2利用公式北航信号与系统实验报告计算系统的频率响应。

三、三、实验内容实验内容1已知系统的差分方程为：

y（n）-1.3y（n-1）+0.4y（n-2）=x（n-1）利用迭代法求得系统的单位样值响应，令N=102利用公式，其中的取值范围是0-2,计算系统的频率响应，计算时令的步长为0.1。

四、四、实验程序源代码、流程图实验程序源代码、流程图和实验结果和实验结果#include#include#definePI3.14159#defineN10doubleh11=0;

doublex11=0;

doubley11=0;

doublefind（intn）if（n0）return0;

elsereturn（xn-1+1.3*find（n-1）-0.4*find（n-2）;

intmain（）inti;

doublere,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