示例代码(GNU Octave(版本9.1.0))

以下可用于求解线性系统A*x=b使用枢轴LU因子分解:

[L,U,P]=lu(A);##现在L*U=P*Ax=U\(L\P)*b;

这是规范矩数组的一种方法十、单位标准:

s=范数(X,“列”);X/=diag;

广播也可以做到这一点(详见广播):

s=范数(X,“列”);十、s

下面的表达式是一种有效计算从置换向量给出的开环符号的方法</p>。它也适用于Octave的早期版本,但速度较慢。

det(eye(长度(p))(p,:))

最后,以下是如何求解线性系统A*x=b使用SVD(骨架非线性)的Tikhonov正则化(岭回归):

m=行(A);n=列(A);[U,S,V]=svd(A);##确定正则化因子alpha##alpha=…##变换到正交基b=U'*b;##使用标准公式,将A替换为S。##S是对角的,因此以下内容将非常快速和准确。x=(S'*S+α^2*eye(n))\(S'*b);##变换到解的基x=V*x;