这种有初始条件的情况不太适合使用拉氏变换,而可以考虑直接使用积分器搭建。
把方程表示成
y′′′′= -(3y′′′+4y′′+5y) + e﹣3t+e﹣5tsin(4t+π/3)
实现如图所示(模型另外作为附件上传)。
注意:
1、你的提问中,初始条件写的有问题:有两个y(0),但没有y'''(0),这里暂且把第一个y(0)当成y'''(0)来用,你根据实际情况再改。
2、本系统是发散的。方程右侧为输入,对左侧进行拉氏变换可得特征多项式,其系数缺项(y'的系数为0),可以由系统稳定的必要条件直接判定系统不稳定。这是由方程左端的系数决定的,与初始条件以及右侧的输入无关。仿真结果也是如此,如下: