CN114169174A - 一种基于正弦激励的频响测量用两点采样优化方法及系统 - Google Patents

Abstract

本发明公开了一种基于正弦激励的频响测量用两点采样优化方法及系统，设置采样起始频率，终止频率和总采样点数；在起始频率和终止频率范围内等间距采样若干个点作为初始信息；根据已有采样信息，估算各子频段的插值误差；挑选出插值误差最大的子频段；在该子频段内合适的位置新增两个采样点；不断循环直至采样点数达到用户设定的总数；询问用户是否需要新增采样点，如果是，则在用户指定新的采样点数后，重新估算各子频段的插值误差并继续采样，否则采样结束。本发明不依赖特定硬件，仅需算法实现，可以直接嵌入现有频响分析仪，其计算负荷低，易用性高，精度高，稳定性强，且具备很好的数据继承性，从而大大提高了基于正弦激励的频响测量方法的实用性。

Description

一种基于正弦激励的频响测量用两点采样优化方法及系统

技术领域

本发明属于频响测量技术领域，具体涉及一种基于正弦激励的频响测量用两点采样优化方法及系统。

背景技术

目前，在各类复杂系统的分析和调试中，不同的频响测量方法的激励信号仅有两类：复合激励信号和正弦激励信号。复合激励信号主要适用于线性系统，否则会引发频谱混叠现象，降低测量精度；而正弦激励信号对线性和非线性系统都有着良好的普适性，因此，它在频响测量中仍然发挥着不可替代的作用。

在使用正弦激励的频响测量方法中，最传统的扫频法很容易引发测量速度或测量精度的问题。例如，如果测量的步长选得太大，便无法保证测量精度；而步长选得太小，又会大大延长测量时间。同时，扫频法的数据继承性也不好，这意味着它无法很好地重复使用已有的采样信息。这也会降低测量效率。

发明内容

本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种基于正弦激励的频响测量用两点采样优化方法及系统，可直接嵌入现有频响分析仪中，不依赖特定硬件，实现总体插值误差的快速下降。

本发明采用以下技术方案：

一种基于正弦激励的频响测量用两点采样优化方法，根据已有采样信息估算由已有采样点划分的各子频段的插值误差；挑选估算的各子频段中插值误差最大的子频段；在插值误差最大的子频段内新增两个采样点；重复以上步骤，直至采样点数达到总采样点数时，本轮采样结束；增加新的采样点，重新开始估算直至采样结束。

具体的，已有采样信息包括初始采样点及后续迭代计算所得的新增采样点。

进一步的，首先采集N_s个点作为初始信息，然后开始迭代，且后续每个采样点都是基于已采样信息所计算出的全局最优点。

进一步的，总采样点数大于等于4。

具体的，估算各子频段的插值误差具体为：将所有已采样点逐一作为起始点构建三角形，并计算对应的三角形面积，估算得到所有子频段的插值误差。

具体的，挑选估算的各子频段中插值误差最大的子频段具体为：将增益和相位的全部2N-4个插值误差一同进行比较，选出最大的值，并确定对应的子频段[f_i，f_i+1，f_i+2]。

进一步的，新增的两个采样点位于子频段[f_i，f_i+1，f_i+2]的[f_i，f_i+1]之中，或者位于[f_i+1，f_i+2]之中。

进一步的，新增的两个采样点选取[f_i，f_i+1]和[f_i+1，f_i+2]范围内的中点或任一等分点。

具体的，用户新增的采样点数为偶数。

本发明的另一个技术方案是，一种基于正弦激励的频响测量用两点采样优化系统，其特征在于，包括：

估算模块，根据已有采样信息估算由已有采样点划分的各子频段的插值误差；

选择模块，挑选估算的各子频段中插值误差最大的子频段；

优化模块，在插值误差最大的子频段内新增两个采样点；

采样模块，重复估算模块至优化模块，直至采样点数达到总采样点数时，本轮采样结束；

循环模块，增加新的采样点，转入估算模块，直至采样结束。

与现有技术相比，本发明至少具有以下有益效果：

本发明一种基于正弦激励的频响测量用两点采样优化方法，可优化频响测量中每个采样点的位置，使采样点位置的分布更为合理。从频响测量的数学实质出发，将基于正弦激励的频响测量视为对被测系统的频域采样，通过估算已有采样值的分段线性插值模型和被测系统理论模型之间的插值误差，在每个迭代内把两个新增采样点放置于插值误差最大的子频段内，从而实现总体插值误差的最速下降。

进一步的，首先采集N_s个点作为初始信息，然后开始迭代，且后续每个采样点都是基于已采样信息所计算出的全局最优点。其优势在于该方法总是高效地使用每个采样点，使得每个点都是全局最优的结果。不同于商用的扫频法，该方法不是在一开始就已经算出了所有待采样的位置(均匀采样)，而是首先采样少量初始点，再根据初始点的信息逐步迭代增加新的采样点。每个迭代都会利用所有已采样信息进行计算，以得到新增采样点的最佳位置。

进一步的，总采样点数最少需要四个采样点，使得后续的迭代计算可以正常进行。后续迭代中，会根据梯形法则进行插值误差的比较，而这个比较的过程至少需要两个子频段参与，每个子频段包含三个采样点。假设这两个子频段是相邻且部分重合的，即第一个子频段为[f_i，f_i+1，f_i+2]，第二个子频段为[f_i+1，f_i+2，f_i+3]，则总采样点数最少为四。

进一步的，估算各子频段的插值误差e_i具体为：将所有已采样点逐一作为起始点构建三角形，并计算对应的三角形面积，估算得到所有子频段的插值误差。其目的是在系统理论模型未知的前提下，直接利用梯形法则估算采样值和系统理论模型之间的插值误差。这一估算之所以可以实现，是因为本发明将基于正弦激励的频响测量视为一个数学插值问题，因此，插值理论中插值误差估算的方法，如梯形法则，就可以直接应用于频响测量。这将有助于频响分析仪自适应地将新增采样点置于插值误差(也即实际采样误差)最大的位置，使得采样精度不断提升。

进一步的，将增益和相位的全部2N-4个插值误差一同进行比较，选出最大的值，并确定对应的子频段[f_i，f_i+1，f_i+2]。这是为了选出增益和相位中插值误差(也即实际采样误差)最大的位置。此外，将增益和相位一同进行比较，也是为了实现全局比较，即该比较并非仅针对增益或者相位，而是把增益与相位的采样信息同时比较。显然，全局比较能够更为有效地选出插值误差最大的位置。

进一步的，新增的两个采样点位于确定的子频段[f_i，f_i+1，f_i+2]的[f_i，f_i+1]之中，或者位于[f_i+1，f_i+2]之中。这是由于子频段[f_i，f_i+1，f_i+2]本身就包含[f_i，f_i+1]和[f_i+1，f_i+2]两个更小的子频段，因此，单次迭代中直接新增两个采样点以分别采样[f_i，f_i+1]和[f_i+1，f_i+2]就是一个自然而然的选择。同时，这种方法也可以尽可能使两个新增采样点分散开，以采样更多的信息。反例则是将两个采样点放置得过于靠近，如两个点都置于[f_i，f_i+1]之中。尽管这也是算法允许的，但这不是效率最高的方式。

进一步的，新增的两个采样点选取[f_i，f_i+1]和[f_i+1，f_i+2]的中点或其他等分点，这是因为这些位置是固定的，易确定的，且为常见的选点位置。但是，在这些选点策略中，仅有中点是唯一确定的，因此具备最好的数据继承性；其他选择亦可，但总体效率有所降低。实际应用中，应以中点为优先选择。

进一步的，新的采样点N′_m为偶数，这是因为该方法每个迭代都会新增两个采样点，因此，假如用户在采样任务结束后，仍要增加新的采样点，那么，新增采样点数N′_m也必须是偶数，否则算法将出错。

综上所述，本发明具有计算负荷低、用性强、精度高、稳定性强，且数据继承性好的特点。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

图1为本发明流程图；

图2为数值测试中待测系统传递函数的波特图；

图3为本发明中根据已有采样信息构建三角形的示意图；

图4为本发明中两点采样优化方法的采样过程示例图；

图5为数值测试中三种基于正弦激励的频响测量方法的采样误差对比图。

具体实施方式

本发明提供了一种基于正弦激励的频响测量用两点采样优化方法，从频响测量的数学实质出发，将基于正弦激励的频响测量视为对被测系统的频域采样，通过估算已有采样值的分段线性插值模型和被测系统理论模型之间的插值误差，在每个迭代内把两个新增采样点放置于插值误差最大的子频段内，从而实现总体插值误差的最速下降。

请参阅图1，本发明一种基于正弦激励的频响测量用两点采样优化方法，包括以下步骤：

S1、设置采样起始频率，终止频率和总采样点数，在起始频率和终止频率范围内等间距采样N_s个点，，此时采样的N_s个点将作为迭代启动的初始信息；

采样起始频率f_start，终止频率f_end和总采样点数N_m是采样任务必须具备的三个预设参数，其中，起始频率f_start和终止频率f_end指定了采样的频段范围，而总采样点数N_m指定了采样目标系统所需要的采样次数。这三个参数的设置完全取决于测量要求，无需任何经验性法则，因此对用户是非常友好的。同时，总采样点数N_m必须是一个偶数，且满足N_m≥4。

值得注意的是，本发明将采用迭代采样方法，即每个采样点都是基于已采样信息所计算出的全局最优点。本发明将首先采集一些点作为初始信息，也即N_s个点，然后再开始迭代；其同时要求初始采样的N_s个点必须为等间距采样，其中，等间距采样是为了更均匀地采样系统，而N_s常取值为4。实际上，这些初始采样点无需严格满足等间距采样(等间距采样只是最优的情况)；而初始采样点数N_s也只需满足N_s≥4。本发明设置为4，即是设置了理论上允许的最小值。这也是总采样点数N_m必须满足N_m≥4的缘故。

S2、根据已有采样信息，估算各子频段的插值误差e_i，i＝1，2，...2N-4，其中N表示已经采样的点的数量；

使用梯形法则估算插值误差，即针对任意连续三个采样点构成的子频段[f_i，f_i+1，f_i+2]，将其采样值相连构成三角形后，三角形面积约等于该子频段插值误差的三倍。该法则是数学原理中相当经典的插值误差估计方法，因此其对应的数学证明在此处略去。

在应用梯形法则时，只需将所有已采样点逐一作为起始点构建三角形，并计算对应的三角形面积，就可以估算对应子频段的插值误差。注意，插值误差本质就是真实的采样误差的良好估计。由于一共有N个采样点，因此对增益(幅值)和相位分别可以构建N-2个三角形(最后两个采样点无法再继续构建新的三角形)，所以对应的插值误差一共有2N-4个。由于这些三角形面积的计算值可能为正，也可能为负，因此所有计算值都应该取绝对值，以体现采样值和系统理论模型之间的真实差异；同时，由于增益(幅值)和相位的三角形会放在一起比较，因此，也必须对这些三角形的面积取相对值，即相对于增益(幅值)和相位中的三角形面积均值分别取相对值，以实现数量级上的公平比较。而如果这些三角形面积既不取绝对值计算，也不转换为相对值，那么，在绝大多数应用中，相位的插值误差将占据绝对主导，从而无法实现增益(幅值)和相位的公平比较。这是因为相位往往在[-180°，180°]之间，而增益(幅值)却可能有着相当小的数量级。因此，对三角形面积的绝对值计算和相对值计算是本发明中绝不可以忽略的步骤，这也是实现插值误差的全局比较的必要条件。

值得说明的是，尽管本方法在估算插值误差时使用了梯形法则，但其他相似的插值误差估计方法也应涵盖在保护范围内。

S3、挑选出插值误差最大的子频段[f_i，f_i+1，f_i+2]，其对应误差为max(e_i)；

增益(幅值)和相位的全部2N-4个插值误差将一同进行比较，以选出其中最大的值，并找到对应的子频段[f_i，f_i+1，f_i+2]。显然，该子频段的误差就是max(e_i)。需要再次强调的是，子频段的插值误差实际就是采样误差的良好估计，因此，搜寻最大插值误差对应的子频段，即是搜寻采样误差最大的子频段。

此外，目前已有许多筛选最大值的算法，本方法在筛选最大误差时并不局限为某一特定算法。常见算法如逐个比对法，或是分治算法，都可以实现寻找最大值的目的。由于这些算法往往有相当成熟的商业工具包或是开源代码，此处不再详述。

S4、在该子频段内合适的位置新增两个采样点；

新增的两个采样点位于[f_i，f_i+1]之中，或者位于[f_i+1，f_i+2]之中。

两个采样点的最佳位置是[f_i，f_i+1]和[f_i+1，f_i+2]的中点，即一个点采样[f_i，f_i+1]的中点，另一个点采样[f_i+1，f_i+2]的中点。

当然，两个采样点也并不一定严格设置为中点，也可以是其他固定的位置，如三分点，四分点或其他等分点；此处选取为中点，是最稳妥的选择。这是因为系统本身是未知的，含大量未采样信息，因此，中点将是一个保守但稳定的策略。

S5、不断循环步骤S2至步骤S4，直至采样点数N达到用户设定的总数，也即当N＝＝N_m时，本轮采样结束；

N是对已采样点数的计数，每经过一次迭代，都会使N＝N+2，并判断是否已经满足N＝＝N_m。如果还未满足该条件，则转入步骤S2，继续采样新点；而如果已经满足该条件，则本轮采样结束。需要指出的是，由于每个迭代新增两个采样点，因此，本发明只需进行(N_m-4)/2次迭代，即可完成一轮采样任务，这将显著降低计算负荷。

S6、询问用户是否需要新增采样点，如果是，则在用户指定新的采样点数N′_m后，转入步骤S2，否则采样结束。

新增采样点数N′_m必须是一个偶数。本发明允许用户在一轮采样结束后，再次输入一个新的采样点数N′_m，并在全部已有采样值的基础上继续进行采样。这也就是数据继承性，即针对同一个系统，它总是可以基于已有采样值进行新的采样任务，并且，先采样N_m个点，再接着采样N′_m个点，和一开始就直接采样N_m+N′_m个点的结果是完全一致的。在实际应用中，用户在完成一次采样任务后，可能对采样精度不满足，应用本发明则可以在已有采样数据的基础上继续进行采样，不断提升采样精度；并且，新增采样点也仍然具备数据继承性，换言之，这个增加采样点的过程可以持续进行下去，直到采样精度达到了用户的需求。而传统商用扫频法则只能重新开始采样，这将大大延长采样时间，降低采样效率。

请参阅图2，为待测系统传递函数进行补充说明的波特图，也即待测系统的频域特性曲线。为了更好地体现本发明的优势与适用性，此处所选波特图中的曲线特性变化较为复杂，既有平直的部分，也有弯曲陡峭的部分。

本发明再一个实施例中，提供一种基于正弦激励的频响测量用两点采样优化系统，该系统能够用于实现上述基于正弦激励的频响测量用两点采样优化方法，具体的，该基于正弦激励的频响测量用两点采样优化系统包括估算模块、选择模块、优化模块、采样模块以及循环模块。

其中，估算模块，根据已有采样信息估算由已有采样点划分的各子频段的插值误差e_i，i＝1，2，...2N-4，N为采样点数；

选择模块，挑选估算的各子频段中插值误差最大的子频段[f_i，f_i+1，f_i+2]；

优化模块，在插值误差最大的子频段[f_i，f_i+1，f_i+2]内新增两个采样点；

采样模块，重复估算模块至优化模块，直至采样点数N达到总采样点数时，本轮采样结束；

循环模块，增加新的采样点N′_m，转入估算模块，直至采样结束。

在使用本方法时，将依据梯形法则估算所有子频段的插值误差，即对所有已采样点，以每一个点为起点构建三角形，并计算对应三角形的面积，如图3所示。之后，从所有估算的误差中，挑选出相对值最大的一个，并找到相应的子频段，然后在该子频段中新增两个采样点。这一采样过程如图4所示。针对图一所示系统进行数值测试时，选用了扫频法，自适应频率注入法和本发明共同进行测试。从图5可见，随着采样点数的不断增多，三种频响测量方法的总体误差随之下降，并且，本发明总是有着最低的误差。这就证明了本方法的有效性。

综上所述，本发明一种基于正弦激励的频响测量用两点采样优化方法及系统，具有以下优点：

1、计算负荷低，本方法每次迭代将采样两个点，因此对于N_m个采样点，仅需估算(N_m-4)/2次插值误差；

2、易用性强，本方法仅需用户指定采样起始频率，终止频率和总采样点数，无需指定任何其他经验参数；

3、精度高，本方法基于梯形法则估算各子频段插值误差，并总是在插值误差最大的子频段内新增采样点，因此，随着采样点增多，总体插值误差会以几乎最快的速度下降；

4、稳定性强，本方法的采样机制确保了更多的采样点总是有更高的采样精度，因此不会存在稳定性问题；

5、数据继承性好，本方法总是可以基于所有已采样信息进行新的采样任务，从而提高采样效率，体现了良好的数据继承性。

以上内容仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明权利要求书的保护范围之内。

Claims (10)

1.一种基于正弦激励的频响测量用两点采样优化方法，其特征在于，根据已有采样信息估算由已有采样点划分的各子频段的插值误差；挑选估算的各子频段中插值误差最大的子频段；在插值误差最大的子频段内新增两个采样点；重复以上步骤，直至采样点数达到总采样点数时，本轮采样结束；增加新的采样点，重新开始估算直至采样结束。

2.根据权利要求1所述的基于正弦激励的频响测量用两点采样优化方法，其特征在于，已有采样信息包括初始采样点及后续迭代计算所得的新增采样点。

3.根据权利要求2所述的基于正弦激励的频响测量用两点采样优化方法，其特征在于，首先采集N_s个点作为初始信息，然后开始迭代，且后续每个采样点都是基于已采样信息所计算出的全局最优点。

4.根据权利要求2所述的基于正弦激励的频响测量用两点采样优化方法，其特征在于，总采样点数大于等于4。

5.根据权利要求1所述的基于正弦激励的频响测量用两点采样优化方法，其特征在于，估算各子频段的插值误差具体为：将所有已采样点逐一作为起始点构建三角形，并计算对应的三角形面积，估算得到所有子频段的插值误差。

6.根据权利要求1所述的基于正弦激励的频响测量用两点采样优化方法，其特征在于，挑选估算的各子频段中插值误差最大的子频段具体为：将增益和相位的全部2N-4个插值误差一同进行比较，选出最大的值，并确定对应的子频段[f_i，f_i+1，f_i+2]。

7.根据权利要求6所述的基于正弦激励的频响测量用两点采样优化方法，其特征在于，新增的两个采样点位于子频段[f_i，f_i+1，f_i+2]的[f_i，f_i+1]之中，或者位于[f_i+1，f_i+2]之中。

8.根据权利要求6所述的基于正弦激励的频响测量用两点采样优化方法，其特征在于，新增的两个采样点选取[f_i，f_i+1]和[f_i+1，f_i+2]范围内的中点或任一等分点。

9.根据权利要求1所述的基于正弦激励的频响测量用两点采样优化方法，其特征在于，新增的采样点数为偶数。

10.一种基于正弦激励的频响测量用两点采样优化系统，其特征在于，包括：

估算模块，根据已有采样信息估算由已有采样点划分的各子频段的插值误差；

选择模块，挑选估算的各子频段中插值误差最大的子频段；

优化模块，在插值误差最大的子频段内新增两个采样点；

采样模块，重复估算模块至优化模块，直至采样点数达到总采样点数时，本轮采样结束；

循环模块，增加新的采样点，转入估算模块，直至采样结束。

Images