前言
某校的拉普拉斯变换主要的考试方法是解线性微分方程(组),解这样的方程组的过程是
- 根据初值条件,把给定的微分方程左边和右边都转换到s域
- 做除法或者消元,算得Y=L(y)在s域的表达式
- 进行部分因式分解
- 把结果转换回时域
其中第一和第四步都是拉普拉斯变换的内容,在上篇文章已经详细说明。而解方程组和因式分解……懂的都懂。不过根据某数学菜鸡的经验,拉普拉斯变换的部分花的时间少但容易错,因式分解的部分虽然计算艰难但可以验算,复变内容背得熟练的话基本上不会出错。
为避免反复抄几个字母无数次展现本人仅有的线性代数水平,本文中所有的线性方程组都用矩阵形式表示。
二次部分分式分解
矩阵……真的存在吗
分母因式分解后,最容易解决的是实数线性项(的最高次项),利用“遮盖法”即可快速算出结果。
而二次项有两个办法,一个是解出对应的两个虚根,算出结果后再合并;另一个办法便是解方程组了。
在没有计算器的情况下,复数乘除法非常难算,且似乎不易优化。而方程组的方法则比较整齐,适合强迫症,且比较容易验算。
例:
(s+1)(s+3)(s2+4)130s=s+1A+s+3B+s2+4Cs+D
A和B很容易求出
AB=(3−1)(12+4)130×(−1)=−13=(−3+1)(32+4)130×(−3)=15
而C和D列方程组求解
|
A |
B |
C |
D |
|
s3 |
1 |
1 |
1 |
0 |
0 |
s2 |
3 |
1 |
4 |
1 |
0 |
s |
4 |
4 |
3 |
4 |
130 |
1 |
12 |
4 |
0 |
3 |
0 |
即
⎣⎡13412114414300143⎦⎤⎣⎡ABCD⎦⎤=⎣⎡001300⎦⎤
左边两列有很明显的周期性,若把第一行的-4倍加到第三行,第二行的-4倍加到第四行,省略前两行
[0000−1−164−1]⎣⎡ABCD⎦⎤=[1300]
直接就可以算出C和D。这是s2+a2类型的简单之处。若是之前没算出A和B,这里也可以带入再计算。
如果有时间验算,除了可以计算这个矩阵方程本身,还可以带入s为特定值直接检验分式分解的正确性。
解符号表达式方程组
某人认为解二元带符号的方程组,最好的办法是直接背诵
[ACBD]−1=AD−BC1[D−C−BA]
比如解
(s+6)X−4Y−2X+(s+4)Y=A=0
先求出
[s+6−2−4s+4]−1=s2+10+161[s+424s+6]
再左乘右边的向量就好了。这样不需要反复地抄写和化简符号表达式,一直都是最简的。
[XY]=s2+10+161[s+424s+6][A0]=s2+10+16A[s+42]