分解.txt

来自「脑电信号EEG的信号波形分析」· 文本 代码 · 共 27 行

x=lvboxue';
s=x(1,160001:185000);
%%脑电有五个频率段，所以用小波四层分解
%%开始小波多尺度分解


% 计算第一组值
[c1,l1]=wavedec(s,4,'db4');
x1_d1=wrcoef('d',c1,l1,'db5',1);
x1_d2=wrcoef('d',c1,l1,'db5',2);
x1_d3=wrcoef('d',c1,l1,'db5',3);
x1_d4=wrcoef('d',c1,l1,'db5',4);
x1_a4=wrcoef('a',c1,l1,'db5',4);

figure(1)
subplot(611)
plot(s)
subplot(612)
plot(x1_a4)
subplot(613)
plot(x1_d4)
subplot(614)
plot(x1_d3)
subplot(615)
plot(x1_d2)
subplot(616)
plot(x1_d1)