gs.m

当计算混沌系统Lyapunov指数时

% 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