[发明专利]一种刻画振动系统周期振幅与系统参数的关系的方法

说明书

技术领域

本发明公开了一种刻画振动系统周期振幅与系统参数的关系的方法，属于振动控制领域。

背景技术

狭义的振动又称机械振动，是指物体在平衡位置附近作的往复运动。广义的振动是指描述系统的状态变量在其基准值上下交替变化的过程。振动可分为线性振动和非线性振动。线性振动通常是对非线性振动简化和近似，自然界中真实的振动都是非线性的。周期振动是常见的振动形式。许多自然现象都具有周期性，大到天体运行，小到电子旋转。社会现象也其周期特性，如社会经济的涨跌，甚至民族的兴衰。研究这些周期现象可以预测至控制其运动演化规律。

振动系统通常用微分方程来描述，形式如下：

x·=g(t,x,ϵ).---(1)]]>

其中，t是标量，表示时间。x通常是向量，用来描述系统的状态变量，如物体的位移、电流的大小、声音的强弱，种群的数量，经济的总量等。表示状态变量对时间的导数，通常描述系统状态变量的变化速率。ε通常是向量，用来描述系统的参数变量，如惯性、弹性、激励的幅值或频率。g(t,x,ε)是关于时间、状态变量以及参数变量的非线性向量函数。

利用摄动法得到原系统的一阶近似解的调制方程具有如下形式：

X′＝G(X,ε).            (2)

其中X描述原系统状态变量的一阶近似解的幅值和相位，X′表示X对一阶时间尺度的导数，G(X,ε)是非线性向量函数。

方程(2)的右端等于零时，调制方程有平衡解，即原系统状态变量的一阶近似解的幅值和相位都是常量，这对应原系统出现周期解。求解方程(2)的平衡解，等价于求解

G(X,ε)＝0.             (3)

方程(3)称为调制方程的平衡解等价方程，是原系统的周期振动幅值和相位X关于系统参数变量ε的隐式函数。由于G(X,ε)通常是非线性函数，X很难表示成ε的显式函数。因此，数值求解变系数非线性方程(3)成为得到原系统周期振动幅值和参数变量的关系的重要手段。有些软件含有求解非线性方程组或带参数非线性方程组的函数，如MATLAB的fsolve函数。该函数需要事先给出初始值，而且只能得到由该初始值求出的单个参数点处的一个解，但是振动系统周期振幅是参数变量的多值函数，而且刻画周期振幅与参数变量的关系需要在参数变量的整个取值范围进行。

根据一些科技文献，国外有些机构研发了一些用于非线性分析的软件，可以用来计算微分振动系统的周期解，如AUTO-07P，AnT，XPP-AUT等。这些软件通常运行在Linux环境下，或者需要在Windows系统中虚拟出Linux环境，而且其安装运行都很繁杂。国内有些文献中出现过类似的计算，但是其中的计算方法都是针对简单的系统，例如单自由度系统，或者单参数系统。

本发明针对以上情况给出一种刻画振动系统周期振幅与系统参数的关系的方法。本发明所述的方法利用嵌套循环结构遍历所有参数变量的取值范围，在每个参数点处计算周期振幅，从而刻画振动系统周期振幅与系统参数的关系。本发明可以在参数变量的取值范围内求出振动系统的多个周期解，并假设每个解在各自参数点处的邻域是连续的，由此以数值解的收敛性作为退出条件向该邻域的上下边界扩展，从而提高计算效率。本发明可以应用于任意可数维参数变量和任意自由度的振动系统，为振动系统的研究提供了有效的手段，为周期现象的研究提供了可靠的方法。

发明内容

本发明的目的在于，为克服现有方法的不足，给出一种刻画振动系统周期振幅与系统参数的关系的方法。针对具体振动系统的调制方程X′＝G(X,ε)，X＝[X_m]表示m维状态变量，ε＝[ε_n]表示n维参数变量，本发明可以应用于任意可数维参数变量和任意可数维状态变量的调制方程。