![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=1739039767-KeAlHouuOii8mPtLNCrmzP5Fo4ZzYCGK-0-fe2bcd98bf30abd148a8d5e69e698ff7)
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=1739039767-2BtusEiCTQuRrI556MGiDNkpJR5iU3tG-0-8de041d3f8bf0ae53e01fbd9c9590a71)
运行结果如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-68.jpg?sign=1739039767-CxiW09ZwGFFQ1pmxqW81aFDV1ExoiTds-0-72edc6b698413c3d0b3885f2886ddaae)
通解表示的程序如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-69.jpg?sign=1739039767-kanQ7CgqbjKZDBI6mRmVCb6cGwYopVRA-0-9cad2482707cc3a23780ee57fc6d3d79)
运行结果如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-70.jpg?sign=1739039767-d7aVRqTU5N5tOfNOa6Bh1ihbDUNwSLlC-0-365baa5b49be5d337a5d736013dfe0a3)
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=1739039767-hgDqua1XcWpAijSkDldyiY7LzqCfrQVx-0-8657122b91d18e7a06e19cd927f6d90c)
运行结果如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-73.jpg?sign=1739039767-xyHUM1zDjkcp1To0FjqXyKIFjIkn9WKq-0-b1452f4c8ab6e150ed8a4542cc809c38)
说明该方程组无解。
【例2-20】 求方程组的通解。
解:在MATLAB中建立M文件如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-75.jpg?sign=1739039767-ySovp14h99iILDJYHXQLZCTcl5setBme-0-b34a3b17535f76cc895c504a3345aa59)
运行结果如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-77.jpg?sign=1739039767-PFZp2jyhLsm056ADf5ZoWzwqVFG1kYqw-0-582d4c4d900e5ffb6fdb79e074c5a745)
所以原方程组的通解为:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-78.jpg?sign=1739039767-EBhF9ORB4kFifNneJQ9hoe2IXZ9ENJFG-0-320593a70e87cf881083761862939440)
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=1739039767-gja7UlbN7ZYcq6xjwNFUqh5iTZrpRrOd-0-169d34354920ed24856bc8d1a2a2a50a)
运行结果如下:
![](https://epubservercos.yuewen.com/0A09AD/16948915905056606/epubprivate/OEBPS/Images/978-7-111-53210-1-Chapter02-81.jpg?sign=1739039767-wu54HjKfTx40BfDnutHElOGVnqpPKMR8-0-1421194007a60e9e469b16e68f65894d)