golub.m

来自「Based on Matlab,Gauss Iteration Method」· M 代码 · 共 13 行

function A = golub(n,s);
%   被交互实验trylu调用
%   GOLUB  生成高度病态的整数测试矩阵
%   GOLUB(n) 给出两个n*n随机整数矩阵的乘积，其中一个是单位下三角的
%   另外一个是单位上三角的。不选主元的 LU 分解无法还原如此坏条件的
%   矩阵。
%   

if nargin < 2, s = 10; end
L = tril(round(s*randn(n)),-1)+eye(n);
U = triu(round(s*randn(n)),1)+eye(n);
A = L*U;