天天看點

《Experiment with MATLAB》讀書筆記(五)

讀書筆記(五)

這是第五部分線性方程求解

%% 前除

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