t;titleImpulseRespo
se;H,WfreqsBum,de
;figure3;plotW,absH;xlabela
g.freq.\omegarad/s;ylabelHj\omega;titleMag
itudeRespo
e;grido
程序运行的结果如图一,图二,图三所示:
图一系统函数零极点分布图形
2
f图二系统的单位冲激响应
图三系统的幅频响应
5、MATLAB系统时域求解的基本函数求取系统单位阶跃响应:step求取系统的冲激响应:impulseystep
umde
t:其中
um和de
分别为系统传递函数描述中的分子和分母
3
f多项式系数,t为选定的仿真时间向量,一般可以由t0stepe
d等步长地产生出来。该函数返回值y为系统在仿真时刻各个输出所组成的矩阵
yxtstep
umde
:此时时间向量t由系统模型的特性自动生成状态变量x返回为空矩阵。
yxtstepABCDiu:其中ABCD为系统的状态空间描述矩阵,iu用来指明输入变量的序号。x为系统返回的状态轨迹。
impulse函数的用法求取脉冲激励响应的调用方法与step函数基本一致。yimpulse
umde
t;yxtimpulse
umde
;yxtimpulseABCDiutcovar:连续系统对白噪声的方差响应i
itial:连续系统的零输入响应lsim:连续系统对任意输入的响应对于离散系统只需在连续系统对应函数前加d就可以,如dstep,dimpulse等。
4
f三、实训过程与内容
1、采用的软件及开发平台
MATLAB70计算机
2、系统的详细设计
(一)、MATLAB中系统S平面分析的基本函数(a)Laplace拉普拉斯变换语法:laplaceFlaplaceFtfourierFwz简介:LlaplaceF表示对F函数关于默认自变量t做拉普拉斯变换。默认结果是关
于s的函数。拉普拉斯变换完成了时域函数到频域函数的转换。LlaplaceFt会令L为t的函数,而非默认的s。LlaplaceFwz规定F是关于w的函数,而L是关于z的函数。
bilaplace拉普拉斯逆变换用法:FilaplaceLFilaplaceLyFilaplaceLyx简介:FilaplaceL将Ls变换成FtFilaplaceLy将Ls变换为Fy,而非FtFilaplaceLyx将Ly变换为Fx
(二)分析与实现(1)由系统函数零、极点分布决定时域特性当Hs极点(一阶)pj中和的取值为006003006,
050250时,画出对应于ht波形,并分析波形。解:
5
f首先,根据s域的函数HS,做拉普拉斯逆变换求出时域函数ht,绘制ht的图形。
代码及波形图如下:
(1)clearclcclearclc
symsss
H1s00605jH1s006
hrealilaplaceHhrealilaplaceH
ezplothezploth
(2)clearclcsymssH1s006025jhrealilaplaceHezploth
(3)syms
(4)clearclcsymss
(5)clearclcsymss
(6)clearclcsymss
6
fH1s00305jH1s0r