原创力文档- 1、本文档共7页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
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
文档评论（0）