Mh抽样 Matlab程序求助

hujf1987

我对MCMC----MH抽样理解不够透彻，拿位能帮帮我？？？？？？


% 系统重采样算法(已知xhat为1*N的行向量，存放待采样粒子，wx 1*N为对应的归一化了的权值)
    Px = cumsum(wx);                                                   
    ut(1)=rand(1)/N;            
    k = 1;                                                       %k记录的是需要进行重采样的粒子
    ix = zeros(1,N);
    for j = 1:N
        ut(j)=ut(1)+(j-1)/N;
        while(Px(k)<ut(j));                                 %比较采样尺度，进行重采样
              k = k + 1;
        end;
        ix(j) = k;
        wx(j)=1/N;                                         %wx(j)记录粒子权值， 每个权值都变为1/N
    end;
                        
    xhat = xhat(:,ix);                                                                        
    xp(t) = mean(xhat);     

以上程序是 用系统重采样算法实现的对粒子的重采样  ,谁能提供 用‘MCMC-----Metropolis Hastings’算法进行重采样的程序。其实没几句，费您点儿时间帮帮我吧，谢谢了