资源简介
矩阵论作业,让写的矩阵QR分解程序,按施密特正交化、归一化,得到酉阵Q和上三角R
代码片段和文件信息
[nm]= size(A); %取A的行数n和列数m
Q = zeros(nm); %Q为n阶单位阵
Q1 = zeros(nm);
R = zeros(mm);
Q(:1) = A(:1);
Q1(:1) = Q(:1)/norm(Q(:1));
R(11) = norm(Q(:1));
for i = 2:m
% Q(:i) = A(:i);
B = zeros(n1);
for j = 1:i-1
B = B + A(:i)
评论
共有 条评论