软件三离散.docxVIP

离散系统频率特性分析实验目的通过本实验学习和利用Matlab软件进行离散系统的频率特性的分析二、实验原理1．离散系统的频率响应对于某因果稳定离散系统，如果激励序列为正弦序列：则，根据《信号与系统》课程给出的结果有，系统的稳态响应为：定义离散系统的频率响应为其中，——称为离散系统的幅频特性；——称为离散系统的相频特性；是以为周期的周期函数，只要分析在范围内的情况，便可分析出系统的整个频率特性。2．用MATLAB实现离散系统的频率特性分析方法（1）直接法设某因果稳定系统的系统函数，则系统的频响特性为：MATLAB提供了专门用于求离散系统频响特性的函数freqz()，调用freqz()的格式有以下两种：[H,w]=freqz(B,A,N) B和A分别为离散系统的系统函数分子、分母多项式的系数向量，N为正整数，返回量H则包含了离散系统频响在范围内N个频率等分点的值，向量w则包含范围内N个频率等分点。调用中若N默认，默认值为512。[H,w]=freqz(B,A,N,’whole’)该调用格式将计算离散系统在范围内N个频率等分点的频率响应的值。因此，可以先调用freqz()函数计算系统的频率响应，然后利用abs()和angle()函数及plot()函数，即可绘制出系统在或范围内的频响曲线。三、实验内容绘制如下系统的频响曲线MATLAB命令如下：B=[1 -0.5];A =[1 0];[H,w]=freqz(B,A,400,whole);Hf=abs(H);Hx=angle(H);clffigure(1)plot(w,Hf)title(离散系统幅频特性曲线)figure(2)plot(w,Hx) title(离散系统相频特性曲线)运行结果如图8-3所示。（2）几何矢量法利用几何矢量求解示意图如图8-4所示。有：则系统的幅频特性和相频特性分别为：（8-7）（8-8）根据式（8-7）和（8-8），利用MATLAB来求解频率响应的过程如下：根据系统函数定义分子、分母多项式系数向量和；调用前述的ljdt()函数求出的零极点，并绘出零极点图；定义Z平面单位圆上的个频率分点；求出所有的零点和极点到这些等分点的距离；求出所有的零点和极点到这些等分点矢量的相角；根据式（8-7）和（8-8）求出系统的和；绘制指定范围内系统的幅频曲线和相频曲线；下面是实现上述过程的实用函数dplxy()。有四个参数：k为用户定义的频率等分点数目；B和A分别为系统函数分子、分母多项式系数向量；r为程序绘制的频率特性曲线的频率范围（）。function dplxy(k,r,A,B)%The function to draw the frequency response of discrete systemp=roots(A); %求极点q=roots(B);%求零点figure(1)ljdt(A,B)%画零极点图w=0:r*pi/k:r*pi;y=exp(i*w);%定义单位圆上的k个频率等分点N=length(p);%求极点个数M=length(q);%求零点个数yp=ones(N,1)*y;%定义行数为极点个数的单位圆向量yq=ones(M,1)*y;%定义行数为零点个数的单位圆向量vp=yp-p*ones(1,k+1);%定义极点到单位圆上各点的向量vq=yq-q*ones(1,k+1);%定义零点到单位圆上各点的向量Ai=abs(vp);%求出极点到单位圆上各点的向量的模Bj=abs(vq);%求出零点到单位圆上各点的向量的模Ci=angle(vp);%求出极点到单位圆上各点的向量的相角Dj=angle(vq);%求出零点到单位圆上各点的向量的相角fai=sum(Dj,1)-sum(Ci,1);%求系统相频响应H=prod(Bj,1)./prod(Ai,1);%求系统幅频响应figure(2)plot(w,H);%绘制幅频特性曲线title(离散系统幅频特性曲线)xlabel(角频率)ylabel(幅度)figure(3)plot(w,fai)title(离散系统的相频特性曲线)xlabel(角频率)ylabel(相位)四、实验总结本实验进行了使用Matlab软件分析离散系统频谱特性的编程，并且通过两种不同的方法进行对比，思路虽不同但结果一致。矢量法接触较少但仍是一个比较不错的方法。