压缩感知CS最全matlab程序，二维三维图像处理，恢复算法

压缩感知CS最全matlab程序，二维三维图像处理，恢复算法，CSDN中最全的压缩感知程序包，12M的程序，也包含陶哲轩的关于压缩感知的PPT及文章

%   感谢亲亲使用此代码，此代码解决您的问题了吗~（@^_^@）~
%   没解决的话告诉亲亲一个好消息，我这里提供分钱成品代码（′▽‘〃）哦~登录淘宝店铺“大成软件工作室”便可领取
%   是的，亲亲真的没有看错，挠破头皮的问题真的1分钱就可以解决了\（^o^）/YES!
%   小的这就把传送门给您，记得要收藏好哦（づ￣3￣）づ╭～
%   传送门：https://shop305228693.taobao.com/?spm=a230r.7195193.1997079397.2.yVBY8j
%   如果传送门失效，亲亲可以来店铺讨要，客服MM等亲亲来骚扰哦~（*/ω╲*）

function [theta posterior bound]=BCSvb（Phi v hyperpara niter tol plotflag）
%BCSvb: Bayesian CS inversion implemented by variational Bayes.
%USAGE: [theta posterior bound]=BCSvb（Phi v hyperpara niter tol plotflag）
%INPUT （number in [] means default value）:  
%   Phi: N x M CS projection matrix
%   v: N x 1 CS observation
%   hyperpara: VB hyperparameters
%       hyperpara.a0: scalar hyperparameter 1 for noise precision [1e-6]
%       hyperpara.b0: scalar hyperparameter 2 for noise precision [1e-6]
%       hyperpara.c0: scalar hyperparameter 1 for non-zero coefficients precision [1e-6] 
%       hyperpara.d0: scalar hyperparameter 2 for non-zero coefficients precision [1e-6] 
%   niter: scalar maximum number of VB iterations [100]
%   tol: scalar tolarance of relative change of VB lower bound for two consecutive steps [1e-4] 
%   plotflag: indicator one means plotting result for each iteration zero means no plot [0] 
%OUTPUT:
%   theta: M x 1 mean of reconstructed signal
%   posterior: hyperparameter for posterior distribution
%       posterior.MUtheta: M x 1 mean of inverted theta 
%       posterior.VARtheta_diag: M x 1 variance of inverted theta （for each coefficient）
%       posterior.a: scalar posterior hyperparameter 1 for noise precision
%       posterior.b: scalar posterior hyperparameter 2 for noise precision
%       posterior.c: M x 1 posterior hyperparameter 1 for coefficient precision
%       posterior.d: M x 1 posterior hyperparameter 2 for coefficient precision
%   bound: 1 x （number of VB iterations） VB lower bound for each iteration

%--------------------------------------------------------------------------
% References:
% S.Ji Y.Xue and L.Carin “Bayesian Compressive Sesning“ （2007）.
% L.He and L.Carin “Exploiting Structure in Wavelet-based Bayesian Compressive Sensing“ （2008）  
%
% Lihan He ECE Duke University
% Created: Nov. 21 2008
% Last change: Mar. 3 2009 allowing [] for input arguments （using default value）
%--------------------------------------------------------------------------

% ---------------------
% check input arguments
% ---------------------

if nargin<6 plotflag=0; end
if nargin<5 tol=1e-4; end
if nargin<4 niter=100; end
if nargin<3
    hyperpara.a0=1e-6;
    hyperpara.b0=1e-6;
    hyperpara.c0=1e-6;
    hyperpara.d0=1e-6;
end

if isempty（hyperpara）
    hyperpara.a0=1e-6;
    hyperpara.b0=1e-6;
    hyperpara.c0=1e-6;
    hyperpara.d0=1e-6;
end
if isempty（niter） niter=100; end
if isempty（tol） tol=1e-4; end
if isempty（plotflag） plotflag=0; end

a0=hyperpara.a0;
b0=hyperpara.b0;
c0=hyperpara.c0;
d0=hyperpara.d0;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5269  2016-09-24 14:10  BCS implemented\BCSvb.m

     文件      65551  2002-12-03 07:45  BCS implemented\cameraman.pgm

     文件       2629  2016-09-24 14:10  BCS implemented\demo_BCSvb_1d.m

     文件       2099  2016-09-24 14:10  BCS implemented\demo_BCSvb_2d.m

     文件       2112  2016-09-24 14:11  BCS implemented\demo_BCSvb_2d_small.m

     文件      65610  2002-12-03 04:04  BCS implemented\indor2.pgm

     文件     107698  2012-07-15 11:18  compressive sensing basic knowledge\anziam_pics.pdf

     文件    8367240  2012-07-15 11:20  compressive sensing basic knowledge\baraniuk-IMA-CScamera-june07.pdf

     文件    1028107  2012-07-15 11:19  compressive sensing basic knowledge\compressed-sensing1.pdf

     文件     284802  2012-07-15 11:18  compressive sensing basic knowledge\Compressive_Sensing.pdf

     文件     194495  2012-07-15 11:18  compressive sensing basic knowledge\TerryTao.pdf

     文件      80896  2013-02-27 13:30  CS（matlab）\分块压缩感知程序（matlab）.doc

     文件       1557  2014-04-27 20:48  CS_Examples\CS_OMP.m

     文件       1220  2014-04-27 20:45  CS_Examples\CS_SPGL1.m

     文件       3861  2009-05-20 19:22  CS_Examples\spgl1_1_7\ChangeLog

     文件        697  2009-05-20 19:22  CS_Examples\spgl1_1_7\Contents.m

     文件      26436  2009-05-20 19:22  CS_Examples\spgl1_1_7\COPYING

     文件        184  2009-05-20 19:22  CS_Examples\spgl1_1_7\NormGroupL2_dual.m

     文件        176  2009-05-20 19:22  CS_Examples\spgl1_1_7\NormGroupL2_primal.m

     文件        375  2009-05-20 19:22  CS_Examples\spgl1_1_7\NormGroupL2_project.m

     文件        221  2009-05-20 19:22  CS_Examples\spgl1_1_7\NormL12_dual.m

     文件        209  2009-05-20 19:22  CS_Examples\spgl1_1_7\NormL12_primal.m

     文件        463  2009-05-20 19:22  CS_Examples\spgl1_1_7\NormL12_project.m

     文件         63  2009-05-20 19:22  CS_Examples\spgl1_1_7\NormL1_dual.m

     文件         63  2009-05-20 19:22  CS_Examples\spgl1_1_7\NormL1_primal.m

     文件        227  2009-05-20 19:22  CS_Examples\spgl1_1_7\NormL1_project.m

     文件       1561  2009-05-20 19:22  CS_Examples\spgl1_1_7\private\ensure.m

     文件       6398  2009-05-20 19:22  CS_Examples\spgl1_1_7\private\heap.c

     文件       3660  2009-05-20 19:22  CS_Examples\spgl1_1_7\private\heap.h

     文件      11849  2009-05-20 19:22  CS_Examples\spgl1_1_7\private\lsqr.m

............此处省略102个文件信息