☰
系统辨识三要素耦合与M序列试验设计:从数据到可用模型的工程实践
2026/10/29 12:33:20 网站建设 项目流程

简介:这份系统辨识教案面向自动化专业学生与相关课程学习者,围绕从输入输出数据中建立系统模型这一核心问题,梳理辨识理论、实验设计与编程计算的学习路径。压缩包内为1个doc文档,约1.42MB,内容涵盖课程参考书目、学习要求与考核办法,并系统讲解实体与模型的概念、模型分类及数学模型的获取方法。文档以Zadeh提出的辨识定义切入,展开数据、模型类与准则三要素,结合长网造纸过程的质量控制实例,完整演示明确辨识目的、收集先验知识、设计试验、数据采集与预处理、结构辨识与参数估计、模型检验及转换应用等步骤,便于读者建立从理论到工程实践的完整认知。目前已有150人学习,适合希望减少课堂依赖、通过自学掌握辨识技术并完成数据分析报告的学习者参考。

1. 系统辨识:为什么“给数据就能建模”这件事没那么简单

很多人第一次听到系统辨识,直觉反应是“不就是拿输入输出数据拟合一条曲线吗”。真上手就会发现,同一组数据,换一个准则、换一个采样间隔、换一个模型类,拟合出来的东西可能天差地别,甚至阶次越高残差越小、模型却越不可用。这背后不是玄学,而是系统辨识从一开始就把“数据、模型类、准则”三件事绑在一起定义:Zadeh 在 1962 年给出的经典表述是,辨识是在输入输出数据基础上,从一组给定模型类中确定一个与所观测系统等价的模型。等价与否由准则裁决,而准则又由辨识目的决定。这份自动化专业的《系统辨识与建模》教案,价值就在于它没有停在公式罗列,而是把“目的如何反向决定实验设计、模型结构和误差定义”这条链路完整摊开,让人看清数据拟合方法为什么不能脱离三要素单独比优劣。

2. 三要素的耦合:为什么同一组数据能拟合出完全不同的模型

2.1 数据不是中立的:观测时间、手段、目的都在改写它

教案里对“数据”的定义很克制:由观测实体而得,不唯一,受观测时间、观测目的、观测手段影响。这句话看着平淡,实际是整门课的第一道分水岭。同一个被控对象,你在负荷稳定时采一段、在启停阶段采一段,得到的动态信息完全不同;你用 1 秒采样和用 10 秒采样,能辨识出的最高阶次也完全不同。所以数据本身不是“客观事实”,而是被实验设计塑造过的观测切片。

这也解释了为什么教案把“设计辨识试验”放在“数据采集”之前,而不是反过来。试验设计要回答的是:系统边界划到哪里、哪些输入叠加扰动、哪些输入保持稳定、扰动信号频带多宽、采样频率多高、试验持续几个周期。这些决定一旦定下,后面能辨识出什么模型基本就被框住了。常见的翻车是先把现场数据抓回来再想怎么建模,结果发现激励不充分、噪声淹没信号、采样间隔和过渡过程严重不匹配,只能返工重采。

2.2 模型类是假设空间,不是答案本身

模型类规定了模型的形式,同样不唯一,受辨识目的和辨识方法影响。教案把数学模型细分成图表模型、解析模型、程序模型、语言模型,又把动态模型按集中/分布、序列/方程、线性/非线性、时变/时不变、确定/随机、单变量/多变量、离散/连续、有时延/无时延逐层切开。这一串分类不是名词堆砌,而是在告诉你:选模型类本质上是在选假设空间的大小和形状。

假设空间太小,真实系统不在里面,怎么估都欠拟合;假设空间太大,参数多到数据撑不住,就会过估计。教案里有一句很关键的判断:模型阶次越高,准则函数 J 值越小,导致模型过估计。这就是为什么不能拿 J 值单独当模型好坏的标准——它天然偏爱复杂模型。真正要做的,是在辨识目的约束下选一个“够用且可分析”的模型类,比如为了设计常规控制器,需要零极点、时延和动态响应准确,且模型便于分析;为了监视过程参数做故障诊断,则要求参数准确、模型能直观体现被监视量。目的不同,模型类的取舍就不同。

2.3 准则定义“等价”,而等价是相对的

准则规定了模型与实体等价的评判标准,同样不唯一。教案把误差分成参数误差、方程误差、输出误差、白化误差加权准则几类,差别不只是公式长相,而是优化问题的性质完全不同。

方程误差定义为 e(k)=y(k)+a·y(k-1)-b·u(k-1)-c·u(k-2),J 是参数的线性函数,好处是能直接套最小二乘,坏处是只能作一步预测,适合 AY=B1U1+B2U2 这类模型。输出误差则是 e(k)=y(k)-ŷ(k),ŷ(k) 由模型递推得到,J 变成参数的非线性函数,但输入序列已知时可以作多步预测,适合 B/A 结构。这两者的选择不是口味问题:你要一步预测还是多步预测,决定了你该用哪个准则。

更隐蔽的是采样间隔对准则的影响。教案明确指出,采样间隔增大一倍,模型阶次大致减少一半;要做不同采样步长下模型的比较,误差公式必须改成 J=∑e(k)ᵀW e(k)·T₀²/N,T₀ 为采样步长。不加这个 T₀² 因子,拿不同采样率下的 J 值直接比大小,结论就是错的。这类细节在教材里往往一笔带过,但恰恰是实际报告里最容易出错的地方。

提示:三要素只有在相同的前提下才可比较方法优劣。换数据、换模型类、换准则之后,原来的“最优方法”结论会失效,这不是方法退化了,而是比较基准变了。

3. 辨识试验设计:M 序列为什么能当“可控的白噪声”用

3.1 白噪声的理想性与物理不可实现性

教案对白噪声的定义很干净:零均值、平稳、谱密度为常数的随机过程。它有几个对辨识极其友好的性质:作为系统输入时,输入与输出的互相关函数正比于系统的单位脉冲响应;能激发系统所有模态,实现充分激励;能防止数据病态,保证辨识精度;还能用输出估计误差是否具有白色性来判断辨识方法和模型结构是否合适。

但教案同时点破了一个现实:工业生产中白噪声在物理上是不存在的,常见的是有色噪声。有色噪声的表示定理说,任何平稳噪声序列都可以看作白噪声经过一个渐近稳定的线性环节(成形滤波器)的输出。这条定理的意义在于,它把“有色”问题转化成了“白噪声 + 已知或可估的滤波器”问题,后面的预测误差类方法基本都建立在这个视角上。

3.2 M 序列的构造与自相关特性

既然理想白噪声造不出来,工程上就需要一个“近似白”的可控替代品,M 序列(二位式最大长度伪随机序列)就是教案重点展开的对象。它由 n 级移位寄存器产生,周期 N=2ⁿ-1,是一个周期内能表示的最多非零状态数。教案给的 4 阶例子是 1 1 1 1 0 1 0 1 1 0 0 1 0 0 0 1…,把 0 换成 -1 后,自相关函数呈现很规整的形状:

R_M(τ) = 1 τ = i·N, i = 0,1,2,… R_M(τ) = -1/N τ = i·N + k, k = 1,2,…,N-1

当 N 充分大时,这个 -1/N 的旁瓣趋近于 0,自相关函数近似为 δ 函数。于是以 M 序列为输入的线性系统,其输入输出互相关函数序列就近似等于脉冲响应序列——这就是“用相关分析法直接测脉冲响应”的底层依据。教案还给了频谱层面的解释:M 序列的频谱是线条谱而非光滑曲线,直流分量与 N² 成正比,加大 N 可减少直流分量;频带约为 1/(N·Δt),减小步宽 Δt 可增加带宽;谱线密度与 N 成正比。

3.3 多输入正交化与逆重复 M 序列

多输入单输出系统有个麻烦:几个输入同时扰动,输出里混着各输入的影响,怎么分离?教案给的方案是用同一 M 序列的不同移位序列作为各输入的扰动,比如 {M(k)}、{M(k+J)}、{M(k+2J)}。当输出对各输入的过渡过程时间都小于 J 时,这些移位序列之间可认为互相正交,从而把多输入问题拆成多个单输入问题分别处理。这个技巧的代价是试验时间变长,因为移位间隔 J 必须大于最长过渡过程。

逆重复 M 序列则是把 M 序列与方波序列相乘,周期变成 2N,一个周期内 1 和 -1 的个数相等,自相关函数在 τ=i·2N 时为 1、τ=i·2N+1 时为 -1、其余为 0。它保留了白噪声式的激励优点,同时幅值完全可控,均值严格为零,对含有直流漂移的现场更友好。

3.4 参数怎么定:n 和 Ts 的两套设计路径

教案给了两套 M 序列参数设计方法,都围绕系统过渡过程时间 T 和最高工作频率 f_max 展开,通常取 f_max=(3~15)/T。

第一套从频带约束出发:先由 1/Ts>3·f_max 得 Ts=T/(3k)(k 为 3~15 中的取值),再由 1.5T<(2ⁿ-1)·Ts 反解 n,近似得到 4.5k≈2ⁿ-1。第二套更直接:先定 n,再由 n·Ts=(0.1~0.8)T 算 Ts。两套方法的共同约束是:扰动信号频带要宽于过程工作频带,持续时间要覆盖 3~5 个扰动周期,幅值由安全工况决定。对于多输入系统,还要求 (2ⁿ-1-1)·Ts>1.5·∑(输出对各输入的过渡过程时间),保证移位序列之间真正错开。

采样间隔则单独由香农定理约束:1/T₀>2·f_max;有扰动信号时还要求 T₀<Ts 且 Ts 是 T₀ 的整数倍。这条“整数倍”约束经常被忽略,导致重抽样时相位对不齐,相关分析结果出现偏移。

注意:M 序列幅值不是越大越好。教案明确要求幅度以不引起断纸等生产异常为限,安全工况优先于激励充分性,这是工程辨识和实验室辨识最本质的区别。

4. 从原始数据到可辨识模型:预处理与模型转换的隐藏工作量

4.1 数据预处理四件事:校正、滤波、重抽样、标准化

教案把预处理拆成四步,每一步都有明确的失败模式。传感器非线性校正依据厂家静态非线性曲线编程修正;滤波要求对输入输出加相同的滤波器,这样输入输出关系不变,若只滤一边就会引入额外动态;重抽样针对过渡过程较长的子系统,用重新选采样间隔的办法降低模型阶次;标准化做零均值、单位方差,目的是降低舍入误差、防止病态方程、提高计算精度。

标准化的推导值得细看。对模型 y(k)+a·y(k-1)=b·u(k-1)+c·u(k-2)+d 先做零均值化,得到 Y(k)+a·Y(k-1)=b·U(k-1)+c·U(k-2),再除以标准差 σ_y、σ_u,得到 Yσ(k)+a·Yσ(k-1)=B·Uσ(k-1)+C·Uσ(k-2),其中 B=b·σ_u/σ_y、C=c·σ_u/σ_y。辨识完成后必须用反演公式 b=B·σ_y/σ_u、c=C·σ_y/σ_u 还原真实参数。漏掉反演这一步,报告里的增益会整体偏一个 σ 比例,而且因为量纲被抹掉,肉眼很难发现。

4.2 参考模型到应用模型的转换

教案把模型转换单列一步,是因为结构辨识和参数估计得到的“参考模型”往往是为了算法方便而选的参数化形式,不一定适合直接使用。比如辨识得到的是离散差分方程,但控制器设计需要连续传递函数;或者得到的是带特定噪声结构的 ARMAX 形式,而仿真训练系统只关心动态响应。转换过程要保证零极点、时延、直流增益这些关键特征不丢,否则前面所有精度都白费。

4.3 模型检验不通过时往哪返

教案的流程图里,模型检验不通过要“分析问题所在,并返回到前期各相应步骤”,应用评价不成功同样要回退。这句话实际操作起来最考验判断力:残差有偏可能是模型类选错,也可能是数据预处理滤波不当;参数置信区间过宽可能是激励不充分,也可能是采样间隔与过渡过程不匹配。教案没有给万能诊断表,但把“返回到相应步骤”而不是“一律重来”写进流程,本身就是一种工程态度——定位到哪一层出错,就只回退到那一层,避免全流程推倒重来。

5. 避坑与排查:五条来自教案细节的踩坑记录

5.1 现象:阶次越加越高,J 值一直降,模型却不能用

原因:准则函数 J 对模型复杂度天然偏袒,阶次越高拟合残差越小,这是教案明确指出的“过估计”现象。解决:把模型检验和应用评价作为独立关卡,用输出误差的多步预测表现、零极点物理合理性、时延是否符合先验知识来交叉判断,而不是只看 J 值。

5.2 现象:不同采样率下比较两个模型,结论完全相反

原因:误差公式没有按采样步长归一化。解决:跨采样步长比较时统一用 J=∑e(k)ᵀW e(k)·T₀²/N,把 T₀² 因子带上,否则采样密的那组天然占便宜。

5.3 现象:多输入系统辨识出的通道互相“串扰”

原因:各输入扰动信号没有做移位正交化,或者移位间隔 J 小于最长过渡过程时间。解决:按 (2ⁿ-1-1)·Ts>1.5·∑T_i 校验移位间隔,确保输出对各输入的过渡过程都短于 J。

5.4 现象:标准化之后参数还原不回去,增益整体偏移

原因:只做了正向标准化,忘了反演公式 b=B·σ_y/σ_u、c=C·σ_y/σ_u。解决:把反演作为标准化流程的固定收尾步骤,并在报告里同时记录标准化前后的参数对照。

5.5 现象:滤波之后输入输出关系变了,辨识出虚假动态

原因:只对输出滤波或对输入输出用了不同滤波器。解决:严格对输入输出施加相同滤波器,保持输入输出关系不变;若必须单边滤波,需在模型类里显式补偿滤波器动态。

6. 设计取舍:为什么“够用的模型”比“精确的模型”更难选

系统辨识里最难的从来不是推最小二乘公式,而是在目的、数据、模型类、准则之间做取舍。教案把辨识目的列在最前面,并且逐条对应模型要求:验证理论模型要零极点、结构、参数都准,模型类还得同理论模型一致;设计常规控制器要动态响应、零极点、时延准确,且模型便于分析;设计数字控制器要便于计算机运算;预报预测要动态响应和时延准确;监视过程参数做故障诊断只要参数准确、模型直观。这些要求彼此并不兼容——追求参数无偏可能牺牲多步预测能力,追求结构简洁可能丢掉高频动态。

我一般会这么理解这套取舍:辨识不是逼近真实系统的全部,而是在给定观测预算和安全约束下,构造一个在特定用途上“等价”的简化描述。教案里造纸过程的例子把这点体现得很清楚:车速因同步困难被设为恒定,流量到定量的延迟约 60 秒而响应只有约 2 秒,浓度到定量的延迟约 120 秒而响应约 80 秒,蒸汽压力到水份的延迟约 45 秒而响应约 60 秒。这些先验知识直接决定了哪些输入叠加扰动、采样频率取多少、试验做多长。先验知识越扎实,假设空间就能收得越紧,需要的试验数据越少,模型也越稳。

验证一套辨识结果是否可信,我习惯做三件事:一是看残差是否接近白色,白色性说明模型把可解释的动态基本吃干净了;二是看参数置信区间和物理量纲是否合理,增益符号反了、时间常数数量级离谱,基本就是模型类或预处理出了问题;三是拿模型去做它本来要服务的那个任务,控制器设计就看闭环响应,故障诊断就看参数漂移能不能被稳定跟踪。教案把“应用评价”放在流程末端并允许回退,本质上就是承认模型的价值只能由用途来裁决。

这份教案真正值得细读的地方,是它把“为什么这么设计”写进了每一步的约束条件里:M 序列的 n 和 Ts 不是拍脑袋定的,是被过渡过程和安全工况夹出来的;标准化不是形式主义,是为了防病态方程;模型转换不是多余环节,是因为参考模型和应用模型服务的目标不同。从那以后我每次做辨识,都强制先把辨识目的和验收准则写成一页纸,再动数据采集设备——目的没写清就开测,后面返工的概率高得离谱。希望帮到你。

本文还有配套的精品资源,点击获取

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询