H923(Sat) 科學計算方法與實務探討

H923(Sat) 科學計算方法與實務探討

 MatLab 

A = zeros(8)
A = zeros(8,8)

% 歸零 8x8 矩陣

A = zeros(1,8)
 0 0 0 0 0 0 0 0

A = zeros(8,1)
 0
 0
 ...
 0


A = rand(8) 
%  亂數產生 (0~1)

A = 100*rand(8) 
%  亂數產生 (0~100)

A = 100*rand(8) - 50
%  亂數產生 (-50 ~50)


A = [ A11 A12 A13 ]    %一維陣列
A = [ A11 A12; A21 A22 ]    %二維陣列



  Anxn * Xnx1 = Bnx1

  A- [ A X ] = A-  B
            X   = A-  B
Note: A-  inv(A)
A- A = 1

-----------------------------------------
A- 與 B ,(反過)來求解 X

A = rand(8);
X = rand(8,1);
B = A * X;

Y = inv(A) B;     % A- => inv(A)

比較 Y與 X 差異
norm(X-Y)
ans =  1.1161e-15
結果為一個很小的誤差(數據沒什麼差異)


註: norm 有三個實現方式
 1. sum | X-Y |
  X-Y 絶對值的總合 2. max | X-Y |
  X-Y 絶對值的最大值
 3. sum (X-Y)^2
   X-Y 平方的總合

-----------------------------------------------

 t = [ 0:0.01:2*pi];
 a = sin(t);
 plot(t,a)

 sin wave 繪圖

如果要控制繪圖細節動作,請參考 
H923(Sat) 科學計算方法與實務探討 課後議題