1 简介
2 部分代码
function A = construct_operator(T, rho, y)
if iscolumn(y) % 1d
L = length(y) - (T+1);
A = zeros(L+1,T+1);
for i = 1:T+1
% take the i-th cannonical basis vector
e = zeros(T+1, 1);
e(i) = 1;
% apply the operator
A(:,i) = direct_operator(e, rho, y);
end
elseif ismatrix(y) && prod(size(y) == size(y')) % 2d
L = length(y) - (T+1);
A = zeros((L+1)^2, (T+1)^2);
for i = 1:T+1
for j = 1:T+1
% take the (i,j)-th cannonical basis matrix
e = zeros(T+1, T+1);
e(i,j) = 1;
% apply the operator
Ae = direct_operator(e, rho, y);
A(:,i+(T+1)*(j-1)) = Ae(:);
end
end
else
error('Wrong input format');
end
end
3 仿真结果
4 参考文献
-
Adaptive Recovery of Signals by Convex Optimization Z. Harchaoui, A. Juditsky, A. Nemirovski, D. Ostrovskii. COLT 2015.
-
Structure-Blind Signal Recovery D. Ostrovskii, Z. Harchaoui, A. Judistky, A. Nemirovski. NIPS 2016.
-
Efficient First-Order Algorithms for Adaptive Signal Denoising D. Ostrovskii, Z. Harchaoui. ICML 2018.
博主简介:擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真,相关matlab代码问题可私信交流。
部分理论引用网络文献,若有侵权联系博主删除。