MATLAB如何求解线性代数方程(组):例子例题
若原方程组写为AX=B(A为系数矩阵,B为右端列向量),可用矩阵左除 \ 做运算,即是 X=A\B。如果将原方程组改写成 XA=B。X 可用矩阵右除 / 求解,即是 X=B/A。若以反矩阵运算求解 AX=B,X=A -1 B,即是 X=inv(A)*B,或是改写成 XA=B,X=BA -1 ,即是X=B*inv(A)。
我们直接以下面的例子来说明这三个运算的用法:
>> A=[3 2 -1;-1 3 2;1 -1 -1];% 将等式的左边系数键入
>> B=[10 5 -1]';% 将等式右边的已知项键入,B要做转置
>> X=A\B % 先以左除运算求解
X = % 注意X为列向量
-2
5
6
>> C=A*X % 验算解是否正确
C = % C=B
10
5
-1
>> A=A';% 将A先做转置
>> B=[10 5 -1];
>> X=B/A % 以右除运算求解的结果亦同
X = % 注意X为列向量
-2
5
6