一个matlab程序有问题
来源:学生作业帮 编辑:大师作文网作业帮 分类:综合作业 时间:2024/11/11 14:22:58
一个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
一直到最后,
这是方程原型
这是一个二阶方程的求解,给你下面这个图片,希望对你有帮助.这里需要注意方程的初值.