生成混沌图像

function gl
%生成混沌序列，并生成图像
clear;
clc;
% o=imread('image/lena.bmp');
%%%%%%%%%%%%读入混沌初始值
prompt={'混沌初始值 '};
mytitle='initial value';
lines=1;
def={'98'};
mya=inputdlg(prompt,mytitle,lines,def,'on');
key11=['0.',mya{1}];
key=str2num(key11);
%%%%%%%%%%%%
[filename pathname]=uiputfile({'*.bmp';'*.tiff';'*.tif'}, '保存图像');
os=[pathname filename];
m=512;
n=512;
l=zeros(m,n);
l(1)=key;
for i=2:m*n
l(i)=1-2*l(i-1)*l(i-1);
end
%二值化的时候需要注意，此处的临界值是0，而不是0.5
key=0;
for i=1:m*n
if (l(i)>=key)
l(i)=0;
else
l(i)=1;
end
% waitbar((i/(m*n)),h);
end
l=~l;
% imwrite(l,'image/wm.bmp');
%进度条
h= waitbar(0,'程序处理中，请耐心等待。。。');
for i=1:100, % computation here %
waitbar(i/100);
end
close(h) ;
imwrite(l,os);
% handles.figure;
subplot(1,1,1),imshow(l),title('混沌图像');
% test=imread(os);
% figure,imshow(test,[]);

生成图像如下：