Question

matlab第二题

Answer 1

用ode45求解常微分方程组的初值问题，应按下列步骤来求解。

1、首先，建立自定义函数，f= func(t,x)

2、其二，用ode45（）函数求解t和x值。说明：x为向量，即x=[x(1) x(2)]。求解格式

[t,x] = ode45(@func,[0 10],y0);

用dsolve（）函数可以得到常微分方程的解析值。求解格式

[x1,y1]=dsolve(Dx+5*x+y==exp(t),Dy-x-3*y==0,x(0)==1,y(0)==0)   

3、将得到t在【0，10】时的x、x1、y1值，用plot（）函数绘制出解析值和数值解对比图形。

具体实现代码和运行结果如下