Lb的维度不对,L是3*3的,b是1*3的不能相除,若原来的方程为Ax=b,那么b应该是3*1的矩阵才行,也就是b要进行一下转置,这样Ax=L*U*x=b',推出x=U(Lb')
>> x=U(Lb')
x =
2.7879
17.5152
16.6667
Lb的维度不对,L是3*3的,b是1*3的不能相除,若原来的方程为Ax=b,那么b应该是3*1的矩阵才行,也就是b要进行一下转置,这样Ax=L*U*x=b',推出x=U(Lb')
>> x=U(Lb')
x =
2.7879
17.5152
16.6667