连续LTI时间系统的频域分析及S域分析.doc

实验报告 连续LTI时间系统的频域分析及S域分析 实验目的 1. 了解傅里叶变换、傅里叶逆变换的实现方式，以及傅立叶变换的时移特性、傅立叶变换的频移特性的实现方式； 2. 了解拉普拉斯变换的相关分析及实现方式，了解连续系统零极点图的绘制方式及利用零极点图判断系统的稳定性。 实验内容 1. 了解如何采用函数fourier()对时域信号进行傅里叶变换，以及如何采用函数ifourier()对频域信号进行傅里叶反变换； 2. 了解傅里叶变换的数值计算方法，以及信号频谱图； 3. 了解傅里叶变换的时移特性、频移特性的实现方法； 4. 了解连续信号拉普拉斯变换及拉普拉斯逆变换的实现形式； 5. 了解连续LTI系统的系统函数零极点图的画法，并从零极点图判断系统的稳定性。 实验验证 实验题目 1.试求下列信号的傅里叶变换的数学表达式。 2.试画出信号的频谱图，并画出信号以及信号的频谱图。 程序清单及频谱图如下图所示 3.已知某连续LTI系统的系统函数为：，试用Matlab求出系统的零极点，并绘出零极点分布图，同时判断系统的稳定性。 程序如下，结果如下图。 a=[8,2,3,1,5]; b=[1,3,2]; lxljdt(a,b) 系统零极点为：-2，-1。 4.绘出下列系统函数的零极点图： 程序如下，结果如下图。 a=[2,1]; b=[1,2,2,1]; lxljdt(a,b) 5.分别求下列信号的拉普拉斯正变换，并绘制在平面的三维曲面图。 取正变换： syms t positive laplace(4*exp(-5*t)*heaviside(t)) ans=4/(s+5) 图像绘制： x=-5.1:0.015:-4.9; y=-0.1:0.015:0.1; [x,y]=meshgrid(x,y); s=x+i*y; L=abs(4./(s+5)); mesh(x,y,L) surf(x,y,L) xlabel(x),ylabel(y),zlabel(z) title(L=4/(s+5)) axis([-5.1,-4.9,-0.1,0.1,0,500]) 取正变换： syms t positive laplace(3*cos(5*t)*heaviside(t)) ans=3*s/(s^2+25) 图像绘制： x=-5:0.6:5; y=-10:0.6:10; [x,y]=meshgrid(x,y); s=x+i*y; L=abs(3.*s./(s.^2+25)); mesh(x,y,L) surf(x,y,L) xlabel(x),ylabel(y),zlabel(z) title(L=3*s/(s^2+25)) 取正变换： syms t positive laplace(8*dirac(t-5)) ans=8*exp(-5*s) 图像绘制： x=0:0.05:1; y=-5:0.6:5; [x,y]=meshgrid(x,y); s=x+i*y; L=abs(8.*exp(-5.*s)); mesh(x,y,L) surf(x,y,L) xlabel(x),ylabel(y),zlabel(z) title(L=8*exp(-5*s)) view([5,5,5]) 取正变换： syms t positive laplace(f4=3*exp(-6*t)*cos(2*t)*heaviside(t)) ans=3/4*(s+6)/(1/4*(s+6)^2+1) 图像绘制： x=-8:0.3:-4; y=-4:0.15:4; [x,y]=meshgrid(x,y); s=x+i*y; L=abs(3./4.*(s+6)./(1./4.*(s+6).^2+1)); mesh(x,y,L) surf(x,y,L) xlabel(x),ylabel(y),zlabel(z) title(L=3/4*(s+6)/(1/4*(s+6)^2+1)) M文件说明： lxljdt .m文件： function [p,q]=lxljdt(A,B) %A:系统函数分母多项式系数向量 %B:系统函数分子多项式系数向量 %p:函数返回的系统函数极点位置行向量 %q:函数返回的系统函数零点位置行向量 p=roots(A) %求系统极点 q=roots(B) %求系统零点 p=p %极点列向量转置为行向量 q=q %零点列向量转置为行向量 x=max(abs([p,q]))%确定坐标范围 x=x+0.1 y=x clf hold on axis([-x,x,-y,y]) %确定坐标轴显示范围 axis(square) plot([-x,x],[0,0]) %画横坐标轴 plot([0,0],[-y,y])