比较使用for循环和响向量化代码计算二维函数图像

matlab

%for循环执行
function [rt,f,g]=twodsin(A,u0,v0,M,N)
%函数的输入值是A,u0,v0,M,N
A=1;u0=1/(4*pi);v0=1/(4*pi);M=50;N=50;
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(r,c);
g=A*sin(u0*R+v0*C);
t2=toc;  %结束时间
%计算两种时间的比率
rt=t1/(t2+eps)
g=mat2gray(g);
imshow(g)