MATLAB 向量化编程实践

mytomorrow

3 h" `2 h/ A9 ~. L& @* WMATLAB的特色就是向量化的编程，所以有必要练习如何使用向量化编程，用一个小例子开始，例如下面的正弦函数求和运算，要求画出图像。



' M5 A) b! S0 y  U毋庸置疑，我们肯定可以使用循环等和其他语言一样的方法实现这个函数的求和运算，但这体现不出MATLAB的特色，所以呢？我们来一起推导下，如何使用向量化的编程：

1 L  j8 ]: y7 C- k) H
" I4 G& `+ M5 F1 Q2 a$ P
, V. ~, C  _; E) V- }; Y) j+ E: r: M上面的手稿也可以看到了，式子（1）就很适合直接用MATLAB编程了，那么我们写出这个代码来，最后我们需要对结果进行转置下，因为我想得到的结果值是一个行向量：

• %try to use vector style to programming by using Matlab
• clc
• clear
• close all
• t = 0:0.01:1;
• k = 1:2:5;
• xt = sin(2*pi*t'*k)*(1./k');
• xt = xt';
• plot(t,xt);
• xlabel('t');
• ylabel('amplititude');
  3 i6 q7 D& D7 c1 o) i- [" }

  
; B  a) M+ `) [) r- ^, i
5 e' D( V5 g4 `% d  J
% a# T! S1 A! M同样，我们也可以对式（1）两边同时转置：

! w4 R5 V$ ?# l" |* ~; l+ C9 S
+ u# e0 J: ?# ^. |% y
3 h8 A$ G  a6 I这样对应的MATLAB代码为：
0 e: m7 ?. x& I- Z

• %try to use vector style to programming by using Matlab
• clc
• clear
• close all
• t = 0:0.01:1;
• k = 1:2:5;
• xt = (1./k)*sin(2*pi*k'*t);
• plot(t,xt);
• xlabel('t');
• ylabel('amplititude');
  5 R, A: A2 M4 d

  



5 Z, R+ X5 r. D7 ~/ j2 B" [1 N* W" l最后，我给出使用循环语句实现这个例子的程序：
* o( z$ A! t+ @

• clc
• clear
• close all
• t = 0:0.01:1;
• xt = zeros(1,length(t));
• for k = 1:2:5
•     xt = xt + (1/k)*sin(2*pi*k*t);
• end
• plot(t,xt);
• xlabel('t');
• ylabel('amplititude');
  6 X5 c' y) W1 s/ T

   
2 E# J- O  k. e( p* \% G

& d' U8 u: ?# Z* \0 O结果都是一样的。
6 N" p; p3 D, l  `& Y5 M
3 x( v3 w; v0 T3 W6 E6 `' F2 N本身暂时不涉及那种方式的优劣以及对于高效性的讨论。
3 p$ U+ s1 m$ Y* s6 a: h
通过上述例子，我希望自己能一步一个脚印的去学会使用向量化编程。
- Z3 ]8 u2 ?  q5 I; u, R

xixihahaheihei

谢谢分享