声发射熵值参数的Matlab处理1/2

熵值是声发射数据处理的特征参数，对表征声发射活动具有重要意义。
Santis ¹研究了地震熵与b值的关系，给出了它们之间的表达式，并将熵的方法应用于震例分析，表明地震孕育的过程是熵增加的过程。本文采用熵处理过程如下：首先将声发射震级序列按照最大和最小震级划分等级，取1000个声发射数据作为统计窗口进行分组，计算每组数据中相同震级的声发射数目，进而计算每一震级的概率，然后按照式1.1和式1.2计算熵值。同b值的计算过程一样，也采用滑动窗口的方式取样计算熵值，滑动窗口设为100个声发射事件，以窗口末尾时间作为计算熵值的时刻，得到声发射熵值随时间（应变）的变化过程。

1. 声发射数据格式

打开hghk-1.txt文本文件，数据格式见下图。左列为时间（秒），右列为声发射能量。

2. 在包含entropy.m的文件夹下运行main.m，结果如下。

其中第一列为时间（秒），第2列为熵值。

3. 熵值函数部分源代码

function [t entro]=entropy(data_in,level)
  
  [M N]=size(data_in);
  %记录输入矩阵末尾的时间
  t=data_in(M,1);
%每个等级里计数
 ******

1.SANTIS A D．The Gutenberg–Richter Law and Entropy of Earthquakes：Two Case Studies in Central Italy[J]．Bulletin of the Seismological Society of America，2011，101(3)：1386-1395．

---

1. 1 ↩︎