压缩感知BP LASSO OMP STOMP算法内含完整代码，包括DWT.m、fdrthresh.m|FWT_PO.m、SolveBP.m等M文件，有详细注释

压缩感知BP LASSO OMP STOMP算法（内含完整代码，包括DWT.m、fdrthresh.m|FWT_PO.m、SolveBP.m等M文件，有详细注释）

%  构造正交小波变换矩阵，图像大小N*N，N=2^P，P是整数。

function ww=DWT（N）

[hg]= wfilters（‘sym8‘‘d‘）;       %  分解低通和高通滤波器‘d‘表示的是分解decomposition滤波器h是LO_D（低通滤波器组）g=HI_D（），

% N=256;                           %  矩阵维数（大小为2的整数幂次）
L=length（h）;                       %  滤波器长度
rank_max=log2（N）;                  %  最大层数
rank_min=double（int8（log2（L）））+1;  %  最小层数
ww=1;   %  预处理矩阵

%  矩阵构造
for jj=rank_min:rank_max
    
    nn=2^jj;
    
    %  构造向量
    p1_0=sparse（[hzeros（1nn-L）]）;
    p2_0=sparse（[gzeros（1nn-L）]）;
    
    %  向量圆周移位
    for ii=1:nn/2
        p1（ii:）=circshift（p1_0‘2*（ii-1））‘;
        p2（ii:）=circshift（p2_0‘2*（ii-1））‘;
    end
    
    %  构造正交矩阵
    w1=[p1;p2];
    mm=2^rank_max-length（w1）;
    w=sparse（[w1zeros（length（w1）mm）;zeros（mmlength（w1））eye（mmmm）]）;
    ww=ww*w;
    
    clear p1;clear p2;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-04-03 16:49  2 频域与空域高斯随机采样对比\
     文件         949  2020-03-21 14:28  2 频域与空域高斯随机采样对比\DWT.m
     文件         525  2006-07-30 14:12  2 频域与空域高斯随机采样对比\fdrthresh.m
     文件         993  2006-03-28 17:01  2 频域与空域高斯随机采样对比\FWT_PO.m
     文件       66614  2009-06-08 18:00  2 频域与空域高斯随机采样对比\lena256.bmp
     文件        3931  2014-04-22 15:13  2 频域与空域高斯随机采样对比\main_FD.m
     文件        4009  2014-04-03 16:25  2 频域与空域高斯随机采样对比\main_SD.m
     文件       54596  2006-07-30 14:13  2 频域与空域高斯随机采样对比\pdco.m
     文件       10186  2006-07-30 14:13  2 频域与空域高斯随机采样对比\pdcoSet.m
     文件         475  2006-07-30 14:18  2 频域与空域高斯随机采样对比\reverse.m
     文件        4483  2007-08-18 20:54  2 频域与空域高斯随机采样对比\SolveBP.m
     文件       10874  2006-08-22 10:54  2 频域与空域高斯随机采样对比\SolveLasso.m
     文件        5451  2006-07-30 14:11  2 频域与空域高斯随机采样对比\SolveOMP.m
     文件        4889  2006-07-30 14:11  2 频域与空域高斯随机采样对比\SolveStOMP.m
     文件        2759  2014-04-02 16:09  2 频域与空域高斯随机采样对比\Wavelet_OMP_WD.m