用findpeaks函数编写一个关于动态阈值法的程序

动态阈值法是一种常用的信号处理技术，通过对信号峰值的动态调整来实现峰值检测。以下是一个使用Matlab的`findpeaks`函数实现动态阈值法的简单示例程序：

% 创建一个测试信号
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间向量
x = sin(2*pi*50*t) + 0.5*sin(2*pi*120*t) + 0.2*cos(2*pi*300*t) + randn(size(t)); % 信号

% 初始化动态阈值
threshold = 0.5;
min_peak_distance = 0.1;
peak_prominence = 0.3;

% 遍历信号，查找峰值
for i = 1:length(x)
    % 计算当前阈值
    threshold = 0.5*max(x(max(1, i-100):min(length(x), i+100)));
    
    % 查找峰值
    [peaks,locs] = findpeaks(x(max(1, i-100):min(length(x), i+100)), 'MinPeakDistance', min_peak_distance, 'MinPeakProminence', peak_prominence, 'Threshold', threshold);
    
    % 将峰值位置转换为全局位置
    locs = locs + i - 101;
    
    % 输出结果
    if ~isempty(peaks)
        disp(['Found ' num2str(length(peaks)) ' peaks at positions: ' num2str(locs)]);
    end
end

在这个程序中，我们首先创建了一个测试信号`x`。然后，我们使用`findpeaks`函数遍历信号，查找峰值。在每个时间步骤中，我们计算当前阈值，然后使用`findpeaks`函数查找在该时间步骤内信号的所有峰值。为了确保峰值的准确性，我们设置了最小峰间距（`min_peak_distance`）和峰值突出度（`peak_prominence`）参数。最后，我们将峰值位置转换为全局位置，并输出结果。

这个程序只是一个简单的示例，您可以根据自己的需要调整阈值、峰间距和峰值突出度等参数，以适应不同的信号处理任务。