对脑电信号使用小波分解技术重建

%利用小波变换的EEG信号的多分辨率分析
%读入EEG信号
%load eeg.dat;
t=0:1/512:1;
t=t‘;
x=sin（t）;
%x=eeg（21:512）‘;
N=length（x）;
%采样频率
fs=512;
%小波分解
[cl]=wavedec（x5‘db5‘）;
%重构第5层逼近信号
a5=wrcoef（‘a‘cl‘db5‘5）;
%重构第1-5层细节信号
d5=wrcoef（‘d‘cl‘db5‘5）;
d4=wrcoef（‘d‘cl‘db5‘4）;
d3=wrcoef（‘d‘cl‘db5‘3）;
d2=wrcoef（‘d‘cl‘db5‘2）;
d1=wrcoef（‘d‘cl‘db5‘1）;
%显示各层信号
figure（1）;
subplot（611）;plot（a5‘LineWidth‘2）;ylabel（‘a5‘）;
subplot（612）;plot（d5‘LineWidth‘2）;ylabel（‘d