CN104808488A - 一种适用于不确定性系统的非线性度量方法 - Google Patents
一种适用于不确定性系统的非线性度量方法 Download PDFInfo
- Publication number
- CN104808488A CN104808488A CN201510098749.9A CN201510098749A CN104808488A CN 104808488 A CN104808488 A CN 104808488A CN 201510098749 A CN201510098749 A CN 201510098749A CN 104808488 A CN104808488 A CN 104808488A
- Authority
- CN
- China
- Prior art keywords
- function
- nonlinear
- linear
- nonlinearity
- linear function
- 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.)
- Granted
Links
Landscapes
- Filters That Use Time-Delay Elements (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种适用于不确定性系统的非线性度量方法,包括以下步骤:设函数空间F表示具有特定分布的随机变量x构成的所有函数集合,将所述函数空间F分为线性函数集合L和非线性函数集合G,设非线性函数gk的非线性度量MoN为非线性函数gk与线性函数集合L的偏差;得有非线性函数gk和线性函数集合L之间的接近程度,然后设定非线性度量vk,在求解离非线性函数gk最近的线性函数,然后根据离非线性函数gk最近的线性函数得不确定性系统的非线性度量v。本发明可以获得不确定性系统的非线性度量。
Description
技术领域
本发明涉及一种非线性度量方法,具体涉及一种适用于不确定性系统的非线性度量方法。
背景技术
任何领域的非线性问题通常都要比线性的问题更难处理。处理的难度伴随着非线性问题的非线性度(Degree of Nonlinearity,DoN)的增加而增加。尽管确定一个系统是不是非线性的并不难,但是仅仅知道系统是非线性的并不够,更希望知道的是这个系统有多么非线性。也就是说,怎样去量化一个问题的非线性。关于此类问题的量化信息正是处理此问题时的困难所在,尤其是在比较不同的问题时。以非线性滤波举例,目前已开发出许多不同的非线性滤波算法,如扩展卡尔曼滤波(EKF),无味卡尔曼滤波(UF),粒子滤波(PF)等等,这些非线性滤波算法在适用性和计算复杂性上各有不同。如果知道了系统的非线性度(DoN),就可以帮助用户明智的选择非线性滤波器。另外,用一个线性化的系统来近似非线性的系统是比较常用的办法,这样可以显著地简化分析。然而,这种方法只有当非线性较弱时效果才会较好,因此,必须寻找一种合适的量化措施,也就是要对非线性进行度量。
非线性度量(Measure of Nolinearity,MoN)的首次研究出现在六十年代的文献[1],[2]中。后来,大量的MoNs被提出用于不同的应用中。这些度量方法可分为两类,第一类是按照非线性函数和与它最接近 的线性函数之间的偏差来度量非线性;Beale在文献[1]中,研究回归分析问题时对非线性度量问题进行的研究,也是第一位对非线性度量的较为严谨的研究。具体方法为:测量一个非线性函数g和基于泰勒级数展开(TSE)的与g最接近的线性函数之间的距离。它是一种局部测量,由于g是在某个点x0处,通过一阶TSE进行线性化的。MoN被定义为标准化的“总距离”,即g和它的线性化近似之间的“总距离”,该线性化近似通过x0点的小领域内的多个采样点计算得到。这个局部性度量使用非线性函数和线性近似之间的距离的基本思想构成了这一类问题所有度量的基础。文献[2]中提出非线性控制系统的非线性度量的方法,它运用了“最优”线性近似和范数的概念,比Beale的基于泰勒级数展开的方法要更优越。文献[3]中提出的方法是运用两个线性系统来捕获一个单输入-单输出(SISO)系统的非线性。文中把MoN定义为从g到它的最大下界和最小上界线性函数的距离中更大的那个距离。这个度量的定义主要用于SISO系统,并且如果用于一般的多输入多输出(MIMO)系统,那么它的扩展也很麻烦;另外,即使系统是有界输入有界输出稳定的,两个线性上下界函数很难找。最近文献[4]中提出一种用于确定性控制系统的相对非线性度量。本质上,它是系统的非线性输入输出映射g和它的最优线性近似之间的标准化差。在以上这一类度量方法中,除了Beale的度量方法外,其他几种度量方法都不用计算函数的导数,这对于那些导数很难计算,甚至不存在的应用中是很理想的。然而,它们各有缺点,都不适用于随机系统。
第二类度量方法是,运用函数在某一点的曲率作为非线性度量。文 献[5][6]提出一种基于曲率的非线性度量用于一回归模型。对于函数z=g(x),MoN由它在x处沿某一方向l的一阶和二阶导数确定的。一阶和二阶导数分别是曲线g(x+cl)在点x处的瞬时“速度”和“加速度”向量。其中,c是独立标量变量。在x处基于曲率的MoN定义为N(x)=maxlNl(x)。这种度量被运用到了非线性动态模型[7]或者非线性量测模型[8]的目标跟踪,地面一定目标指示器雷达跟踪[9]和视频跟踪[10],文献[11]通过仿真研究了关于MoN的滤波器性能。这种非线性度量方法有如下特点:容易计算给定导数;有清楚的物理和几何解释;它们是局部度量而不是总体度量;对于导数不存在的许多应用就不能用;它运用非线性函数g的一阶和二阶项目,而忽略了高阶项,不好评判;此外,这种非线性度量已经应用于要么动态模型是非线性,要么量测模型是非线性的目标跟踪中,但对于动态模型和量测模型都是非线性时如何度量其非线性仍不清楚。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供了一种适用于不确定性系统的非线性度量方法,该方法可以获得不确定性系统的非线性度量。
为达到上述目的,本发明所述的适用于不确定性系统的非线性度量方法包括以下步骤:
设函数空间F表示随机变量x构成的所有函数集合,将所述函数空间F分为线性函数集合L和非线性函数集合G,设非线性函数gk的非线性度量MoN为非线性函数gk与线性函数集合L的偏差;
则有非线性函数gk和线性函数集合L之间的接近程度Jk为:
其中,E为关于随机变量xk的期望,L为所有线性函数L(x)=Ax+b的集合,
Lk为线性函数集合L中的第k个线性函数,gk(x)为非线性函数;
L(x)=Ax+b和gk具有相同的维数,则设定非线性度量vk为:
其中,Cgk为gk(x)的协方差矩阵,tr(x)为矩阵的迹,A为斜率,且A≠0,而b为L(x)在y轴上的截距;
设为离gk最近的线性函数,则有
其中,Cx=cov(x),Cgx=cov(g,x)为协方差矩阵,为Cx的逆;由式(5)、(6)及(7),得不确定性系统的非线性度量v,其中
设J(g1,g2)为非线性函数集合G中的两个非线性函数g1和g2间的接近程度,则有
由式(4)得非线性函数gk和线性函数集合L之间的接近程度Jk为:
离非线性函数gk最近的线性函数通过以下方程求解:
得其中,E[Ax+b-g(x)]为Ax+b-g(x)的均值,Cx=cov(x),则有
本发明具有以下有益效果:
本发明所述的适用于不确定性系统的非线性度量方法用函数空间F表示随机变量x构成的所有函数的集合,把该函数空间F分成线性函数集合L和非线性函数集合G,给定一非线性函数gk,它的非线性度定义为gk与线性函数集合L的偏差,也就是说一个点距离子空间L有多远,然后运用最大下界的概念得到非线性度量,本发明适用于任何形式的不确定性系统,不需要求导数,它是一个全局性的度量方法。
具体实施方式
下面将对本发明做进一步详细描述:
本发明所述的适用于不确定性系统的非线性度量方法包括以下步骤:
设离散时间非线性随机系统为
xk+1=fk(xk)+uk+wk (1)
zk=hk(xk)+vk (2)
其中,xk为随机状态,uk为已知的线性控制输入,wk和vk分别为零均值白色过程噪声和量测噪声;
将式(1)和(2)堆叠得:
yk=gk(xk)+Uk+ηk (3)
其中,
yk=[x′k+1,z′k]′,gk(xk)=[fk(xk)′,hk(xk)′]′
Uk=[u′k,0′]′,ηk=[w′k,v′k]′
yk和xk之间为非线性关系,yk在线性控制输入Uk和噪声hk的作用下是线性的。
设函数空间F表示随机变量x构成的所有函数集合,然后把所述函数空间F分成线性函数集合L和非线性函数集合G,设非线性函数gk,则非线性函数gk的非线性度量MoN为非线性函数gk与线性函数集合L的偏差;
设J(g1,g2)表示非线性函数集合G中的两个非线性函数g1和g2间的接近程度,则有
由式(4)得非线性函数gk和线性函数集合L之间的接近程度Jk为:
其中,E为期望E是关于随机变量xk的期望,L是所有线性函数(实际上是仿射)L(x)=Ax+b的集合,Lk为线性函数集合L中的第k个线性函数,gk(x)为非线性函数,L(x)=Ax+b和gk具有相同的维数非线性度量vk为:
其中,Cgk为gk(x)的协方差矩阵,tr(x)为矩阵的迹,A为斜率,且A≠0,而b为L(x)在y轴上的截距;
为离gk最近的线性函数,通过一下方程求解,
这两个等式解为其中,E[Ax+b-g(x)]为Ax+b-g(x)的均值,Cx=cov(x),Cgx=cov(g,x)是协方差矩阵,故有,
其中,Cx=cov(x),Cgx=cov(g,x)为协方差矩阵,为Cx的逆;
结合(5)、(6)及(7),得不确定性系统的非线性度量v,其中
Claims (3)
1.一种适用于不确定性系统的非线性度量方法,其特征在于,包括以下步骤:
设函数空间F表示随机变量x构成的所有函数集合,将所述函数空间F分为线性函数集合L和非线性函数集合G,设非线性函数gk的非线性度量MoN为非线性函数gk与线性函数集合L的偏差,则有非线性函数gk和线性函数集合L之间的接近程度Jk为:
其中,E为关于随机变量xk的期望,L为所有线性函数L(x)=Ax+b的集合,Lk为线性函数集合L中的第k个线性函数,gk(x)为非线性函数;
L(x)=Ax+b和gk具有相同的维数,则设定非线性度量vk为:
其中,Cgk为gk(x)的协方差矩阵,tr(x)为矩阵的迹,A为斜率,且A≠0,而b为L(x)在y轴上的截距;
设为离gk最近的线性函数,则有
其中,Cx=cov(x),Cgx=cov(g,x)为协方差矩阵,为Cx的逆;
由式(5)、(6)及(7),得不确定性系统的非线性度量v,其中
2.根据权利要求1所述的适用于不确定性系统的非线性度量方法,其特征在于,设J(g1,g2)为非线性函数集合G中的两个非线性函数g1和g2间 的接近程度,则有
由式(4)得非线性函数gk和线性函数集合L之间的接近程度Jk为:
3.根据权利要求1所述的适用于不确定性系统的非线性度量方法,其特征在于,离非线性函数gk最近的线性函数通过以下方程求解:
得其中,E[Ax+b-g(x)]为Ax+b-g(x)的均值,Cx=cov(x),则有
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510098749.9A CN104808488B (zh) | 2015-03-05 | 2015-03-05 | 一种适用于不确定性系统的非线性度量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510098749.9A CN104808488B (zh) | 2015-03-05 | 2015-03-05 | 一种适用于不确定性系统的非线性度量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104808488A true CN104808488A (zh) | 2015-07-29 |
CN104808488B CN104808488B (zh) | 2017-07-07 |
Family
ID=53693423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510098749.9A Expired - Fee Related CN104808488B (zh) | 2015-03-05 | 2015-03-05 | 一种适用于不确定性系统的非线性度量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104808488B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002032104A (ja) * | 2000-07-14 | 2002-01-31 | Nagoya Industrial Science Research Inst | 遺伝的アルゴリズムによる非線形システム同定方法 |
CN103759742A (zh) * | 2014-01-22 | 2014-04-30 | 东南大学 | 基于模糊自适应控制技术的捷联惯导非线性对准方法 |
-
2015
- 2015-03-05 CN CN201510098749.9A patent/CN104808488B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002032104A (ja) * | 2000-07-14 | 2002-01-31 | Nagoya Industrial Science Research Inst | 遺伝的アルゴリズムによる非線形システム同定方法 |
CN103759742A (zh) * | 2014-01-22 | 2014-04-30 | 东南大学 | 基于模糊自适应控制技术的捷联惯导非线性对准方法 |
Non-Patent Citations (2)
Title |
---|
杜静静: "基于非线性度量和MLD_MPC的多模型方法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
翟德明: "多视度量和回归学习方法及应用研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN104808488B (zh) | 2017-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shen et al. | The maximum principle for a jump-diffusion mean-field model and its application to the mean–variance problem | |
CN101452575B (zh) | 一种基于神经网络的图像自适应增强方法 | |
CN103902812B (zh) | 一种粒子滤波方法、装置及目标跟踪方法、装置 | |
Mishura et al. | Mixed stochastic differential equations with long-range dependence: Existence, uniqueness and convergence of solutions | |
Vilgis et al. | Mean-field theory of concentrated polyelectrolyte solutions: Statics and dynamics | |
CN103822598B (zh) | 地基sar在时间去相关严重区域的形变监测方法 | |
Ak et al. | Numerical study of Rosenau-KdV equation using finite element method based on collocation approach | |
Tian et al. | Finite-time stability for impulsive switched delay systems with nonlinear disturbances | |
Liu et al. | New stability criterion of uncertain systems with time-varying delay | |
Cao et al. | Effect of weather on agricultural futures markets on the basis of DCCA cross-correlation coefficient analysis | |
CN104808488A (zh) | 一种适用于不确定性系统的非线性度量方法 | |
CN105095153B (zh) | 一种广义数据融合方法 | |
CN104573216A (zh) | 一种天线性能的优化方法及装置 | |
Saini et al. | Ratio estimators for the finite population mean under simple random sampling and rank set sampling | |
Karabutov | Structural identification of systems with distributed lag | |
CN107995684A (zh) | 提升位置指纹的wlan室内定位精度方法及系统 | |
CN104459654B (zh) | 基于逆变换法的空间相关极化sar杂波图像仿真方法 | |
CN103761238B (zh) | 一种折线滑动窗口累积差值对比分析相似程度的方法 | |
Ghysels et al. | Simple Granger causality tests for mixed frequency data | |
Harlim | Numerical strategies for filtering partially observed stiff stochastic differential equations | |
Kiselev et al. | Investigating a two-sector model of economic growth with the Cobb-Douglas production function | |
Forde | On the Markovian projection in the Brunick–Shreve mimicking result | |
Ran et al. | Information fusion multi-stage identification method for multisensor multi-channel ARMA models | |
Gorban | Transformation of hyper-random quantities and processes | |
Dillschneider | Characteristic Function-Based Estimation of Affine Option Pricing Models |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170707 Termination date: 20180305 |