gs.m

来自「Lorenz系统的李雅谱诺夫指数」· M 代码 · 共 17 行

% G-S正交化
% A = GS(V)
% V 为n*n向量
% A 为n*n正交向量
% Example:
%   A=GS([1 0 1;2 2 0;3 1 1])
% Author：yujunjie
% Author's email: [e-mail]
function A = GS(V)
n=length(V);
for i=1:n
    k=0;  
    for j=1:i-1
        k=k+(A(:,j)'*V(:,i)/(A(:,j)'*A(:,j)))*A(:,j);
    end
    A(:,i)=V(:,i)-k;
    end