matlab 样本熵代码

根据引用提供的信息，该作者对一系列熵特征提取进行了封装，其中包括样本熵。以下是样本熵的MATLAB代码实现：

function [SampEn]=sampen(X,M,r)
% X: time series;
% M: embedding dimension;
% r: tolerance;
N=length(X);
% create matrix
XMat=zeros(N-M+1,M);
for i=1:N-M+1
XMat(i,:)=X(i:i+M-1);
end
% calculate Euclidean distance
DMat=zeros(N-M+1,N-M+1);
for i=1:N-M+1
for j=i:N-M+1
DMat(i,j)=max(abs(XMat(i,:)-XMat(j,:)));
DMat(j,i)=DMat(i,j);
end
end
% count similar patterns
B=zeros(1,N-M);
for i=1:N-M
B(i)=length(find(DMat(i+1:N-M+1,i)<=r))-1;
end
% calculate SampEn
SampEn=-log(sum(B)/length(B));
end

其中，X为时间序列，M为嵌入维度，r为容差。函数返回样本熵SampEn的值。