matlab画相轨迹,非线性振动系统相轨迹对应波形图的绘制方法

非线性振动系统的研究方法目前主要有:定性分析法、定量分析法、数值分析法.定性分析方法是在相平面上研究解或平衡点的性质和相图性质从而定性地确定解的性态;是利用相平面内的相轨迹直观地显示出非线性系统运动的主要性质和特征.但是,在很多非线性振动问题中,希望通过相轨迹能得到其波形图,即通过相轨迹就能定性知道系统的运动情况[1-4].到目前为此,这方面的文献很少.尽管在非线性电路问题研究中有人提出过根据相轨迹方程绘制对应波形图的理论[5],但对一般的相轨迹方程,真正操作时难度太大,很难画出对应的波形图.在非线性振动问题中,更少见到关于该问题的研究报道.本文对这一问题进行了研究,得出了一种通过非线性振动系统相轨迹方程,利用matlab软件绘制其对应波形图的方法.利用该方法能绘制任何形状的相轨迹对应的波形图.1相轨迹上各点时间的标定根据非线性电路波形图绘制理论,可以得出非线性振动系统相轨迹上各点时间的标定方法.设非线性振动系统的状态方程为x觶=y,y觶=(fx,y)!.(1)如果要给出x=X(t)和y=Y(t)的波形图,必须对相轨迹上各点标定出时间.在如图1所示的相平面上假设有一条通过点P(0x0,y0)的相轨迹.设沿相轨迹从P0到达P1,P2,P3,……,Pn等点所需要的时间为t1,t2,t3,……,tn.假定这些点都相距很近,因相轨迹是x、y的函数,若把y表示为x的函数,则各段相轨迹可用函数y=C(jx)来表示(j=1,2,3,4,……).例如:求P1点所对应的时间,相轨迹P0P1可表示为y=C(1x),于是式(1)可写为dxdt=C(1x),dydt=(fx,y)"$$$$$#$$$$$%.(2)图1相轨迹上各点时间的标定Fig.1Differenttimegivenbyphase根据式(2),可得:dt=dxC(1x),对其两边积分,则有t1t0乙dt=x(x(tt0)1)乙C(d1xx).(3)由式(3)可以求得t1如下:t1=t0+x(t1)x(t0)乙C(d1xx).(4)对于任意的相轨迹,有的部分曲率不大,有的部分曲率较大.曲率不太大的部分,二相邻状态点pj和pj+1间的距离若足够小,则可用一直线段来表示其间的相轨迹,如图2所示.曲率较大的部分,可以用一段圆弧来表示,如图3所示.在图2所示的直线段相轨迹情况下,由于y=y(t)j+m[x-x(t)j],其中,m表示直线线段的斜率.根据式(4),tj+1可用下式计算:tj+1=tj+x(tj+1)x(t)j乙y(t)j+md[xx-x(t)j]=tj+1mlny(y(ttj+)j1)(m0)(5)当m=0时(相轨迹平行于水平轴),根据y=y(tj+1)和y=dxdt可求得tj+1:tj+1=tj+x(tj+1)-x(t)jy(tj+1).(6)在图3所示的圆弧线情况下,由于y=姨2-(x-xQ)2,根据式(4),有tj+1=tj+x(tj+1)x(t)j乙姨2-(dxx-xQ)2=tj+cos-1x(t)j-xQ-cos-1x(tj+1)-xQ=tj+i-i+1=tj-.(7)其中Q为圆弧圆心,为弧线的曲率半径,为圆弧所张的圆心角.可见,通过一段圆弧所需的时间与该圆弧所张的圆心角有关,与弧长无关.设此段圆弧为BA,如图4所示,根据式(7),则tBA可用圆心角AOB表示.SA是圆周上A点的切线,过B点作BMOA,则:AOB=2SAB=2ABM,又ABM=tanAMBM,所以AOB=2arctanABMM,因此tBA=2arctanAMBM[6].因此,相轨迹为圆弧时,tj+1又可用下式计算,即tj+