relaxzzy.m

来自「MATLAB 6.5 从基础讲解其功能」· M 代码 · 共 32 行

function elapt=relaxzzy(iter)
% 本程序供试验 JIT 和加速器的作用
sz=102;
plate=magic(sz)*64/(sz*sz);
newplate=plate;
im=image(plate);
axis off
set(gcf,'DoubleBuffer','on')
pause
shg
tic
for i=1:iter
    for j=2:(sz-1)
        jm1=j-1;
        jp1=j+1;
        for k=2:(sz-1)
            km1=k-1;
            kp1=k+1;
            newplate(j,k)=(plate(jm1,km1)/2+plate(jm1,k)+...
            plate(jm1,kp1)/2+plate(j,km1)+plate(j,kp1)+...
            plate(jp1,km1)/2+plate(jp1,k)+plate(jp1,kp1)/2)/6;
        end
    end
    plate=newplate;
    if (0==rem(i,5))
        set(im,'cdata',plate)
        drawnow
    end
end
elapt=toc;
close