《Mathematical and Computational Applications》：A Unified Framework for Determining Post-Jump Initial Conditions for Biproper LTI Systems Subjected to Analytic Causal and Dirac Delta Inputs
编辑推荐：
确定跳变后初始条件对于受不连续输入作用的因果线性时不变（LTI）系统的精确分析与数值仿真至关重要。尽管传统时域方法在切换事件后立即需要这些初始条件，但尚无通用解析流程可直接从控制方程导出。本文提出一种统一解析框架，用于确定由常系数非齐次常微分方程（ODE）描述
确定跳变后初始条件对于受不连续输入作用的因果线性时不变（LTI）系统的精确分析与数值仿真至关重要。尽管传统时域方法在切换事件后立即需要这些初始条件，但尚无通用解析流程可直接从控制方程导出。本文提出一种统一解析框架，用于确定由常系数非齐次常微分方程（ODE）描述的因果单输入单输出LTI系统在已知跳变前初始条件且承受解析因果或狄拉克δ（Dirac delta）输入时的跳变后初始条件。研究人员发展了三种互补方法：基于常微分方程（ODE）的方法导出显式递归闭式表达式；基于状态空间模型（SSM）的方法建立等价向量形式；基于拉普拉斯变换初值定理（IVT）的方法不仅确定输出值在t=0+的值，还确定其高阶导数初始条件。三种方法被证明数学等价，若干算例验证了所提框架及其在LTI系统分析与数值仿真中的适用性。
研究背景方面，线性时不变（LTI）系统因可套用强有力分析工具而备受关注，其中真实物理因果系统由真分式（proper）传递函数表示，可分为严格真分式（strictly proper）与双真分式（biproper，即相对阶r=0）。当输入在t=0处施加且不连续时，含输入导数的双真分式系统输出及其导数会出现跳变，传统教材通常假定跳变后初始条件已知，未给出由控制方程、跳变前条件与因果输入直接推导完整跳变后初始条件向量{y(0+), y′(0+), …, y(n-1)(0+)}的通用步骤。分布论与经典跳变条件法虽严谨，但未组织成任意阶系统的递归闭式。为此，研究人员在《Mathematical and Computational Applications》发表本文，建立覆盖解析因果与狄拉克δ两类输入的跳变后初始化统一框架。
关键技术方法上，研究人员采用三条互补路线。其一为ODE逐次积分法：对n阶常系数非齐次ODE在t∈(-ε, ε)做n次分布意义下反微分，利用引理区分局部有界正则项（积分为零）与δ及其导数项（按阶贡献单位质量），对解析因果输入u(t)=f(t)·1(t)用莱布尼茨型导数公式分离Heaviside跳变产生的δ项，对狄拉克δ输入u(t)=δ(t)则直接分布配对，分别推得递归式(26)(41)。其二为状态空间法：取可观测实现，借状态向量x(t)在解析因果输入下连续（x(0+)=x(0-)）这一性质，由可观测矩阵O与马尔可夫参数矩阵M将ODE初值向量表为x(0+)与输入导数向量的线性组合(73)。其三为初值定理扩展法：将经典IVT推广至y(k)(0+)=lims→∞sk+1Y(s)-已确定低阶项，从单边拉氏域代数方程(80)按s大幂次配系数得同一递归关系。
研究结果各节要点如下。
第2节常微分方程公式化中，2.1款针对解析因果输入证得定理1：跳变Δy(j)(0)=y(j)(0+)-y(j)(0-)满足anΔy(j) = -Σi=0j-1aiΔy(i) + bjf(j)(0) + Σq=0j-1bq·0（经索引整理为式27），由此顺次解出y(j)(0+)。2.2款针对狄拉克δ输入证得定理2：双真分式系统拆为直接透传项bnδ(t)加严格真分式响应，递归式(41)(42)仅由b0…n-1与δ系数决定，推论1给出严格真分式令bn=0的特例。
第3节状态空间分析公式化提出定理3：解析因果输入下y(j)(0+) = cTAjx(0-) + Σq=0j djq u(q)(0+)，其中矩阵M在相似变换下不变（引理5），故跳变后初值属输入—输出系统本身属性而非实现依赖；注记明确本SSM公式不推广至δ输入。
第4节初值定理公式化证明定理4：若yreg(t)在0+邻域解析且指数阶，则y(k)(0+)=lims→∞sk+1Y(s)-Σr=0k-1sk-ry(r)(0+)，推论2将其用于严格真分式冲激响应。
第5节等价性中，定理5比对SSM式(108)与ODE式(27)经马尔可夫参数D, cTb, cTAb…与传递函数大s展开系数一致而等价；定理6比对IVT配系数与ODE逐次积分递归式同源。三部分在解析因果输入域完全等价，δ输入则仅ODE与IVT两路适用。
讨论与结论部分指出，所称“统一框架”意指三法同解一个初始化问题而非每法通吃所有输入类；解析因果输入下三者等价且可互验，δ输入限用ODE与IVT。研究人员结论为：所提递归闭式(27)(42)可直接由ODE系数、跳变前初值与输入在原点各阶导数算出完整跳变后初值向量，算法1给出计算流程，解决了双真分式及含输入导数严格真分式LTI系统在切换点继续时域求解的前置条件缺失问题，为电路、控制与数值仿真提供不依赖先猜后校的初始化依据。
