控制工程理论的创新性数学补充及创新性概念

李春明，尹晓丽，张晓玲

(1.中国石油大学(华东)机电工程学院，山东 青岛，266580；2.山东石油化工学院，山东 东营，257061)

在工业生产上，往往需要控制一些关键参数，如温度、压力、气体含量等,这属于机械控制工程[1-3]。机械控制工程已形成了较固定而完备的理论[4]，但是在一些计算公式的推导过程中，由于学生都学习过高等数学、线性代数等数学课程，往往忽略一些重要的数学推导步骤,这给初学者带来负担。在一般院校，尤其是职业教育院校，学生的数学基础往往不够扎实，难以理解被省略的数学推导步骤，导致该课程的学习难度高。本文研究理论性强的公式、系统建模及其必要的数学推导过程。

1 拉氏变换微积分步骤

由于微积分是高等数学的基本知识，专业著作[5]中往往省略关键步骤。

1.1 指数函数的拉氏变换

补充复合函数的定积分及以下求导结果：

则该式的拉氏变换为

1.2 一次幂函数的拉氏变换

补充如下分部积分的详细计算及求不定式的洛必达法则(L′Hopital′s rule)计算式：

则该式的拉氏变换等于以下两项：

上式的第一项为

则该式为

1.3 正弦函数和余弦函数的拉氏变换

补充复变函数的欧拉公式及复数定义：

正弦函数的拉氏变换推导过程为

余弦函数的拉氏变换推导过程为

2 代数运算

2.1 传递函数的部分分式展开法

该展开式是根据二次方程的求根公式和数学建模的待定系数法而获得的。将函数的分子和分母分别因式分解，二次项因子由求根公式获得两个一次项相乘的形式，假设该函数可展开为

式中ri为待定系数，称为留数。

留数由特殊情况确定，在s=pi时，可以将(s-pi)乘以等号两边，左边能算出来，右边为留数。由于等式仍然成立，所以可求得该待定系数。

分子的幂次如果大于分母的，则须用分母去除分子而得到合理的分式。如可补充多项式相除的竖式计算方式而获得：

该竖式为

2.2 劳斯表及其两个特例

计算二阶行列式可采用以下简易计算式：

劳斯表是判断系统稳定性的，有时会出现两个特殊情况：

1)在新行中，首列元素等于零。为避免被零除，该项须用很小的正数代替，如

2)新行的所有元素均为零。这导致以下所有元素均为零。将上一行所对应的多项式求导，由求导结果确定该新行。所依据的理论是洛必达法则。如

分别对应于多项式：

2.3 基于二次假设将系统微分方程转化为差分方程

设采样周期为T，用差分代替微分，根据后(左)向差分的定义及二次假设，变量x的一阶和二阶差分为[6-7]：

Δx(k)=0.5[3x(k)-4x(k-1)+x(k-2)]，

Δ2x(k)=x(k)-2x(k-1)+x(k-2)。

将以下微分方程离散化：

(m+1.5cT+kT2)x(k)+(-2m-2cT)x(k-1)+

(m+0.5c)x(k-2)=0。

经典的转化算法是基于两层线性假设的，必然没有上述基于二次假设的准确。

3 数学建模

3.1 弹簧-质点-干摩擦系统的数学模型

忽略最大静摩擦力大于动摩擦力的特性，图1所示系统可建立如下分段方程：

图1 弹簧-质点-摩擦振动系统模型

k[xi(t)-xo(t)]=0，

该方程不能用来建立传递函数。数值算法可处理分段形式的微分方程，而传递函数算法不能。

3.2 受到阻尼的转子系统模型

设转子受到阻尼力和外力矩M，则力矩平衡方程为

式中：J为转子的转动惯量；ω为转子的角速度；fM为黏性力矩阻尼系数，如果半径为r，阻尼系数为c，则该系数为cr2，目前将其定义为摩擦系数不确切。

4 结束语

S(科学)T(技术)E(工程背景)M(数学)教学理念引领教学内容的改进[8]。其中，数学基础往往是各专业课所忽视的。本文补充的数学推导表面上创新性平常，实际上使高难度知识点简易化，比一般的理论创新更有实践价值。本文的理论创新可直接沉淀到教材，因此，也是山东省本科教育改革研究重点项目(Z2020057)的研究内容。