北航信号与系统实验2.docx

信号与系统实验名称：连续时间系统的频率响应一、实验目的：1 加深对连续时间系统频率响应理解； 2 掌握借助计算机计算任意连续时间系统频率响应的方法。二、实验原理： 连续时间系统的频率响应可以直接通过所得表达式计算，也可以通过零极点图通过用几何的方法来计算，而且通过零极点图可以迅速地判断系统的滤波特性。 根据系统函数H(s)在s平面的零、极点分布可以绘制频响特性曲线，包括幅频特性 H(jw) 曲线和相频特性H(w)曲线。这种方法的原理如下：假定，系统函数H(s)的表达式为当收敛域含虚轴时，取s = jw，也即在s平面中，s沿虚轴从- j∞移动到+ j∞时，得到容易看出，频率特性取决于零、极点的分布，即取决于Zj 、Pi 的位置，而式中K 是系数，对于频率特性的研究无关紧要。分母中任一因子(jw- Pi )相当于由极点 p 引向虚轴上某点 jw的一个矢量；分子中任一因子(jw-Zj)相当于由零点Zj引至虚轴上某点 jw的一个矢量。在右图示意画出由零点Zj和极点 Pi 与 jw点连接构成的两个矢量，图中Nj、Mi 分别表示矢量的模，ψj、θi 表示矢量的辐角（矢量与正实轴的夹角，逆时针为正）。 对于任意零点Zj 、极点Pi ，相应的复数因子（矢量）都可表示为：于是，系统函数可以改写为式中当ω延虚轴移动时，各复数因子（矢量）的模和辐角都随之改变，于是得出幅频特性曲线和相频特性曲线。这种方法称为s平面几何分析。通过零极点图进行计算的方法是： 1 在S平面上标出系统的零、极点位置； 2 选择S平面的坐标原点为起始点，沿虚轴向上移动，计算此时各极点和零点与该点的膜模和夹角； 3 将所有零点的模相乘，再除以各极点的模，得到对应频率处的幅频特性的值； 4 将所有零点的幅角相加，减去各极点的幅角，得到对应频率处的相角。三、实验内容用 C 语言编制相应的计算程序进行计算，要求程序具有零极点输入模块，可以手工输入不同数目的零极点。1.假设系统1的零极点为：零点Z1=0；极点 p1=-1-j,p2=-1+j.2.假设系统2的零极点为：零点Z1=2+1.5j,Z2=2-1.5j;极点p1=-0.5-0.7j,p2=-0.5+0.7j;p3=-2.5-2j,-2.5+2j;计算频率从0～5频段的频谱，计算步长为0.1，分别计算上面两个系统的幅频特性和相频特性，将所得结果用表格列出，并画出相应的幅频特性曲线和相频特性曲线。 判断所给系统的滤波特性，对于带通滤波器，计算出 3dB 带宽的起始频点和结束频点；对于低通或高通滤波器，计算出3dB带宽的截止或开始的频率。四、画出系统一和系统二的零极点图 系统一 系统二五、程序流程图和程序代码程序流程图如下：程序代码如下：#include stdio.h#include math.hint m,n,u=0;float a[100][2],b[100][2];void laf(){int i,o;printf(the number of zero point:\n);scanf(%d,m);printf(Laf the zero point.(x,y)\n);for(i=1;i=m;i++){printf(%d:,i);scanf(%f,%f,a[i][0],a[i][1]);}printf(the number of pole point:\n);scanf(%d:,n);printf(input the pole points.(x,y)\n);for(o=1;o=n;o++){printf(%d:,o);scanf(%f,%f,b[o][0],b[o][1]);}}float r(float x,float y,float z){return sqrt(x*x+(y-z)*(y-z));}float p(float x,float y,float z) {const float pi=3.1415927; float t; if(x==0){t=pi/2;}else t=atan((y-z)/x);if(x0) t=pi+t;return t;}void main(){float f,Ts=0.1,h,w,mo,ao;int e;laf();for(f=0;f=5;f=f+Ts){h=1;w=0;mo=1;ao=0;for(e=1;e=m;e++) {h=h*r(a[e][0],a[e][1],f); w+=p(a[e][0],a[e][1],f); }for(e=1;e=n;e++){mo=mo*r(b[e][0],b[e][1],f);ao+=p(b[e][0],b[e][1],f);}h=h/mo;w=w-ao;printf(