有输出L和U矩阵此外用max函数求第一列的最大值比for循环要快但其他列并没有求最大值而是比较了对角元素和eps的大小关系决定是否换行所以理论上此方法未完全达到选取列主元素的功能。
b比较二
ffPLU1
要比较长若将第二个程序中第二行PLU1的调用改为MATLAB自带函数lu同样可以完成此功能得到PLU3则耗时比较短。
总结MATLAB中编写函数需要规范调用格式如xlupiplzyluAb一共输出五个量在调用时必须全部写出若不写则只能输出x结果得不得其他结果。
本题目需要记录主行交换的信息输出数列IP因此直接调用lu函数不合适。
参考文献
1易大义沈云宝李有法计算方法第2版浙江大学出版社pp175181
2廖波数值计算与最优化原理MATLAB实现北京邮电大学出版社
fr