clear;clc;?
%?close?all;?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

img_name=‘.\sourcePicture\frame179.bmp‘;
I=imread（img_name）;

I1?=?imresize（I[240320]）;
%??subplot（221）;
figure（1）;
imshow（I1）;
imagesc（I1）;
I1=double（I1）/255;

I?=?I1;?
%?figure;imagesc（I）;??%??

Im_dealt?=?I;
[rowcolz]?=?size（Im_dealt）;

landline?=?0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Iw?=?I;?
A?=?0.8;
m=100;
n=300;

tidx?=?2;?
%?for?beta?=?0.05:0.05:0.15
for?beta?=?0.20:0.05:0.30
????for?i=1:3
?????????for?j=landline+1:row
?????????????for?l=1:col????
??????????????d（jl）?=?1/（（j-landline）^.05?+?0.0001）;
?????????????????d2（jl）?=?d（jl）*8;
????????????????if?j??????????????????????d（jl）?=?-0.04*landline?+?18;????
???????????????