matlab写自适应峰值检测算法

自适应峰值检测算法是一种用于信号处理和数据分析的常见算法。以下是在MATLAB中实现自适应峰值检测算法的一般步骤：

1. 对输入信号进行预处理：例如，使用高通滤波器去除低频噪声。

2. 计算信号的局部极大值：使用MATLAB中的 findpeaks 函数查找信号中的峰值，并记录它们的幅度、位置和宽度。

3. 计算信号的峰值阈值：根据信号的峰值统计信息，使用自适应算法计算信号的峰值阈值。

4. 选择真正的峰值：将局部极大值与峰值阈值比较，选择真正的峰值。

以下是一个示例代码，演示了如何在MATLAB中实现自适应峰值检测算法：

% 生成测试信号
fs = 1000; % 采样率
t = linspace(0,1,fs);
x = sin(2*pi*50*t) + randn(size(t));

% 预处理：使用高通滤波器去除低频噪声
fc = 40; % 截止频率
[b,a] = butter(6,fc/(fs/2),'high');
x = filter(b,a,x);

% 计算局部极大值
[pks,locs,w] = findpeaks(x,'MinPeakDistance',50);

% 计算峰值阈值
p = 0.5; % 阈值系数
n = length(x); % 信号长度
m = round(n/10); % 统计窗口长度
h = zeros(size(x)); % 峰值阈值
for i = 1:length(x)
    if i <= m
        h(i) = p * max(x(1:i));
    elseif i >= n-m
        h(i) = p * max(x(i:end));
    else
        h(i) = p * max(x(i-m:i+m));
    end
end

% 选择真正的峰值
idx = pks > h(locs);
pks = pks(idx);
locs = locs(idx);
w = w(idx);

% 绘制结果
plot(t,x)
hold on
plot(t(locs),pks,'ro')
hold off

在这个示例中，我们首先生成了一个包含噪声的正弦信号，并使用高通滤波器去除低频噪声。然后，我们使用MATLAB中的 findpeaks 函数计算信号的局部极大值，记录它们的幅度、位置和宽度。接下来，我们使用自适应算法计算信号的峰值阈值，并将局部极大值与峰值阈值比较，选择真正的峰值。最后，我们绘制了信号和检测到的峰值。