核递推最小二乘辨识算法仿真研究
最小二乘参数辨识的matlab仿真程序注释与剖析

最小二乘一次完成算法的MATLAB 仿真 例2-1 考虑仿真对象)()2(5.0)1()2(7.0)1(5.1)(k v k u k u k z k z k z +-+-=-+-- (2-1)其中,)(k v 是服从正态分布的白噪声N )1,0(。
输入信号采用4阶M 序列,幅度为1。
选择如下形式的辨识模型)()2()1()2()1()(2121k v k u b k u b k z a k z a k z +-+-=-+-+ (2-2)设输入信号的取值是从k =1到k =16的M 序列,则待辨识参数LS θˆ为:LSθˆ=L τL 1L τL z H )H H -( (2-3)其中,被辨识参数LSθˆ、观测矩阵z L 、H L 的表达式为 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=2121ˆb b a a LSθ,⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=)16()4()3(z z z L z ,⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=)14()2()1()15()3()2()14()2()1()15()3()2(u u u u u u z z z z z z L H(2-4)例2-1程序框图如图2.1所示:例2-1Matlab仿真程序如下:%二阶系统的最小二乘一次完成算法辨识程序,文件名:FLch3LSeg1.mu=[-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,1]; %系统辨识的输入信号为一个周期的M序列z=zeros(1,16); %定义输出观测值的长度for k=3:16z(k)=1.5*z(k-1)-0.7*z(k-2)+u(k-1)+0.5*u(k-2); %用理想输出值作为观测值endsubplot(3,1,1) %画三行一列图形窗口中的第一个图形stem(u) %画输入信号u的径线图形subplot(3,1,2) %画三行一列图形窗口中的第二个图形i=1:1:16; %横坐标范围是1到16,步长为1plot(i,z) %图形的横坐标是采样时刻i, 纵坐标是输出观测值z, 图形格式为连续曲线subplot(3,1,3) %画三行一列图形窗口中的第三个图形stem(z),grid on %画出输出观测值z的径线图形,并显示坐标网格u,z %显示输入信号和输出观测信号%L=14 %数据长度HL=[-z(2) -z(1) u(2) u(1);-z(3) -z(2) u(3) u(2);-z(4) -z(3) u(4) u(3);-z(5) -z(4) u(5) u(4);-z(6) -z(5) u(6) u(5);-z(7) -z(6) u(7) u(6);-z(8) -z(7) u(8) u(7);-z(9) -z(8) u(9)u(8);-z(10) -z(9) u(10) u(9);-z(11) -z(10) u(11) u(10);-z(12) -z(11) u(12) u(11);-z(13) -z(12) u(13) u(12);-z(14) -z(13) u(14) u(13);-z(15) -z(14) u(15) u(14)] %给样本矩阵H L赋值ZL=[z(3);z(4);z(5);z(6);z(7);z(8);z(9);z(10);z(11);z(12);z(13);z(14);z(15);z(16)] % 给样本矩阵z L赋值%Calculating Parametersc1=HL'*HL; c2=inv(c1); c3=HL'*ZL; c=c2*c3 %计算并显示θˆLS%Display Parametersa1=c(1), a2=c(2), b1=c(3),b2=c(4) %从θˆ中分离出并显示a1、a2、b1、b2LS%End例2-1程序运行结果:u =[ -1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,1]z =[ 0,0,0.5000,0.2500,0.5250,2.1125, 4.3012,6.4731,6.1988,3.2670,-0.9386, -3.1949,-4.6352,6.2165,-5.5800,-2.5185] HL =1.0000-1.0000-0.5000 0 -1.0000 1.0000 -0.2500 -0.5000 1.0000-1.0000-0.5250 -0.2500 1.0000 1.0000 -2.1125 -0.5250 1.0000 1.0000 -4.3012 -2.1125 1.0000 1.0000 -6.4731-4.3012 -1.0000 1.0000-6.1988-6.4731 -1.0000 -1.0000-3.2670-6.1988 -1.0000 -1.00000.9386-3.2670 1.0000 -1.00003.19490.9386 -1.0000 1.00004.63523.1949 -1.0000 -1.00006.21654.6352 1.0000 -1.00005.58006.2165 1.0000 1.0000(14*4)ZL =[ 0.5000,0.2500,0.5250,2.1125,4.3012,6.4731,6.1988,3.2670,-0.9386,-3.1949, -4.6352,-6.2165,-5.5800,-2.5185]T (14*1)c =[ -1.5000,0.7000,1.0000,0.5000]Ta1 = -1.5000 a2 = 0.7000 b1 = 1.0000 b2 =0.5000-101-10010-10010对比:)()2(5.0)1()2(7.0)1(5.1)(k v k u k u k z k z k z +-+-=-+-- (2-1) 可以看出,由于所用的输出观测值没有任何噪声成分,所以辨识结果无任何误差。
递推阻尼最小二乘法辨识算法公式的详细推导与说明

控制理论与控制工程学位课程《系统辨识》考试报告递推阻尼最小二乘法公式详细推导专业:控制理论与控制工程班级:2011双控(研)学生姓名:江南学号:20110201016任课教师:蔡启仲老师2012年06月29 日摘要在参数辨识中,递推最小二乘法是用得最多的一种算法。
但是,最小二乘法存在一些缺点,如随着协方差矩阵的减小,易产生参数爆发现象;参数向量和协方差矩阵的处置选择不当会使得辨识过程在参数收敛之前结束;在存在随机噪声的情况下,参数易产生漂移,出现不稳定等。
为了防止参数爆发现象,Levenberg 提出在参数优化算法中增加一个阻尼项,以增加算法的稳定性。
本文在一般的最小二乘法中增加了阻尼因子,构成了阻尼最小二乘法。
又根据实时控制的要求,详细推到了递推阻尼最小二乘公式,实现在线辨识。
关键字:系统辨识,最小二乘法,递推算法正文1.题目的基本要求已知单入单出系统的差分方程以及噪声,在应用最小二乘法进行辨识的时候,在性能指标中加入阻尼因子,详细推导阻尼最小二乘法的递推公式。
2.输入辨识信号和系统噪声的产生方法和理论依据 2.1系统辩识信号输入选择准则(1)输入信号的功率或副度不宜过大,以免使系统工作在非线性区,但也不应过小,以致信噪比太小,直接影响辩识精度;(2)输入信号对系统的“净扰动”要小,即应使正负向扰动机会几乎均等; (3)工程上要便于实现,成本低。
2.2白噪声及其产生方法 (1) 白噪声过程(2)白噪声是一种均值为0、谱密度为非0常数的平稳随机过程。
(3)白噪声过程定义:如果随机过程()t ω的均值为0,自相关函数为()()2R t t ωσδ= (2.2.1)式中()t δ 为狄拉克(Dirac) 分布函数,即(){(),00,01t t t dt δδ∞∞=≠∞==⎰-且t (2.2.2)则称该随机过程为白燥声过程。
2.3白噪声序列 (1) 定义 如果随机序列{()}w t 均值为0,并且是两两不相关的,对应的自相关函数为()2,0,1,2w l R l l σδ==±± 式中{1,00,0l l l δ=≠=则称这种随机序列{()}w t 为白噪声序列。
相关_最小二乘两步法在辨识中的应用与改进

图 2 Figure.2 the comparison of R 梯 and R矩
(11)
这种形式便于计算机求解,是目前普 遍使用的求解方法。式中N △ t 要大于系统 的调整时间。同样,自相关函数 Rxx(τ) 和互相关函数 R xv(τ)都是用有限和来 近似计算积分的,引起的截断误差正比于 △ t 2。这样,在辨识中存在的误差源主要 有随机干扰,量化误差,有限和代替积分产 生的误差。其中人为误差只有有限和代替 积分产生的截断误差,结合仿真结果,我们 得出:截断误差仅仅影响到对传递函数零 点的估计,对极点估计没有影响。这样我们 可以快速,高精度地估计出系统传递函数
Байду номын сангаас
的极点,从而判定系统的稳定性。
6 互相关函数计算的梯形法改进
在前面提到,数字计算机中积分是化 为有限和近似计算的。由分析可以得知, 这种有限和计算实际上是积分的矩形近 似,引起的截断误差正比于△ t 2。我们可 以用梯形法作一改进,以减小截断误差。
将梯形近似公式
(12) 运用在互相关函数中,可得:
(2) 则(1 )式可以表示为:
(3) 所以只要能正确的利用输入、输出序 列估计出θ就能得到被辨识系统的传递函 数, 最小二乘法是一种估计θ值的方法。 3.2 辨识方法 3.2.1 互相关函数法 相关性可用互相关函数来表示:
(4) τ是时间间隔。 当 x(t)=y(t)时,互相关
函数就变成了自相关函数。 设 t=0 时,x(t)=0,而且 t
若想得到系统的精确模型可以采取减小截断误差的方法如梯形近似法甚至四阶龙10可以把r看作是系统在rxvxx入下的响应这样要获得该系统的参数模型只需把分别看作xxxv系统的输入和输出然后按最小二乘法估计参数
递推最小二乘算法 simulink

递推最小二乘算法simulink 递推最小二乘算法是一种常用的数学算法,用于估计数据的最佳拟合曲线。
Simulink是一款强大的系统建模和仿真工具,它可以方便地进行算法的模拟和验证。
本文档将介绍如何在Simulink中使用递推最小二乘算法来实现数据的最佳拟合。
首先,我们需要明确递推最小二乘算法的基本原理。
递推最小二乘算法是一种迭代算法,根据已知的数据点,通过不断更新拟合曲线的参数,使得拟合曲线与实际数据的残差平方和最小化。
其核心思想是通过迭代计算,不断调整参数,逐步优化拟合效果。
在Simulink中,我们可以使用递归结构来实现递推最小二乘算法。
首先,我们需要建立一个递推模型框架,包括数据输入、参数更新、残差计算和拟合曲线输出等组件。
其次,我们需要确定递推算法的初始参数值,通常可以使用数据的平均值作为初始参数。
然后,我们可以通过迭代计算的方式,不断更新参数值,直到拟合曲线达到最佳效果。
在模型的具体实现中,我们可以使用Simulink中的线性系统模块、运算模块和数据存储模块等进行建模。
通过连接和配置不同的模块,我们可以构建一个完整的递推最小二乘算法模型。
在模型的验证过程中,我们可以使用Simulink提供的仿真功能,输入实际数据,并观察拟合曲线与实际数据的吻合程度。
需要注意的是,在使用递推最小二乘算法之前,我们需要确定一些参数,如迭代次数、收敛准则和阈值等。
这些参数的选择将直接影响到递推算法的收敛速度和拟合效果。
因此,我们需要根据具体的应用场景和数据特点进行合理选择和调整。
总结来说,递推最小二乘算法是一种有效的数据拟合算法,可以在Simulink中得到方便的实现和验证。
通过合理配置模型和参数,我们可以得到与实际数据最佳拟合的拟合曲线。
在实际应用中,递推最小二乘算法具有广泛的应用领域,如信号处理、系统辨识和机器学习等。
希望本文档能够帮助您更好地理解和应用递推最小二乘算法。
递推最小二乘估计及模型阶次辨识

实验二 递推最小二乘估计(RLS)及模型阶次辨识(F-Test )1 实验方案设计1.1 生成输入数据和噪声用M 序列作为辨识的输入信号,噪声采用标准正态分布的白噪声。
生成白噪声时,首先利用乘同余法生成U[0,1]均匀分布的随机数,再利用U[0,1]均匀分布的随机数生成标准正态分布的白噪声。
1.2 过程仿真辨识模型的形式取)()()()()(11k e k u z B k z z A +=--,为方便起见,取n n n b a == 即nn n n zb a b z b z B z a a a z a z A ------++++=++++=...1)(...1)(22112211用M 序列作为辨识的输入信号。
1.3 递推遗忘因子法数据长度L 取534,初值⎪⎪⎪⎩⎪⎪⎪⎨⎧⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡==1000010000100001)0(001.0)0(P θ 1.4 计算损失函数、噪声标准差损失函数⎥⎦⎤⎢⎣⎡+---+-=μθμττ)()1()()]1(ˆ)()([)1()(2k h k P k h k k h k z k J k J噪声标准差θλdim )(ˆ-=L L J1.6 F-Test 定阶法计算模型阶次统计量t)22,2(~222)1()1()()1,(----++-=+n L F n L n J n J n J n n t其中,)(∙J 为相应阶次下的损失函数值,L 为所用的数据长度,n 为模型的估计阶次。
若a t n n t >+)1,(,拒绝00:n n H >,若a t n n t <+)1,(,接受00:n n H >,其中αt 为风险水平α下的阀值。
这时模型的阶次估计值可取1+n 。
1.6 计算噪信比和性能指标噪信比22ye σση= 参数估计平方相对偏差i i i ni i i θθθθθδˆ~,~1221-=⎪⎪⎭⎫ ⎝⎛=∑= 参数估计平方根偏差ii i n i ini iθθθθθδˆ~,)()~(2122122-==∑∑== 2 编程说明M 序列中,M 序列循环周期取15124=-=p N ,时钟节拍t ∆=1Sec ,幅度1=a ,特征多项式为1)(56⊕⊕=s s s F 。
系统辨识实验二

《系统辨识》实验二要点实验二 递推最小二乘估计(RLS)及模型阶次辨识(F-Test )一、实验目的① 通过实验,掌握递推最小二乘参数辨识方法 ② 通过实验,掌握F-Test 模型阶次辨识方法二、实验内容1、仿真模型实验所用的仿真模型如下: 框图表示模型表示)()2(5.0)1()2(7.0)1(5.1)(k v k u k u k z k z k z λ+-+-=-+-- 其中u (k )和z (k )分别为模型的输入和输出变量;v (k )为零均值、方差为1、服从正态分布的白噪声;λ为噪声的标准差(实验时,可取0.0、0.1、0.5、1.0);输入变量u (k )采用M 序列,其特征多项式取1)(4⊕⊕=s s s F ,幅度取1.0。
2、辨识模型辨识模型的形式取)()()()()(11k e k u z B k z z A +=--为方便起见,取n n n b a ==,即nn nn zb z b z b z B z a z a z a z A ------+++=++++= 22112211)(1)(根据仿真模型生成的数据{}L k k u ,,1),( =和{}L k k z ,,1),( =,辨识模型的参数n n b b b a a a ,,,,,,2121 和;并确定模型阶次n ,同时估计出模型误差)(k e 的方差(应近似等于模型噪声)(k v 的方差,即为2λ)和模型的静态增益K 。
3、辨识算法① 采用递推遗忘因子法:[][][]⎪⎪⎩⎪⎪⎨⎧--=+--=--+-=-)1()()(1)()()1()()()1()()1()()()()1()(1k k k μk μk k k k k k k k k z k k k P h K I P h P h h P K h K τττθθθ 其中,遗忘因子10≤<μ(具体值根据情况自已确定);数据长度L 可取100、300、500;初始值⎩⎨⎧==IP 2)0()0(a εθ。
递推最小二乘法原理

递推最小二乘法原理递推最小二乘法(Recursive Least Squares, 简称RLS)是一种经典的自适应滤波算法,它在信号处理、通信系统、控制系统等领域得到了广泛的应用。
本文将介绍递推最小二乘法的原理及其在实际应用中的一些特点。
首先,让我们来了解一下最小二乘法。
最小二乘法是一种数学优化方法,用于寻找一组参数,使得给定的模型与观测数据之间的误差平方和最小。
在线性回归问题中,最小二乘法可以用来拟合一个线性模型,以最小化观测数据与模型预测值之间的差异。
最小二乘法的基本思想是通过最小化误差的平方和来寻找最优的参数。
递推最小二乘法是最小二乘法的一种变种,它的特点在于可以实时地更新参数估计,适用于需要动态调整的系统。
在实际应用中,由于系统参数可能随时间变化,传统的最小二乘法在每次参数更新时都需要重新计算整个数据集,计算复杂度较高,不适合实时性要求高的场景。
而递推最小二乘法则可以通过递推的方式,实时地更新参数估计,适用于动态环境下的参数估计问题。
递推最小二乘法的原理可以用数学公式来描述。
假设我们有一个线性模型,\[y_k = \theta^T x_k + e_k\]其中\(y_k\)是观测数据,\(x_k\)是输入向量,\(\theta\)是待估计的参数,\(e_k\)是噪声。
我们的目标是通过观测数据\(y_k\)和输入向量\(x_k\)来估计参数\(\theta\)。
递推最小二乘法的核心思想是通过递推的方式,实时地更新参数\(\theta\)的估计值。
具体来说,我们可以通过以下递推公式来更新参数\(\theta\)的估计值,\[\theta_k =\theta_{k-1} + \frac{P_{k-1}x_k}{1 + x_k^T P_{k-1} x_k}(y_k x_k^T \theta_{k-1})\]其中\(\theta_k\)是第\(k\)次的参数估计值,\(\theta_{k-1}\)是第\(k-1\)次的参数估计值,\(P_{k-1}\)是第\(k-1\)次的参数估计误差的协方差矩阵。
递推最小二乘法原理

递推最小二乘法原理递推最小二乘法(Recursive Least Squares, RLS)是一种经典的自适应滤波算法,广泛应用于信号处理、通信系统、控制系统等领域。
它通过不断地更新参数估计,实现对信号的实时跟踪和预测。
本文将介绍递推最小二乘法的原理及其应用。
首先,我们来了解一下最小二乘法(Least Squares, LS)的基本原理。
最小二乘法是一种常见的参数估计方法,它通过最小化观测数据的残差平方和来确定参数的估计值。
对于线性模型,最小二乘法的估计值可以通过求解正规方程或者利用矩阵运算来得到。
然而,在实际应用中,数据通常是逐步到达的,因此需要一种能够实时更新参数估计的方法,这就是递推最小二乘法。
递推最小二乘法的核心思想是利用递推的方式不断更新参数估计,以适应数据的动态变化。
在每次新的数据到达时,根据当前的参数估计和新的数据,通过一定的计算方法得到更新后的参数估计。
这样,参数估计可以随着时间的推移而不断优化,从而更好地适应信号的变化。
在递推最小二乘法中,通常会引入一个遗忘因子(ForgettingFactor),用于衰减历史数据的影响,使算法更加关注最近的数据。
遗忘因子的选择对于算法的性能有着重要的影响,合理的遗忘因子可以平衡对历史数据和最新数据的重视程度,从而提高参数估计的准确性和稳定性。
除了参数估计,递推最小二乘法还可以用于信号预测和滤波。
通过不断更新参数估计,可以实现对信号的实时跟踪和预测,从而在控制系统和通信系统中发挥重要作用。
在实际应用中,递推最小二乘法通常与自适应滤波器相结合,构成自适应滤波系统,用于抑制噪声、提取信号等。
总之,递推最小二乘法是一种重要的自适应参数估计方法,具有实时性强、适应性好的特点,适用于动态环境下的信号处理和控制。
通过不断更新参数估计,可以实现对信号的准确跟踪和预测,为实际应用提供了有力的支持。
在未来的研究和工程实践中,递推最小二乘法仍将发挥重要作用,为各种领域的应用提供有效的解决方案。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
qa K L ) e ne…于 20 0 4年 提 出 能 实现 在线 实 时 辨 识 , 因而 得 到 了 广 泛 的应 用 。然 su r , R S 是 由法 国学者 E gl
而, 该方 法仅 是一 种 线性 辨识 方法 , 于非 线性 系 统 的 , 是 将 传 统 R S算 法 在 特 征 空 间 中 的直 接 应 对 它 L
KRL F —KR S,i ic se .T e KRL lo t m n a i o a S ag rt m o a e rn n i e r y tm d n f ai n a S, B L sd s u s d h S ag r h a d t d t n l i r i RL o h i c mp r d f o l a s l i s o n s e ie t c t . i i o b l y b h i l t ne a l .T e s lt n r s l h w a e KRL e t c t n meh d i b i u l u e o ot e t dt n i t y t e smu ai x mp e h i ai e ut s o t tt i o mu o s h h S i ni ai t o so vo sys p r r h r i o ・ d i f o i t a i l RL n f ai t o h te n t d n f a i c u a y,s b l y o n i n ef r n e a S i e i c t n meh d w eh ri e ie i c t n a c r c dti o h ti o t ii ra —i tre e c .Ho v ri e tr e KRL e t a t t we e n t h e h S i ni d —
方 法是 一种 非 常经 典 的且 应 用最 广泛 的 系统辨 识 方 扰 严 重 的情 况 下 , 识 的稳定 性也 不令 人满 意 。 辨 法 。由于该 方法 不 仅 可 减 少 计 算 量 和存 储 量 , 且 而 核 递 推 最 小 二 乘 法 ( e e R crv es kr l eus elat— n i
i e t c t n e e t h n t e AL —K S a d F d n i ai f c a h D i f o t RL B—KRL t o .t e S —KRL d n f ain meh d i a t u a l ut l i n S me d h W h S i e t c t t o sp ri l y s i e t i i o c r b a me
核递 推 最 小 二 乘 辨识 算 法 仿 真研 究
李 志农
( . 州大学 , 南 1郑 河
赵 匡 何
南 昌 30 6 ) 3 0 3
况
郑州 4 0 0 ;. 昌航空大学 , 50 12 南 江西
[ 关键词 ] 核方法 ; 非线性 回归 ; 递推最小 二乘 ; 稀化表示 [ 摘 要] 论述 了核递 推最小二乘辨识思想 和三种典型算法 即 A D—K K 、W — R S 和 F L R S S K L , B—K L , R S 通过仿真研究 , 比较 了
K L R S法 比其他两种 K L R S辨识算法获得 了更好 的辨识效果 。S —K K W R S法特别适用 于时变非线性系统辨识 。
[ 中图分类号 ]T 7 T 15 3 T 9 1 H1 ; H 6 . ;N 1
[ 文献标志码 】A
[ 文章 编号】 10 — 9 6 2 1 )2— 0 1— 6 0 1 4 2 (0 1 0 0 0 0
Ke r s k r e t o y wo d : en lmeh d;n n i e ri e t c t n;r c r ie la t q a e;s a s e r s na in o l a d n f ai n i i o e u sv e s —s u r p e r p e e tt . r o
2
21年6 01 月
第2 5卷 第 2期
核递推最小二乘辨识算法仿真研究
用 , 过 核 函数 将 原 始 的非线 性 数 据 映 射 到 高 维特 通
由上式 可看 出 , 在每次 迭代 中 , 了求估 计 权 值 为
征空 间 , 其线 性 可分 , 使 然后 在特征 空 间 中实施 递推 向量 & , 要计 算核 矩 阵 K 的逆 , 核 矩 阵 的维 数 都 且 这样 , 矩阵 的求逆 涉 及 的 最小二乘算法 , 巧妙 的将非线性 问题近似为线性 问 随着数 据点 的增加 而增 大 ,
mn( ) l l=∑l l i W = I 一T J e d
() 度 过 高 , 数据 点 将 被 丢 掉 ; 2 则该 如果 6 过 一 个 给 超
其 中 , 是 一 个 向量 , 含 当前 所 有 数 据 的期 望 输 定 的阀值 , 认 为 该 数 据 点 不 能 被 已有 的数 据 点 d 包 则 出d =[ d】o 是更 新 向量 权 值 。每 次迭 代 线 性表 出 , 么 已有 数 据 集 应 该 扩 展 , 纳 该 数 据 d , T , 那 接 中通过 新 到来 的数 据点 和期望 输 出来更 新权 值 向量 点 。由此可 见 , 技 术 相 当于 只使 用 了在 给 定 精 度 该 直到 满足迭 代 终 止 条件 , 最终 输 出 的 即为 所 求 下 样本数 据 中线 性 无 关 的点 , 而最 终 选 取 的那 些数
2 1 年 6月 0 1 第2 5卷 第 2期
南昌航空大学学报 ( 自然科学版 )
J UR O NALOFNA HAN HAN NC G GKO I RST NA UR CE E NG UNVE IY( T ALS INC S)
du 11 n 20 Vo . 5. 12 7 修 20 — 6 0
[ 基金 项 目] 国家自 然科学 基金项 目 505 25750)湖南省机械设备健康维护重点实验室开 (1 3 ; 728 ; 77 0 放基金项目 200) ( 94 0
[ 作者 简介 ] 李志农(96 16一 ) 男, , 南昌航空大学无损检测技术教育部重点实验室教授 , 博士。主要研究方向: 智能检测与信号处理。
S m u a in Re e r h o r e c r i e Le s q a e i l t s a c fKe n lRe u sv a tS u r o
I e i c to g rt m d ntf a i n Al o ih i
L i—n n IZh o g, ZHAO a g Ku n HE a 2 Ku ng
( . hn zo n e i , hnzo , ea 50 1 hn 1Z eghu U i nt Z eghu H n n4 0 0 ,C i v y a; 2 N nh n ag og U i  ̄ t, aca g J nx 3 0 3 C i ) . a cagH n kn n e i N nh n ,i g i 0 6 , n v y a 3 h a
6 =m n【 a ( ( )l i l i c)一 I
‘= l
() 5
其 中 , 一 个 向量 , ( 代 表 已有 数据 的任 意 a是 0 c)
线性组合。如果 6 不能超过一个给定的阀值 , 则
认 为该新 到来 的数 据点 与 已有数 据点 的线 性依 赖 程
—
v ri g n n i e r s se i e t c t n a yn o n a y t m d n i ai . l i f o
递推 最 d-乘 法 ( eus elat qae R S x R crv s—sur, L ) i e
的辨识 , 往往 很 难 得 到令 人 满 意 的结 果 , 另外 , 干 在
的权值 向量 , 或称 之 为滤 波器 系数 。
据 点所组 成 的集合 即是 给定 精度 下样 本数据 的极 大
核 递 推最小 二 乘 ( R S) 法 则 是 在 特 征 空 间 线 性无 关组 。这 样 一 来 , 与计 算 的相 当于 只 有 极 KL 算 参 中实 施 标 准 最 小 二 乘 算 法 。在 迭 代 次 数 为 n时 , 大 线性 无关 组 中 的数 据 , 么 核矩 阵 的维 数 也 就 被 那
AbtatT ekre rcri es sur K L )iet ct n ter ad t e pclagr m,ie L —K L ,W — src:h enl eus el t q a v a e( R S d n f ao h o n h et i l i i i i y r y a ot h . .A D R SS
1. ALD —KRLS 1
1 核 递推最小二乘辨识算法
给定 系统 的输 入输 出模 型为 Y n H( ) )+ ( ) ( )= n R( e n
近似线 性依 赖 技 术 ( pr iael er eed A po m t i a pn— x n d
e c , L 技 术 由法 国学 者 E gl 2 0 n y A D) ne 于 0 4年 提 出 。 其 基本 原理 为 当一 个 新 数 据 点 到来 时 , 要 考 虑 它 需 ( ) 与 现有 数据 集 的线 性依 赖程 度 , 过计 算误 差 : 1 通
i t to f ain meh d,b s d o h W —KR S meh d h st e u i u l i gwi d w p r r n e,t e S —KRL t o b mn et r c o ae n teS L t o a h nq e si n n o e o ma c h W d f S me d o t sb t h e
传统递推 最小二乘 ( L ) R S 辨识算法 和核递 推最小二乘 ( R S 辨识算法对非线性 系统的辨识能力 。仿 真研究表 明 , K L) 不论是 在 辨识精度 , 稳定性还是抗 干扰性方面 ,K K R S辨识算 法 明显 优 于传 统 R S辨 识法 。在这 三种典 型的 K L K R S辨识 算法 ,S — W