Create the program to return the orthogonal basis and orthonormal basis
76 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Nguyet Nguyen
am 13 Aug. 2022
Bearbeitet: Torsten
am 14 Aug. 2022
Write a program to input any number of vectors in R^n and return return the orthogonal basis and orthonormal basis of the subspace spanned by these vectors using Gram Schmidt process. Can someone check if I do right
for j= 1:n
v = A (: , j) ;
for i = 1: j-1
R(i,j) = Q (:, i )' * A ( :, j) ;
v = v - R ( i ,j ) * Q ( :, i);
end
R ( j, j ) = norm ( v );
Q (:, j) = v/R ( j, j ) ;
end
0 Kommentare
Akzeptierte Antwort
Torsten
am 13 Aug. 2022
Bearbeitet: Torsten
am 13 Aug. 2022
As far as I can see, you only return the orthonormal basis from the Gram Schmidt process in the matrix Q.
And saving the norms is superfluous if you have to return the orthonormal basis, too.
n = 2;
A = [1 3; 4 -7; -1 -12];
Q = zeros(size(A));
QN = zeros(size(A));
for j= 1:n
v = A (:,j) ;
for i = 1: j-1
rij = Q(:,i).' * A(:,j) / (Q(:,i).' * Q(:,i));
v = v - rij * Q(:,i);
end
Q(:,j) = v;
QN(:,j) = v/norm(v);
end
sqrt(Q.'*Q)
QN.'*QN
2 Kommentare
Torsten
am 14 Aug. 2022
Bearbeitet: Torsten
am 14 Aug. 2022
If the resulting vector v is the zero vector, you shouldn't include it as column in Q. This will happen if vi is already in the span of v1,...,v_i-1. And this will definitely happen if more than n vectors are supplied - as is possible according to the formulation of the assignment.
Weitere Antworten (1)
Siehe auch
Kategorien
Mehr zu Parallel Computing Fundamentals finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!