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