2. 一维系统数学描述

要对一维系统做参数辨识与校准,第一步不是急着选算法,而是把“系统到底用什么数学语言描述”这件事定下来。同一台设备、同一个热工或机电对象,可以用常微分方程、传递函数、状态空间和离散差分模型四种方式表达,它们描述的是同一个物理对象,只是观察角度和适用场景不同。辨识算法、采样方式、参数含义都直接取决于你选了哪种描述,因此这一节要先把四种模型讲清楚,并明确输入输出变量与参数各代表什么。

一、常微分方程:最贴近物理机理的描述

一维系统的常微分方程(ODE)通常写成输入 u(t) 与输出 y(t) 及其各阶导数之间的关系。以一阶惯性系统为例:

τ · dy(t)/dt + y(t) = K · u(t)

其中 τ 是时间常数(单位 s),反映输出跟随输入的快慢;K 是静态增益,反映稳态时输出与输入的比值。若系统含二阶动态,则会出现 d²y/dt²,并引入阻尼比 ζ 与自然频率 ωn。ODE 的优点是物理意义明确,参数可直接对应到质量、阻尼、热阻、热容等真实量;缺点是求解与辨识都需要处理导数,对噪声敏感。

二、传递函数:把微分方程变成代数关系

在零初始条件下对 ODE 做拉普拉斯变换,就得到传递函数 G(s) = Y(s)/U(s)。一阶惯性系统对应:

G(s) = K / (τ·s + 1)

二阶系统对应 G(s) = K·ωn² / (s² + 2ζωn·s + ωn²)。传递函数把时域的微分运算变成了复频域的代数运算,便于分析极点位置、稳定性、频率响应。对辨识而言,传递函数形式直接决定了待估参数集合:一阶只需估 K 和 τ,二阶则需估 K、ζ、ωn。需要注意的是,传递函数默认系统线性、时不变、零初始条件,实际对象若存在明显非线性或大初始偏差,直接套用会带来系统误差。

三、状态空间:把高阶方程拆成一阶方程组

状态空间模型用一组一阶微分方程描述系统内部状态 x 与输入输出之间的关系:

dx/dt = A·x + B·u
y     = C·x + D·u

对一维二阶系统,状态可取 x1 = y、x2 = dy/dt,则 A、B、C、D 为相应维度的矩阵。状态空间的优势是能同时描述内部状态和外部输出,便于处理多输入多输出、初始条件非零、以及现代控制与卡尔曼滤波类辨识方法。对一维系统而言,状态空间与传递函数可以互相转换,参数含义一致,只是表达形式不同。

四、离散差分模型:面向数字采集与计算机辨识

实际辨识几乎都在计算机上完成,采集到的是离散序列 u(k)、y(k),采样周期为 Ts。此时更自然的是差分方程,例如一阶 ARX 形式:

y(k) = a·y(k-1) + b·u(k-1) + e(k)

其中 a、b 是离散参数,e(k) 是残差。离散参数与连续参数之间存在映射关系,例如一阶系统在零阶保持下 a = exp(-Ts/τ)、b = K·(1 - a)。这意味着辨识出 a、b 后,可以反推 τ 和 K。离散模型直接对应最小二乘等算法,是参数辨识最常用的落地形式。

一维系统四种数学描述的等价关系与参数映射 常微分方程 ODE τ·dy/dt + y = K·u 参数:τ (s)、K 物理意义最直接 传递函数 G(s) G(s)=K/(τs+1) 参数:K、τ(或 ζ、ωn) 频域分析、极点判稳 状态空间 ẋ=Ax+Bu, y=Cx+Du 参数:A、B、C、D 内部状态、多变量 拉氏变换 零初始条件 选状态变量 x1=y, x2=ẏ 离散差分模型 y(k)=a·y(k-1)+b·u(k-1) 参数:a、b(采样 Ts) 直接用于最小二乘辨识 a=exp(-Ts/τ), b=K(1-a) ZOH 离散化 输入 / 输出变量 u(t):激励(电压/流量/功率) y(t):响应(温度/转速/位移) Ts:采样周期 采集为序列 辨识落地选择 连续参数 → 物理校准 离散参数 → 算法辨识 两者可互相换算 参数映射

五、四种描述之间的转换与选择

四种模型并非彼此独立。ODE 经拉普拉斯变换得传递函数,传递函数经状态选取得状态空间,连续模型经零阶保持离散化得差分模型;反过来,差分模型辨识出 a、b 后,也可反算 τ、K 回到物理参数。选择哪种描述,取决于任务:若目标是理解物理机理、做参数校准,优先用 ODE 或传递函数;若要做数字辨识、在线估计,优先用差分模型;若系统有内部状态需要观测或涉及多变量,则用状态空间。

六、参数含义与辨识对应关系

模型形式典型表达式待辨识参数参数物理含义适用辨识方法
常微分方程τ·dy/dt + y = K·uτ、K时间常数、静态增益非线性最小二乘、阶跃响应法
传递函数G(s)=K/(τs+1)K、τ(或 ζ、ωn)增益、惯性、阻尼、频率频域拟合、最小二乘
状态空间ẋ=Ax+Bu, y=Cx+DuA、B、C、D状态转移、输入输出耦合子空间辨识、卡尔曼滤波
离散差分y(k)=a·y(k-1)+b·u(k-1)a、b离散极点、输入增益最小二乘、递推最小二乘
教学示例(假设数据):设某一维温控对象真实参数为 τ=20 s、K=2.5,采样周期 Ts=1 s。按 ZOH 离散化可得 a=exp(-1/20)≈0.9512,b=2.5×(1-0.9512)≈0.1220。若用最小二乘从 u(k)、y(k) 序列辨识出 a≈0.95、b≈0.12,则可反算 τ≈-Ts/ln(a)≈19.5 s、K=b/(1-a)≈2.4,与真值接近,说明离散辨识结果可映射回连续物理参数。以上数值为示例假设,用于说明换算过程,并非实测结果。
注意:离散参数 a、b 与连续参数 τ、K 的换算依赖采样周期 Ts 和离散化方式(ZOH、双线性等)。若采样过快,a 接近 1,辨识矩阵易病态;若采样过慢,动态信息丢失,反算误差增大。一般建议每个时间常数采样 5~10 个点。