資源簡介
實現對符合廣義高斯分布樣本的廣義高斯參數的估計,利用Newton–Raphson迭代方法實現參數的數值解。
代碼片段和文件信息
function?[theta?sigma]?=?ggd_fit(x?theta0)
%
%[theta?sigma]?=?ggd_fit(x?theta0);
%?根據樣本估計單維廣義高斯分布的參數值
%?x是樣本向量
%?theta0是形狀參數的初始值
%?theta是迭代完成后的形狀參數值
%?sigma是迭代完成后與方差有關的參數值
%?該程序是根據文獻????A?globally?convergent?and?consistent?method?for?estimating?the?shape?parameter?of?a?generalized?Gaussian?distribution
%?編寫而成,在原作者Lingyun?Zhang?基礎上修改某些部分?增強了程序的魯棒性。
%?number?of?sample?points
x?=?x(find(abs(x)>0.001));
n?=?length(x);
%?estimate?the?shape?parameter
theta?=?theta0;
i=1;
T(1)=theta0;
while?(1)
????Y1?=?mean(abs(x).^theta);
????Y2?=?mean(abs(x).^(2*theta));
????Z?=?Y2?/?Y1?/?Y1?-?(theta+
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????1330??2014-12-04?14:59??ggd_fit.m
-----------?---------??----------?-----??----
?????????????????1330????????????????????1
- 上一篇:一種簡易發射機電路圖
- 下一篇:Qt實現notepad記事本源碼工程(版本三
評論
共有 條評論