function?stevenson?=?stevenson_ht（imagefilename）
%?誤差擴散算法——stevenson半調算法。將灰度圖像很好的轉換為二值圖像。
%?對輸入的image圖像進行stevenson半色調處理并將生成的半色調圖像保存為文件
%?image:待半色調處理的圖像
%?filename:保存為圖像文件的文件名
%?stevenson半色調算法是有半色調算法中誤差擴散的一種。

[mn]=size（image）;
t=ones（mn）/2;
e=zeros（mn）;???????%e為量化誤差矩陣，初始值為零
image=double（image）;

for?i=1:m
????for?j=1:n
????????lu1?=?LU1（ije）;
????????lu2?=?LU2（ije）;
????????lu3?=?LU3（ije）;
????????lu4?=?LU4（ije）;
????????lu5?=?LU5（ije）;
????????u??=?U（ije）;
????????ru1?=?RU1（ije）;
????????ru2?=?RU2（ije）;
????????ru3?=?RU3（ije）;
????????ru4?=?RU4（ije）;
????????ru5?=?RU5（ije）;
????????l??=?L（ije）;
????????aij?=?lu1?+?lu2?+?lu3?+?lu4?+?lu5?+?u?+?ru1?+?ru2?+?ru3?+?ru4?+ru5?+?l;
????????fij?=?image（ij）/255?-?aij;
????????if?fij?????????????stevenson（ij）?=?0;?????%?0為黑
????????else
????????????stevenson（ij）?=?1;
????????end
????????e（ij）?=?stevenson（ij）?-?fij;
????end
end

%imshow（stevenson）;
imwrite（stevensonstrcat（filename‘_stevenson_halftone.bmp‘）‘bmp‘）;

%|-------------------------------|???|-------------------------------------|?
%|??????????????x???????32???????|???|?lu5???????lu4??????ru4?????????ru5??|
%|?12???????26??????30??????16???|???|??????lu3???????u?????????ru3????????|
%|??????12??????26??????12???????|???|?lu2???????lu1??????ru1?????????ru2??|
%|??5???????12??????12??????5????|???|???????l????????x????????????????????|
%|-------------------------------|??