gensol.m

来自「这是在网上下的一个东东」· M 代码 · 共 13 行

%
% soln=gensol(G,lower,upper,a,beta)
%
% This routine generates a random solution according to the posterior
% pdf.  
%
function soln=gensol(G,lower,upper,a,beta)
[m n]=size(G);
soln=zeros(n,1);
for i=1:n,
  soln(i)=lower(i)+qinv(rand,a(i),upper(i)-lower(i),beta(i));
end;