讀書筆記(五)
這是第五部分線性方程求解
%% 前除
format bank
A = [3 12 1; 12 0 2; 0 2 3]
b = [2.36 5.26 2.77]'
x = A\b
% X = A\B 就是方程 A*X = B 的解
% 先求逆再左乘
%% 後除
x = b'/A'
% X = B/A 就是方程 X*A = B 的解
% B/A 就是 B*INV(A) 也就是B*A的逆
%% 除數為零或者出現無窮大
A(3,:) = [6 0 1]
A\b
%% 一種避開除數為零的求解方法
b(3) = 2.63
x = A(1:2,1:2)\b(1:2);
x(3) = 0
A*x
%% 另一種避開除數為零的求解方法
z = null(A)
A*z % A*z=0
t = rand % 随機參數不影響結果
y = x + t*z
A*y % A*z=b