如何把脑电信号转换成功率谱?matlab脑电信号特征提取程序|天天速讯
编写一程序。。驾驶员清醒和疲劳状态时脑电信号特征提取,只需要把脑电信号转换成功率谱
(资料图片)
脑电信号已采集
t=0.001:0.001:1;
x=load("D:\EEGdata_data"); %读取文件
y=load("D:EEGdata_newdata\");
xx={}; %每个导联的数据存储
yy={};
n=1000; %数据数目
sc=7; %小波包的分解尺度
for i=1:1:8 %导联的数据分离
xx{i}=x(:,i);
yy{i}=y(:,i);
end
for i=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("rest close 原始信号")
ylabel("幅值")
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("audio close 原始信号")
ylabel("幅值")
%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("功率谱幅值(mv^2)");
title("rest_close功率谱")
subplot(2,2,4);
plot(f,pyy1(1:n/2));
ylabel("功率谱幅值(mv^2)");
title("audio_close的功率谱")
%分解信号选择分解尺度为7,同时重构信号
wpt=wpdec(xx{i},sc,"db7","shannon"); %小波包分解信号
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,"db7","shannon");
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("audio close delta");
ylabel("幅值");
subplot(2,1,1);plot(xx80);
title("rest close delta");
ylabel("幅值");
figure
subplot(2,1,1);plot(xx81);
title("rest close theta");
ylabel("幅值");
subplot(2,1,2);plot(yy81);
title("audio close theta");
ylabel("幅值");
figure
subplot(2,1,1);plot(xx82);
title("rest close alpha");
ylabel("幅值");
subplot(2,1,2);plot(yy82);
title("audio close alphta");
ylabel("幅值");
figure
subplot(2,1,1);plot(xx83);
title("rest close beta");
ylabel("幅值")
subplot(2,1,2);plot(yy83);
title("audio close beta");
ylabel("幅值");
n=1000;
%fft_重构信号的频谱分析
xx180=fft(xx80,n);
pxx180=xx180.*conj(xx180)/n;
xx181=fft(xx81,n);
pxx181=xx181.*conj(xx181)/n;
xx182=fft(xx82,n);
pxx182=xx182.*conj(xx182)/n;
xx183=fft(xx83,n);
pxx183=xx183.*conj(xx183)/n;
yy180=fft(yy80,n);
pyy180=yy180.*conj(yy180)/n;
yy181=fft(yy81,n);
pyy181=yy181.*conj(yy181)/n;
yy182=fft(yy82,n);
pyy182=yy182.*conj(yy182)/n;
yy183=fft(yy83,n);
pyy183=yy183.*conj(yy183)/n;
%画出重构信号的功率谱图
n=60;
f=1:n/2;
figure
subplot(2,1,1);
plot(f,pyy180(1:n/2));
ylabel("功率谱幅值(mv^2)");
title("rest close delta的功率谱");
subplot(2,1,2);
plot(f,pyy180(1:n/2));
ylabel("功率谱幅值(mv^2)");
title("audio close detta的功率谱");
figure
subplot(2,1,1);
plot(f,pxx181(1:n/2));
ylabel("功率谱幅值(mv^2)");
title("rest close theta的功率谱");
subplot(2,1,2);
plot(f,pyy181(1:n/2));
ylabel("功率谱幅值(mv^2)");
title("audio close theta的功率谱");
figure
subplot(2,1,1);
plot(f,pxx182(1:n/2));
ylabel("功率谱幅值(mv^2)");
title("rest close alpha的功率谱");
subplot(2,1,2);
plot(f,pyy182(1:n/2));
ylabel("功率谱幅值(mv^2)");
title("audio close alpha的功率谱");
figure
subplot(2,1,1);
plot(f,pxx183(1:n/2));
ylabel("功率谱幅值(mv^2)");
title("rest close beta的功率谱");
subplot(2,1,2);
plot(f,pyy183(1:n/2));
ylabel("功率谱幅值(mv^2)");
title("audio close beta的功率谱");
end;
运行不粗来,得不到功率谱图.........TAT,求帮助!
标签:
相关推荐:
最新新闻:
- 如何把脑电信号转换成功率谱?matlab脑电信号特征提取程序|天天速讯
- 全球热文:如何把真实的文件隐藏到一张图背后?内涵图制作
- 【技术】基于SIP协议的语音通信技术
- 有没有MAC的WPS?金山公司明年将有全平台的WPS:全球热门
- 今日看点:安卓手机如何设置便携式WLAN热点?设置教程
- 利勃海尔冰箱怎么样?利勃海尔冰箱有哪些优势?_全球速看料
- Trinity是什么?AMDTrinityAPU最新发布
- 怎么把Word转PDF格式?转换方法步骤
- 三星S6358怎么刷机?三星S6358移动定制版完整刷机教程|实时焦点
- 华为g9怎么样?华为g9参数配置评测
- 渲染农场:分布式并行集群计算系统
- win10系统显示不了桌面怎么回事?解决步骤如下_今日精选
- excel/word如何快速输入当前的时间和日期?输入方法技巧
- SCM供应链管理系统为企业带来哪些益处?SCM供应链管理系统介绍
- office怎么卸载?office卸载步骤-天天时讯
- 微信聊天记录怎么导出?微信聊天记录导出方法 环球观天下
- 天天头条:dr.dre耳机怎么样?dr.dre耳机简介及价格
- 数据透视表是什么?excel中数据透视表的基本概念 世界今热点
- 全球消息!新倩女幽魂哪个区人最多?2020倩女幽魂哪个区人多?
- 电脑扫描文件怎么弄?微信一物一码怎么开通? 环球关注
- Word转PDF用什么软件?全能PDF转换助手-全球观点
- 热门看点:Win8如何开启802.1x网络认证?Win8设置和修改系统电源方案的方法
- 数控铣床法兰克系统代码完整版 数控铣床法兰克系统代码大全
- 每日精选:网盘哪个好?网盘性能推荐
- 环球观天下!华为P50有哪些配置?华为P50参数配置介绍
- 【全球报资讯】荣耀MagicBook14用什么处理器?MagicBook14处理器介绍
- 天天安卓模拟器是什么?安卓模拟器电脑版下载
- n95手机怎么样?n95手机详细测评:今日讯
- 即时焦点:一文读懂程序化交易、算法交易、量化投资、高频交易
- 重磅!2022年第四届全国高校计算机能力挑战赛决赛晋级公示:世界看热讯
- 魔兽世界怀旧服:炼金术(350)技术介绍
- 环球报道:java基础知识:标签和target的属性
- 焦点关注:编辑好正文和引用内容的文本如何添加?
- 什么是沙盒游戏?沙盒类游戏哪十大排行?
- C++/MFC串口通信——光源控制器控制
- Linux下常见的锁机制-自旋锁 自旋锁spinlock见解
- typec数据线是什么?Type-C数据线有哪些优势?:每日播报
- 奖学金名称都有哪些?奖学金名称中英文对照表:焦点短讯
- 世界观察:苹果care是什么意思?苹果care的含义
- 华为matebook和荣耀magicbook有什么区别?两者区别介绍
- 小米手机如何刷机?小米手机刷机教程步骤
- 打印机硒鼓的寿命一般多久?硒鼓多久换一次?
- 开机后桌面显示慢怎么办? 电脑桌面刷新很慢解决方法
- PSAM部分的硬件设计 PSAM部分的硬件设计图
- 双通道内存有什么好处?双通道内存好处介绍
- 天天滚动:梅西C罗姆巴佩等众星发文缅怀球王:安息吧!贝利
- 环球快看点丨什么是静电耳机? 静电耳机和普通耳机什么区别?
- 库尔勒经济开发区:电子政务“一站式”服务平台 焦点速看
- cmake-3.19.6-Linux自解压程序 简单的自解压程序|天天动态
- 被覆盖的文件如何恢复?不小心把文件覆盖了怎么恢复?
- 图省事、紧裤腰,65亿用户规模,影视OST怎么失灵了?|全球观点
- 被覆盖的文件如何恢复?不小心把文件覆盖了怎么恢复?
- 防火墙ddos防御设置 windows防火墙ddos
- IE浏览器在电脑哪里?查找IE浏览器的方法
- 清空浏览器缓存怎么操作?清理浏览器历史浏览记录方法
- 两寸证件照尺寸是多少?二寸证件照快速制作步骤
- 如何设置开机u盘启动?开机优先启动u盘
- kdh是什么文件?kdh文件怎么转换成word文件
- 什么是dns服务器?设置DNS提升网速方法
- 云手机是什么东西?云手机用途
- 手机没有信号无服务是怎么回事?手机没信号解决思路
- 电脑虚拟内存不足怎么解决?教你彻底清理C盘垃圾
- 打印机墨水怎么加?佳能型号打印机加墨水教程
- 固态硬盘怎么安装?固态硬盘安装后需要怎么设置?
- 摄像机像素和分辨率区别是什么?摄像机像素和分辨率区别介绍
- 主板温度多少算正常? 主板温度高怎么办?
- 联想最新笔记本系统如何恢复?一键还原的使用方法
- 来了!《潜行者2:切尔诺贝利之心》首个实机演示公开 2023年发售-焦点速讯
- “球王”贝利去世 终年82岁 当前热闻
- Xbox金会员1月会免游戏公布 《彩虹坠入》等入库|今日观点
- 今日精选:疫情下的农村老人
- 傲意科技完成近亿元Pre-B轮融资,智慧互联产业基金、中原前海基金领投
- 封杀半年后 俄罗斯用户可以下载Windows 11了
- 前沿热点:三星Galaxy Book2 Pro 360发布 搭载高通骁龙 8cx Gen 3
- 无需越狱 大神利用iOS16漏洞成功修改系统字体 环球消息