随着控制过程复杂性的提高和对非线性问题认识的深入，非线性系统辨识越来越成为控制问题的关键。建立描述非线性现象的模型是研究非线性问题的基础，如Bilinear模型、Hammerstein模型、Wiener模型和输入输出仿射模型等^[1-3]。但这些模型由于受到自身结构限制，只能描述部分特殊的非线性系统，所以一些普遍意义上的非线性模型相继被提出，如Volterra模型、NARMAX模型、神经网络模型、模糊逻辑模型、支持向量学习机模型和遗传规划(GP)模型等^[4-11]。综观上述非线性模型，各有优势和不足，但其均可用非线性的线参数多项式组合模型来表示。国内外针对此类模型结构和参数的一体化辨识方法开展了一系列研究。文献[12]采用经典Gram-Schmidt(CGS)算法和改进Gram-Schmidt(MGS)算法对NARMAX模型进行子模型集合的选取；文献[13]采用新模型选项准则对NARMAX模型进行结构确定和参数估计研究；文献[14]采用前向正交最小二乘法(OLS)和误差衰减(ERR)准则对线参数非线性多项式模型进行结构预选；文献[15]提出自适应遗传算法与相关度评估结合方法对组合非线性系统进行结构和参数一体化辨识研究。文献[16]采用正交回归与嵌套最优正则化对一类非线性模型进行非线性辨识。

本文提出一种非线性多项式模型结构与参数一体化辨识策略。首先，采用递归改进Gram-Schmidt (RMGS)算法从向量空间的全集中执行“择优过程”，即寻找出一个能够描述输入输出映射关系的集合，该集合是包含核心结构项和少量冗余项的优选集。其次，在优选集合里执行“劣汰过程”，根据基于改进正交化次序的模型结构劣汰策略逐个删除对实际输出贡献相对较小的结构项，以系统完备性指标为约束，确认结构与参数。二者结合即为结构辨识与参数辨识的一体化方法，此方法避免了正交最小二乘法和前向逐步回归法在非线性系统辨识过程会漏选关联性结构项的问题。为方便对比分析，本文称该一体化方法为MTGS (Modified Total Gram-Schmidt)算法。

1 线参数多项式组合模型描述

普适性的非线性系统有一大类可以用线参数多项式组合模型^[15]来描述。这类模型通用表达式为

式中：f(x)为模型输出, x为训练样本的输入向量；θ₀为系统广义截距；θ_i(i=1, 2, …, m)为模型系数；ϕ_i(x)(i=1, 2, …, m)为多项式结构项；e(x)为干扰噪声。

若ϕ_i(x)∈u(k-j)，j=1, 2, …, m，则模型为脉冲响应模型，u为与输入有关的结构项，k为采样时刻；若ϕ_i(x)∈{u(k-i), y(k-j)}，i=0, 1, …, m₁，j=1, 2, …, m₂，则模型为AMAX模型，y为与输出有关的结构项；若ϕ_i(x)∈{y(k-j)}，j=1, 2, …, m₂，则模型为自回归(AR)模型；若ϕ_i(x)∈{u(j)，u(j)u(k)，u(j)u(k)u(l), …}，j, k, l=1, 2, …, n，则模型为Volterra级数模型；若ϕ_i(x)∈{u(j), y(j), u(j)u(k), u(j)y(k), y(j)y(k), …}，j, k=1, 2, …, n，则模型为非线性NARMAX模型；若ϕ_i(x)∈{K(x, x₁), K(x, x₂), …, K(x, x_l), b}，其中x₁, x₂, …, x_l为l个训练样本，K(·)为样本输入的非线性映射，则模型为基于核函数映射的非线性模型，b为常数；若ϕ_i(x)∈{(x, cos x, e^x, ln x, …)，[+, /, ∑·, …]}，其中(·)为对输入向量非线性描述，[·]为非线性运算符，则模型为非线性多项式组合模型。

在辨识过程中，首先对数据进行中心化处理，消除广义截距对结构项的共线性影响。根据实际N个采样点数据，式(1)转化为

式中：Y=[y(x₁), y(x₂), …, y(x_N)]^T；θ=[θ₁, θ₂, …, θ_m]^T；e=[e(x₁), e(x₂), …, e(x_N)]^T；Φ=[ϕ₁, ϕ₂, …, ϕ_m]_N×m, i=1, 2, …, N, N≥m。

任何系统在其寿命范围内，严格非时变情况是不存在的，因此多项式系数也是时变的。为有效观测或控制一个系统，本文以非线性时变多项式模型为研究对象，开展模型结构与参数的一体化辨识研究。

2 基于向量空间投影理论的正交化分析

非线性多项式模型结构辨识问题的关键在于如何从一个满覆盖的空间集合V^m中寻找一个次优且满足精度的子空间集合V^r，r≤m是有界但未知。

引理1^[17]  Y(V)为向量Y在空间V的投影，其与向量Y夹角为β，投影残差为δ_Y^V，且满足：。

定理1  假定由张成的满空间为V^m=span{ϕ₁, ϕ₂, …, ϕ_m}，任意的V^r=span{ϕ₁, ϕ₂, …, ϕ_r}，V^r-1=span{ϕ₁, ϕ₂, …, ϕ_i-1, ϕ_i+1, …, ϕ_r}。记向量Y在满空间的投影为Y(V^m)，在V^r空间与V^r-1空间的投影分别为Y(V^r)与Y(V^r-1)，投影残差分别为δ_Y^Vr与δ_Y^Vr-1，且分别满足δ_Y^Vr⊥V^r，δ_Y^Vr-1⊥V^r-1。记向量ϕ_i在V^r-1空间的投影为ϕ_i(V^r-1)，在V^r-1空间的投影残差δ_ϕi^Vr-1=ϕ_i-ϕ_i(V^r-1)，且满足δ_ϕi^Vr-1⊥V^r-1；则有如下结论：

1) Y(V^r)=Y(V^r-1)+δ_Y^Vr-1δ_ϕi^Vr-1。

2) ||Y(V^r)||₂=||Y(V^r-1)||₂+||δ_Y^Vr-1δ_ϕi^Vr-1||₂。

证明  结论1)的证明:

由V^r=V^r-1+span(δ_ϕi^Vr-1), 得Y(V^r)=Y(V^r-1+span(δ_ϕi^Vr-1))=Y(V^r-1)+Yδ_ϕi^Vr-1；又由Y(V^r-1)⊥span(δ_ϕi^Vr-1)，span(δ_ϕi^Vr-1)⊥V^r-1，V^r-1⊥δ_ϕi^Vr-1得Y(V^r-1)⊥δ_ϕi^Vr-1；即Yδ_ϕi^Vr-1=Yδ_ϕi^Vr-1-Y(V^r-1)·δ_ϕi^Vr-1=δ_Y^Vr-1δ_ϕi^Vr-1。故Y(V^r)=Y(V^r-1)+δ_Y^Vr-1δ_ϕi^Vr-1。

结论2)的证明:

由V^r-1⊥δ_ϕi^Vr-1，得||Y(V^r)||₂=||Y(V^r-1)||₂+||Yδ_ϕi^Vr-1||₂=||Y(V^r-1)||₂+||δ_Y^Vr-1δ_ϕi^Vr-1||₂。  证毕

定理1揭示了向量在空间的投影与残差在空间的投影之间的对应关系。由定理1可得以下2条推论：

推论1  δ_Y^Vr=δ_Y^Vr-1-δ_Y^Vr-1δ_ϕi^Vr-1。

推论2  ||δ_Y^Vr||₂=||δ_Y^Vr-1||₂(1-cos²(δ_ϕi^Vr-1, δ_Y^Vr-1))。

证明  推论1的证明:

由δ_Y^Vr=Y-Y(V^r)=Y-(Y(V^r-1)+Y·δ_ϕi^Vr-1)=δ_Y^Vr-1-Yδ_ϕi^Vr-1；又Y(V^r-1)⊥δ_ϕi^Vr-1得δ_Y^Vr=δ_Y^Vr-1-Yδ_ϕi^Vr-1+Y(V^r-1)δ_ϕi^Vr-1=δ_Y^Vr-1-δ_Y^Vr-1δ_ϕi^Vr-1。

推论2的证明:

由||δ_Y^Vr||₂=||δ_Y^Vr-1-δ_Y^Vr-1δ_ϕi^Vr-1||₂=(δ_Y^Vr-1)²-2δ_Y^Vr-1(δ_Y^Vr-1δ_ϕi^Vr-1)+(δ_Y^Vr-1δ_ϕi^Vr-1)²和Y-Y(V^r)=δ_Y^Vr, Y-Y(V^r-1)=δ_Y^Vr-1，得δ_Y^Vr-1=Y(V^r)+δ_Y^Vr-Y(V^r-1)=δ_Y^Vr-1δ_ϕi^Vr-1+δ_Y^Vr；又由δ_Y^Vr-1δ_ϕi^Vr-1∈V^r, δ_Y^Vr⊥V^r得δ_Y^Vr⊥δ_Y^Vr-1δ_ϕi^Vr-1；所以δ_Y^Vr-1(δ_Y^Vr-1δ_ϕi^Vr-1)=(δ_Y^Vr+δ_Y^Vr-1δ_ϕi^Vr-1)(δ_Y^Vr-1δ_ϕi^Vr-1)=(δ_Y^Vr-1δ_ϕi^Vr-1)²。故||δ_Y^Vr||₂=||δ_Y^Vr-1||₂-||δ_Y^Vr-1δ_ϕi^Vr-1||₂=||δ_Y^Vr-1||₂(1-||δ_Y^Vr-1δ_ϕi^Vr-1||₂/||δ_Y^Vr-1||₂)=δ_Y^Vr-1||₂(1-cos²(δ_Y^Vr-1, δ_ϕi^Vr-1))。  证毕

推论描述了向量在空间的投影残差之间的递推关系。

正交化算法有2个功能^[18]：①将变量集合中的信息进行正交分解；②排除ϕ₁，ϕ₂，…，ϕ_m中的冗余变量(即被变换成0的变量)。由于这m个向量非正交，为求解模型系数θ，需对其正交化得到正交空间及向量Y在正交空间的投影。而误差就产生在上述过程中。

由推论可知，正交化次数越多，正交化向量范数越小，舍入误差对正交化进程的影响越大。因此CGS算法^[12]每一次迭代都是原始向量ϕ_i相对前面已正交化小向量的正交化，即“大向量面向小向量正交化”的过程，故误差较大，如图 1所示，其中x、x′、x″为正交空间的正交向量。

MGS算法^[12]由于每一次迭代选择的向量是ϕ_i^k而非原始向量。而ϕ_i^k的获得是基于原始向量已正交化后的向量，可理解为“小向量面向小向量正交化”的过程，故误差较小，如图 2所示。

综上，正交化过程选择不同，直接导致计算误差的不同传递，造成计算结果变差甚至错误，此现象在大矩阵计算中尤为明显。以上2种算法未对正交化向量选择次序进行处理，计算中仍可能出现“大向量面向小向量正交化”的过程。考虑到正交化次序对输出向量残差化的影响，本文提出非线性多项式模型结构和参数一体化辨识方法。

3 多项式模型的结构和参数一体化辨识 3.1 模型结构优选分析与停止条件

由推论1和2可知：①随着正交化次数的增多，输出向量的投影残差模越来越小，最终趋于噪声化；②在正交化过程中，输出向量残差模在某些项时其二阶范数迅速递减。分析可知，下降较大时对应为模型结构项。若正交化过程中执行一种残差监视线程，则可区分这些结构项，实现多项式模型的结构辨识。

改善这种投影有2种思路：①在每一步投影时，向与输出残差相关度最大的向量投影；②采用向二阶范数最大的基向量投影。思路①的风险在于：若范数小，正交化残差偏差较大；思路②的风险在于：若范数大，但有可能不是结构项，且增加了正交化次数。因此，本文折中选择某一项使得输出残差向量范数下降最快，这样迭代过程就较快实现输出向量残差的噪声化。

针对线参数模型式(2)，根据QR分解有e_N×1=Y-Φθ=Y-QRθ=Y-QG，其中：Q=[q₁, q₂, …, q_m]_N×m，R_m×m为对角线为1的上三角矩阵。数学上为方便求解，设D=Q^TQ=diag[d₁, d₂, …, d_m]，经过对输出向量Y的处理可得m维实向量，即G=D^-1Q^TY；然后针对Rθ=G回代计算求解。对Y=QG+e_N×1两边取二阶范数的平方得

将式(3)变形可得

因此，若在选项中选择对应最大的项，将得到较大的贡献项，这就构成了RMGS算法的结构优选策略。因此，基于输出向量残差化的结构停止条件归结为：

1) 监视残差比，设置阈值，当残差比小于阈值时停止。

2) 监视残差化输出，同样通过设置阈值项进行判别。

3.2 基于RMGS算法的模型结构优选策略

在MGS算法基础上提出RMGS算法，即采用选择主元策略，调整正交化次序，使得正交化过程的误差影响最小。

原先的正交化待选项为Φ=[ϕ₁, ϕ₂, …, ϕ_m]和Q=[q₁, q₂, …, q_m]，对应的解为θ，原先解的序号对应为λ(0, :)=[1, 2, …, m]。

记重新选择的正交化矩阵为Φ_new=[ϕ_{1, new}, ϕ_{2, new}, …, ϕ_{m, new}]和Q_new=[q_{1, new}, q_{2, new}, …, q_{m, new}]，对应的解为θ_new，选主元后的序号为λ(m, 1), λ(m, 2), …, λ(m, m)。

RMGS算法的迭代流程如下：

步骤1  记ϕ_j⁰=ϕ_j, j=1, 2, …, m；记Y⁽⁰⁾=Y；记q₁ⁱ=ϕ_i⁰, i=1, 2, …, m。计算：d₁ⁱ= < q₁ⁱ, q₁ⁱ>；。

假设d₁^i₁=max{γ₁ⁱ, 1≤i≤m}，记q_{1, new}=q₁^i₁并作为Q_new的第1列。分别计算g_{1, new}=；λ(m, 1)=λ(0, i₁)；θ_{1, new}=θ_i1。以g_{1, new}作为G_new第1个元素，以ϕ_i1作为Φ_new第1列，以θ_i1作为θ_new的第1个元素。

执行一个提取过程：for i=1:k-1，有R_{i, 1}^new⇐α^new_{i, i1}。即将前面i=1~(k-1)步骤中得到的α_i^new第i₁列元素作为R^new的第1列元素，以确保回归迭代过程中结构向量与系数是对应的。

计算且j∉{i₁}；更新特征矩阵：ϕ_j¹=ϕ_j⁰-α_{1, j}q_{1, new}，j=1, 2, …, m且j∉{i₁}；更新输出向量：Y⁽¹⁾=Y⁽⁰⁾-g_{1, new}q_{1, new}；计算残差比：。

步骤k  1 < k≤m：此时i=1, 2, …, m，j∉{i₁, i₂, …, i_(k-1)}；记q_kⁱ=ϕ_i^k-1，计算d_kⁱ= < q_kⁱ, q_kⁱ>，。

假设d_k^i_k=max{γ_kⁱ, k≤i≤m}，记q_{k, new}=q_k^i_k并作为Q_new的第k列。分别计算g_{k, new}=，λ(k)=i_k，θ_{k, new}=θ_ik。以g_{k, new}作为G_new第k个元素，以ϕ_ik作为Φ_new的第k列，以θ_ik作为θ_new的第k个元素。

执行一个提取过程：for i=1:k-1，有R_{i, k}^new⇐α^new_{i, ik}。即将步骤i(i=1, 2, …, k-1)中得到的α_i^new第i_k列元素作为R^new的第k列元素，确保回归迭代过程中结构向量与系数是对应的。

计算且j∉{i₁, i₂, …, i_(k-1)}；更新特征矩阵：ϕ_j^k=ϕ_j^k-1-α_{k, j}q_{k, new}，j=1, 2, …, m且j∉{i₁, i₂, …, i_(k-1)}；更新输出向量：Y^(k)=Y^(k-1)-g_{k, new}q_{k, new}；计算残差比：, k=k+1。

如果k≥m或者γ_k < ε，ε为阈值，则跳出循环，给出模型重要结构项系数的序号λ=[i₁, i₂, …, i_k]，进入步骤Last；否则迭代循环。

步骤Last  至此已得到G_new和R^new，其中：

对R^newθ_new=G_new回代，回代结果对应序号{λ(1), λ(2), …, λ(m-1)}的原始结构项系数。

3.3 基于改进正交化次序模型结构劣汰策略

由第2节可知，后期正交化过程会产生较大误差。在此过程中，选项个数仍记为m个，原始输入向量矩阵的重新排列为Φ_new=[ϕ_λ(1), ϕ_λ(2), …, ϕ_λ(m)]，其正交化矩阵Q_new=[q_{1, new}, q_{2, new}, …, q_{m, new}]中第i个向量q_{i, new}相对于q_{1, new}, q_{2, new}, …, q_{i-1, new}是正交的，但相对于后面向量的正交化程度是不均等的。为平等对待正交化向量，将结构项顺序重新排列，每一向量分别最后正交化，并计算其贡献，再根据贡献大小排序，实现“平等”劣汰过程。

针对ϕ_λ(i)，重新排列Φ_newⁱ=[ϕ_λ(1), ϕ_λ(2), …, ϕ_λ(i-1), ϕ_λ(i+1), …, ϕ_λ(m)]，记Q_newⁱ=[q_{1, new}, q_{2, new}, …, q_{i-1, new}, q′_{i+1, new}, …, q′_{m, new}, q″_{i, new}]，故ϕ_λ(i)相对于ϕ_λ(1), ϕ_λ(2), …, ϕ_λ(i-1), ϕ_λ(i+1), …, ϕ_λ(m)对描述Y的贡献为C(λ(i))=，对应贡献因子为。

对Q_newⁱ的求解可重新执行一次正交化过程，但会占用大量计算资源。考虑到上述迭代过程中已进行正交化，因此通过适当的变形即可求出Q_newⁱ。下面给出推理过程：

在已正交化过程中，以第i个向量q_{i, new}为例，其后续向量为q_{i+1, new}；进行重新排列后，得到的q′_{i+1, new}相对q_{i+1, new}少执行了一次针对q_{i, new}的正交化过程。因此：q′_{i+1, new}-·q_{i, new}=q_{i+1, new}。

依次类推，其后面的向量q′_{i+2, new}, q′_{i+3, new}, …, q′_{m, new}相对q_{i+2, new}, q_{i+3, new}, …, q_{m, new}都少执行了一次针对q_{i, new}的正交化过程。因此，q′_{j, new}=q_{i, new}+q_{j, new}，j=i+2, i+3, …, m。

同理，由于第i个向量在重排后放置于最后位置的q″_{i, new}，其相对于原始向量增加了针对q′_{i+1, new}, q′_{i+2, new}, …, q′_{m, new}的正交化过程。故可得其推理过程：

综上，上述推理方法可减少求解Q_newⁱ时，再次正交化的全部计算量。分别对每个向量重新排列，使每个候选项在正交化过程中受到的误差影响是一致的，然后在这种一致的基础上，删除贡献较小的项，谨慎地保护了结构项不易漏选。单纯的劣汰过程也可针对特征项的全集进行，但需要注意某些项在计算过程会出现 < q′_{j, new}, q′_{j, new}>接近于零，此时需将该项删除，避免数值计算病态问题。

4 仿真算例

为说明MTGS算法的有效性，以2个典型算例为研究对象，将MTGS算法与CGS算法、MGS算法进行仿真对比。给出了MTGS算法详细仿真步骤，限于篇幅对CGS算法和MGS算法只给出仿真对比结果。

4.1 算例1

非线性多项式模型——Volterra模型设定三阶脉冲响应输入输出表达式

式中:e(k)为白噪声干扰信号，噪信比为0.1；一、二、三阶时延分别为d₁=5、d₂=5和d₃=2。

在对上述模型进行辨识前，认为对系统的输入输出关系有一定了解，预先假定系统为一个阶次为3的非线性多项式模型，其中一阶时延为10，二阶时延为7，, 三阶时延为5，共含有C₁₀¹+(C₇²+C₇¹)+(C₅³+C₅¹C₄¹)+C₅¹=73个结构项(包含了系统式(6)的核心结构项)。此时问题转化为如何从这些结构项集合中选择结构并辨识参数。

由排列顺序可假定系统的结构项为：ϕ₁=u(k-1), …, ϕ₁₀=u(k-d₁), ϕ₁₁=u(k-1)², …, ϕ₁₇=u(k-1)u(k-d₂), …, ϕ₃₈=u(k-d₂)², …, ϕ₃₉=u(k-1)³, ϕ₄₀=u(k-1)²u(k-2), …, ϕ₇₃=u(k-d₃)³。

由采样点数N(N>73)，可计算73个结构项分别对应的N个采样数据，并组成矩阵x_ϕ1=[u(k), u(k+1), …, u(k+N-1)]^T，x_ϕ2=[u(k-1), u(k), …, u(k+N-2)]^T，…，x_ϕ73=[u(k-d₃+1)³, u(k-d₃+2), …, u(k-d₃+N)³]^T，Y=[y(k+1), y(k+2), …, y(k+N)]^T。根据排列可知：系统的实际结构项对应为ϕ₁, ϕ₂, ϕ₃, ϕ₅, ϕ₁₁, ϕ₁₂, ϕ₁₃, ϕ₁₉, ϕ₂₄, ϕ₂₆, ϕ₄₀。

辨识步骤如下：

步骤1  择优过程残差减小比率如图 3所示，共筛选出19项，图中序号1~19对应为2、1、24、3、40、11、12、13、26、19、5、28、21、66、39、17、34、68、8；基本上前面11项对应系统式(6)实际结构项，但由于择优过程需要保持一定的冗余度，故多选择了部分项。

步骤2   劣汰过程贡献因子变化如图 4所示。图 4(a)、(b)、(c)依次表示算法每一步迭代后淘汰贡献因子小的结构项的动态过程。可见每一项的贡献在正交化顺序不同时发生了变化，因此劣汰过程中采取的方式对每一个向量都是公平的；在劣汰过程，通过2步迭代就达到了满意度。系统值与3种算法辨识的结构项及对应参数值的对比如表 1所示。

设输入u(k)=sin (k/20)+sin (k/60)，k=1, 2, …, 500，设置初始时刻值为y(1)=y(2)=y(3)=0.01。图 5所示为实际模型输出与3种算法辨识模型输出的仿真对比。

4.2 算例2

非线性多项式模型——NARMAX模型系统的结构选择如下模型

式中: e₁(k)为输入信号强度的1/50随机信号。为了验证方法，预先假定系统为一个阶次为2的非线性多项式模型，其中u(k)与y(k)的时延均为3，共含有C₃₊₃¹+C₃₊₃¹+C₃₊₃²=27个结构项(包含了系统式(7)的核心结构项)。由排列顺序可假定系统的结构项为：ϕ₁=u(k-1), …, ϕ₃=u(k-3), ϕ₄=y(k-1), …, ϕ₆=y(k-3)，…，ϕ₇=u(k-1)², …, ϕ₁₀=u(k-1)u(k-3), …, ϕ₁₃=u(k-3)², …, ϕ₁₄=y(k-1)², …, ϕ₂₇=y(k-3)²。根据采样点数N(N>27)，可计算27个结构项分别对应的N个采样数据，原理同算例1，不在赘述。此时问题转化为如何从这些结构项集合中选择结构并辨识参数。

步骤1  择优过程残差减小比率如图 6所示，共筛选出7项，图中选项次序1~7分别为5、11、4、1、10、2、22项，包含了系统结构项。

步骤2  劣汰过程贡献因子变化如图 7所示。2个分图依次表示算法每一步迭代淘汰贡献因子小的结构项的动态过程。删除了第22结构项，剩下前6项，即正确的结构项，同时得到结构项对应系数，见表 2。上述结构辨识，虽然结构辨识正确，但辨识的参数仍存在不足(如第1项)，见表 2。分析原因，笔者认为是采样数据奇异或误差的吸收造成的；重新对数据进行处理后产生第2次仿真数据，进行系统结构与参数辨识，其过程如下：

步骤1  择优过程残差减小比率如图 8所示。该步骤中选出了7项，图中选项次序1~7分别为4、11、5、1、10、2、18项，包含了系统结构项。

步骤2  劣汰过程贡献因子变化如图 9所示。图 9(a)、(b)依次表示算法每一步迭代淘汰贡献因子小的结构项的动态过程。删除了第18项后，只剩下6项，即正确的结构项，对应参数如表 3所示。这个过程的结构辨识正确，参数偏差相对前次较小，可看出正交化次序的不同会对辨识结果有影响。

设输入为u(k)=0.001(rand-0.5)+0.01·(sin(k/60)+sin(k/20))，k=1, 2, …, 500，初始时刻值y(1)=y(2)=y(3)=0.000 1。图 10为实际模型与3种算法辨识模型的仿真输出对比。

5 结论

本文提出了一种非线性多项式模型结构和参数辨识方法，该方法包含基于贡献项的“择优过程”与基于冗余项的“劣汰过程”，研究了线参数多项式组合模型的结构与参数一体化辨识：

1) 采用RMGS算法可以从一个满覆盖的向量空间集合V^m中寻找出一个次优且满足精度需求的子空间集合V^r。

2) 基于改进正交化次序的模型结构劣汰策略，在包含核心结构项和少量冗余项的优选集中筛选出对实际输出贡献相对较小的结构项，以系统完备性指标为约束，确认了系统模型的结构与参数，稳妥地保护了随正交化次数增加可能会漏选的非显著模型项。

3) 对2个典型的非线性多项式模型进行仿真验证。无论从辨识模型的模型系数还是基于辨识模型的再次仿真对比结果，都说明MTGS算法比CGS算法和MGS算法更具有效性且计算时间短。

4) 需要说明的是，合理地定量选择阈值仍需进一步研究。若设置过大，则会将某些微弱贡献项漏选；若设置过小，则又可能无端增加迭代次数，增加计算时间。