fid = fopen('test_001.img','w'); r=random('Normal',100,0,64,64); z=random('Uniform',0,5,64,64); %%%%%%%%%%%%%%%%%%%%%%%%%%%%% a=14; b=16; for i=15:18 for j=a:b r(i,j)=random('Normal',102,0,1,1); end a = a-1; b = b+1; end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%