灵鹊做题网一个matlab程序有问题
来源：学生作业帮 编辑：灵鹊做题网作业帮 分类：综合作业 时间：2024/04/29 03:00:37
一个matlab程序有问题
解一个微分方程组,大家看看我的程序有何问题?
m文件：function dv=eq(t,y)
dy=[y(2);(-3.986e+014*y(1))/((y(1))^2+(y(3))^2)^1.5;y(4);(-3.986e+014*y(3))/((y(1))^2+(y(3))^2)^1.5];
下面是程序：
[t,y]=ode45('eq',[0 2],[0 1e4 6.371e6 0]);
运行结果：
Warning:Mixture of single and double data for 't0','y0',and 'f(t0,y0)' in call to ode45.
> In funfun\private\odearguments at 139
In ode45 at 173
然后就是一堆数据,但结果都是这样的,
0 0.0100 6.3710 0
0.0000 0.0100 6.3710 0.0000
0.0000 0.0100 6.3710 0.0000
一直到最后,
这是方程原型
解一个微分方程组,大家看看我的程序有何问题?
m文件：function dv=eq(t,y)
dy=[y(2);(-3.986e+014*y(1))/((y(1))^2+(y(3))^2)^1.5;y(4);(-3.986e+014*y(3))/((y(1))^2+(y(3))^2)^1.5];
下面是程序：
[t,y]=ode45('eq',[0 2],[0 1e4 6.371e6 0]);
运行结果：
Warning:Mixture of single and double data for 't0','y0',and 'f(t0,y0)' in call to ode45.
> In funfun\private\odearguments at 139
In ode45 at 173
然后就是一堆数据,但结果都是这样的,
0 0.0100 6.3710 0
0.0000 0.0100 6.3710 0.0000
0.0000 0.0100 6.3710 0.0000
一直到最后,
这是方程原型
这是一个二阶方程的求解,给你下面这个图片,希望对你有帮助.这里需要注意方程的初值.