[求助]哪位有杜芬方程和范德波方程的matlab实现程序

河边草

rt，急求，请好心人发到我油箱里：sun20022727@56.com
小生在此先谢过了！

happy

参考ode45的例子，自己写一下，很方便的

wenyue

可以采用保辛的摄动方法，参考钟万勰院士最近一两年的文章。计算精度非常高

河边草

可以采用保辛的摄动方法，参考钟万勰院士最近一两年的文章。计算精度非常高

这一方面我实在是个外行，请问能否给指点一下文章名字，先谢过了！

河边草

继续求助，哪位能把杜芬方程的M文件上传一下？？<BR>多谢了

happy

1. % Duffing方程仿真图==相轨迹图
   
2. % 方程在duffing.m文件中定义
   
3. [t,x]=ode45(@duffing,[0,2800],[0,1.5]);
   
4. x1=x(:,1);
   
5. x2=x(:,2);
   
6. x1=mod(x1,2*pi);
   
7. x1(x1>pi)=x1(x1>pi)-2*pi;
   
8. figure
   
9. h=plot(x1,x2)
   
10. 
    
11. %%%%%%%%%% 以下部分定义Duffing方程，此文件另存为duffing.m
    
12. function dx=duffing(t,x)
    
13. % 定义达芬Duffing方程：
    
14. % diff(x,2)+mu*diff(x,1)-x+x^3=F*cos(w*t)
    
15. mu=0.1;
    
16. F=2.09;
    
17. w=0.5;
    
18. dx=[x(2);F*cos(w*t)-mu*x(2)+x(1)-x(1).^3];

复制代码

jumpwolf

上面的程序中ode45是不是要改成ode15s,达分方程是刚性方程,而解刚性方程要用ode15s。不知道对不对！

wenyue

你可以看看《时间有限元与保辛》这篇文章，是钟院士2005年发表在《机械强度》上的一篇文章。里面专门讲到求解duffing方程的例子了。

happy

1. 上面的程序中ode45是不是要改成ode15s,达分方程是刚性方程,而解刚性方程要用ode15s。不知道对不对！

复制代码

就是ode45，上述方程并不是刚性的

jumpwolf

恩,是的。
但是我用ode45和ode15s得到的结果是一样的。是不是ode15s的适用面更广一点，他的精度是几阶的？

happy

恩,是的。
但是我用ode45和ode15s得到的结果是一样的。是不是ode15s的适用面更广一点，他的精度是几阶的？

一般用的都是ode45，如果ode45失败可以用ode15s参试求解
ode15s采用的是多步法是一种变阶法

河边草

谢谢楼上各位的帮忙，谢谢！