极大值法去噪的matlab代码

以下是MATLAB中使用极大值法进行去噪的示例代码：

% 读取原始信号
original_signal = load('original_signal.mat');

% 添加高斯白噪声
noisy_signal = original_signal + 0.1*randn(size(original_signal));

% 设置窗口大小和极大值阈值
win_size = 20;
max_threshold = 0.5;

% 极大值去噪
denoised_signal = zeros(size(noisy_signal));
for i = 1:length(noisy_signal)
    % 窗口边界处理
    start_idx = max(i-win_size,1);
    end_idx = min(i+win_size,length(noisy_signal));
    
    % 获取窗口内的最大值
    local_max = max(noisy_signal(start_idx:end_idx));
    
    % 如果当前点大于等于最大值减去阈值，则保留，否则置为0
    if noisy_signal(i) >= local_max - max_threshold
        denoised_signal(i) = noisy_signal(i);
    else
        denoised_signal(i) = 0;
    end
end

% 绘制去噪结果
figure;
subplot(2,1,1); plot(original_signal); title('Original Signal');
subplot(2,1,2); plot(denoised_signal); title('Denoised Signal');

这个示例代码使用一个固定大小的窗口和一个固定的极大值阈值进行去噪。在每个窗口内，寻找最大值并将窗口内小于最大值减去阈值的点置为0，其余点保留。