跳至主要內容

数学基础

大约 15 分钟

数学基础

拉普拉斯变换

对于详细部分也可参考复变函数与积分变换笔记的傅里叶变换与拉普拉斯变换部分

拉普拉斯变换的定义

f(t)f(t) 为定义在 [0,+∞)[0,+\infty) 的实函数, s=δ+jωs=\delta+j\omega 为复参数, 其中 β\beta 为一个足够大的实数, 象函数 F(s)F(s) 在复平面上某一区域内收敛, 定义拉普拉斯变换 (Laplace)

L[f(t)]=F(s)=∫0+∞f(t)e−stdt \mathscr{L}[f(t)]=F(s)=\int_{0}^{+\infty}f(t)e^{-st}dt

基本函数的拉普拉斯变换

使用 f(t)f(t) 表示原函数, 即表示时域上的信号, 使用 F(s)F(s) 表示变换后的函数, 称为象函数, 其函数变量 ss 为一个复数
如果没有说明, 通常认为原函数 f(t)=0,t<0f(t)=0,t<0. 部分情况下需要使用 f(t)H(t)f(t)H(t) 强调此特性

有拉氏变换关系 F(s)=L[f(t)]F(s)=L[f(t)]F(s)=L[f(t)]F(s)=\mathscr{L}[f(t)] 或者拉氏反变换 L−1[F(s)]=f(t)\mathscr{L}^{-1}[F(s)]=f(t)

单位阶跃函数

通常使用 H(t)H(t) 或 u(t)u(t) 表示单位阶跃函数, 满足

u(t)={1,t>00,t<0 u(t)=\begin{cases} 1,&t>0\\ 0,&t<0 \end{cases}

易得阶跃函数有性质

H(t−a)=−H(a−t) H(t-a)=-H(a-t)

其有拉氏变换

L[H(t)]=1s \mathscr{L}[H(t)]=\frac{1}{s}

由于拉式变换中规定了信号函数 f(t)f(t) 必定有 t≥0t\ge 0, 因此常数 11 作为信号时, 也可直接视为单位阶跃函数, 因此有

L[1]=1s \mathscr{L}[1]=\frac{1}{s}

阶跃函数变形

对于其他含阶跃的函数均可使用单位阶跃函数表示, 以图中函数举例

f1(t)f_1(t) 可表示为 f1(t)=H(t+2)f_1(t)=H(t+2)
f2(t)f_2(t) 可表示为 f2(t)=H(t+2)−H(t−3)f_2(t)=H(t+2)-H(t-3)

对于方波, 矩形波等也可使用单位阶跃函数表示, 如图三角波有

f(t)=t[H(t)−H(t−1)] f(t)=t[H(t)-H(t-1)]

其拉氏变换可根据延迟定律得出

L[t(H(t)−H(t−1))]=L[tH(t)−(t−1)H(t−1)−H(t−1)]=1s2−e−ss2−e−ss=1−e−s(1+s)s2 \begin{aligned}\mathscr{L}[t(H(t)-H(t-1))]&=\mathscr{L}[tH(t)-(t-1)H(t-1)-H(t-1)]\\ &=\frac{1}{s^2}-\frac{e^{-s}}{s^2}-\frac{e^{-s}}{s}\\ &=\frac{1-e^{-s}(1+s)}{s^2}\end{aligned}

单位冲激函数

使用 δ(t)\delta(t) 表示单位冲激函数.
单位冲激函数在 t≠0t\neq 0 的位置值为 00, 在 t=0t=0 处为一条长为 ∞\infty, 宽为 00, 面积为 11 的直线

其有拉氏变换

L[δ(t)]=1 \mathscr{L}[\delta(t)]=1

幂函数

对于幂函数的拉氏变换有通式

L[tn]=n!sn+1 \mathscr{L}[t^n]=\frac{n!}{s^{n+1}}

例如

单位斜坡函数 f(t)=tf(t)=t, 其拉氏变换为 1s2\frac{1}{s^2}

根据公式有拉氏反变换

L−1(3s8)=L−1(37!7!s7+1)=37! \mathscr{L}^{-1}(\frac{3}{s^{8}})=\mathscr{L}^{-1}(\frac{3}{7!}\frac{7!}{s^{7+1}})=\frac{3}{7!}

指数函数

对于指数函数有拉氏变换通式

L[eat]=1s−a \mathscr{L}[e^{at}]=\frac{1}{s-a}

其中 aa 可取负数, 例如

L[e−5t]=1s+5 \mathscr{L}[e^{-5t}]=\frac{1}{s+5}

三角函数

正弦函数的拉氏变换

L[sin⁡(ωt)]=ωs2+ω2 \mathscr{L}[\sin(\omega t)]=\frac{\omega}{s^2+\omega^2}

余弦函数的拉氏变换

L[cos⁡(ωt)]=ss2+ω2 \mathscr{L}[\cos(\omega t)]=\frac{s}{s^2+\omega^2}

可根据函数名首字母有 ss 变换后无 ss 记忆

例如

L−1(5s2+29)=L−1(52929s2+29)=529sin⁡(29t) \mathscr{L}^{-1}(\frac{5}{s^{2}+29})=\mathscr{L}^{-1}(\frac{5}{\sqrt{29}}\frac{\sqrt{29}}{s^{2}+29})=\frac{5}{\sqrt{29}}\sin(\sqrt{29}t)

复合函数的拉氏变换

根据拉氏变换的线性性, 因此对于拉氏变换与逆变换可先将复合函数拆分为相加形式.

对于复杂分式的拆分方法可见部分分式分解

例如拉氏变换

L[2e−4t+3cos⁡(5t)+4sin⁡(3t)+t5+5δ(t)]=2s+4+3ss2+25+12s2+9+5!s6+5 \mathscr{L}[2e^{-4t}+3\cos(5t)+4\sin(3t)+t^5+5\delta(t)]\\ =\frac{2}{s+4}+\frac{3s}{s^2+25}+\frac{12}{s^2+9}+\frac{5!}{s^{6}}+5

拉氏逆变换

L−1[3s+5+2s+5s2+16+6+7s6]=3e−5t+2cos⁡(4t)+54sin⁡(4t)+6δ(t)+75!t5 \mathscr{L}^{-1}[\frac{3}{s+5}+\frac{2s+5}{s^2+16}+6+\frac{7}{s^6}]\\ =3e^{-5t}+2\cos(4t)+\frac{5}{4}\sin(4t)+6\delta(t)+\frac{7}{5!}t^5

注意在逆变换配凑系数时的方法, 应先写出基本变换公式的形式, 再写出系数使等式左右相等

L−1[5s2+16]=L−1[4s2+1654]=54sin⁡(4t) \mathscr{L}^{-1}[\frac{5}{s^2+16}]=\mathscr{L}^{-1}[\frac{4}{s^2+16}\frac{5}{4}]=\frac{5}{4}\sin(4t)

L−1[7s6]=L−1[5!s675!]=75!t5 \mathscr{L}^{-1}[\frac{7}{s^6}]=\mathscr{L}^{-1}[\frac{5!}{s^6}\frac{7}{5!}]=\frac{7}{5!}t^5

拉普拉斯变换的基本特性

线性性

拉氏变换及其逆变换具有线性性, 即和的变换等于变换的和, 表达式为

L[Af(t)+Bg(t)]=AL[f(t)]+BL[g(t)]=AF(s)+BG(s)L−1[AF(s)+BG(s)]=AL−1[F(s)]+BL−1[G(s)]=Af(t)+Bf(t) \mathscr{L}[Af(t)+Bg(t)]=A\mathscr{L}[f(t)]+B\mathscr{L}[g(t)]=AF(s)+BG(s)\\ \mathscr{L}^{-1}[AF(s)+BG(s)]=A\mathscr{L}^{-1}[F(s)]+B\mathscr{L}^{-1}[G(s)]=Af(t)+Bf(t)

第一移位定理

对于拉氏变换 L[f(t)]=F(s)\mathscr{L}[f(t)]=F(s) 当信号 f(t)f(t) 乘上一个收敛因子 eate^{at} 时, 有拉氏变换

L[eatf(t)]=F(s−a) \mathscr{L}[e^{at}f(t)]=F(s-a)

正变换中, 当存在收敛因子时, 可先忽略, 最后再带入 s−as-a 例如

L[e−3tcos⁡(2t)]=s+3(s+3)2+4=s+3s2+6s+13 \begin{aligned}\mathscr{L}[e^{-3t}\cos(2t)]&=\frac{s+3}{(s+3)^2+4}\\ &=\frac{s+3}{s^2+6s+13}\end{aligned}

逆变换中, 当无法直接套用公式时, 可先尝试提出 s−as-a, 有例子

L−1[2s+3s2+2s+10]=L−1[2(s+1)+3−2(s+1)2+10−1]=L−1[2(s+1)(s+1)2+9+3(s+1)2+913]=e−t[2cos⁡(3t)+13sin⁡(3t)] \begin{aligned}\mathscr{L}^{-1}[\frac{2s+3}{s^2+2s+10}]&=\mathscr{L}^{-1}[\frac{2(s+1)+3-2}{(s+1)^2+10-1}]\\ &=\mathscr{L}^{-1}[\frac{2(s+1)}{(s+1)^2+9}+\frac{3}{(s+1)^2+9}\frac{1}{3}]\\ &=e^{-t}[2\cos(3t)+\frac{1}{3}\sin(3t)]\end{aligned}

第二移位定理 (延迟定理)

对于拉氏变换 L[f(t)]=F(s)\mathscr{L}[f(t)]=F(s)

如图所示, 如果要让信号整体向右平移 aa, 不能仅变换信号函数的自变量为 f(t−a)f(t-a), 因为还会将函数 <0<0 的部分一起引入. 为了消除额外引入的部分, 还需要乘上阶跃函数, 因此整体平移后的信号为

f(t−a)H(t−a) f(t-a)H(t-a)

对于这个整体平移的信号, 有拉氏变换 (注意, 第二移位定理中 aa 的符号没有改变且变换后变量保持为 ss, 与第一移位定理不同)

L[f(t−a)H(t−a)]=F(s)e−as \mathscr{L}[f(t-a)H(t-a)]=F(s)e^{-as}

其中称 e−ase^{-as} 为延迟环节

对于正变换, 可先根据 H(t−a)H(t-a) 部分配凑, 其中对常数的延迟变换满足(见阶跃函数变换说明) L[mH(t−a)]=mse−as\mathscr{L}[mH(t-a)]=\frac{m}{s}e^{-as}, 有例子

L[t⋅H(t−5)]=L[(t−5)⋅H(t−5)+5H(t−5)]=1s2e−5s+5se−5s \begin{aligned}\mathscr{L}[t\cdot H(t-5)]&=\mathscr{L}[(t-5)\cdot H(t-5)+5H(t-5)]\\&=\frac{1}{s^2}e^{-5s}+\frac{5}{s}e^{-5s}\end{aligned}

对于反变换, 可先忽视延迟环节变换, 再对结果代入 t→t−at\to t-a 并乘上 H(t−a)H(t-a) , 有例子

L−1[e−3sss2+16]=cos⁡[4(t−3)]H(t−3) \mathscr{L}^{-1}[e^{-3s}\frac{s}{s^2+16}]=\cos[4(t-3)]H(t-3)

微分特性定理

当对原函数求导后, 其导数的象函数与原来的象函数与原函数存在关系

L[f(n)(t)]=snF(s)−sn−1f(0)−sn−2f′(0)−...−f(n−1)(0)=snF(s)−∑i=1nsn−if(i−1)(0) \begin{aligned}\mathscr{L}[f^{(n)}(t)]&=s^{n}F(s)-s^{n-1}f(0)-s^{n-2}f'(0)-...-f^{(n-1)}(0)\\ &=s^{n}F(s)-\sum_{i=1}^{n}s^{n-i}f^{(i-1)}(0)\end{aligned}

可根据 ss 的次数从第一项开始递减, 且 ss 的次数与 f(0)f(0) 的导数之和始终为 n−1n-1 记忆

注意区分 snF(s)s^nF(s) 部分为象函数, saf(b)(0)s^af^{(b)}(0) 部分为以 00 处原函数导数的值为系数的多项式例如

L[f′(t)]=sF(s)−f(0)L[f′′(t)]=s2F(s)−sf(0)−f′(0) \begin{aligned}\mathscr{L}[f'(t)]&=sF(s)-f(0)\\ \mathscr{L}[f''(t)]&=s^2F(s)-sf(0)-f'(0)\end{aligned}

以 yy 表示原函数, YY 表示象函数, 以常微分方程为例

y′′+4y′+3y=3cos⁡(2t),y(0)=0,y′(0)=1 y''+4y'+3y=3\cos(2t),y(0)=0,y'(0)=1

首先先求出各个原函数导数对应的象函数

L[y′]=sY−y(0)L[y′′]=s2Y−sy(0)−y′(0) \begin{aligned}\mathscr{L}[y']&=sY-y(0)\\ \mathscr{L}[y'']&=s^2Y-sy(0)-y'(0)\end{aligned}

对方程两侧同时变换有

(s2+4s+3)Y−1=3ss2+4Y=s2+3s+4(s2+4)(s+1)(s+3)Y=165−3s+48s2+4+151s+1+−2131s+3 \begin{aligned}(s^2+4s+3)Y-1&=\frac{3s}{s^2+4}\\ Y&=\frac{s^2+3s+4}{(s^2+4)(s+1)(s+3)}\\ Y&=\frac{1}{65}\frac{-3s+48}{s^2+4}+\frac{1}{5}\frac{1}{s+1}+\frac{-2}{13}\frac{1}{s+3} \end{aligned}

通过求 YY 的逆变换解出

y(t)=L−1[Y(s)]=−365cos⁡(2t)+2465sin⁡(2t)+e−t5−2e−3t13 y(t)=\mathscr{L}^{-1}[Y(s)]=\frac{-3}{65}\cos(2t)+\frac{24}{65}\sin(2t)+\frac{e^{-t}}{5}-\frac{2e^{-3t}}{13}

积分特性定理

当对原函数积分后, 其导数的象函数与原来的象函数与原函数存在关系

L[∫0tf(λ)dλ]=1sF(s) \mathscr{L}[\int_0^t f(\lambda)d\lambda]=\frac{1}{s}F(s)

同理对于多重积分有

L[∫0t∫0tf(λ)dλ]=1s2F(s) \mathscr{L}[\int_0^t\int_0^t f(\lambda)d\lambda]=\frac{1}{s^2}F(s)

终值与初值定理

在对信号函数求极限时有如下定理

求信号在 00 处的极限有如下公式, 称为初值定理

lim⁡t→0f(t)=lim⁡s→∞sF(s) \lim_{t\to 0}f(t)=\lim_{s\to \infty}sF(s)

求信号在 ∞\infty 处的极限有如下公式, 称为终值定理

lim⁡t→∞f(t)=lim⁡s→0sF(s) \lim_{t\to \infty}f(t)=\lim_{s\to 0}sF(s)

卷积定理

设原函数 g(t)g(t) 的象函数为 G(t)G(t)

定义卷积运算有 (注意卷积运算后得到的结果为一个函数), 注意与傅里叶变换不同, t<0t<0 时, f(t)=g(t)=0f(t)=g(t)=0, 被积函数在区间 [0,t][0,t] 外的值均为 00, 因此拉普拉斯变换下的卷积中, 要以 tt 为积分上限 (如果参与运算的函数严格复合信号函数的条件, 则以 ∞\infty 作为积分上限时与 tt 等价)

(f∗g)(t)=∫0tf(τ)g(t−τ)dτ=∫0tf(t−τ)g(τ)dτ (f*g)(t)=\int_0^t f(\tau)g(t-\tau)d\tau=\int_0^t f(t-\tau)g(\tau)d\tau

卷积的拉氏变换满足 (此定理重要性较低)

L[(f∗g)(t)]=F(s)⋅G(s) \mathscr{L}[(f*g)(t)]=F(s)\cdot G(s)

有例题 (注意, 积分变量为 λ\lambda, 因此 e2te^{2t} 对于被积函数来说相当于常量)

L[e2t∫0te−2λsin⁡(3λ)dλ]=L[∫0te2(t−λ)sin⁡(3λ)dλ]=L[e2t∗sin⁡(3t)]=1s−2⋅3s2+9 \begin{aligned}\mathscr{L}[e^{2t}\int^{t}_0 e^{-2\lambda}\sin(3\lambda)d\lambda]&=\mathscr{L}[\int^{t}_0 e^{2(t-\lambda)}\sin(3\lambda)d\lambda]\\ &=\mathscr{L}[e^{2t}*\sin(3t)]\\ &=\frac{1}{s-2}\cdot\frac{3}{s^2+9} \end{aligned}

微分与积分关系定理

原函数与象函数之间有微分与积分关系 (此定理重要性较低)

L[tf(t)]=−dF(s)ds \mathscr{L}[tf(t)]=-\frac{dF(s)}{ds}

L[1tf(t)]=∫s+∞F(s)ds \mathscr{L}[\frac{1}{t}f(t)]=\int_{s}^{+\infty}F(s)ds

例如

L[tsin⁡2t]=−dds[2s2+4]=4s(s2+4)2 \begin{aligned}\mathscr{L}[t\sin 2t]&=-\frac{d}{ds}[\frac{2}{s^2+4}]\\ &=\frac{4s}{(s^2+4)^2}\end{aligned}

Z 变换

参考自 https://www.bilibili.com/video/BV1nT4y167Zoopen in new window

用于离散系统分析, 处理离散时间信号
Z 变换可以将离散信号转换为复频域信号

冲激序列

假设采样间隔为 TT, 有冲激序列

δT(t)=∑n=0∞δ(t−nT) \delta_T(t)=\sum_{n=0}^{\infty}\delta(t-nT)

每隔 TT 间隔有一个冲激函数

信号采样

将信号采样表示为冲激序列与信号相乘的形式, 有

fs(t)=δT(t)⋅f(t) f_s(t)=\delta_T(t)\cdot f(t)

虽然 fs(t)f_s(t) 名为对信号的采样, 但 fs(nT)=f(nT)⋅δ(nT−nT)≠f(nT)f_s(nT)=f(nT)\cdot \delta(nT-nT) \neq f(nT)

信号采样的拉普拉斯变换

显然 fs(t)f_s(t) 在时间上是连续的, 因此可以进行拉普拉斯变换
根据冲激函数 δ(t)\delta(t) 的采样定律有

Fs(s)=∫−∞∞δT(t)⋅f(t)e−stdt=∑n=−∞∞∫−∞∞δ(t−nT)⋅f(t)e−stdt=∑n=−∞∞f(nT)e−snT \begin{split} F_s(s)&=\int_{-\infty}^{\infty}\delta_T(t)\cdot f(t)e^{-st}\mathrm{d}t\\ &=\sum_{n=-\infty}^{\infty}\int_{-\infty}^{\infty}\delta(t-nT)\cdot f(t)e^{-st}\mathrm{d}t\\ &=\sum_{n=-\infty}^{\infty}f(nT)e^{-snT} \end{split}

其中 f(nT)f(nT) 为一个间隔 TT 的序列, 取值仅与顺序 nn 有关, 因此记为 f[n]f[n]

变量代换

定义变量代换 (注意 zz 仍然为复数)

z=esT z=e^{sT}

将此时信号采样的拉普拉斯变换即为 ZZ 变换

Z[fs(t)]=Fs(z)=∑n=0∞f[n]z−n Z[f_s(t)]=F_s(z)=\sum_{n=0}^{\infty}f[n]z^{-n}

与拉普拉斯变换的关系

即 ZZ 变换的一般形式, 其本质为一个复数域的幂级数

对于拉普拉斯变换 ss 为一个复平面上的点
经过变换 z=esTz=e^{sT}, ss 的左半平面映射到 zz 的单位圆内, 右半部分则映射到单位圆外, 虚轴映射为单位圆

Z 变换举例

单位阶跃函数 1(t)1(t) 的 ZZ 变换
对于 f(t)=1(t)f(t)=1(t) 显然有 f[n]=1f[n]=1
根据此类幂级数可得

Z[Fs(z)]=∑n=0∞z−n=∑n=0∞(1z)n=zz−1 Z[F_s(z)]=\sum_{n=0}^{\infty}z^{-n}=\sum_{n=0}^{\infty}(\frac{1}{z})^{n}=\frac{z}{z-1}

单位冲激函数 δ(t−kT)\delta(t-kT) 的 ZZ 变换

根据 L[δ(t−kT)]=e−kTs\mathscr{L}[\delta(t-kT)]=e^{-kTs}, 以及变量代换 z=esTz=e^{sT} 可得

Z[δ(t−kT)]=z−k Z[\delta(t-kT)]=z^{-k}

代数知识

多项式分式

基本形式

对于任意实数多项式分式具有以下通常形式

F(x)=P(x)Q(x)=∏i=0naixi∏i=0mbixi F(x)=\frac{P(x)}{Q(x)}=\frac{\prod_{i=0}^{n}{a_ix^i}}{\prod_{i=0}^{m}{b_ix^i}}

多项式方程的解

多项式方程 P(x)=0P(x)=0 的解有以下的特性

  1. 方程的解必定为实数 x=mx=m 或一对共轭复数 x=a+bi,x=a−bix=a+bi,x=a-bi
  2. 当多项式的最高次幂为 nn, 则方程包括重根共有 nn 个解
  3. 根据方程的解, 可将多项式因式分解为 K∏(x−xi)K\prod(x-x_i) 的形式, 对于其中的复数根有 (x−a−bi)(x−a+bi)=x2−2ax+a2+b2(x-a-bi)(x-a+bi)=x^2-2ax+a^2+b^2

零极点形式

对于多项式分式
当 P(z)=0P(z)=0 时, F(z)=0F(z)=0, 因此称 zz 为 F(x)F(x) 的零点
当 Q(p)=0Q(p)=0 时, F(p)=∞F(p)=\infty, 因此称 pp 为 F(x)F(x) 的极点

根据上述结论, 可将多项式分式化简为如下形式

F(x)=K∏(x−zi)∏[x2−2aix+(ai2+bi2)]∏(x−pi)∏[x2−2cix+(ci2+di2)] F(x)=\frac{K\prod(x-z_i)\prod[x^2-2a_ix+(a_i^2+b_i^2)]}{\prod(x-p_i)\prod[x^2-2c_ix+(c_i^2+d_i^2)]}

其中
ziz_i, ai±biia_i\pm b_ii 为方程 P(x)=0P(x)=0 的解, 即 F(x)F(x) 的零点
pip_i, ci±diic_i\pm d_ii 为方程 Q(x)=0Q(x)=0 的解, 即 F(x)F(x) 的极点

复数

注意

该部分笔记尚未完成

设 u=a+bi=Aϕ,v=c+di=Bθu=a+bi=A\phase\phi,v=c+di=B\phase\theta 其中

  • AA 为 uu 的幅值
  • ϕ\phi 为 uu 的相位
  • aa 为 uu 的实部 ReRe
  • bb 为 uu 的虚部 ImIm

分母有理化

定义 uˉ=a−bi\bar{u}=a-bi 为 uu 的共轭复数, 共轭复数与复数间的相乘满足

uuˉ=a2+b2 u\bar{u}=a^2+b^2

因此可将分母乘上其共轭复数, 将其转为实数

a+bic+di=(a+bi)(c−di)(c+di)(c−di)=ac+bd−adi+bcic2+d2 \frac{a+bi}{c+di}=\frac{(a+bi)(c-di)}{(c+di)(c-di)}=\frac{ac+bd-adi+bci}{c^2+d^2}

复数乘除法

根据复数的几何意义可得

复数乘法为幅值相乘, 相位相加

u×v=(A×B)(ϕ+θ) u\times v=(A\times B)\phase{(\phi+\theta)}

复数除法为幅值相除, 相位相减

u÷v=(A÷B)(ϕ−θ) u\div v=(A\div B)\phase{(\phi-\theta)}

复数幅值

对于单个复数, 其幅值及其实部与虚部平方和开根号

∣a+bi∣=a2+b2 |a+bi|=\sqrt{a^2+b^2}

根据复数乘除法
对于多个复数组成的分式, 其幅值为分子因式幅值之积除以分母因式幅值之积

∣∏(am+bmi)∏(cn+dni)∣=∏∣am+bmi∣∏∣cn+dni∣=∏am2+bm2∏cn2+dn2 \begin{split}|\frac{\prod (a_m+b_mi)}{\prod (c_n+d_ni)}| &=\frac{\prod |a_m+b_mi|}{\prod |c_n+d_ni|}\\ &=\frac{\prod\sqrt{a_m^2+b_m^2}}{\prod\sqrt{c_n^2+d_n^2}}\end{split}

复数辐角

对于单个复数即关于其实部与虚部的四象限 arctan⁡\arctan

a+bi=Arctan(a,b)⁡ \phase{a+bi}=\operatorname{Arctan(a,b)}

根据复数乘除法 对于多个复数组成的分式, 其相位为分子因式相位之和减去分母因式幅值之和

∏(am+bmi)∏(cn+dni)=∑am+bmi−∑cn+dni=∑Arctan⁡(am,bm)−∑Arctan⁡(cn,dn) \begin{split}\phase{\frac{\prod (a_m+b_mi)}{\prod (c_n+d_ni)}} &=\sum \phase{a_m+b_mi}-\sum \phase{c_n+d_ni}\\ &=\sum \operatorname{Arctan}(a_m,b_m)-\sum\operatorname{Arctan}(c_n,d_n)\end{split}

部分分式分解

部分分式分解将分式转变为几个简单分式之和, 此技巧在使用拉普拉斯变换求解微分方程与分式不定积分中有重要作用

根据拉氏反变换中, F(s)F(s) 经常为复杂的分式形式, 为了完成拉氏反变换, 需要分解这些分式才能更好地套用基本函数的变换公式.

部分分式分解即将分式由左侧分解为右侧形式

(x)=∑i=0maixi∑j=0nbjxj=∑k=0pdk1x+dk2ck2x2+ck1x+ck0 (x)=\frac{\sum_{i=0}^{m} a_ix^i}{\sum_{j=0}^{n} b_jx^j}=\sum_{k=0}^{p}\frac{d_{k1}x+d_{k2}}{c_{k2}x^2+c_{k1}x+c_{k0}}

定义其中

∑i=0maixi=P(x) \sum_{i=0}^{m} a_ix^i=P(x)

∑j=0nbjxj=Q(x) \sum_{j=0}^{n} b_jx^j=Q(x)

部分分式分解通常需要以下步骤

假分式转为真分式

当 P(x)P(x) 的最高次项 mm 大于等于 Q(x)Q(x) 的最高次项 nn, 则称 f(x)f(x) 为假分式, 此时可以从高次项到第次项逐步分解为真分式与多项式之和

P(x)=∑i=0m−ngixi+P′(x)Q(x) P(x)=\sum^{m-n}_{i=0}g_{i}x^{i}+\frac{P'(x)}{Q(x)}

此时真分式的分子 P′(x)P'(x) 的最高次项 m′<nm'<n, 之后对真分式部分 f′(x)=P′(x)Q(x)f'(x)=\frac{P'(x)}{Q(x)} 进行分解.

例如 (也可使用多项式除法的方法, 留下余数)

2x2+x−4x+1=2x+2x2+x−4−2x(x+1)x+1=2x−1+−x−4+(x+1)x+1=2x−1+−3x+1 \begin{aligned}\frac{2x^2+x-4}{x+1}&=2x+\frac{2x^2+x-4-2x(x+1)}{x+1}\\ &=2x-1+\frac{-x-4+(x+1)}{x+1} \\&=2x-1+\frac{-3}{x+1}\end{aligned}

分母因式分解

对分母部分 Q(x)Q(x) 因式分解 (使用试根法, 十字相乘法与多项式除法等), 得到

Q(x)=∏i=0uqi(x) Q(x)=\prod_{i=0}^{u}q_i(x)

注意对于因式 ax2+bx+cax^2+bx+c 还需要检查是否能通过十字相乘法进一步分解

待定系数分解

对于因式 (实根部分) qt(x)=(ax+b)kq_t(x)=(ax+b)^k 可转变为如下分解分式

ht(x)=∑i=1kAi(ax+b)i h_t(x)=\sum_{i=1}^k\frac{A_i}{(ax+b)^i}

对于因式 (虚根部分) qt(x)=(ax2+bx+c)kq_t(x)=(ax^2+bx+c)^k 可转变为如下分解因式

ht(x)=∑i=1kBix+Ci(ax2+bx+c)i h_t(x)=\sum_{i=1}^k\frac{B_ix+C_i}{(ax^2+bx+c)^i}

其中 Aj,Bj,CjA_j,B_j,C_j 为待定系数

求解待定系数

经过待定系数分解后的形式为

P′(x)Q(x)=∑t=0uht(x) \frac{P'(x)}{Q(x)}=\sum_{t=0}^{u}h_{t}(x)

两侧同时乘以 Q(x)Q(x), 可将两侧同时化为一般多项式

P′(x)=Q(x)∑t=0uht(x) P'(x)=Q(x)\sum_{t=0}^{u}h_{t}(x)

根据一般多项式同次项系数相同的条件, 可列出关于待定系数的线性方程, 将待定系数解出便完成分解

分解示例

部分分式分解

f(x)=x2+x+1x4+5x3+x2−9x−6 f(x)=\frac{x^2+x+1}{x^4+5x^3+x^2-9x-6}

题目为真分式, 因此不需要假分式分解, 使用试根法可得 Q(−1)=0Q(-1)=0, 因式分解得到待定系数式

f(x)=x2+x+1x4+5x3+x2−9x−6=x2+x+1(x+1)2(x2+3x−6)=A1(x+1)2+A2x+1+A3x+A4x2+3x−6 \begin{aligned}f(x)&=\frac{x^2+x+1}{x^4+5x^3+x^2-9x-6}\\ &=\frac{x^2+x+1}{(x+1)^2(x^2+3x-6)}\\ &=\frac{A_1}{(x+1)^2}+\frac{A_2}{x+1}+\frac{A_3x+A_4}{x^2+3x-6} \end{aligned}

两边同时乘以分母

P(x)=Q(x)[A1(x+1)2+A2x+1+A3x+A4x2+3x−6]=A1(x2+3x−6)+A2(x+1)(x2+3x−6)+(A3x+A4)(x+1)2=A1(x2+3x−6)+A2(x3+4x2−3x−6)+A3(x3+2x2+x)+A4(x2+2x+1) \begin{aligned}P(x)&=Q(x)[\frac{A_1}{(x+1)^2}+\frac{A_2}{x+1}+\frac{A_3x+A_4}{x^2+3x-6}]\\ &=A_1(x^2+3x-6)+A_2(x+1)(x^2+3x-6)+(A_3x+A_4)(x+1)^2\\ &=A_1(x^2+3x-6)+A_2(x^3+4x^2-3x-6)+A_3(x^3+2x^2+x)+A_4(x^2+2x+1)\end{aligned}

整理式子, 根据同次项系数相等的条件得到线性方程组

[−6−6013−31214210110][A1A2A3A4]=[1110] \begin{bmatrix} -6&-6&0&1\\ 3&-3&1&2\\ 1&4&2&1\\ 0&1&1&0\end{bmatrix} \begin{bmatrix}A_1\\A_2\\A_3\\A_4\end{bmatrix}= \begin{bmatrix}1\\1\\1\\0\end{bmatrix}

求解得到

{A1=−18A2=764A3=−764A4=2932 \begin{cases}A_1=\frac{-1}{8}\\A_2=\frac{7}{64}\\A_3=\frac{-7}{64}\\A_4=\frac{29}{32}\end{cases}

代入待定系数得到分解结果

f(x)=−18(x+1)2+764(x+1)+−7x+5864(x2+3x−6) f(x)=\frac{-1}{8(x+1)^2}+\frac{7}{64(x+1)}+\frac{-7x+58}{64(x^2+3x-6)}