如何将一个信号的频率“强制”到另一个信号上；两个不同频率的信号

Question

我需要测量两个不同频率信号之间的某些参数。即，滞后/相位差。

我知道我不能使用 xcorr，例如，因为频率不同。

数据示例：

如果需要，我可以附上数据。

我可以使用以下方法测量两个信号的频率：

%calculate frequencies 
[maxValue1,indexMax1] = max(abs(fft(sig1-mean(sig1))));
f1 = indexMax1 * 100 / 100;

[maxValue2,indexMax2] = max(abs(fft(sig2-mean(sig2))));
f2 = indexMax2 * 100 / 100;


%which one is higher?
maxF = max (f1, f2);

如何强制/更改任一信号的频率以使其与另一个信号相同？例如，两个信号的频率都应该是 maxF。

Answer 1

FFT 可能不是完成您想做的事情的最佳方式。如果您的两个信号中的每一个都是真正的单分量（即，它们不是多个正弦波的混合），那么您可以使用希尔伯特变换获得更好的性能。希尔伯特将您的实信号转换为复信号。它还允许您非常轻松地直接评估相位、频率和幅度。

%convert to complex domain
h_sig1 = hilbert(sig1);
h_sig2 = hilbert(sig2);

%plot instantaneous phase of the signals
phase_sig1_rad = angle(h_sig1);  %radians
phase_sig2_rad = angle(h_sig2);  %radians
dphase_rad = wrapTo2Pi(phase_sig1_rad - phase_sig2_rad);
plot([phase_sig1_rad(:) phase_sig2_rad(:) dphase_rad(:)]);

%other info: compute the frequency of the signals
dt_sec = 1/sample_rate_Hz;  %what is the time between samples
freq_sig1_Hz = diff(unwrap(phase_sig1_rad))/(2*pi)/dt_sec;
freq_sig2_Hz = diff(unwrap(phase_sig2_rad))/(2*pi)/dt_sec;

%other info: compute the amplitude of the signals
amp_sig1 = abs(h_sig1);
amp_sig2 = abs(h_sig2);

希望这会有所帮助！