Matlab 广义高斯分布建模

供参考，用matlab实现广义高斯分布建模，附这一方法的入门文献

% 一种广义高斯分布的参数快速估计法（http://wenku.baidu.com/view/02bde6e9856a561252d36f9f.html）
%代码实现
clear;close all;clc;

x=-20:0.4:20;

%1.画GGD图象
%u代表均值，o代表方差，a形状参数，b尺度参数
u=0;
 %a 的值决定了图象的形状
a=0.7;o=10; 

b=o.*（gamma（1/a）/gamma（3/a））.^0.5;
f=（a/（2*b*gamma（1/a）））.*exp（-（abs（x-u）/b）.^a）;
plot（xf）;

% 2.画R（x）图象
for n=1:length（x）
    r1（n）=gamma（2/x（n））.*gamma（2/x（n））;
    r2（n）=gamma（1/x（n））.*gamma（3/x（n））; 
end
r=r1./r2;
plot（xr）;

%3.拟合R（x）
cftool（xr）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      334225  2011-05-06 15:46  一种广义高斯分布的参数快速估计法.pdf
     文件         527  2011-05-06 15:39  ggd.m