%利用小波變換的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