clear?all

disp（‘Data?loading?and?generation...‘）;
ima?=?double（imread（‘lena.png‘））;
L?=?3;
s?=?zeros（size（ima））;
for?k?=?1:L
????s?=?s?+?abs（randn（size（ima））?+?i?*?randn（size（ima）））.^2?/?2;
end
ima_nse?=?ima?.*?sqrt（s?/?L）;

disp（‘Non-it?PPB?filter...‘）;
hW?=?10;
hD?=?3;
alpha?=?0.88;
T?=?inf;
nbit?=?1;
ima_fil_1?=?ppb_nakagami（ima_nse?L?...
?????????????????????????hW?hD?...
?????????????????????????alpha?T?...
?????????????????????????nbit）;

disp（‘It?PPB?filter...‘）;
hW?=?[1?3?5?10];
hD?=?[0?1?2?3];
alpha?=?0.92;
T?=?0.2;
nbit?=?[1?2?3?4];
ima_fil_2?=?ppb_nakagami（ima_nse?L?...
?????????????????????????hW?hD?...
?????????????????????????alpha?T?...
?????????????????????????nbit）;

figure
subplot（221）?plotimage_sar（ima）;
title（‘Noise-free?image‘）;
subplot（222）?plotimage_sar（ima_nse?ima）;
title（sprintf（‘Noisy?image?（PSNR?%.3f）‘?snr（ima_nse?ima）））;
subplot（223）?plotimage_sar（ima_fil_1?ima）;
title（sprintf（‘Non-it?PPB?filter?（PSNR?%.3f）‘?snr（ima_fil_1?ima）））;
subplot（224）?plotimage_sar（ima_fil_2?ima）;
title（sprintf（‘It?PPB?filter?（PSNR?%.3f）‘?snr（ima_fil_2?ima）））;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1114??2010-03-23?13:03??ppbNakagami\example_nakagami.m
?????文件????????1114??2010-03-23?13:03??ppbNakagami\exemple_nakagami.m
?????文件??????151199??2010-03-23?13:03??ppbNakagami\lena.png
?????文件?????????187??2010-03-23?13:03??ppbNakagami\plotimage.m
?????文件?????????326??2010-03-23?13:03??ppbNakagami\plotimage_sar.m
?????文件????????7029??2013-05-17?08:21??ppbNakagami\ppb_nakagami.m
?????文件??????637537??2010-03-23?13:03??ppbNakagami\ppbNakagami.mexa64
?????文件??????618644??2010-03-23?13:03??ppbNakagami\ppbNakagami.mexglx
?????文件???????77824??2010-03-23?13:03??ppbNakagami\ppbNakagami.mexw32
?????文件??????263168??2012-05-26?13:12??ppbNakagami\ppbNakagami.mexw64
?????文件?????????137??2010-03-23?13:03??ppbNakagami\psnr.m
?????文件????????2624??2012-05-26?13:29??ppbNakagami\README
?????文件?????????129??2010-03-23?13:03??ppbNakagami\snr.m