9. Matlab仿真二阶弹簧阻尼系统

前言

本篇文章介绍了如何使用MATLAB分析二阶系统

被控系统

$$\ddot{x}+5\dot{x}+6x=F$$
可以对上述系统进行拉式变化，得到传递函数为
$$G=\frac{1}{s^2+5s+6}$$
可以分析系统的伯德图，具体代码和图像如下：

s=tf('s');
G=1/(s^2+5*s+6);
bode(G)

输入正弦得到输出

t=0:0.05:10;
xt=1+5*cos(t);
lsim(G,xt,t)

已知输出得到输入

首先对输出进行Laplace变换而后进行反变换得到输入，但是在matlab中的输入会包含狄拉克函数项

syms s t;
y=0.1+0.5*cos(t);
Ys=laplace(y);
G=1/(s^2+5*s+6);
Xs=Ys/G;
xt=ilaplace(Xs)

输出为：
$$xt=(5\ast cos(t))/2+3\ast dirac(t)-(5\ast sin(t))/2+(3\ast dirac(1,t))/5+3/5$$
同样也可以手动计算输入的幅值、相位差等，从而手动得到输入信号的表达式。