• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

matlab计算样本熵

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

计算14通道得脑电数据吗,将得出的样本熵插入Excel表格

a = zeros(1,14);
b = a';
  for i =1:14
      b(i) = SampEn(d1_1(i,1:3000),2,0.2*std(d1_1(i,1:3000)));
  end
xlswrite('C:\Users\25626\Desktop\滤波后数据\14\样本熵\样本熵.xlsx',b,'Sheet1','J');%数据写入A列

  调用的SampEn函数

function SampEnVal = SampEn(data, m, r)
% SampEn  计算时间序列data的样本熵
% 输入:data是数据一维行向量
%      m重构维数,一般选择1或2,优先选择2,一般不取m>2
%      r 阈值大小,一般选择r=0.1~0.25*Std(data) 数据偏差
% 输出:SampEnVal样本熵值大小

data = data(:)';
N = length(data);
Nkx1 = 0;
Nkx2 = 0;

for k = N - m:-1:1
    x1(k, :) = data(k:k + m - 1);
    x2(k, :) = data(k:k + m);
end

for k = N - m:-1:1
    x1temprow = x1(k, :);
    x1temp    = ones(N - m, 1)*x1temprow;   
    dx1(k, :) = max(abs(x1temp - x1), [], 2)';   
    Nkx1 = Nkx1 + (sum(dx1(k, :) < r) - 1)/(N - m - 1);    
    x2temprow = x2(k, :);
    x2temp    = ones(N - m, 1)*x2temprow;
    dx2(k, :) = max(abs(x2temp - x2), [], 2)';
    Nkx2      = Nkx2 + (sum(dx2(k, :) < r) - 1)/(N - m - 1);
end
Bmx1 = Nkx1/(N - m);
Bmx2 = Nkx2/(N - m);
SampEnVal = -log(Bmx2/Bmx1);
end

  


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
MATLAB描绘极坐标图像——polar发布时间:2022-07-22
下一篇:
matlab图形handle发布时间:2022-07-22
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap