outofplane.m

matlab模拟产生中心应力集中光学条纹图

function y=outofplane
M=512;
N=512;

angle=zeros(M, N);
I=zeros(M, N);
for i=1:M
    for j=1:N
        angle(i,j)=0.188*(j-256)*log(((i-256)*(i-256)+(j-256)*(j-256))/246/246);
        I(i,j)=220.0*abs(sin(0.5*angle(i,j)));
    end
end

H=uint8(I);
y=H;