python解二阶微分方程组_python：利用多种方式解微分方程(以二阶微分系统零状态响应为例)...

1.问题：求系统的零状态响应

image.png

2.引入

首先用高数知识求解非齐次常系数微分方程

image.png

再利用信号与系统中冲激响应求解验证

image.png

利用MATLAB求解验证

y=dsolve('D2y+3*Dy+2*y=exp(-t)','y(0)=1','Dy(0)=2','t')

得出结果：

y =

(t - 2 exp(-t) + 3) exp(-t)

根据结果检验，上述手动计算与实际计算机得出结果一致。

t=0:0.1:20;

y = (t - 2 .*exp(-t) + 3) .*exp(-t);

y1=-exp(-t) .*(t - 2 .*exp(-t) + 3) + exp(-t).* (1 + 2.* exp(-t));

plot(t,y,'r-',t,y1,'b-'),legend('y','y’')

用MATLAB模拟图像结果：

image.png

3.利用Python求解该方程

通过上述计算，我们利用Python求解系统