gs.m

来自「基于定义法计算Lyapunov指数。只要在matlab命令窗口中直接输入文件中函」· M 代码 · 共 22 行

% 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