function?thr?=?Birge_Massart（c?s）
%用Birge-Massart策略估計閾值
%?給定分解層數j，對j+1及更高層，所有系數保留
%?對第i層（1<=i<=j），保留絕對值最大的n_i個系數，n_i由下式確定
%???n_i?=?M/（j+2-i）^alpha
%?其中M和alpha為經驗系數，缺省情況下M取第一層分解后的系數長度L（1）
%?一般情況下，M滿足?L（1）<=M<=2L（1）
%?alpha的取值因用途不同，壓縮情況下一般取alpha=1.5，降噪一般取alpha=3

%設定alpha、M和j的值
alpha?=?3;
M?=?prod（s（1?:））;
j?=?size（s?1）?-?2;

%計算每層的閾值
for?i?=?1:j
????n?=?round（?M/（j+2-i）^alpha?）;???????????????????%每層保留的系數個數
????D?=?detcoef2（‘compact‘?c?s?i）;???????????????%提取第i層的所有細節系數?
????if?n?>=?size（D）
????????thr（:?i）?=?zeros（3?1）;
????else
????????D?=?sort（abs（D））;???????????????????????????%排序
????????thr（i）?=?D（end-n）;??????????????????????????%設置第i層的閾值
????end
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????844??2007-07-09?15:05??畢設\Birge_Massart.m

?????文件????????767??2007-07-09?15:55??畢設\den1.m

?????文件???????1113??2007-07-09?15:55??畢設\den1_10.m

?????文件???????1600??2007-07-09?15:54??畢設\den1_5_1.m

?????文件????????796??2007-07-09?15:55??畢設\den1_9.m

?????文件????????756??2007-07-09?15:55??畢設\den2.m

?????文件????????856??2007-06-07?08:28??畢設\den3.m

?????文件???????2808??2007-07-09?15:57??畢設\den4.m

?????文件???????4356??2007-07-09?15:58??畢設\denoise0615.m

?????文件???????5379??2007-07-09?15:22??畢設\denoise0701.m

?????文件????????492??2007-06-28?19:33??畢設\Donoho.m

?????文件?????150923??2000-10-07?15:18??畢設\lena.png

?????文件?????476195??2000-10-07?15:05??畢設\lena_color.png

?????文件????????492??2007-07-09?15:18??畢設\meanfilter.m

?????文件????????581??2007-07-09?15:18??畢設\medfilter.m

?????文件????????189??2007-07-02?09:22??畢設\MSE.m

?????文件????????281??2007-07-09?15:52??畢設\MSE_color.m

?????文件???????1051??2007-07-01?09:57??畢設\noise.m

?????文件????????101??2007-05-23?15:27??畢設\PSNR.m

?????文件????????113??2007-07-02?09:29??畢設\PSNR_color.m

?????文件???????1754??2007-06-28?21:27??畢設\wdenoise.m

?????文件???????1778??2007-07-09?16:03??畢設\文件清單及說明.txt

?????目錄??????????0??2009-09-15?16:45??畢設

-----------?---------??----------?-----??----

???????????????653225????????????????????23