将import pyemd.*; %改成直接用MATLAB编写CEEMDAN

如果你想在MATLAB中直接编写CEEMDAN算法，而不是依赖外部库如PyEMD，你可以参考以下步骤。由于MATLAB有自己的信号处理工具箱，我们可以利用其内置的信号分析功能。下面是不依赖第三方库的CEEMDAN算法的基本实现：

function imfs = ceemdan_matlab(signal, num_realizations, noise_level)
    % 初始化变量
    imfs = [];
    imf_list = {};
    
    % 参数验证
    if nargin < 3 || isnan(noise_level)
        error('Noise level is required.');
    end
    
    % Step 1: 增加随机噪声并执行EMD
    for iteration = 1:num_realizations
        % 生成随机噪声
        noise = randn(size(signal)) * noise_level;
        
        % 加噪声后的信号
        noisy_signal = signal + noise;
        
        % 使用MATLAB的emd函数（需要安装Wavelet Toolbox）
        try
            [~, imf] = wavedec(noisy_signal, 1, 'sym8', 'reflection');
        catch
            warning('EMD not available, using another method for the first IMF.');
            % 如果没有Wavelet Toolbox，可以用其他方法代替，比如基于小波分解的快速emd
            % 这里简化示例，你可以自定义一个替代方法
            % imf = your_alternative_method(noisy_signal);
        end
        % 将第一个IMF添加到列表
        imf_list{end+1} = imf;
    end
    
    % Step 2: 计算第一阶IMF
    imf_1 = mean(cell2mat(imf_list), 2); % 平均所有第一次IMFs
    imfs{1} = imf_1;
    
    % Step 3: 计算剩余残差
    residual = signal - imf_1;
    
    % Step 4: 递归获取后续IMFs
    for k = 2
        next_imfs = cell(length(imf_list), 1);
        
        % 对每个噪声样本的残差执行EMD
        for i = 1:length(imf_list)
            try
                % 使用EMD提取残差的第一IMF
                [~, next_imf] = wavedec(residual + noise(i,:), 1, 'sym8', 'reflection');
                next_imfs{i} = next_imf;
            catch
                % 同样，如果缺少工具箱，处理异常
                % next_imf = alternative_method(residual + noise(i,:));
            end
        end
        
        % 平均下一次IMF
        if ~isempty(next_imfs)
            imf_k = mean(cell2mat(next_imfs), 2);
            imfs{end+1} = imf_k;
            
            % 更新残差
            residual = residual - imf_k;
        else
            break; % 没有找到足够的极点，结束迭代
        end
        
        % 停止条件：残差极点数小于2
        if length(find(diff(find(diff(sign(diff(residual)))) ~= 0))) < 2
            break;
        end
    end
end

% 示例使用
signal = sin(2*pi*0.05*(0:500)); % 示例合成信号
num_realizations = 100; % 实验次数
noise_level = 0.2; % 噪声水平
imfs = ceemdan_matlab(signal, num_realizations, noise_level);

这个版本的CEEMDAN使用了MATLAB内置的`wavedec`函数进行小波分解，并在每次迭代中处理残差。注意，对于没有Wavelet Toolbox的情况，可能需要自定义一种替代方法。另外，上述代码未包含所有的错误处理细节，实际应用中可能需要更完善的错误检查。