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

mojojo

方程：

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

sogooda

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

sogooda

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

1. clear;clc
   
2. B0=0;
   
3. yip=1;
   
4. global B0 yip
   
5. options = odeset('RelTol',1e-4,'AbsTol',[1e-6 1e-6]);
   
6. [T,y]=ode15s('fxy',[0 1],[0 0],options);
   
7. plot(T,y)

复制代码

1. function x1=fxy(t,x);
   
2. global B0 yip
   
3. x1(1)=(x(2)+x(1)*(1-x(1)^2))/yip;
   
4. x1(2)=-x(1)+B0;
   
5. x1=x1';

复制代码

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

mojojo

原帖由 sogooda 于 2008-5-31 15:58 发表

                               
登录/注册后可看大图

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

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

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

ch_j1985

原帖由 mojojo 于 2008-5-31 16:13 发表

                               
登录/注册后可看大图

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

lz有没有方程的一些初始值？有的话，贴出来，OK？

mojojo

原帖由 ch_j1985 于 2008-5-31 16:33 发表

                               
登录/注册后可看大图

lz有没有方程的一些初始值？有的话，贴出来，OK？

第一个方程eps=0.1，B0=0.5694。。。

sogooda

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

mojojo

原帖由 sogooda 于 2008-5-31 16:50 发表

                               
登录/注册后可看大图

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

多谢指教，谢谢！