![MATLAB 8.5基础教程](https://wfqqreader-1252317822.image.myqcloud.com/cover/310/31304310/b_31304310.jpg)
上QQ阅读APP看书,第一时间看更新
2.5 利用矩阵解线性方程组
可将一般的线性方程组使用矩阵形式表示,利用矩阵运算及函数可以很容易解决线性齐次、非齐次方程的有解、无解、无穷多解等问题。
2.5.1 线性方程组的唯一解
线性方程组的矩阵形式为AX=b(A为系数矩阵,b为常数项列向量,X为未知数列向量),其唯一解为X=A−1b。
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-65.jpg?sign=1739186739-fKdF5iyiRysnhl0HmZup0xVUKzuBG9eP-0-05f13a73e8f9ca0303aac3038141a5bf)
2.5.2 齐次线性方程组的通解
齐次线性方程组矩阵形式:AX=0。
格式:Z=null(A,′r′)%Z的列向量是方程AX=0的有理基础解系
【例2-18】 求方程组的通解。
解:求基础解系的程序如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-67.jpg?sign=1739186739-vyhnKfO8Z7hYOA7gtEcUae7y9WERSFds-0-a6faf8053b0308a78c89313b68a6c9c0)
运行结果如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-68.jpg?sign=1739186739-BIRObWpSKPPgIjJ0RMoFlQ6QOynI1yAb-0-045856f24279ba068995b45713a4bc45)
通解表示的程序如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-69.jpg?sign=1739186739-cE4jVqV6m8X1BQtzr25sqSF4lpTswLHx-0-145d1456b51274c57c785fba60c0dbc4)
运行结果如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-70.jpg?sign=1739186739-wZGKs2wLbVlNAbKWVp4ek80RRDAPn6Mg-0-ea84e264752853fdb3f2c379ee93aeac)
2.5.3 非齐次线性方程组的通解
对于非齐次线性方程组,需要先判断方程组是否有解,若有解,则再求通解。其步骤如下:
1)判断AX=b是否有解,若有解则进行第2)步。
2)求AX=b的一个特解。
3)求AX=0的通解。
4)AX=b的通解=“AX=0”的通解+“AX=b”的一个特解。
【例2-19】 求方程组的解。
解:在MATLAB中建立M文件如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-72.jpg?sign=1739186739-IJyRKOo6xBPeE8CSVrRq3v99ADsJtOFK-0-27833d5ceebd592cc5aa728ac17d7586)
运行结果如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-73.jpg?sign=1739186739-jgvq1kYgBnjaic1gwTUwBhKT0Y30ZzzJ-0-42c5aa0e498b3ab5fe08087b2e3d2da1)
说明该方程组无解。
【例2-20】 求方程组的通解。
解:在MATLAB中建立M文件如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-75.jpg?sign=1739186739-gqxQWzgaOQgAxV3HmZkRSS1eRf3ZIUlt-0-b39fd715afca29eda29fecaeb9e29e3a)
运行结果如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-77.jpg?sign=1739186739-b8Q1qgXxXjUQ0fUbV0IuK9cvNVFtNnOm-0-a7c16bd34dc6d7580002529aeb4ff835)
所以原方程组的通解为:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-78.jpg?sign=1739186739-eax8mG8Dho4trx92IljDPWjv06Y7cidA-0-8fcbd0e5d8f1d059e2a4e1772f3f28c3)
2.5.4 超定方程组
超定方程组是指方程的个数大于未知数的个数的线性方程组,通常无精确解,但存在近似的最小二乘解。其解法不需要检查系数矩阵的秩是否小于行数、列数,而直接利用广义逆矩阵函数pinv()计算即可。
格式:X=pinv(A)*b %A为超定方程组的系数矩阵,b为常数项列向量
【例2-21】 求方程组的解。
解:在Matlab中建立M文件如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-80.jpg?sign=1739186739-Uiw9a0gRejr1YFkH7nj3GWzJGX12Lo91-0-b353a1a38e54a5b13083bf36d94c3e52)
运行结果如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-81.jpg?sign=1739186739-UohMM4JbBw2tJKVsty4Frx4P7BmCPIGA-0-5bd902b5bd2398d99af38211e8a653fb)