nrpshdfig1.m

来自「一种基于压缩感知技术的图像重建程序」· M 代码 · 共 94 行

% Plot Rho F (V-S threshold)

rhomax = 0.999;
rhomin = 0.005;
olddelta = 0.005;
FaceNeigh = [];
mesh_density=100;
for delta = linspace(0.005,0.999,mesh_density)
rhomax = min(0.999,rhomax + 5.* (delta - olddelta));  % Assume drho/ddelta < 5
rhomin = max(0.005,rhomin - 5*(delta-olddelta));
olddelta = delta;
errtol = .000001;
psitol = .0000001;
which=3;
refine=1;
refine_steps=1;
while( rhomax - rhomin > errtol),
    newrho = (rhomax + rhomin)/2;
    nu = linspace(delta, min(.9999999, delta+0.1),200);
    Psi = PsiNet(nu,newrho,delta,which);
    if refine==1
      for m=1:refine_steps
        j=1;
        while (Psi(j+1)>=Psi(j) & j+1<length(Psi)),
          j=j+1;
        end
        deltanu=nu(2)-nu(1);
        nu = linspace(max(delta,nu(j)-5*deltanu),min(0.9999,nu(j+1)+5*deltanu),200);
        Psi = PsiNet(nu,newrho,delta,which);
      end
    end
%      plot(nu,Psi)
%      [delta, rhomin, rhomax]
%      pause
    if all(Psi < -psitol),
        rhomin = newrho;
    else
        rhomax = newrho;
    end
%    [delta,rhomin,rhomax]
end
  delta
%  [delta, rhomin]
FaceNeigh = [FaceNeigh; delta rhomin ];
end

plot(FaceNeigh(:,1),FaceNeigh(:,2))
hold on

% Plot Rho N

% Solve for Rho_N

global newdelt
mesh_density=100;
delta = linspace(.0001,0.999,mesh_density);
Neigh = [];
for d = delta,
     newdelt = d;

%zone=linspace(0.01,.99,100);

     if d < .12,
         limits = [1/500,1/15];
     elseif d< .5,
         limits = [1/20,1/6];
     else 
         limits = [1/12,1/2];
     end

%plot(zone,GenRhoNDiff(zone));
%pause

limits=[0.01,0.99];

    rhostar = fzero(@GenRhoNDiff,limits);
    Neigh = [ Neigh ; d rhostar ];
    disp([d,rhostar])
end

plot(Neigh(:,1),Neigh(:,2),'r')

%
% Copyright (c) 2006. David Donoho
%  

%
% Part of SparseLab Version:100
% Created Tuesday March 28, 2006
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%