clear; clc;

%% 参数设置
global N M lambda tau iter_max method trail_times epi extend_target K;

N = 64;
M = 4;
K = 10;
lambda = zeros(M * N, 1);
lambda(:) = 0.3;
tau = 1e-4;
iter_max = 500;
% method = "debiased_LASSO";
% method = "LASSO";
method = "BP";
% method = "VAMP";
trail_times = 1e2;
epi = 0;
extend_target = 1;

%% FAR
P_fa_FAR = FAR_simu(3);
% fprintf("FAR: %f\n", P_fa_FAR);

%% PD
P_fa_PD = PD_simu(1);
fprintf("PD: %f\n", P_fa_PD);