压缩感知BP算法

%一维信号BP重构算法
clc
clear all;
close all;
M = 200;   %观测信号长度
N = 1000;  %稀疏信号长度
K = 30;   %稀疏度
%% -----1.生成稀疏度为K的稀疏信号-----
x0 = zeros（N1）;
p = randperm（N）;
x0（p（1:K）1） = rand（K1）-0.5;
%------  高斯感知矩阵Phi   -------------
Phi = sqrt（1/M） * randn（MN）;
for i = 1:N
    Phi（:i） = Phi（:i） / norm（Phi（:i））;
end
%-------- 测量向量 y  ----------
y  = zeros（M1）;
y = Phi * x0;
%%  -----2. 含高斯白噪声观测 SNR=40dB --------
ye2 = mean（y.^2）;
SNR = 40;
sgmav = sqrt（ ye2*10^（-SNR/10） ）;
noisev = sgmav*randn（M1）;
y = y + noisev;
%% ------3.BP重构算法 ----------------