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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=15:18
for j=31:34
r(i,j)=random('Normal',102,0,1,1);
end
end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=14:19
for j=47:48
r(i,j)=random('Normal',102,0,1,1);
end
end for i=16:17
for j=45:50
r(i,j)=random('Normal',102,0,1,1);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a=13;
b=17;
for i=30:35
for j=a:b
r(i,j)=random('Normal',102,0,1,1);
end a=a-1;
b=b+1; end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=30:35
for j=30:35
r(i,j)=random('Normal',102,0,1,1);
end
end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i =29:35
for j=47:49
r(i,j)=random('Normal',102,0,1,1);
end
end for i=31:33
for j=45:51
r(i,j)=random('Normal',102,0,1,1);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a=13;
b=18;
for i=46:53
for j=a:b
r(i,j)=random('Normal',102,0,1,1);
end
a=a-1;
b=b+1;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=45:53
for j=29:37
r(i,j)=random('Normal',102,0,1,1);
end
end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=44:53
for j=46:49
r(i,j)=random('Normal',102,0,1,1);
end
end for i=47:50
for j=43:52
r(i,j)=random('Normal',102,0,1,1);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
r=r+z;count =fwrite(fid,r);
fclose('all');
fid = fopen('test_002.img','w');
s=random('Normal',100,0,64,64);z=random('Uniform',0,5,64,64);
for i=8:12
for j=8:12
s(i,j)=random('Normal',102,0,1,1);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
s=s+z;count=fwrite(fid,s);
fclose('all');