神经电压缩感知matlab

压缩感知源代码，内含重采样，压缩感知，BSBL 重构。特别好

clear;
ECG=load（‘ECG360.txt‘）;
fs = 360;
%t1 = 0:1/fs:1;
x = ECG（1:1600）;
xout = resample（x1000160）
%t2 = （0:（length（y）-1））*6/（5*fs）;
save（‘ECG.mat‘‘xout‘）;
CF=4;
NNN=CF*100;
zu=25;

load ECG.mat;
for ii=1:zu
xin=xout（NNN*（ii-1）+1:NNN*ii）;
%xin=dct（xin）;
if（xin==zeros（NNN1））
    yout（NNN*（ii-1）+1:NNN*ii）=zeros（NNN1）
else
%A=load（‘matrix.txt‘）;
A=BernoulliMtx（50NNN）;
%y=A*xin;
[rec]=DEMO_nonSparse（xinNNN50A）;
%rec=idct（rec_1）;
%eval（[‘save date_‘num2str（i）]）
yout（NNN*（ii-1）+1:NNN*ii）=rec;
end
end
xout1=xout（1:zu*NNN）;
save （‘recdata.mat‘‘yout‘）;
figure（1）;
subplot（211）
plot（xout1）;
xlabel（‘Time / s‘）; ylabel（‘Voltage / mV‘）;
string=[‘x[n] ‘‘before‘];
title（string）

subplot（212）
 plot（yout）;
xlabel（‘Time / s‘）; ylabel（‘Voltage / mV‘）;
string=[‘x[n] ‘‘after‘];
title（string）
disp（[‘With DCT SNDR     :‘num2str（20*log10（（norm（xout1）^2）/（norm（xout1-yout‘）^2）））]）;
disp（[‘with DCT PRD   :‘num2str（100*（norm（xout1-yout‘）/norm（xout1）））]）;
%disp（[‘with DCT ESR   :‘num2str（sum（（xout1-yout）.^2）/sum（xout1.^2））]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1092  2018-10-12 09:20  neuraltest.m