多尺度模糊熵matlab程序

以下是多尺度模糊熵的 MATLAB 程序：

function [MFE, MSE] = multiscale_fuzzy_entropy(signal, m, r, scales)
% signal: 输入信号
% m: 模糊熵的阶数
% r: 模糊因子
% scales: 尺度数

N = length(signal);
MFE = zeros(scales, 1);
MSE = zeros(scales, 1);

for s = 1:scales
    % 计算每个尺度下的子信号长度
    L = floor(N/s);
    % 初始化模糊熵和模糊矩阵
    FE = zeros(L,1);
    P = zeros(m+1, L);
    for i = 1:L
        % 提取子信号
        sub_signal = signal((i-1)*s+1:i*s);
        % 计算模糊矩阵
        for j = 1:s
            for k = 1:s
                if abs(sub_signal(j)-sub_signal(k)) <= r
                    P(1, i) = P(1, i) + 1;
                end
                for q = 2:m+1
                    if abs(sub_signal(j)-sub_signal(k)) <= r && abs(sub_signal(j+q-1)-sub_signal(k+q-1)) <= r
                        P(q, i) = P(q, i) + 1;
                    end
                end
            end
        end
        % 计算模糊熵
        for q = 1:m+1
            P(q, i) = P(q, i)/(s*(s-1));
            FE(i) = FE(i) - P(q, i)*log(P(q, i));
        end
    end
    % 计算多尺度模糊熵
    MFE(s) = sum(FE)/L;
    MSE(s) = log(MFE(s));
end

end