LU用doolittle***求源自LUP用doolittle***求
问题来源:已知Ax=b① 求x
L是下三角矩阵U是上三角矩阵,P是一个置换矩阵(P将在下一篇博客中写出)
忽略P下面说明LU的求法:
1.参数矩阵A做如下划分
2.对划汾好的做用doolittle***求(这个用doolittle***求证明简单,问题是怎么想到的(我不知道))
然后就递归的求这个舒尔补的LU用doolittle***求
这样就能求出LU用doolittle***求P是用来防止A是奇异矩阵的。
参考:《算法导论》第28章这上面讲的比较啰嗦