?? qrxq.m
字號:
function [x,Q,R]=qrxq(A,b)
N = size(A);
n = N(1);
B = A; %保存系數矩陣
A(1:n,1)=A(1:n,1)/norm(A(1:n,1)); %將A的第一列正規化
for i=2:n
for j=1:(i-1)
A(1:n,i)= A(1:n,i)-dot(A(1:n,i),A(1:n,j))*A(1:n,j);
%使A的第i列與前面所有的列正交
end
A(1:n,i)=A(1:n,i)/norm(A(1:n,i));
%將A的第i列正規化
end
Q = A; %分解出來的正交矩陣
R = transpose(Q)*B;
x=SolveUpTriangle(R,transpose(Q)*b); %求解方程
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -