matlab脑电信号处理.docxVIP

matlab脑电信号处理

matlab脑电信号处理

t=0.001:0.001:1;

x=load(#39;C:\Users\yxzhang\Desktop\rest_close.txt#39;);%读取文件

y=load(#39;C:\Users\yxzhang\Desktop\audio_close.txt#39;);

xx={};%每个导联的数据存储

yy={};

n=1000;%数据数目

sc=7;%小波包的分解尺度

fori=1:1:8%导联的数据分离

xx{i}=x(:,i);

yy{i}=y(:,i);

end

fori=1:1:8

%画出原始信号图像

figure

subplot(2,2,1)

plot(t,xx{i})

axis([min(t)max(t)1.1*floor(min(xx{i}))1.1*ceil(max(xx{i}))])

title(#39;restclose原始信号#39;)

ylabel(#39;幅值#39;)

subplot(2,2,2)

plot(t,yy{i})

axis([min(t)max(t)1.1*floor(min(yy{i}))1.1*ceil(max(yy{i}))])

title(#39;audioclose原始信号#39;)

ylabel(#39;幅值#39;)

%fft_原始信号的频谱分析

xx1=fft(xx{i},n);

pxx1=xx1.*conj(xx1)/n;

yy1=fft(yy{i},n);

pyy1=yy1.*conj(yy1)/n;

%画出0-30hz内的功率谱图像

n=60;

f=1:n/2;

subplot(2,2,3);

plot(f,pxx1(1:n/2));

ylabel(#39;功率谱幅值（mv^2)#39;);

title(#39;rest_close功率谱#39;)

subplot(2,2,4);

plot(f,pyy1(1:n/2));

ylabel(#39;功率谱幅值（mv^2)#39;);

title(#39;audio_close的功率谱#39;)

%分解信号选择分解尺度为7，同时重构信号

wpt=wpdec(xx{i},sc,#39;db7#39;,#39;shannon#39;);%小波包分解信号

xx80=wprcoef(wpt,[sc,0]);%重构信号

xx81=wprcoef(wpt,[sc,1]);

xx82=wprcoef(wpt,[sc,2]);

xx83=wprcoef(wpt,[sc,3]);

wpt=wpdec(yy{i},8,#39;db7#39;,#39;shannon#39;);

yy80=wprcoef(wpt,[sc,0]);

yy81=wprcoef(wpt,[sc,1]);

yy82=wprcoef(wpt,[sc,2]);

yy83=wprcoef(wpt,[sc,3]);

%画出重构信号

figure

subplot(2,1,2);plot(yy80);

title(#39;audioclosedelta#39;);

ylabel(#39;幅值#39;);

subplot(2,1,1);plot(xx80);

title(#39;restclosedelta#39;);

ylabel(#39;幅值#39;);

figure

subplot(2,1,1);plot(xx81);

title(#39;restclosetheta#39;);

ylabel(#39;幅值#39;);

subplot(2,1,2);plot(yy81);

title(#39;audioclosetheta#39;);

ylabel(#39;幅值#39;);

figure

subplot(2,1,1);plot(xx82);

title(#39;restclosealpha#39;);

ylabel(#39;幅值#39;);

subplot(2,1,2);plot(yy82);

title(#39;audioclosealphta#39;);

ylabel(#39;幅值#39;);

figure

subplot(2,1,1);plot(xx83);

title(#39;restclosebeta#39;);

ylabel(#39;幅值#39;)

subplot(2,1,2);plot(yy83);

tit