mojojo 发表于 2008-5-31 14:57

这样的vdp方程能用matlab的工具箱画图吗？

方程：

需要画此方程图，但刚接触matlab，麻烦高手指教一二……

sogooda 发表于 2008-5-31 15:27

版面搜索一下“微分方程”，用matlab求数值解应该没问题。

sogooda 发表于 2008-5-31 15:58

学解微分方程。可是为啥结果看起来好像不太对？是不是因为方程刚性？可是我用的ode15s啊。
不解。

clear;clc
B0=0;
yip=1;
global B0 yip
options = odeset('RelTol',1e-4,'AbsTol',);
=ode15s('fxy',,,options);
plot(T,y)
function x1=fxy(t,x);
global B0 yip
x1(1)=(x(2)+x(1)*(1-x(1)^2))/yip;
x1(2)=-x(1)+B0;
x1=x1';

[ 本帖最后由 sogooda 于 2008-5-31 16:01 编辑 ]

mojojo 发表于 2008-5-31 16:13

原帖由 sogooda 于 2008-5-31 15:58 发表 http://www.chinavib.com/forum/images/common/back.gif
学解微分方程。可是为啥结果看起来好像不太对？是不是因为方程刚性？可是我用的ode15s啊。
不解。

clear;clc
B0=0;
yip=1;
global B0 yip
options = odeset('RelTol',1e-4,'AbsTol',);
=od ...

图好像大致画出是这个样子的。。。

ch_j1985 发表于 2008-5-31 16:33

原帖由 mojojo 于 2008-5-31 16:13 发表 http://www.chinavib.com/forum/images/common/back.gif
图好像大致画出是这个样子的。。。
35264
lz有没有方程的一些初始值？有的话，贴出来，OK？

mojojo 发表于 2008-5-31 16:36

原帖由 ch_j1985 于 2008-5-31 16:33 发表 http://www.chinavib.com/forum/images/common/back.gif

lz有没有方程的一些初始值？有的话，贴出来，OK？
第一个方程eps=0.1，B0=0.5694。。。

sogooda 发表于 2008-5-31 16:50

方程的结果和参数的选择、初始条件有密切的关系，所以楼主还是根据上面的程序自己改改吧

mojojo 发表于 2008-5-31 16:54

原帖由 sogooda 于 2008-5-31 16:50 发表 http://www.chinavib.com/forum/images/common/back.gif
方程的结果和参数的选择、初始条件有密切的关系，所以楼主还是根据上面的程序自己改改吧
多谢指教，谢谢！
页: [1]
查看完整版本: 这样的vdp方程能用matlab的工具箱画图吗？