第四章数值计算
4.1矩阵的LU分解、行列式、逆,线性方程组的解
4.1.1 常用算法
B=A‘ % 求矩阵的转置矩阵
det(A) %求矩阵A对应的行列式值
inv(A) %求矩阵A的逆矩阵
x=A\b %求线性方程Ax=b的解
x=b/A %求线性方程xA=b的解
[L,U,P]=lu(A) %矩阵A的上三角分解,U是上三角矩阵,L*U=P*X
[说明]
第四条指令中得斜杠”\”是”左除”符号, 对Ax=b, A在变量x得左边,所以指令中得A必须在”\”得左边,不可放错. 如要解xC=d形式得方程,则要用”右除”指令, x=d/C
【例1】“求逆”法和“左除”法解线性方程的性能对比(1)
0.2368x1+0.2471x2+0.2568x3+1.2671x4=1
.8471
0.1968x1+0.2071x2+1.2168x3+0.2271x4=1
.7471
0.1582x1+1.1675x2+0.1768x3+0.1871x4=1
.6471
1.1161x1+0.1254x2+0.1397x3+0.1490x4=1
.5471
A=[0.2368 0.2471 0.2568 1.2671
0.1968 0.2071 1.2168 0.2271