样本熵的matlab程序

转自http://www.ilovematlab.cn/thread-89427-1-1.html MATLAB下的动态样本熵计算




  1. function SampEnVal = SampEn(data, m, r)
  2. %SAMPEN  计算时间序列data的样本熵
  3. %        data为输入数据序列
  4. %        m为初始分段,每段的数据长度
  5. %        r为阈值
  6. % $Author: lskyp
  7. % $Date:   2010.6.20
  8. % Orig Version: V1.0--------分开计算长度为m的序列和长度为m+1的序列
  9. %                           这一版的计算有些问题,需要注意两个序列总数都要为N-m
  10. % Modi Version: V1.1--------综合计算,计算距离时通过矩阵减法完成,避免重循环
  11. % V1.1 Modified date: 2010.6.23
  12. data = data(:)';
  13. N = length(data);
  14. Nkx1 = 0;
  15. Nkx2 = 0;
  16. % 分段计算距离,x1为长度为m的序列,x2为长度为m+1的序列
  17. for k = N - m:-1:1
  18.     x1(k, :) = data(k:k + m - 1);
  19.     x2(k, :) = data(k:k + m);
  20. end
  21. for k = N - m:-1:1
  22.     % x1序列计算
  23.     % 统计距离,由于每行都要与其他行做减法,因此可以先将该行复制为N-m的矩阵,然后
  24.     % 与原始x1矩阵做减法,可以避免两重循环,增加效率
  25.     x1temprow = x1(k, :);
  26.     x1temp    = ones(N - m, 1)*x1temprow;
  27.     % 可以使用repmat函数完成上面的语句,即
  28.     % x1temp = repmat(x1temprow, N - m, 1);
  29.     % 但是效率不如上面的矩阵乘法
  30.     % 计算距离,每一行元素相减的最大值为距离
  31.     dx1(k, :) = max(abs(x1temp - x1), [], 2)';
  32.     % 模板匹配数
  33.     Nkx1 = Nkx1 + (sum(dx1(k, :) < r) - 1)/(N - m - 1);
  34.     
  35.     % x2序列计算,和x1同样方法
  36.     x2temprow = x2(k, :);
  37.     x2temp    = ones(N - m, 1)*x2temprow;
  38.     dx2(k, :) = max(abs(x2temp - x2), [], 2)';
  39.     Nkx2      = Nkx2 + (sum(dx2(k, :) < r) - 1)/(N - m - 1);
  40. end
  41. % 平均值
  42. Bmx1 = Nkx1/(N - m);
  43. Bmx2 = Nkx2/(N - m);
  44. % 样本熵
  45. SampEnVal = -log(Bmx2/Bmx1);

  • 9
    点赞
  • 77
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
样本MATLAB实现可以使用以下代码: ```matlab function H = SampleEntropy(N, m, r, input) for i = 1 : length(input) for j = 1 : N-m+1 temp = input(i, j:j+m-1); A = []; for k = 1 : N-m+1 if k ~= i dist = max(abs(temp-input(k, j:j+m-1))); if dist <= r A = [A; k]; end end end C(i, j) = length(A); end end B = sum(C); H = -sum(B .* log(B/size(C, 2))) / size(C, 2); end ``` 这段代码实现了计算样本的函数。其中,N表示样本数据序列的长度,m表示采样的模板长度,r表示匹配窗口的差异阈值,input为输入的样本数据序列。在函数中,首先遍历每个样本,然后遍历每个子序列,计算其与其他子序列的距离,并根据阈值判断是否为匹配点。最后,统计每个窗口内的匹配点数,计算概率并计算样本。需要注意的是,若概率为0或1,需要取特殊值进行计算。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* [样本及其matlab实现方法](https://blog.csdn.net/weixin_44463965/article/details/130672201)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *3* [matlab求信号的样本](https://blog.csdn.net/weixin_44463965/article/details/130555499)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值