hdcpnpdfig11.m

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

% Figure 8.2b


gamma = linspace(.03,.5);
sgamma = SeriesMillsSolve(1-gamma);
ygamma = yFroms(sgamma);
Xi = -sgamma.^2 .* (gamma ./(1-gamma)) ./2  ...
                - log(2/pi)/2 ...
                + log( -sgamma ) - log(1 - gamma) ;
diffxinu = diff(Xi,2) ./(gamma(2)-gamma(1)).^2;
plot(gamma(1:length(diffxinu)),diffxinu);
title('\xi_\gamma''''(y_\gamma)')
xlabel('\gamma')

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