非线性动力学的两个问题

lingfei.kong

请教各位大侠：
本人是初学非线性动力学，有几个问题想请教大家：
1.  不知道响应周期是多少，怎么求非自治系统周期解和响应周期。（书上说用打靶法）
是否利用龙格－库塔法，将激励频率带入迭代N次，然后观察庞加莱相图消除响应周期部分，从而得到周期解？响应周期怎么得到？例如下程序：
function df=dafen(t,x,flag,force)
df=[x(2);force*cos(1.2*t)-x(1)^3+x(1)-0.3*x(2)];
clear
ff=0.222;
options=odeset('RelTol',1e-7);
tt=3*pi/1.2;％带入激励频率
[t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],ff);
plot(x(7000:end,1),x(7000:end,2),'.')
%pojialai
hold on
i=2000:100:3000
plot(x(i,1),x(i,2),'*')

2．上面程序中ff和[]是什么意思？我查了ODE45，options是相关参数例如误差精度

[ 本帖最后由 eight 于 2007-8-15 21:02 编辑 ]