matlab 二元高阶微分方程的求解

[复制链接]
查看11 | 回复1 | 2013-2-2 12:43:09 | 显示全部楼层 |阅读模式
1、你的ode函数写的有问题,应该是function f=tuoyuan(t,x)f=[x(3);
x(4);
(-5*x(1)/sqrt((x(1)^2+x(2)^2)))/(x(1)^2+x(2)^2);
(-5*x(2)/sqrt((x(1)^2+x(2)^2)))/(x(1)^2+x(2)^2)];这样,f(1)-f(4)分别对应x'、y'、x"、y",你原来的写法把顺序都搞乱了。 2、按照改后的函数,求出来的结果也并非椭圆。求得x'和y'经过一段时间后会变为常值,而x和y则会维持固定的速度增大。你确定微分方程没错吗?怎样见得理论上是一个椭圆呢? 我知道下面的微分方程
x"=-...
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

主题

0

回帖

4882万

积分

论坛元老

Rank: 8Rank: 8

积分
48824836
热门排行