fourgs.m

来自「计算连续方程Lyapunov指数的程序,自己编写的」· M 代码 · 共 15 行

% 四维G-S变换
function A = FourGS(V)  % V 为4*4向量
v1 = V(:,1);
v2 = V(:,2);
v3 = V(:,3);
v4 = V(:,4);
a1 = zeros(4,1);
a2 = zeros(4,1);
a3 = zeros(4,1);
a4 = zeros(4,1);
a1 = v1;
a2 = v2-((a1'*v2)/(a1'*a1))*a1;
a3 = v3-((a1'*v3)/(a1'*a1))*a1-((a2'*v3)/(a2'*a2))*a2;
a4 = v4-((a1'*v4)/(a1'*a1))*a1-((a2'*v4)/(a2'*a2))*a2-((a3'*v4)/(a3'*a3))*a3;
A = [a1,a2,a3,a4];