nrm.m

来自「这是一个三面阵列的数学工具软件」· M 代码 · 共 8 行

function N=nrm(A)
% makes columnwise normalized version of A
[n,m]=size(A);
d=sum(A.^2);
w=d>1e-30;
w=logical(w);
N=zeros(n,m);
if sum(w)>0,N(:,w)=A(:,w)./(ones(n,1)*d(w).^.5);end;