gensol.m

基于李雅普诺夫的睡眠脑电动力学分析 核心算法程序。

%
% 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;