Question

用matlab求解一个三元的二阶微分方程组，求大神指点。用不用ode45之类的？好像解析解解不出，数值解也行

最好帮忙画下图

参数k=0.413719,c=0.1,a=0.0015,d=0.12
初值dy/dt(t=0)=0,y(t=0)=314,x(t=0)=0.121
求dy/dt的 图像
真心求代码，在线等，虚心求教
求matlab大神解救

Answer 1

function dx=myfun(t,x)

dx=zeros(3,1);

k=0.413719;

c=0.1;

a=0.0015;

d=0.12;

dx(1)=(a*sqrt(c)*sqrt(x(1)^2-d^2))/(x(1)^(3/2))*x(3);

dx(2)=x(3);

dx(3)=k*(x(1)/c-c^2/x(1)^2)*a*sqrt(c/x(1))*sqrt(x(1)^2-d^2)/x(1);

clc

clear

options=odeset('reltol',1e-8);

x0=[0.121;314;0];

[t,y]=ode45(@myfun,[0,20],x0,options)

plot(t,y(:,2))

help 一下ode45的用法你就知道……