中对曲线进行斜率提取_Chapter9|9.810 揭秘曲线拟合

9.8 揭秘曲线拟合

曲线拟合初看起来像个巫术，但在这儿我想通过几个简单的类推帮助你理解曲线拟合，其实曲线拟合相当直观。下面的例子确实非常简单。在理论部分讨论了系统传递函数和频响函数的一些信息。单自由度系统的传递函数写成部分分式形式如下所示

频响函数方程为

现在这两个方程表明在第一个方程中自变量为“s”，在第二个方程中自变量为“ω”，而函数“h”的值取决于这些自变量。同时也注意到这有两个常量或参数：留数“a”和极点“p”。所以，在一些给定的“ω”值处定义函数“h”值的这些参数，我们称之为模态参数。

现在考察系统传递函数或称为频响函数的系统传递函数切片。唯一要意识到的是单自由度系统的传递函数曲面和频响函数曲线仅仅由两个参数定义，即极点“p”和留数“a”。所以，观察图9-13，意识到只有两个参数定义了曲面和曲线，这非常神奇。

图9-13 系统传递函数和FRF。来源：Vibrant Technology, Inc.

现在让我们后退一步来看一些更简单、更普遍的例子。让我们考虑对一些测量数据进行非常简单的直线拟合。对图9-14中的数据执行最小二乘误差最小化的拟合。可以用任何曲线来拟合该组数据，但似乎一阶拟合最合适。当然，要使用的模型是y=mx+b，定义这条直线需要两个参数，即斜率和y方向的截距。

图9-14 简单的直线拟合例子

在图9-14中，数据的最小二乘拟合得到的两个参数结果为：斜率12.097，y方向截距为-0.019。此外，认识到这两个参数是从一组测量数据中获得的，存在一定的偏差，最小二乘回归分析确定了斜率和y方向截距的最合适的参数，用来表示这些测量数据。

因此，如果将相同的逻辑应用于单自由度系统的频响函数，则可以对图9-15中的数据执行频响函数(如上所述)的二阶模型拟合。从示意图中很容易看到有一组数据和曲线拟合，曲线拟合获得了两个参数，即极点和留数。它与直线拟合完全相同，只是数据是复数，曲线更复杂一些。但原理上，跟直线拟合是相同的方法论。测量的频响函数分布在离散数据点上(是复数)。对频响函数进行曲线拟合，按最小二乘方式能找到描述这组数据最合适的参数。

图9-15 SDOF曲线拟合概念

当然，图9-15中的数据是针对单自由度系统的。这个方法可以扩展到更高阶的方程，如图9-16所示。因此，用这种方式可以对由频响函数测量的离散复数数据进行拟合得到多阶模态(或本质上更高阶的多项式)。所有与模态参数估计过程相关的问题在图9-16中的数据拟合时都会遭遇到。

图9-16 MDOF曲线拟合概念

因此，理解了简单的直线拟合过程，那么，相同的过程也会在模态参数估计过程加以应用(但是，当然数据是复数并且曲线更复杂)。本质上，在这两种情况下，描述函数的参数都是按最小二乘方式提取的。

因此，整个曲线拟合过程真的不是巫术。与我们常用的简单直线回归分析实际上是相同的过程，模态参数估计只是简单数据曲线拟合的延伸。

9.9 曲线拟合不同频带的极点和留数

一般来说，进行模态参数估计时，通常希望用包含少数几阶模态的窄带宽。但是曲线拟合过程分两步，并且有两种不同的方法用于提取极点和留数。在一种方法中，极点可以从几个分离较窄的频带中提取到，留数也要从这些频带中提取到，定义的频带如图9-17所示。但是，在参数估计的第二步获得留数时，有时人们会使用一个包含所有较窄频段的宽频段，如图9-18所示。真正需要关注的是综合的频响函数是否与实际测量数据吻合较好。只要综合的频响函数与测量的频响函数吻合很好，那么，两种方法都是可接受的。

图9-17 相同频带的曲线拟合获得极点和留数

图9-18 包含所有频带的曲线拟合获得留数

9.10 使用多个带宽拼接的模态参数综合FRF

当然，一旦处理了所有的数据，会由提取的极点和留数重新创建或“综合”频响函数，并与原始测量的数据进行比较，以确保数据拟合良好。现在如果所有的模态(极点和留数)都用一个频带来拟合，这是一个直截了当的过程。但是，如果使用多个较窄的频带，则需要将数据拼接在一起以重新创建或综合出频响函数。

假设用三个分离的频带进行数据拟合，如图9-19至9-21所示。每个频带都会拟合得到该频带的模态和这个频带的上下残余项。现在，当三个分离的频带组合在一起拟合时，需要组合每个频带的模态，同时也应包括第一个频带的下残余项和最后一个频带的上残余项，以便生成一个正确的综合频响函数，如图9-22所示。

图9-19 第1个频带的数据拟合

图9-20 第2个频带的数据拟合

图9-21 第3个频带的数据拟合

图9-22 使用来自于第1、2、3个频带的模态(极点和留数)和第1个频带的下残余项以及第3个频带的上残余项综合出频响函数

作者介绍

钱小猛，硕士毕业于浙江大学机械工程学院。具有扎实的模态分析理论功底和丰富的模态测试经验。曾参与开发试验模态分析软件，完成了多个基于Test.Lab的二次开发项目，对振动控制也有较深的研究。钱小猛私人微信：zjuchien001，请注明：单位+姓名

注：翻译自Peter Avitabile《Modal Testing - A Practitioner's Guide》， 

翻译：钱小猛；校对：谭祥军

请你为《模态试验实用技术：实践者指南》的购买册数投票：

觉得不错，请点赞！

正版图书，评论数3300+

扩展阅读

1.交稿了！《模态试验实用技术：实践者指南》交稿了！

2.什么是曲线拟合？

3.各种常见的曲线拟合方法

4.99.02：曲线拟合令我感到迷惑！各种不同的拟合技术意味什么？

5.Chapter9|9.1-2 模态指示工具帮助识别模态

6.Chapter9|9.3-5 SDOF和MDOF拟合，局部拟合和整体拟合对比

7.Chapter9|9.6-7 风机叶片实例与揭秘稳态图