离散时间信号通过线性时不变系统实验报告教材Word格式文档下载.docx

1、其中f为数字频率，w为数字角频率。产生两个信号: N取大于500.设两个离散时间线性时不变系统分别为：系统1：系统2：其中： 2、实验内容 （1）令分别通过1、2（注意：通过系统1可以用线性卷积求解，通过系统2必须用差分方程求解），注意系统1、2输出信号的区别和共性。 #include#includemath.h#define pi 3.1415926void conv（float u,float v,float w, int m, int n）;/void plot（int xlenth,float *y,char *name,int gain,char *xlabel,char *ylab

2、el）;void main（） int i,j; int m,n; int gain=1; int y11,y22; char name=f; char xlabel=X,ylabel=Y float xa650=0.0,xb650=0.0,h650=0.0,y1650=0.0,y2650=0.0; for（i=0;i600;i+） xai=sin（2*pi/55.0*i）; xbi=sin（2*pi/55.0*i）+sin（2*pi*20.0/55.0*i）; for（j=0;j=49;j+） hj=1.0; m=sizeof（xa）/sizeof（xa0）; n=sizeof（xb）/si

3、zeof（xb0）; /a=sizeof（h）/sizeof（h0）; y11=sizeof（y1）/sizeof（y10）; y22=sizeof（y2）/sizeof（y20）; conv（xa,h,y1,m,n）; conv（xb,h,y2,m,n）; /plot（y11,y1,s,gain,xlabel,ylabel）; initgraph（1000,500）; initgraph（960,480,SHOWCONSOLE）;/初始化图像 setbkcolor（WHITE）;/ 用背景色清空屏幕 cleardevice（）; setcolor（BLACK）; outtextxy（480,

4、 10, name）;/图像名 outtextxy（900, 220, xlabel）;/ outtextxy（5, 10, ylabel）; line（0,240,960,240）;/画横坐标轴 line（30,0,30,480）;/画纵坐标轴 setlinestyle（PS_DASHDOT）; setlinestyle（PS_SOLID）; setlinecolor（BLUE）;650; line（100+i,50,100+i,xai*25+50）; line（100+i,150,100+i,xbi*25+150）; line（100+i,280,100+i,y1i*5+280）; lin

5、e（100+i,400,100+i,y2i*5+400）; getchar（）; closegraph; /卷积子程序void conv（float u,float v,float w, int m, int n） int k=m+n-1;k; for（j=max（0,i+1-n）;=min（i,m-1）; j+） wi+=uj*vi-j; /画图函数参数xlength:序列长度 y:序列幅度 s:图像名 gain:y坐标轴增益，用来显示y坐标 xlabel:x轴坐标标识符 ylabel：y轴坐标标识符void plot（int xlenth,float *y,char *name,int g

6、ain,char *xlabel,char *ylabel） int i=0; int j=0; char ss6;xlenth;i+） if（!（i%50） sprintf（ss, %d, i）; outtextxy（i+30, 460, ss）; line（i+30,0,i+30,480）; 480;（i%40） line（0,i,960,i）; setorigin（30, 240）; setaspectratio（1, -1）; /并未将 outtextxy进行反转,outtextxy用的原来的坐标 line（i,yi*gain,i+1,yi+1*gain）;运行结果：结论：单频信号通过

7、系统1和系统2输出结果的频率相同。系统存在过冲，系统2效果更好。 （2）令分别通过1、2，注意系统1、2输出信号的区别和共性与（1）有何不同之处。 int i,j,m;float xa650=0.0,xb650=0.0, y1650=0.0,y2650=0.0,a650=0.0,b650=0.0;for（i=0;xai=sin（2*pi/55.0*i）;xbi=sin（2*pi/55.0*i）+sin（2*pi*20.0/55.0*i）;b0=0.0223,b1=0.01,b2=0.0223,a1=-1.7007;a2=0.7613;y10=b0*xa0;for（m=1;m0）&3; y1m=

8、y1m+bj*xam-i; for（i=1;（m-i） y1m=y1m-ai*y1m-i;y20=b0*xb0; y2m=y2m+bj*xbm-i; y2m=y2m-ai*y2m-i;initgraph（1000,800）; 运行结果: 结论：信号有两个频率分量f1和f2，f1与频率相同，过系统1和系统2后，输出结果与过系统1和系统2类似，但过系统后有毛刺，曲线不光滑，且过系统2曲线比过系统1曲线光滑。说明，这两个系统相当于一个低通滤波器。1.4 实验总结 通过本次实验对离散时间信号通过线性时不变系统的原理有了更深的了解。掌握了使用C语言求解线性卷积的方法。对于单频率信号和双频信号通过相同系统后输出的波形的不同也有了更好的认识。