当前位置：首页 > news >正文

信号与系统/控制理论必备：手把手教你搞定拉普拉斯反变换中的部分分式展开（含MATLAB代码）

news 2026/5/28 1:36:42

信号与系统实战：拉普拉斯反变换中的部分分式展开技术精解

在分析RLC电路瞬态响应时，我们常会遇到这样的场景：当输入一个阶跃电压后，电路输出端的电压表达式在复频域中呈现为复杂的有理分式。此时，部分分式展开就像一把精巧的瑞士军刀，能够将这个"缠绕在一起"的表达式分解成多个简单分式的和，每个分式对应着系统的一个自然响应模式。这种分解不仅揭示了系统的内在特性，更为时域分析打开了方便之门。

1. 工程视角下的部分分式展开原理

部分分式展开的核心思想，是将一个复杂的有理分式分解为若干简单分式的线性组合。这种分解之所以重要，是因为每个简单分式都对应着一个基本的时域函数。当我们进行拉普拉斯反变换时，这些简单分式能够直接转换为指数函数、正弦函数等基本时域信号。

考虑一个典型的二阶系统传递函数：

$$ H(s) = \frac{3s+5}{s^2+4s+13} $$

这个传递函数的分母可以因式分解为$(s+2)^2+3^2$，对应着一对共轭复根。通过部分分式展开，我们可以将其表示为：

num = [3 5]; % 分子多项式系数 den = [1 4 13]; % 分母多项式系数 [r,p,k] = residue(num,den)

执行这段MATLAB代码后，我们会得到：

r = 1.5000 - 0.1667i 1.5000 + 0.1667i p = -2.0000 + 3.0000i -2.0000 - 3.0000i k = []

这表明系统由两个共轭的极点组成，每个极点对应着一个特定的响应模式。在时域中，这将转换为一个衰减的正弦振荡。

2. 三类典型情况的处理方法

2.1 单实数极点情况

单实数极点是最简单的情况，对应着系统的指数响应模式。例如，RC电路的传递函数通常具有这种形式：

$$ H(s) = \frac{1}{RCs+1} = \frac{1/RC}{s+1/RC} $$

在MATLAB中处理这类问题非常直接：

% 单极点示例 num = [1]; den = [3 1]; % 对应时间常数τ=3 [r,p,k] = residue(num,den)

输出结果为：

r = 0.3333 p = -0.3333 k = []

这表示系统的时间常数为3秒，反变换后得到时域函数为$\frac{1}{3}e^{-t/3}$。

2.2 重极点情况

重极点常见于包含多个相同储能元件的系统，如两个相同电容组成的滤波电路。考虑以下传递函数：

$$ H(s) = \frac{s+3}{(s+1)^3}

对应的MATLAB处理： ```matlab num = [1 3]; den = conv([1 1],[1 1]); % (s+1)^2 den = conv(den,[1 1]); % (s+1)^3 [r,p,k] = residue(num,den)

输出将显示三个系数，分别对应$(s+1)^3$、$(s+1)^2$和$(s+1)$的项。在时域中，这将产生$t^2e^{-t}$、$te^{-t}$和$e^{-t}$的组合响应。

2.3 共轭复极点情况

机械振动系统、LC电路等常会出现共轭复极点，对应着系统的振荡响应。例如：

$$ H(s) = \frac{2s+5}{s^2+2s+5}

使用MATLAB分析： ```matlab num = [2 5]; den = [1 2 5]; [r,p,k] = residue(num,den)

输出结果展示了共轭的极点和留数：

r = 1.0000 - 0.5000i 1.0000 + 0.5000i p = -1.0000 + 2.0000i -1.0000 - 2.0000i k = []

这对应着时域中的衰减振荡$e^{-t}(2\cos2t + \sin2t)$。

3. MATLAB实战：从理论到工程应用

3.1 residue函数的高级用法

residue函数不仅能处理标准的展开情况，还能应对各种工程实际问题。例如，当分子次数高于分母时（物理系统中较少见但数学上可能）：

num = [2 3 1 2 0]; % 2s^4 + 3s^3 + s^2 + 2s den = [1 4 3]; % s^2 + 4s + 3 [r,p,k] = residue(num,den)

输出中的k将包含多项式部分的系数：

k = 2 -5 15

这对应于展开式中的$2s^2 -5s +15$多项式项。

3.2 数值稳定性问题处理

当处理高阶系统时，可能会遇到数值不稳定的情况。这时可以采用以下策略：

增加计算精度：

format long num = [1e-6 1]; den = [1 1e6 1e6]; [r,p,k] = residue(num,den)

使用符号计算：

syms s F = (s+1)/(s^2 + 1e6*s + 1e6); partfrac(F)

系统降阶：对于极高阶系统，可先进行模型降阶处理。

4. 工程案例分析：多极点系统的时域响应预测

考虑一个实际工程中的四阶系统：

$$ H(s) = \frac{10(s+2)}{(s+1)(s+3)(s^2+2s+5)}

使用MATLAB进行完整分析： ```matlab % 定义系统 num = [10 20]; den1 = [1 1]; den2 = [1 3]; den3 = [1 2 5]; den = conv(conv(den1,den2),den3); % 部分分式展开 [r,p,k] = residue(num,den); % 绘制极点分布 figure; plot(real(p),imag(p),'rx','MarkerSize',10,'LineWidth',2); grid on; xlabel('实部'); ylabel('虚部'); title('系统极点分布'); hold on; plot(0,0,'bo'); % 原点 axis equal;

这个分析揭示了系统包含：