EEG小波变换

%利用小波变换的EEG信号的多分辨率分析
%读入EEG信号
%load eeg.edf;

Fs=500; %采样频率 
n=0:1/Fs:1;%读取edf信号并显示波形
[hdr eeg] = edfread（‘clo_reduced.edf‘）;%读取三分钟信号
s=eeg（3:）‘;
d=zeros（1120*500）;%选取2分钟信号
for i=1:120*500
    d（i）=s（i）;
end
d1=zeros（160*500）;%选取0-60s
for i=1:60*500
    d1（i）=d（i）;
end
d2=zeros（160*500）;%60-120s
for i=1:60*500
    d2（i）=d（i+60*500）;
end
d3=zeros（130*500）;%60-90s
for i=1:30*500
    d3（i）=d2（i）;
end
d4=zeros（130*500）;%90-120s
for i=1:30*500
   d4（i）=d2（i+30*500）; 
end
x=d4;
N=length（x）;
%采样频率
fs=500;
%小波分解
[cl]=wavedec（x7‘db4‘）;
%重构第6层逼近信号
a7=wrcoef（‘a‘cl‘db4‘7）;
% A=fft（a6）;
% plot（（1:N）/30abs（A））
% axis（[0 50 0 25000]）
%重构第1-6层细节信号
d7=wrcoef（‘d‘cl‘db4‘7）;
d6=wrcoef（‘d‘cl‘db4‘6）;
d5=wrcoef（‘d‘cl‘db4‘5）;
d4=wrcoef（‘d‘cl‘db4‘4）;
d3=wrcoef（‘d‘cl‘db4‘3）;
d