《矩阵计算》chapter 01 在 Octave 中整理运行示例算法

---Algo1---------------------------------------------------------------------------------------------------------------------

算法1.1.1(点积)

mc_1_1_1_algo_dot_product.m

function c = mc_1_1_1_algo_dot_product (x, y)
  c = 0
  for i=1:size(x)(1)
    c=x(i)*y(i)
  endfor

endfunction

应用：

>> x1=[3;2;1;4]
x1 =

   3
   2
   1
   4

>> y1=[3;2;1;4]
y1 =

   3
   2
   1
   4

>> mc_1_1_1_algo_dot_product(y1, x1)
c = 0
c = 9
c = 4
c = 1
c = 16
ans = 16
>>

---Algo2----------------------------------------------------------------------------------------------------------------------

mc_1_1_2_algo_saxpy.m

function y = mc_1_1_2_algo_saxpy (n, a, x, y)
  for i=1:n
    y(i) = a*x(i) + y(i)
  endfor
endfunction

>> x1=[1;2;3]
x1 =

   1
   2
   3

>> y1=[1;2;3]
y1 =

   1
   2
   3

>> y1=mc_1_1_2_algo_saxpy(3,2,x1,y1)
y1 =

   3
   6
   9

>>

---Algo3--------------------------------------------------------------------------------------------------------------------