非线性单摆的混沌分析麻烦各位大侠给调试一下代码

非线性单摆的混沌分析clearglobal a w wp p; %定义全局变量,a = 0.25;w =1;wp =0.6667; %设定阻尼项系数,系统频率,周期性外力的频率p = input(' 输入周期性外力矩的振幅大小 ');tspan = [0 600]; %设定积分时间y0 = [0 0.2 0]; %初时条件t = 0,角度 = 0rad,角速度 = 0.2rad/s% [t,y] = ode23(' xxdb ',tspan,y0); %求解名为“xxdb”的线性微分方程[t,y] = ode23(' fcdo ',tspan,y0); %求解名为“fcdo”的非线性微分方程subplot(2,1,1)plot(t,y(:,1),' b '); %描绘出非线性单摆的振动曲线hold ony0 = [0.01 0.02] ; %对初时条件稍作变动,t = 0,角度 = 0.01rad,角速度 = 0.2rad/s%[t,y] = ode23('xxdb',tspan,y0); %重解线性微分方程[t,y] =ode23('fcdo',tspan,y0); %重解非线性微分方程plot(t,y(:,1),' k '); %在同一图上描绘出初时条件稍作变动后的振动曲线xlabel(' 时间 ');ylabel(' 角位移 ');%非线性单摆的微分方程function yp = fcdo(t,y) global a w wp p; %定义全局变量yp = [y(2)- 2*a*y(2)- w*sin(y(1))+ p*cos(y(3))*wp]; %写入非线性单摆的微分方程%线性单摆的微分方程%function yp = xxdb(1,y) global a w wp p; %定义全局变量%yp = [y(2)- 2*a*y(2)- w*(y(1))+ p*cos(y(3))*wp]; %写入线性单摆的微分方程