物理学毕业论文 - 线性动、非线性受迫振动及混沌效应.doc

线性动、非线性受迫振动及混沌效应 在这里我研究线性与非线性的单摆的运动, 对于存在阻尼及外力情况下, 将线性系统与非线性系统的振动情况作一比较, 从而了解非线性系统在某些临界状态时混沌现象. 1 存在阻尼时的线性线性受迫振动 对于线性系统, 根据牛顿第二定律得系统振动方程 (1) 其中为外驱动力频率,、、为系统参数. 1.1 不存在阻尼及外驱动力 假定, 此时系统作简谐振动.假定初始位移和初始速度分别为, .位移、速度与时间的关系如图1、图2.function Q=funxz(t,x,w0,b,f0,w) Q=[x(2);-w0^2*x(1)-2*b*x(2)+f0*cos(w*t)]; 图1计算程序:clear w0=1*pi; b=0.0; f0=0.0; w=0*pi; [t,x]=ode45(@funxz,[0:0.001:10],[1.0 0.0], [], w0,b,f0,w); plot(t,x(:,1),-b,LineWidth,2) xlabel(\fontsize{18}\fontname{Time New Roman}t, Color,k) ylabel(\fontsize{18}\fontname{Time New Roman}x, Color,k) grid on set(gca,xgrid,on,GridLineStyle,-) set(gca,ygrid,on,GridLineStyle,-) title(Linear vibration displacement versus time,fontsize,18,color,r) %设置标题 图计算程序:clear w0=2*pi; b=0.0; f0=0.0; w=0*pi; [t,x]=ode45(@funxz,[0:0.001:10],[1.0 0.0], [], w0,b,f0,w); plot(t,x(:,1),-b,LineWidth,2) xlabel(\fontsize{18}\fontname{Time New Roman}t, Color,k) ylabel(\fontsize{18}\fontname{Time New Roman}x, Color,k) grid on set(gca,xgrid,on,GridLineStyle,-) set(gca,ygrid,on,GridLineStyle,-) title(Linear vibration displacement versus time,fontsize,18,color,r) %设置标题 1.2 存在阻尼但不存在外驱动力 假定, 此时系统.假定初始位移和初始速度分别为, .位移、速度与时间的关系如图、图. 注意此时振动在时振动基本停止. 图计算程序:clear w0=1*pi; b=0.2; f0=0.0; w=0*pi; [t,x]=ode45(@funxz,[0:0.001:50],[1.0 0.0], [], w0,b,f0,w); plot(t,x(:,1),-b,LineWidth,2) xlabel(\fontsize{18}\fontname{Time New Roman}t, Color,k) ylabel(\fontsize{18}\fontname{Time New Roman}x, Color,k) grid on set(gca,xgrid,on,GridLineStyle,-) set(gca,ygrid,on,GridLineStyle,-) title(Linear vibration displacement versus time,fontsize,18,color,r) %设置标题 图计算程序:clear w0=1*pi; b=0.2; f0=0.0; w=0*pi; [t,x]=ode45(@funxz,[0:0.001:50],[1.0 0.0], [], w0,b,f0,w); plot(t,x(:,2),-b,LineWidth,2) xlabel(\fontsize{18}\fontname{Time New Roman}t, Color,k) ylabel(\fontsize{18}\fontname{Time New Roman}v, Color,k) grid on set(gca,xgrid,on,GridLineStyle,-) set(gca,ygrid,on,GridLineStyle,-) title(Linear vibration speed versus time,fontsize,18,color,r) %设置标题 1.3 存在阻尼外驱动力