Matlab课程设计连续时间系统的复频域分析与仿真设计.doc

《Matlab课程设计连续时间系统的复频域分析与仿真设计.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、下载可编辑郑州航空工业管理学院《电子信息系统仿真》课程设计2013级电子信息工程专业131308143班级题目连续时间系统的复频域分析与仿真姓名学号131308143指导教师二О一五年十二月十日.专业.整理.下载可编辑连续时间系统的复频域分析与仿真一．实验目的1.掌握研究连续时间信号和系统频域分析的理论知识进行。2.绘出典型单边信号的时域波形。3.绘出拉普拉斯变换的曲面图及连续时间系统极零点图。4.能够分析系统的稳定性。二．实验原理1.连续时间系统的复频域描述系统函数的实质就是系统单位冲激响应的拉普拉斯变换。

2、因此，系统函数也可以定义为：所以，系统函数的一些特点是和系统的时域响应的特点相对应的。假设描述一个连续时间系统的线性常系数微分方程为：1对式1两边做拉普拉斯变换，则有即2.专业.整理.下载可编辑式2告诉我们，对于一个能够用线性常系数微分方程描述的连续时间系统，它的系统函数是一个关于复变量s的有理多项式的分式，其分子和分母的多项式系数与系统微分方程左右两端的系数是对应的。根据这一特点，可以很容易的根据微分方程写出系统函数表达式，或者根据系统函数表达式写出系统的微分方程。在MATLAB中，表达系统函数的方法是给出

3、系统函数的分子多项式和分母多项式的系数向量。由于系统函数的分子和分母的多项式系数与系统微分方程左右两端的系数是对应的，因此，用MATLAB表示系统函数，就是用系统函数的两个系统向量表示。2.系统函数的零极点分布图系统函数的零极点图能够直观地表示系统的零点和极点在s平面上的位置，从而比较容易分析系统函数的收敛域和稳定性。3.系统函数的极点分布与系统的稳定性和因果性之间的关系3同时，我们还应该记得，一个信号的傅里叶变换的存在条件就是这个信号满足绝对可积条件，所以，如果系统是稳定的话，那么，该系统的频率响应也必然是

4、存在的。又根据傅里叶变换与拉普拉斯变换之间的关系，可进一步推理出，稳定的系统，其系统函数的收敛域必然包括虚轴。稳定的因果系统，其系统函数的全部极点一定位于s平面的左半平面。所以，对于一个给定的LTI系统，它的稳定性、因果性完全能够从它的零极点分布图上直观地看出。三．实验步骤与实验程序.专业.整理.下载可编辑1.将绘制零极点图的扩展splane为文函数文件splane以件名存盘。%   splane%   Thisfunctionisusedtodrawthezero-poleplotinthes-planef

5、unctionsplane(num,den)p=roots(den);                       %Determinethepolesq=roots(num);                       %Determinethezerosp=p';q=q';x=max(abs([pq]));  %Determinetherangeofreal-axis  x=x+1;y=x;              %Determinetherangeofimaginary-axis       pl

6、ot([-xx],[00],':');holdon;          %Drawthereal-axisplot([00],[-yy],':');holdon;      %Drawtheimaginary-axisplot(real(p),imag(p),'x');holdon;        %Drawthepolesplot(real(q),imag(q),'o');holdon;        %Drawthezerostitle('zero-poleplot');xlabel('RealPart'

7、);ylabel('ImaginalPart')axis([-xx-yy]);               %Determinethedisplay-range2.因果系统函数，绘制出系统的零极点图、系统的单位冲激响应、系统的幅度频率响应和相位频率相应的图形。%Program3_1%ThisProgramisusedtocomputeanddrawtheplotsofthefrequencyresponse%ofacontinuous-timesystemb=[1];%Thecoefficientvector

8、oftherightsideofthedifferentialequationa=[132];%Thecoefficientvectoroftheleftsideofthedifferentialequation[H,w]=freqs(b,a);%Computethefrequencyresponse.专业.整理.下载可编辑HHm=abs(H);%Computethemagnituderespons