几种阻尼比识别的方法3

B．基于多输出时域识别方法 B1 随机衰减

随机衰减方法是一种非常典型的当输入未知识别模态参数方法。由于识别结果，这种方法实际上是一种无参数识别方法，即随机衰减符号差，是对特定的初始条件的自由衰减响应。得到的随机衰减图形可以用来识别系统模态参数。去相关是这一方法的基本理论，一个简单的导数如下：

对于一个单输入单输出的线性系统，任何力输入的系统响应可以这么解释 (0)V(t)x(t)x(0)D(t)xth(t)f()d (B-1)

0其中D(t)是对单位初始位移的响应，V（t）是对单位初始电压的响应，h（t）是脉冲响

应，f（t）是外部输入的力，假设外部输入力f（t）是一个定常的零均值的随机过程，可以证实x（t）也是一个定常的零均值过程，也证明了x（t）的初始条件为0，考虑到系统响应x(t-ti)中的x(ti)要满足以下条件：

Ax(ti)A (B-2)

由于系统假设是线性的，整个系统的响应包含了3部分： 1. x(ti)的系统响应

(ti)的系统响应 2. x3.f（t）的系统响应，其中f（t）假设是随机的并且是定常的，即：

(ti)V(tti)x(tti)x(ti)D(tti)xttih(t)f()d (B-3)

假设X是x(t-ti)的随机过程，F是f(t-ti)的随机过程， x（t）的平均值为：

E[X(t)]Ex(0)|Ax(0)A(0)|AExx(0)At0h(t)E[F()]d (B-4)

(ti)也是一个平均值为0的定常随机由于x（t）是一个平均值为0的定常随机过程，x系统并且与x（t）是独立的，因此：

(0)]E[x(0)|AxA]0 (B-5) E[x假设

AE[x(0)|Ax(t)A]A (B-6)

且

b(t)h(t)E[F()]d (B-7)

0tX（t）的期望值为：

E[x(t)]AD(t)b(t) (B-8)

如果f（t）是零均值、定常、白噪声随机过程，它与x（t）是相互独立的，因此输入的

力是一个白噪声随机过程：

E[F(t)]E[f(t)]0 (B-9)

且E[x(t)]AD(t) (B-10)

理论上来讲，如果输入的力f（t）不是一个白噪声随机过程，b（t）不为0，已经证实了如果输入信号时零-意义的定常随机过程，由随机衰减方法产生的自由衰减系统响应的误差在允许的工程限制之内，由于实际的样本数不可能是无穷的，可用数学平均值来产生随机衰减图。

(t)1NNx(ti) (B-11)

i1它是E[X(t)]的一个近似值。作为拇指的一个条件（？？）建议平均数N取值要大，τ至少要比最低系统频率长3倍。由于模态有关的因素由相应模态下的相关幅值所决定，通常是振动模态，模态参与因素越高就越准确，后来的识别也越准确。

多输入多输出的线性系统的随机衰减方法也相同，除了一个测试点要作为一个参照点，平均时间ti，是：

AxR(ti)A (B-12)

其中xR(t)是系统在参照点的响应，得到的随机衰减图，被证明含有对特定初始条件甚至是其他测试点的确定性系统响应。

随机衰减方法并不仅局限于位移，电压和加速度的测量也可以使用，对于转子轴承系统，由于自然的输入信号，随机衰减方法只用于当一个随机外部输入的力是已知的，这可以通过振动器或者电磁轴承作为输入来得到。

B2 yule-walker 方程

yule-walker 等式是一种直接基于假设输入信号是独立的随机过程的一种方法，这个性质直接被用来消除ARMA模型中的输入信号部分，用最小二乘法来求AR参数。

在离散时间上输入和输出的关系可表述为：

x(k)a1x(k1)apx(kp)b0f(k)b1f(k1)bpf(kp) (B-13)

假设输入的力是一个独立的随机过程:

E[f(i)f(j)]0 ij (B-14)

T从误差公式可以看出，没有参照源：

E[f(kps1)x(kp1s2)]0 (B-15)

T当 s10,s20时,

Ex(k)pi1Taix(ki)x(kp1s)0 (B-16)

当s>0，结合公式16，可以得到yule-walker等式：

E[Rk]E[Tk] (B-17)

其中

Ta1I,a2I,apI (B-18)

xkp1TRkxk1xkppTxkp (B-19)

xkp1TTkxk (B-20)

xkpp公式17中的期望值可以通过求和来近似得到，yule-walker等式可表示为：

RkTk (B-21)

其中

kRkRi2p1ki (B-22)

TkTi2p1i (B-23)

典型的，用最小二乘法来求θ，递归图来求时间变化的AR参数。

事实上，输入的力不可能是单纯的白噪声或者一个独立的随机系统，这些条件可以通过假设输入的力是一个过滤噪声来缓减，这个滤波器假设是稳定地，典型的系统阶次冗余可用来解释这些寄生极点，这些寄生极点的估计是比较不稳定的具有较高的阻尼比。

几种参数识别方法总结 方法 后向ARMA模型 应用条件 1 SISO信号测量 2 非常适合单输入输出的转子轴承系统模型参数识别 1 单输出信号测量 2 或者自由响应输出 随即衰减方法 1 仅仅输出测量 2 输入信号必须是白噪声 3 不能直接应用到转子轴承系统 4 如是电磁轴承且滤波后可用 Yule-Walker 方程 1 仅输出测量 2 输入信号必须是独立的随机过程 1 AR参数 2 用于识别AR参数的模型参数识别 1 随即衰减特征信号或自由响应信号 2 识别结果可以用于后续的其他方法中 识别结果 模型参数 后向AR模型 模型参数 3 不能直接应用到转子轴承系统 4 如是电磁轴承且滤波后可用于转子轴承系统 改进的Yule-Walker 方程 1 仅输出测量 2 输入信号可以是白噪声信号或非白噪声信号 3 电磁轴承且滤波后可用此方法 4 后向AR和其他方法使得AR参数的识别更加准确 改进的随机脉冲输入方法 1 仅输出测量 2 锤击或电磁轴承施加随机脉冲信号 3 降低平均数值 4 后向AR方法和其他改进直接的AR模型方法 方法能提高参数识别精度 1 仅输出测量 2 输入信号是正弦信号 3 通过正弦信号的输入产生大量瞬态响应 4 后向AR识别方法 一般可逆特征值问题的1 仅输出测量 New-Raphson修正方法 2 通过其他方法识别的几个模型参数 3 数值上更严格 4 可能产生复数结果 转换特征值问题的迭代算法 1 仅输入测量 2 通过其他方法识别模型参数 3 计算量更少 4大的收敛范围

几个物理参数 几个物理参数 1 AR参数 2 采用AR参数识别的模型参数识别 1 自由响应信号 2 采用后向AR方法识别模型参数 1 AR参数 2用于识别AR参数的模型参数识别