70/CH3/EX3.4.5/3_4_5.sci


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

//page 166
clear;
close;
clc;
A=[1 0 1;1 0 0;2 1 0];//independent vectors stored in columns of A
disp(A,'A=');
[m,n]=size(A);
for k=1:n
    V(:,k)=A(:,k);
    for j=1:k-1
        R(j,k)=V(:,j)'*A(:,k);
        V(:,k)=V(:,k)-R(j,k)*V(:,j);
    end
    R(k,k)=norm(V(:,k));
    V(:,k)=V(:,k)/R(k,k);
end
disp(V,'Q=')