部分分式展开法

部分分式展开法

引入

从高中学习数列的时候,我们曾经学过裂项法来将一个多项式分式裂开成多个一次分式相加的形式。

例如这样:
\[
\dfrac{1}{(x - 1)(x + 1)} =\dfrac{1}{2}( \dfrac{1}{x - 1} - \dfrac{1}{x + 1})
\]
当分母的次数是二次的时候,我们能比较容易地猜出来答案是怎么样的。但如果式子是这样的:
\[
f(x) = \dfrac{1}{(x - a_0)(x - a_1)\cdots(x - a_n)}
\]
分母次数大于 \(2\) ,或者不给出具体的 \(a_i\) 的时候,我们就很难通过猜根的方法来得到答案了。

观察

为了得到一个通用的裂项方法,我们先来观察一下式子:
\[
f(x) = \dfrac{1}{(x - a_0)(x - a_1)\cdots(x - a_n)} = \dfrac{b_0}{x-a_0} + \dfrac{b_1}{x-a_1} + \cdots+\dfrac{b_n}{x-a_n}
\]
其中 \(b_i\) 都是常数。

显然,我们的目标是得到一系列的系数 \(b_0, b_1, \cdots b_n\) ,那我们要怎么去得到系数呢?

如果我们把系数 \(b_i\) 当成一系列的未知数,那么其实这是不是一个 \(n\) 元方程?分别代入 \(n\) 个不同的 \(x\) 的值,就可以得到 \(n\) 个不同的方程,组成一个 \(n\) 元一次方程组。

不过,难道我们要用高斯消元去解这 \(n\) 个方程???我们再仔细思考一下,什么样的方程组比较好解?

假如,当每个方程都只含一个变量 \(b_i\) 的时候,就是一个普通的一元一次方程,通过移项,我们就可以得到 \(b_i\) 的值了。

那么接下来的关键,就是找到一个特殊的值 \(x_i\) ,去消去 \(b_i\) 外的其他的 \(b\) 的值了。

解决

为了在代入 \(x_i\) 的值的时候消去除 \(b_i\) 以外的 \(b\) 的值,我们可以先保证 \(b_i\) 的系数里没有 \(x\) 。

那么这就启发我们,可以两边乘以 \(b_i\) 的系数的分母 \(x- b_i\) 。

当 \(i = 0\) 时,式子两边同时乘以 \(x - b_0\) ,那么便有如下的式子:
\[
\begin{align}
\dfrac{1}{(x - a_0)(x - a_1)\cdots(x - a_n)} &= \dfrac{b_0}{x-a_0} + \dfrac{b_1}{x-a_1} + \cdots+\dfrac{b_n}{x-a_n} \\dfrac{1}{(x - a_1)(x - a_2)\cdots(x - a_n)} &= b_0 + \dfrac{b_1}{x-a_1}(x-a_0)+ \dfrac{b_2}{x-a_2}(x-a_0) + \cdots+\dfrac{b_n}{x-a_n}(x-a_0)
\end{align}
\]
我们再令 \(x = a_0\),就可以得到其实 \(b_0\) 的值:
\[
b_0 = \dfrac{1}{(a_0 - a_1) (a_0 - a_2) \cdots (a_0 - a_n)} = \prod_{i = 1}^n\dfrac{1}{a_0 - a_i}
\]
同理,当我们要求 \(b_i\) 的时候,就有:
\[
b_i = \prod_{j = 0, j \ne i}^n\dfrac{1}{a_i - a_j}
\]
至此,我们就可以很轻易地得到任意 \(b_i\) 的值了。这就是部分分式展开式法。

注意

分子系数不为1

有时候,我们会遇到要裂项的式子里,分子不是 \(1\) ,也就是要裂项的式子如下:
\[
f(x) = \dfrac{g(x)}{(x - a_0)(x - a_1) \cdots (x - a_n)}
\]
其中 \(g(x)\) 为一个多项式( \(g(x)\) 次数要低于 \(n\) )。

其实通过上面的手法,我们可以知道,通过一样的手法,来得到答案。但为什么我要单独拿出来说呢?

因为这里涉及到一个可不可以用部分分式展开的点,也就是我上面加粗的内容。

试想一下,我们要得到一个这样的答案:
\[
\dfrac{b_0}{x-a_0} + \dfrac{b_1}{x-a_1} + \cdots+\dfrac{b_n}{x-a_n}
\]
通过通分,我们可以很明显地看到,分子的次数最高就是 \(n - 1\),也就是当 \(g(x)\) 的次数高于 \(n - 1\) 的时候,我们就不可用这个方法来得到 \(b_i\) 了。

分母拆不成一次项相乘

由代数基本定理我们可以知道,任意一个多项式,都可以拆分为多个一次或者二次项相乘。所以我们只需要考虑分母为二次的时候的情况就行了。

对于二次项来说,分为三种情况:

  1. 有两个不同实数解
  2. 有两个共轭复数解
  3. 有两个相同实数解

对于1,2两种来说,我们可以用同样的方法来得到上面的答案。但是对于第三种来说,就拆不成两个一次的形式相加了。我们需要转变一下裂项的结果。

假如二次项如下:
\[
\dfrac{p(x - a)+q}{(x - a)^2}
\]
那么通过拆分,我们可以得到:
\[
\dfrac{p}{x - a} + \dfrac{q}{(x - a)^2}
\]
也就是说,当分母里含有一个有两个相同实数解的二次项的时候,需要将二次项拆开成如上的形式。

练习

  1. \(\dfrac{x + 4}{x^3 + 3x^2+2x}\)
  2. \(\dfrac{x + 3}{(x + 1) ^3(x + 2)}\)

小技巧:求得各个系数后,可以令 \(x\) 等于一个还没算过的值进行校验等式两边是否相等。

(习题出自信号与《线性系统分析第四版》(吴大正主编))

答案

  1. \(\dfrac{2}{x}-\dfrac{3}{x + 1} + \dfrac{1}{x + 2}\)
  2. \(\dfrac{2}{(x + 1) ^3} - \dfrac{1}{x + 1}+\dfrac{1}{x + 1}-\dfrac{1}{x + 2}\)

------------恢复内容结束------------

原文地址:https://www.cnblogs.com/ljxtt/p/12182295.html

时间: 2024-10-15 03:08:41

部分分式展开法的相关文章

用Matlab进行部分分式展开

[r p k]=residue[num,den] 例如H(s)=(2s3+5s2+3s+6)/(s3+6s2+11s+6) num=[2 5 3 6]; den=[1 6 11 6]; [r p k]=residue[num,den] r=[-6 -4 3]   %留数 p=[-3 -2 1]     %极点 k=2                %余项

【信号与系统】多项式化简方法

多项式长除法 简介 ??? ? ? ? 多项式长除法?是代数中的一种算法,用一个同次或低次的多项式去除另一个多项式.是常见算数技巧长除法的一个推广版本.它可以很容易地手算,因为它将一个相对复杂的除法问题分解成更小的一些问题. 示例 计算 ? 把被除式.除式按某个字母作降幂排列,并把所缺的项用零补齐,写成以下这种形式: ? 然后商和余数可以这样计算: 将分子的第一项除以分母的最高次项(即次数最高的项,此处为x).结果写在横线之上(x3?÷?x?=?x2). 将分母乘以刚得到结果(最终商的第一项),

学习笔记(信号与系统)

学习笔记(信号与系统) 来源:网络 第一章 信号和系统 信号的概念.描述和分类 信号的基本运算 典型信号 系统的概念和分类 1.常常把来自外界的各种报道统称为消息: 信息是消息中有意义的内容: 信号是反映信息的各种物理量,是系统直接进行加工.变换以实现通信的对象. 信号是信息的表现形式,信息是信号的具体内容:信号是信息的载体,通过信号传递信息. 2.系统(system):是指若干相互关联的事物组合而成具有特定功能的整体. 3.信号的描述--数学描述,波形描述. 信号的分类: 1)确定信号(规则信

DSP 第三章

变换作为一种工具来分析信号和LTI系统. Z变换是处理离散信号的,拉普拉斯是处理连续信号的. Z域就是Z的复平面. Z变换只是一种信号的表达形式. 某个时间序列和X(z)以及收敛域相对应. 表3.1就是对前面的完全总结: 有限时间的因果,反因果,和双边信号.对应收敛域:除掉0,除掉正无穷,除掉0和正无穷. 无限时间的因果,反因果,和双边信号.对应收敛域:圈外,圈内,环. 单边的Z变换是为因果信号准备的.公式化的逆Z变换. 表3.2又是对3.2节的完全总结: 每一种Z变换分两块考虑,Z域变换以及收

[离散时间信号处理学习笔记] 8. z逆变换

z逆变换的计算为下面的复数闭合曲线积分: $x[n] = \displaystyle{\frac{1}{2\pi j}}\oint_{C}X(z)z^{n-1}dz$ 式中$C$表示的是收敛域内的一条闭合曲线.该积分表达式可以利用复数变量理论下的柯西积分定理推导得到.不过本门课程用不上这条式子,因为在离散LTI系统分析中所遇到的典型序列和z变换,有如下更简单的z逆变换求解办法. 观察法(查表) 下面是一个常见序列的z变换表格,通过查表可以由z变换所得的函数反过来求得原序列 Sequence Tr

matlab 常用函数汇总

1. 特殊变量与常数 主题词 意义 主题词 意义 ans 计算结果的变量名 computer 确定运行的计算机 eps 浮点相对精度 Inf 无穷大 I 虚数单位 inputname 输入参数名 NaN 非数 nargin 输入参数个数 nargout 输出参数的数目 pi 圆周率 nargoutchk 有效的输出参数数目 realmax 最大正浮点数 realmin 最小正浮点数 varargin   实际输入的参量 varargout 实际返回的参量     2. 操作符与特殊字符 主题词

MATLAB函数表(转自:http://bbs.06climate.com/forum.php?mod=viewthread&tid=16041&extra=page%3D4)

MATLAB函数表 4.1.1特殊变量与常数 ans 计算结果的变量名 computer 确定运行的计算机 eps 浮点相对精度 Inf 无穷大 I 虚数单位 inputname 输入参数名 NaN 非数 nargin 输入参数个数 nargout 输出参数的数目 pi 圆周率 nargoutchk 有效的输出参数数目 realmax 最大正浮点数 realmin 最小正浮点数 varargin 实际输入 的参量 varargout 实际返回的参量     4.1.2操作符与特殊字符 + 加 -

Matlab信号处理工具箱函数

波形产生和绘图chirp 产生扫描频率余弦diric 产生Dirichlet函数或周期Sinc函数gauspuls 产生高斯调制正弦脉冲pulstran 产生脉冲串rectpuls 产生非周期矩形信号sawtooth 产生锯齿波或三角波sinc 产生sinc函数square 产生方波strips 产生条图tripuls 产生非周期三角波 滤波器分析和实现abs 绝对值(幅值)angle 相位角conv 卷积和多项式乘法conv2 二维卷积fftfilt 基于FFT重叠加法的数据滤波filter

MATLAB命令大全

一.常用对象操作:除了一般windows窗口的常用功能键外.1.!dir 可以查看当前工作目录的文件. !dir& 可以在dos状态下查看.2.who 可以查看当前工作空间变量名, whos 可以查看变量名细节.3.功能键:功能键 快捷键 说明方向上键 Ctrl+P 返回前一行输入方向下键 Ctrl+N 返回下一行输入方向左键 Ctrl+B 光标向后移一个字符方向右键 Ctrl+F 光标向前移一个字符Ctrl+方向右键 Ctrl+R 光标向右移一个字符Ctrl+方向左键 Ctrl+L 光标向左移