CN107688554B - 基于自适应傅里叶分解的频域辨识方法 - Google Patents
基于自适应傅里叶分解的频域辨识方法 Download PDFInfo
- Publication number
- CN107688554B CN107688554B CN201710780066.0A CN201710780066A CN107688554B CN 107688554 B CN107688554 B CN 107688554B CN 201710780066 A CN201710780066 A CN 201710780066A CN 107688554 B CN107688554 B CN 107688554B
- Authority
- CN
- China
- Prior art keywords
- frequency
- data
- function
- transfer function
- spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 19
- 238000012546 transfer Methods 0.000 claims abstract description 45
- 230000004044 response Effects 0.000 claims abstract description 41
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 6
- 238000007781 pre-processing Methods 0.000 claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims description 30
- 230000003595 spectral effect Effects 0.000 claims description 17
- 230000003044 adaptive effect Effects 0.000 claims description 11
- 238000012795 verification Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000001427 coherent effect Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000001268 conjugating effect Effects 0.000 claims description 2
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 238000011084 recovery Methods 0.000 claims description 2
- 230000000717 retained effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000005259 measurement Methods 0.000 description 15
- 230000008901 benefit Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000009191 jumping Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Discrete Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Feedback Control In General (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种基于自适应傅里叶分解的频域辨识方法,主要包含两阶段,第一阶段是估计频率响应:首先采集系统的输入和输出数据,然后对数据进行预处理和谱分析,得到系统的频率响应估计值;第二阶段是辨识传递函数模型:主要运用AFD算法,根据第一阶段得到的频率响应估计值将模型表示为加权Blaschke乘积,从而得到离散传递函数。本发明的频域辨识方法通用性强,计算速度快,辨识精度高。
Description
技术领域
本发明属于系统辨识领域,特别是一种基于自适应傅里叶分解的频域辨识方法。
背景技术
建模的方法主要有两种:机理建模与系统辨识方法建模。理论上,基于基本原理建立的数学模型物理意义明确,能够充分反映系统的内在运动机理。然而对于部分复杂的系统,无法对其内部进行精确分析,而系统辨识方法建模是利用系统的输入输出数据,通过辨识算法和误差准则去得到一个具有等价输入输出特性的模型,其优势在于不需要进行详细受力分析和精密测量。
系统辨识可以分为时域辨识和频域辨识。与时域辨识相比,频域辨识无需对模型的结构做任何假设即可得出非参数模型,可以提供系统的初步信息,直观了解系统的特性。通常,在外界干扰比较大的情况下如果采用时域辨识则不能通过一个简单的滤波进行直接处理,还要通过复杂的计算在原有模型基础上再建立相对繁琐的噪声模型。而在频域辨识中采用合适的方式,可以去除指定频段内的噪声,使得最终辨识结果集中于所关注的频段,使得辨识结果准确性得到很大提高。
传统的辨识方法经过多年发展已经变得比较成熟,但是对于实际应用中的一些复杂系统,难以得到精确的模型或令人满意的结果,其主要存在以下局限和不足:1)部分辨识方法存在着不能同时确定系统的结构与参数的问题,需根据假设与尝试来分别进行辨识;2)部分辨识方法要求输入信号已知,且必须具有较丰富的变化,这一条件在许多普通闭环控制系统是可以满足的,但在某些系统中对输入有严格的限制,因此这些方法不便直接应用;3)部分辨识方法运行过程复杂,计算时间耗费大、效率低,而且可能得到的是局部最优解,因此辨识精度相对较低。
发明内容
本发明的目的在于提供一种基于自适应傅里叶分解的频域辨识方法,能够更加精确估计频率响应以及提高辨识模型的精度。
实现本发明目的的技术解决方案为:一种基于自适应傅里叶分解的频域辨识方法,包括以下步骤:
步骤1、采集实验数据:确定采样率与数据长度、选择试验信号、选择采集输入输出数据的方式;
步骤2、数据预处理:对输入输出数据剔除野值,去除均值和趋势项,进行数字滤波以及数据串接;
步骤3、谱分析:采用层叠分窗法得到输入功率谱密度、输出功率谱密度、互谱密度;
步骤4、频率响应辨识:计算频率响应估计值;
步骤5、传递函数辨识:将频率响应估计值作为AFD算法的输入,选择需要拟合的阶数,经过算法的迭代与分解输出一系列参数,得到离散传递函数;
步骤6、传递函数转化:将离散传递函数经过双线性变换得出连续传递函数;
步骤7、使用评估拟合精度高低的代价函数和反映模型预测能力和可靠性的时域验证,对步骤6中连续传递函数进行评价。
本发明与现有技术相比,其显著优点为:
(1)本发明使用层叠分窗法,通过对多段数据上得到的粗略谱估计求平均,可以大幅减小谱估计中的随机噪声,进而得到平滑的谱估计结果;而且在对频率响应进行估计时,同时考虑了输入输出测量噪声的影响,使之更接近于真实值;
(2)本发明运用了一种高效的有理逼近算法-AFD,由于其自适应特性和在能量角度的快速收敛性,辨识得出的模型具有更高的精度,并可以满足后续对系统设计控制器的需求;
(3)本发明的辨识方法通用性强,计算效率快,辨识精度高。
附图说明
图1是本发明基于自适应傅里叶分解的频域辨识方法流程图。
图2是本发明中的频率响应辨识的总流程图。
图3(a)和图3(b)是输入时间历程数据和输出时间历程数据的图像。
图4(a)和图4(b)是本发明中根据四种频率响应估计公式计算得出的结果的幅频特性曲线和相频特性曲线波特图。
图5(a)和图5(b)是本发明中使用AFD算法拟合得出的传递函数与频率响应估计值的对比图,图5(a)为幅频特性曲线图,图5(b)为相频特性曲线图。
图6(a)和图6(b)是本发明中使用CIFER软件得出的连续传递函数与频率响应数据的对比图,其中图6(a)为幅频特性曲线图,图6(b)为相频特性曲线图。
图7是本发明中的AFD辨识方法与CIFER软件进行比较的方案示意图。
图8是本发明中使用测试信号1来进行两种辨识方法时域验证的结果对比图。
具体实施方式
结合图1、图2,一种基于自适应傅里叶分解的频域辨识方法,包括以下步骤:
步骤1、采集实验数据:确定采样率与数据长度、选择试验信号、选择采集输入输出数据的方式;
步骤2、数据预处理:对输入输出数据剔除野值,去除均值和趋势项,进行数字滤波以及数据串接;
步骤3、谱分析:采用层叠分窗法得到输入功率谱密度、输出功率谱密度、互谱密度;
步骤4、频率响应辨识:计算频率响应估计值;
步骤5、传递函数辨识:将频率响应估计值作为自适应傅里叶分解算法(AFD)的输入,选择需要拟合的阶数,经过算法的迭代与分解输出一系列参数,得到离散传递函数;
步骤6、传递函数转化:将离散传递函数经过双线性变换得出连续传递函数;
步骤7、使用评估拟合精度高低的代价函数和反映模型预测能力和可靠性的时域验证,对步骤6中连续传递函数进行评价。
进一步的,步骤2中的剔除野值运用的是拉依达准则,首先计算测量数据列的均值以及标准偏差,以置信概率99.7%为标准,经判断若序列值大于三倍标准偏差,则属于粗大误差然后剔除,若小于则予以保留;
具体公式为:
进一步的,步骤2中的去除趋势项具体是指先用最小二乘拟合出数据的趋势线,然后减去该趋势线;其中的数字滤波是指先根据实验对象的要求设计IIR滤波器或者FIR滤波器,然后将数据通过滤波器进行滤波。
进一步的,步骤3中的层叠分窗法是一种谱估计方法,其具体方法包括以下步骤:
步骤3-1、计算时域数据的傅里叶变换;
步骤3-2、计算谱函数的粗略估计值;
步骤3-3、计算谱函数的平滑估计;
输入自谱的平滑估计:
输出自谱的平滑估计:
互谱的平滑估计:
其中,U是窗函数功率恢复系数,用于修正谱密度的幅值;m是窗口片段的数量;
相干函数表示了在频率点f,由输入自谱线性地引起的输出自谱的部分占整体输出自谱的比例,其取值范围是0~1;若相干函数值满足大于等于0.6且不振荡,则认为辨识精度可接受。
进一步的,步骤4中的频率响应估计公式需要根据具体情况来选择,具体分为以下几种:
(1)H1估计
若只考虑输出测量噪声,由H1估计得到的频率响应值能保证对所有非相关的输出测量噪声情况是无偏的,但结果会导致估计值相对于真实响应偏小;
(2)H2估计
如果假定噪声仅存在于输入测量中,而输出的测量中不含噪声,则H2估计可获得无偏的频率响应估计值。但在实际中,会导致估计值相对于真实响应偏大;
在一般的实际测量中,输入输出噪声是同时存在的,考虑到H1估计和H2估计分别为欠估计和过估计,为了减小误差影响,提出了以下折中算法;
(3)H3估计
(4)H4估计
进一步的,步骤5中的AFD算法是一种有理逼近算法,其具体方法包括以下步骤:
步骤5-2、根据极大选择原理,计算出系数aj,j=1,2…n′;
其中涉及到的公式有:
步骤5-3、根据已有的系数aq和gq,更新下一项gq+1,并计算部分分量的和;
其中涉及到的公式有:
步骤5-4、将gq+1返回到步骤5-2,求出下一个系数aq+1,如此循环,直到求出n′个系数aj,最后得到离散传递函数公式;
式中,mj表示第j项前面的系数;
进一步的,步骤7中的时域验证是一种通过比较在同样激励下模型的输出和测量的系统输出来实现对模型的评价;具体使用了均方误差(Mean Squared Error,MSE)这一指标,公式如下:
而代价函数是用来描述所得到的传递函数的拟合精度,当其值小于100时,大体上反映了飞行动态建模可接受的精度水平,公式如下:
其中,||符号表示每个频率ω处的幅值(dB),∠符号表示每个频率ω处的相位(°),表示频率响应估计值,T表示传递函数模型,nω表示频率采样点的数量,ω1和表示拟合的起始频率值和结束频率值,Wg和Wp分别是幅值和相位平方误差的相对权重,通常取Wg=1.0,Wp=0.01745;
Wγ是加权函数,它依赖于每个频率ω处的相干函数值rxy,公式如下:
下面结合附图和实施例对本发明进行详细说明。
实施例
四旋翼无人飞行器具有造价低、飞行特性良好、独特的位姿控制方式等多方面优点。然而,作为一个多输入多输出、非线性时变系统,四旋翼飞行器的输入量和输出量之间存在耦合,再加上系统本身的不确定性和外部干扰的引入,使得系统的控制问题变得十分复杂。若要实现对四旋翼飞行器的精确控制,必须先建立精准的模型,而系统辨识就是利用系统的输入输出数据,通过辨识算法和误差准则去得到一个可以描述系统等价输入输出特性的模型。
四旋翼飞行器主要包括四个输入通道,分别是横滚通道、俯仰通道、油门通道、偏航通道。而输出响应包括轴向的速度分量、角速度分量(横滚角速率、俯仰角速率、偏航角速率)、机体欧拉角(横滚角、俯仰角、偏航角)等。而本实施例将需要辨识的对象简化为四旋翼飞行器悬停状态点附近的线性化模型,然后运用本发明提出的频域辨识方法来辨识。具体辨识的是以遥控器横滚通道c1为输入,飞行器横滚角roll为输出的单输入单输出传递函数模型。
结合图1、图2,下面将按照步骤进行详细说明:
步骤1,首先进行试验获取四旋翼飞行器的输入输出数据。采集和记录飞行数据的平台是四旋翼飞行器的飞控,飞控中含有惯性测量单元,包括三轴加速度计、三轴陀螺和磁力计等。频域辨识的推荐输入信号是扫频输入,因为扫频信号可在感兴趣的频率范围内提供很均匀的谱激励,从而利于得到好的辨识结果;另外对系统动态中先验知识的不确定性也具有鲁棒性。
本实验扫频输入的频率覆盖范围在0.1Hz到3Hz,采样率设置为10Hz。具体操作如下:四旋翼飞行器首先在高度1.5-2米处实现悬停,同时保持其他通道相对稳定,悬停3s后在横滚通道施加扫频输入,单次采集的数据长度约100s左右,最后回到平衡点悬停一段时间后降落;重复采集5组数据,从通道记录的数据中挑选几组较好的数据作为待辨识数据,并将输入信号做归一化处理。最后得到的一组数据图像如图3(a)和图3(b)所示。
步骤2,在Matlab软件中,编写数据预处理程序,对采集得到的数据进行初步处理。首先运用根据拉依达准则编写的程序,把由于测试仪器工作不稳定而使信号发生偶然跳动或因外界异常干扰的影响导致测量数据中出现的不合理跳点去除;接着运用detrend函数,把数据测量中引入的不必要的直流分量和由于传感器测量时产生的零点漂移、低频性能的不稳定以及周围的环境干扰等导致的趋势项成分去除;然后为了增加信息量,将几组处理后的实验数据串接成一个数据记录;最后根据需要关注的频率范围,选择合适的截止频率和窗函数种类设计FIR低通滤波器,来滤除信号中的高频噪声。
步骤3,对预处理后的信号进行谱分析。Welch法也叫层叠分窗法,是一种在实际谱分析中用于大幅减小谱估计中随机噪声的关键技术。其基本原理是:将原始的时间历程数据记录分割为短的、相互重叠的片段,每个数据片段用窗口函数进行加权,然后对每一段进行傅里叶变换得到谱密度,最后取平均。
在分窗法中,有几个关键点:1、数据重叠在很大程度上能抵消窗口边缘缩减所造成的误差增大,增加重叠比可以减小随机误差,采用50%到80%的重叠比。2、所用的窗函数一般选择汉宁窗,因为渐变的拱形窗口可以减小旁瓣泄露引起的谱估计误差。其中窗口大小的选择需要特别注意,增大窗口尺寸,频率分辨率会提高,能改善低频段辨识结果,但同时会降低高频段的辨识精度,也会减少数据段的数量,从而使随机误差增加;选择小的窗口则正好与此相反。
步骤4,根据求得的谱函数来计算频率响应估计值。在实际辨识应用中,输入输出的测量值中通常包含噪声。输入端噪声来源,如对驾驶杆位移,作动器位移等测量时引入的扰动;输出端噪声来源,如飞行器的动态响应测量量,如陀螺仪、加速度计,受到来自机械、电子、大气扰动等方面的干扰。
选择何种频率响应估计公式(有)更适合于当前应用,取决于测量噪声的主要来源有关。在本实施例中,考虑到输入输出噪声都存在,所以选择了公式得出的频率响应估计值一般以波特图的形式直观显示,如图4(a)和图4(b)所示,分别为四种公式对应的幅频特性曲线和相频特性曲线。
步骤5,运用AFD算法对得到的频率响应估计值进行拟合逼近,得出离散传递函数。首先将频率响应估计值作为算法的初值,然后根据先验知识选择需要拟合的阶数,本实施例选择3阶或者4阶皆可。AFD算法拟合的结果同样以波特图来显示,如图5(a)和图5(b)所示,同时列举得出的4阶离散传递函数如下:
步骤6,由于后续对比的需要,在Matlab中,将离散传递函数经过双线性变换得到连续传递函数。列举得出的4阶连续传递函数如下:
步骤7,得出传递函数模型后,需要进行时域验证来对模型进行评价。模型的辨识是用频域法完成的,但是对于模型进行检验需要在时域中进行,以验证其是否具有较好的预测能力和对输入形式的鲁棒性。具体方法:先用某种测试输入信号激励系统,再测量系统在这种信号作用下的真实输出,最后通过比较在同样激励下模型的输出与测量的真实输出来实现。本实施例将双极方波作为测试信号,因为这种信号可以很好地显示关键的飞行器动态以及模型性能,而且该机动的瞬态特性是常见机动飞行中有代表性的瞬态阶段,是一种确保模型不被过调的很好的检测信号。
最后,为了验证本发明的辨识精度高,在频域辨识领域选取了一款专门针对飞行器的软件CIFER来进行比较。CIFER软件是由艾姆斯研究中心开发的,在众多的系统辨识软件包中,它是第一个完成对端到端的频率响应辨识的集成软件包,实践证明它是解决旋翼飞行器系统辨识中困难问题的有效工具。CIFER在它的用户界面里引入了较多非常实用的有关于飞行器系统辨识的准则,并将许多高效的辨识功能进行了高度集成,目前已被广泛应用于飞行器项目中。
选取由CIFER软件辨识得出的4阶连续传递函数模型,图6(a)和图6(b)用波特图的形式展示了其拟合的效果,而传递函数公式如下所示:
两种方法的比较方案如图7所示,图中的“FRESPID模块”指CIFER中对时域数据进行处理,最终得到频率响应数据的一个模块;而“NAVFIT模块”指CIFER中通过拟合频率响应数据,然后得到零极点式传递函数模型的一个模块。
在得到本发明辨识方法和CIFER软件计算出的传递函数模型后,分别选取多组试验信号对其进行时域验证,然后计算均方误差,结果列表如下:
表1 4阶情况下两种方法的均方误差比较结果
表2 3阶情况下两种方法的均方误差比较结果
从两个表格的数据上来看,本发明方法的均方误差指标比CIFER软件的更小,表明本方法辨识得出的模型精度更高。同时,仅选取4阶模型下的时域验证结果,如图8所示,两种方法所得模型的输出很接近且与真实输出相差不大。
综上所述,本发明首先通过一系列预处理将辨识所需的数据变得优质,接着使用层叠分窗法减小了谱估计的随机误差,然后根据实际情况选择频率响应估计公式来减小测量噪声的影响,得到更加接近真实的系统频率响应。之后运用高效的AFD算法来逼近频率响应估计值,得到离散传递函数模型,最后通过时域验证等方式进行评价,还与专业的CIFER软件进行比较,证明了本发明辨识得出的模型具有很高的精度,可以在一定程度上真实反映系统的输入输出情况,具有可靠的实用性。
Claims (5)
1.一种基于自适应傅里叶分解的频域辨识方法,其特征如下,包括以下步骤:
步骤1、采集实验数据:确定采样率与数据长度、选择试验信号、选择采集输入输出数据的方式;
步骤2、数据预处理:对输入输出数据剔除野值,去除均值和趋势项,进行数字滤波以及数据串接;
步骤3、谱分析:采用层叠分窗法得到输入功率谱密度、输出功率谱密度、互谱密度;
步骤4、频率响应辨识:计算频率响应估计值;
步骤5、传递函数辨识:将频率响应估计值作为自适应傅里叶分解算法的输入,选择需要拟合的阶数,经过算法的迭代与分解输出一系列参数,得到离散传递函数;自适应傅里叶分解算法是一种有理逼近算法,其具体方法包括以下步骤:
步骤5-2、根据极大选择原理,计算出系数aj,j=1,2…n′;
其中涉及到的公式有:
aj=arg max{|<gi,e{a}>|2:a∈D} (18)
式中,D={z∈C:|z|<1}表示单位圆盘,e{a}(z)是算法中定义的一种字典,<>符号表示内积;
步骤5-3、根据已有的系数aq和gq,更新下一项gq+1,并计算部分分量的和;
其中涉及到的公式有:
式中,Bn′(z)是算法中选用的一种有理正交系,也叫Takenaka-Malmquist system;
步骤5-4、将gq+1返回到步骤5-2,求出下一个系数aq+1,如此循环,直到求出n′个系数aj,最后得到离散传递函数公式;
式中,mj表示第j项前面的系数;
步骤6、传递函数转化:将离散传递函数经过双线性变换得出连续传递函数;
步骤7、使用评估拟合精度高低的代价函数和反映模型预测能力和可靠性的时域验证,对步骤6中连续传递函数进行评价。
3.根据权利要求1所述的基于自适应傅里叶分解的频域辨识方法,其特征在于:步骤2中的去除趋势项具体是指先用最小二乘拟合出数据的趋势线,然后减去该趋势线;其中的数字滤波是指先根据实验对象的要求设计IIR滤波器或者FIR滤波器,然后将数据通过滤波器进行滤波。
4.根据权利要求1所述的基于自适应傅里叶分解的频域辨识方法,其特征在于,步骤3中的层叠分窗法是一种谱估计方法,其具体方法包括以下步骤:
步骤3-1、计算时域数据的傅里叶变换;
步骤3-2、计算谱函数的粗略估计值;
步骤3-3、计算谱函数的平滑估计;
输入自谱的平滑估计:
输出自谱的平滑估计:
互谱的平滑估计:
其中,U是窗函数功率恢复系数,用于修正谱密度的幅值;m是窗口片段的数量;
相干函数表示了在频率点f,由输入自谱线性地引起的输出自谱的部分占整体输出自谱的比例,其取值范围是0~1;若相干函数值满足大于等于0.6且不振荡,则认为辨识精度可接受。
5.根据权利要求1所述的基于自适应傅里叶分解的频域辨识方法,其特征在于,步骤7中的时域验证是通过比较在同样激励下模型的输出和测量的系统输出来实现对模型的评价;具体使用均方误差MSE这一指标,公式如下:
代价函数是用来描述所得到的传递函数的拟合精度,公式如下:
其中,||符号表示每个频率ω处的幅值,∠符号表示每个频率ω处的相位,表示频率响应估计值,T表示传递函数模型,nω表示频率采样点的数量,ω1和表示拟合的起始频率值和结束频率值,Wg和Wp分别是幅值和相位平方误差的相对权重;
Wγ是加权函数,公式如下:
式中,rxy是每个频率ω处的相干函数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710780066.0A CN107688554B (zh) | 2017-09-01 | 2017-09-01 | 基于自适应傅里叶分解的频域辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710780066.0A CN107688554B (zh) | 2017-09-01 | 2017-09-01 | 基于自适应傅里叶分解的频域辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107688554A CN107688554A (zh) | 2018-02-13 |
CN107688554B true CN107688554B (zh) | 2021-09-03 |
Family
ID=61155860
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710780066.0A Active CN107688554B (zh) | 2017-09-01 | 2017-09-01 | 基于自适应傅里叶分解的频域辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107688554B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110025308B (zh) * | 2019-04-09 | 2021-09-10 | 澳门大学 | 一种心电特征提取方法、心拍识别方法及装置 |
CN111125885A (zh) * | 2019-12-03 | 2020-05-08 | 杭州电子科技大学 | 一种基于改进克里金插值算法的asf修正表构建方法 |
CN111859606B (zh) * | 2020-04-29 | 2022-12-02 | 北京理工大学 | 一种射手模型参数辨识方法 |
CN111506117A (zh) * | 2020-06-04 | 2020-08-07 | 南京理工大学 | 基于极限环的四旋翼空间圆形编队解耦控制方法 |
CN111695637B (zh) * | 2020-06-16 | 2023-04-25 | 中国科学院国家天文台南京天文光学技术研究所 | 一种机电系统数学模型辨识方法和系统 |
CN112305611A (zh) * | 2020-12-21 | 2021-02-02 | 西南交通大学 | 边坡地震动特性分析方法、装置、设备及可读存储介质 |
CN115420842B (zh) * | 2022-05-10 | 2024-04-12 | 华谱科仪(北京)科技有限公司 | 一种高效液相色谱仪高压泵的精准压力控制方法 |
CN116956650B (zh) * | 2023-09-21 | 2023-11-28 | 沈阳航天新光集团有限公司 | 基于matlab app designer的姿控动力系统力学环境试验条件设计系统 |
CN118250145B (zh) * | 2024-05-28 | 2024-08-09 | 广东车卫士信息科技有限公司 | 车联网通信管理方法、系统、计算机设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106228031A (zh) * | 2016-09-05 | 2016-12-14 | 北京航空航天大学 | 一种基于试飞数据的飞行动力学模型频宽自适应辨识方法 |
CN106524967A (zh) * | 2016-11-07 | 2017-03-22 | 重庆理工大学 | 一种汽车轮心实际行驶位移测量与提取方法 |
-
2017
- 2017-09-01 CN CN201710780066.0A patent/CN107688554B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106228031A (zh) * | 2016-09-05 | 2016-12-14 | 北京航空航天大学 | 一种基于试飞数据的飞行动力学模型频宽自适应辨识方法 |
CN106524967A (zh) * | 2016-11-07 | 2017-03-22 | 重庆理工大学 | 一种汽车轮心实际行驶位移测量与提取方法 |
Non-Patent Citations (4)
Title |
---|
Adaptive Fourier Decomposition and Rational Approximation, Part I Theory;Tao Qian;《International Jounal of Wavelets, Multiresolution and Information Processing》;20141111;第12卷(第5期);第1-13页 * |
一种基于短时机动试飞数据的动力学辨识优化方法;王保印 等;《航空学报》;20170625;第38卷(第6期);第1-12页 * |
基于AFD算法的滚动轴承故障诊断方法;梁瑜 等;《中国铁道科学》;20130131;第34卷(第1期);第95-100页 * |
车辆横摆转动惯量及轮胎侧偏刚度识别方法的研究;白艳 等;《科学技术与工程》;20150331;第15卷(第9期);第125-128页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107688554A (zh) | 2018-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107688554B (zh) | 基于自适应傅里叶分解的频域辨识方法 | |
CN108827299B (zh) | 一种基于改进四元数二阶互补滤波的飞行器姿态解算方法 | |
CN109873586B (zh) | 一种基于高阶滑模观测器的电机机械参数辨识方法及系统 | |
CN105136145A (zh) | 一种基于卡尔曼滤波的四旋翼无人机姿态数据融合的方法 | |
CN105043348A (zh) | 基于卡尔曼滤波的加速度计陀螺仪水平角度测量方法 | |
US8209069B1 (en) | Aircraft sensor anticipator system | |
CN109507706B (zh) | 一种gps信号丢失的预测定位方法 | |
CN110375772B (zh) | 自适应卡尔曼滤波的环形激光器随机误差建模与补偿方法 | |
CN108170127A (zh) | 一种无人机飞行控制系统的故障检测方法 | |
CN106200377A (zh) | 一种飞行器控制参数的估计方法 | |
CN105912013A (zh) | 一种组合体航天器姿态无模型自适应控制方法 | |
Ahmad et al. | Dynamic modelling and open-loop control of a twin rotor multi-input multi-output system | |
CN111538236A (zh) | 基于分数阶近似微分实现阻尼的飞行器纵向过载控制方法 | |
Holzel et al. | Real-time frequency response estimation from flight data | |
CN111427267B (zh) | 一种采用力与力矩自适应估计的高速飞行器攻角跟踪方法 | |
CN109443393B (zh) | 一种基于盲分离算法的捷联惯导信号提取方法及系统 | |
CN110672127B (zh) | 阵列式mems磁传感器实时标定方法 | |
CN116974293A (zh) | 一种编队飞行卫星姿态控制系统可检测性评估方法 | |
CN110888142A (zh) | 基于mems激光雷达测量技术的航天器隐藏目标点测量方法 | |
Hall Jr et al. | System identification for nonlinear aerodynamic flight regimes | |
CN112987054B (zh) | 标定sins/dvl组合导航系统误差的方法和装置 | |
CN111856941B (zh) | 一种基于主动抗扰的自适应终端动态滑模控制方法 | |
CN114091180A (zh) | 基于飞行数据的扰动风定制化建模和大气数据估计方法 | |
CN112129284B (zh) | 一种仅测量角速率的简单过载控制方法 | |
Walker et al. | Real-time flutter analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |