MATLAB小波变换分析ECG信号

%首先加载模拟的纯净ECG信号，并设置相关参数load('ecg.mat'); fs=360;%采样频率t=(0:(length(ecg)-1))/fs;%时间序列N=length(t);%信号长度引入运动伪影（motion artifacct）噪声信号load('ma.mat');去均值并合成带噪声信号ecg=ecg-mean(ecg); q=[ecg ecg]; ma=0.01; d=ecg+ma; %使用“db4”小波对带噪声ECG信号进行7层分解[C,L] = wavedec(d,7,'db4'); %重建每个子频带信号D1 = wrcoef('d',C,L,'db4',1); D2 = wrcoef('d',C,L,'db4',2); D3 = wrcoef('d',C,L,'db4',3); D4 = wrcoef('d',C,L,'db4',4); D5 = wrcoef('d',C,L,'db4',5); D6 = wrcoef('d',C,L,'db4',6); D7 = wrcoef('d',C,L,'db4',7);