MCMC的matlab源代码

MCMC的matlab源代码，非常适合用于多目标跟踪。 希望对大家有所帮助

function x = gengamma（alpha beta）

if （alpha==1）
   x = -log（1-rand（11））/beta;
   return
end
flag=0;       % teste si alpha<1 ou alpha>1
if （alpha<1）
   flag=1;
   alpha=alpha+1;
end
gamma=alpha-1;
eta=sqrt（2.0*alpha-1.0）;
c=.5-atan（gamma/eta）/pi;
aux=-.5;
while（aux<0）
   y=-.5;
   while（y<=0）
      u=rand（11）;  
      y = gamma + eta * tan（pi*（u-c）+c-.5）;
   end
   v=-log（rand（11））;
   aux=v+log（1.0+（（y-gamma）/eta）^2）+gamma*log（y/gamma）-y+gamma;
end;

if （flag==1）
   x = y/beta*（rand（1））^（1.0/（alpha-1））;
else
   x = y/beta;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     目录          0  2007-03-29 20:28  MCMC

     文件        539  2007-03-29 20:27  MCMC\gengamma.m

     文件       2277  2007-03-29 20:26  MCMC\radialBirth.m

     文件       1988  2007-03-29 20:26  MCMC\radialDeath.m

     文件       3570  2007-03-29 20:26  MCMC\radialMerge.m

     文件       2629  2007-03-29 20:26  MCMC\radialRW.m

     文件       3115  2007-03-29 20:27  MCMC\radialSplit.m

     文件       2638  2007-03-29 20:26  MCMC\radialUpdate.m

     文件        301  2007-03-29 20:27  MCMC\rjCubic.m

     文件       4658  2007-03-29 20:28  MCMC\rjdemo1.asv

     文件       4658  2007-03-29 20:28  MCMC\rjdemo1.m

     文件        396  2007-03-29 20:27  MCMC\rjGaussian.m

     文件        315  2007-03-29 20:27  MCMC\rjMultiquadric.m

     文件       8225  2007-03-29 20:27  MCMC\rjnn.m

     文件        300  2007-03-29 20:27  MCMC\rjtpSpline.m

     文件        218  2007-01-22 19:53  www.pudn.com.txt

----------- ---------  ---------- -----  ----

                35827                    16