%一維信號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重構算法?----------------