OLD_MATLAB/InputParameters/grillelambda.m

   1 function[gridLambda] = grillelambda(phiInit,rhoInit,piInit,gamInit,X,Y,gamma,mini,maxi,tau)
   2
   3     n = size(X, 1);
   4     [p,m,k] = size(phiInit);
   5     [phi,rho,pi,~,S] = EMGLLF(phiInit,rhoInit,piInit,gamInit,mini,maxi,1,0,X,Y,tau);
   6
   7     gridLambda = zeros(p,m,k);
   8     for j=1:p
   9         for mm=1:m
  10             gridLambda(j,mm,:) = abs(reshape(S(j,mm,:),[k,1])) ./ (n*pi.^gamma);
  11         end
  12     end
  13
  14     gridLambda = unique(gridLambda);
  15     gridLambda(gridLambda()>1) = [];
  16
  17 end