河边草 发表于 2006-3-12 19:23

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

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

happy 发表于 2006-3-12 19:26

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

wenyue 发表于 2006-3-13 07:23

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

河边草 发表于 2006-3-13 10:40

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

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

河边草 发表于 2006-3-15 10:11

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

happy 发表于 2006-3-16 10:25

% Duffing方程仿真图==相轨迹图
% 方程在duffing.m文件中定义
=ode45(@duffing,,);
x1=x(:,1);
x2=x(:,2);
x1=mod(x1,2*pi);
x1(x1>pi)=x1(x1>pi)-2*pi;
figure
h=plot(x1,x2)

%%%%%%%%%% 以下部分定义Duffing方程，此文件另存为duffing.m
function dx=duffing(t,x)
% 定义达芬Duffing方程：
% diff(x,2)+mu*diff(x,1)-x+x^3=F*cos(w*t)
mu=0.1;
F=2.09;
w=0.5;
dx=;

jumpwolf 发表于 2006-3-16 16:30

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

wenyue 发表于 2006-3-16 21:04

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

happy 发表于 2006-3-17 07:25

回复：（jumpwolf）ode45还是ode15s?

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

jumpwolf 发表于 2006-3-17 15:18

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

happy 发表于 2006-3-17 15:42

回复：（jumpwolf）恩,是的。但是我用ode45和ode15s得...

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

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

河边草 发表于 2006-3-19 12:00

谢谢楼上各位的帮忙，谢谢！
页: [1]
查看完整版本: [求助]哪位有杜芬方程和范德波方程的matlab实现程序