求matlab计算duffing 模型的原文件或代码
http://heaps.cn/cgi-bin/mimetex.cgi?formdata=duffing+model%3A+%5Cfrac%7Bd%5E2+x%7D%7Bd+t%5E2%7D++%2B+g+%5Cfrac%7Bdx%7D%7Bd+t%7D+%2Bx%5E3+%3D+b%5Ccdot+cos+t+求解好像用数值分析 ode 45
[ 本帖最后由 eight 于 2007-6-27 00:30 编辑 ] %达芬振子 m*x''+c*x'+k*x+x^3=f*cos(w*t)
function xdot=duffing(t,x,options,f)
%parameter input
m=1;
c=0.3;
k=1;
w=1.2;
%state functionxdot
xdot=;
%达芬振子 m*x''+c*x'+k*x+x^3=f*cos(w*t) 相图
options=odeset;
f=0;%0.2 0.280.28780.2980.31
tspan=;
x0=;
=ode45('duffing',tspan,x0,options,f);
plot(t,x(:,1));
figure
plot(x(:,1),x(:,2));
% plot(x(end-500:end,1),x(end-500:end,2));
页:
[1]