可以归一化啊：
设：f=y1,f'=y2,f''=y3
设：r=y4,r'=y5
dy1/dt=y2
dy2/dt=y3
dy3/dt=-3*y1*y3+2*y2^2-y4
dy4/dt=y5
dy5/dt=-3*p*y1*y5
就是dy=a*y
a是系数矩阵,尺度5*5
y=[y1;y2;y3;y4;y5]，尺度5*1
dy=[y2;y3;-3*y1*y3+2*y2^2-y4;y5;-3*p*y1*y5],5*1
可以写函数文件了
function
dy=fun(t,y)
dy=[y(2);y(3);-3*y(1)*y(3)+2*y(2)^2-y(4);y(5);-3*p*y(1)*y(5)];
脚本里面可用ode45调用
t=...;
y0=[...;;;];
[t,y]=ode45('fun',t,y0)
其他的自己补充，呵呵
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024