如何构造测量矩阵

最新推荐文章于 2022-04-22 10:35:32 发布

童话WN

最新推荐文章于 2022-04-22 10:35:32 发布

阅读量977

点赞数

分类专栏：压缩感知

本文链接：https://blog.csdn.net/qq_36596540/article/details/105645921

版权

本文分享了由Thong Do编写的MATLAB代码，用于构造测量矩阵，涉及IDWT（逆离散小波变换）和随机块基感知算子。通过示例展示了如何将离散小波变换的稀疏向量转换为样本，并进行块DCT和块Hadamard变换。此外，还提供了相关函数的反向操作，以从测量中重构信号。

摘要由CSDN通过智能技术生成

今天找到一个测量矩阵的构造方法，废话不多说，直接上代码：
% Written by Thong Do, JHU, 2007
% modified, Feb. 2008

function b = A_idwt1d(A,alpha,f0,f1,L,m,n)

% alpha: sparse vector which is taken from Wavelet transform of image x
% A: random projection matrix K x N
% f0, f1: lowpass and highpass filter at the synthesis part
% L: level of decomposition
% m, n: size of image
% per : take per percent of wavelet coefficients
% Return b: vector K x 1

% converting wavelet coefficients (sparse representation) into samples
u = reshape(alpha,m,n);
im = idwt2d(u,f0,f1,L);

% converting samples into measurements
x = im(😃;
if ~isa(A, ‘function_handle’)
b = A*x;
else
b = A(x);
end

% Written by Thong Do, JHU, 2007
% modified, Feb. 2008

function alpha = At_dwt1d(At,b,h0,h1,L,m,n)
% alpha: sparse vector which is taken fr