plotstats.m.txt

这是马尔可夫-蒙特卡罗算法的MATLAB源程序.

function PlotStats(llvec,lpvec,tlP,tlL,count,countmax)
%
% plot the log Likelihood and log Prior

%lag = 1000; % final plot
lag=50; % real-time display

figure(2),subplot(2,1,2)
plot(1:count,lpvec(1:count),'b',[0 countmax],[tlP tlP],'r:'), ylabel('log Prior')
axis([1 countmax min([lpvec(max(1,count-lag):count) tlP])-10 max([lpvec(max(1,count-lag):count) tlP])+10])

figure(2),subplot(2,1,1)
plot(1:count,llvec(1:count),'b',[0 countmax],[tlL tlL],'r:'), ylabel('log Likelihood')
axis([1 countmax min([llvec(max(1,count-lag):count) tlL])-10 max([llvec(max(1,count-lag):count) tlL])+10])

drawnow