前段时间改写老师的一个计算纳米细颗粒物的C语言程序，不过运行的效率实在是太低了，工作站跑都很吃劲，晚上看了冈萨雷斯的数字图像处理意外发现了问题所在。在将C语言的程序改写成Matlab程序的时候我还是保持的写C语言的想法，上来全是for循环而忽视了Matlab的矩阵操作。
这里举的例子是冈萨雷斯的书里面的一个例子，比较的是f(x，y)=A sin(u0x+v0y)，用for循环和向量化代码实现，具体如下：
function [rt,f,g] = twodsin(A,u0,v0,M,N)
tic %开始时间
for r = 1:M
u0x = u0*(r-1);
for c = 1:N
v0y = v0*(c-1);
f(r, c) = A*sin(u0x+v0y);
end
end
t1 = toc; %结束时间 tic %开始时间
r = 0:M-1;
c = 0:N-1;
[C,R] = meshgrid(c,r);
g = A*sin(u0*R+v0*C);
t2 = toc; %结束时间
rt = t1/(t2+eps);
测试代码：
[rt,f,g] = twodsin(1, 1/(4*pi), 1/(4*pi), 512, 512);
输出：
rt = 4.0483 (Intel Core i5 CPU 2.67GHz)
可以看出for循环和向量化操作的差距，冈萨雷斯书中给出的是34.2520，这里是因为电脑配置不同的缘故。
输出图像：
