这个微分方程用matlab怎么解

weiqaz123 · 发表于 2006-10-9 18:39

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

dv/dt=5sinx-0.1*v*v
dx/dt=5cosx/v
初始条件为t=0时，v=150,x=0.8
这个微分方程用matlab怎么解
怎么样得出v和x随着t的变化曲线

[ 本帖最后由 lxq 于 2006-10-9 19:57 编辑 ]

fanghuikeer · 发表于 2006-10-9 21:10

function xdot=ab(t,x)
xdot=[5*sin(x(2))-0.1*x(1)^2;5*cos(x(2))/x(1)]；
[t,x]=ode45(@ab,[0 100],[150 0.8])
plot(t,x)
t的取值范围我设为０到１００，自己可以根据需要而去改变，作出的图形如下：

fenghai222 · 发表于 2006-10-9 22:50

请问怎么分出哪个是v，哪个是x的曲线？
我刚接触这个，谢了！

fanghuikeer · 发表于 2006-10-9 22:56

x(1)代表v，x(2)代表x

fanghuikeer · 发表于 2006-10-9 22:57

用龙格＿库塔法解这个常微分方程，很简单的，多看看数值计算方面的书！！

weiqaz123 · 发表于 2006-10-10 20:03

Input argument 'x' is undefined.

Error in ==> C:\MATLAB6p5\work\ab.m
On line 2 ==> xdot=[5*sin(x(2))-0.1*x(1)^2;5*cos(x(2))/x(1)];

我做的时候出现这种问题是怎么回事啊

sally_mn · 发表于 2006-10-11 11:23

我也是，出现同样错误

jimin · 发表于 2006-10-11 11:32

原帖由 weiqaz123 于 2006-10-10 20:03 发表
Input argument 'x' is undefined.

Error in ==> C:\MATLAB6p5\work\ab.m
On line 2 ==> xdot=;

我做的时候出现这种问题是怎么回事啊

不知道是不是版本问题，我在7.0下fanghuikeer 的程序没问题，图也是和他一样的
只是那个分号要改成英文下的

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 这个微分方程用matlab怎么解

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

评分