logical.m

matlab宝典,电子工业出版社,包含该书的源码4-6章

% Script file logical.m
%
%目的:
%   这段程序代码计算数组中超过6000的平方根 
%  使用下面两种不同的方法：
%  1.使用循环和if结构
%  2.使用逻辑数组

%定义变量
% ii,jj    循环变量
% ave1     第一种方法的平均时间
% ave2     第二种方法的平均时间
% maxcount 循环结算的时间次数

%计算循环结构
maxcount=1;
tic
for jj=1:maxcount
    a=1:10000;
    for ii=1:10000
        if a(ii)>6000
            a(ii)=sqrt(a(ii));
        end
    end
end

%使用逻辑数组的方法计算
ave1=(toc)/maxcount;
maxcount=10;
tic;
for jj=1:maxcount
    a=1:10000;
    b=a>6000;
    a(b)=sqrt(a(b));
end
ave2=(toc)/maxcount;

%显示结果
fprintf('Loop if approach= %9.5f\n',ave1);
fprintf('Logical array approach= %9.5f\n',ave2);