Question

跪求下列微分方程组求解MATLAB的程序

L=l也为正整数

Answer 1

给定的微分方程组，可以用matlab的ode（）函数求解。求解方法：

1、根据方程，自定义微分方程组函数。

2、根据已经条件以及初始条件（因问题没有给出具体数值，本例自行设定），用ode45（）函数命令求得其数值解。

3、用plot函数绘制X(t)—t、Y(t)—t、φ(t)—t的曲线图。

4、按照上述方法，编程后运行可以得到如下图形。

数值解

5、代码

z0=[0.02,0.01,0.03,2,1,2];

t0=0;tf=10; 

tspan=t0:0.1:tf;

options = odeset('RelTol',1e-3,'AbsTol',1e-4);

[t,z]=ode45(@myfun,tspan,z0,options);

disp('            t        x(t)          y(t)          φ(t)')

A=[t z(:,1) z(:,2) z(:,3)];

disp(A)

figure(1)

plot(t,z(:,1),t,z(:,2),t,z(:,3),'-')

legend('x(t)','y(t)','φ(t)')