示波器信号时频分布图脊线提取代码

以下是一个基于值域滤波的示波器信号时频分布图脊线提取的 Matlab 代码示例：

% 读入示波器信号
[x, Fs] = audioread('oscilloscope_signal.wav');

% 计算短时傅里叶变换（STFT）
window_size = round(0.01*Fs);
noverlap = round(0.5*window_size);
nfft = 2^(nextpow2(window_size)+1);
[S, f, t] = spectrogram(x, window_size, noverlap, nfft, Fs);

% 计算能量的对数值
S_log = 10*log10(abs(S));

% 对能量的对数值进行高通滤波
cutoff_freq = 1000;
[b, a] = butter(6, cutoff_freq/(Fs/2), 'high');
S_log_filt = filtfilt(b, a, S_log);

% 计算脊线
thresh = 0.2*max(S_log_filt(:));
[row, col] = find(S_log_filt >= thresh);

% 绘制时频分布图和脊线
figure;
imagesc(t, f, S_log_filt);
set(gca, 'YDir', 'normal');
hold on;
plot(t(col), f(row), 'g', 'LineWidth', 2);
xlabel('Time (s)');
ylabel('Frequency (Hz)');
title('STFT Spectrogram and Ridges');

在这个例子中，我们首先读入一个示波器信号。然后，我们使用短时傅里叶变换（STFT）计算信号的时频分布图，并计算能量的对数值。接着，我们对能量的对数值进行高通滤波，滤除低于某个阈值的能量，保留高于该阈值的能量，从而提取出脊线。最后，我们将时频分布图和脊线绘制出来，以便进行进一步的分析。