信号频域相关提取有用信号偏移位置相关问题

---

前言

今天犯了一个低级错误，费了不少时间去填补自己给自己埋的坑，主要就是频域相关提取信号所在索引位置相关问题，既然犯了错就写个博客记录下，免得下次重蹈覆辙，也算给自己一个教训。

---

一、信号时域和频域相关原理

有关信号时域和频域相关原理请参考我之前的博客内容：信号时域和频域相关原理

二、问题详情

我这里有两段信号 x 和 y，其中 y 信号是我们认为的有用信号，且 x 中包括 y 信号，我现在要做的就是提取出 y 信号在 x 信号中的偏移位置，原理自然就是做信号互相关了，在这里我要用频域互相关的方法解决这个问题：

1、MATLAB 代码

% 假设 x 和 y 是你的两个输入信号
y =[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16];    % 举例的信号 y[n] 
x =[ 9 6 14 0 y 7 6 7 12 5 12 52 64 48 11 -38 -61];          % 举例的信号 x[n]

x_length = length(x);
y_length = length(y);

lags = 1 : x_length;

% 进行频域计算
X = fft(x, x_length );
Y = fft(y, x_length );
R_freq = ifft(X .* conj(Y));

figure(1);
plot(lags,R_freq,'r')
xlabel('时延')
ylabel('相关值幅度')
legend('频域计算结果')

% 计算噪声（x中去除y）
noise = x - [zeros(1, 4), y, zeros(1, length(x) - length(y) - 4)];

% 计算信号和噪声的功率
P_signal = sum(y .^ 2) / length(y);
P_noise = sum(noise .^ 2) / length(noise);

% 计算信噪比
SNR = 10 * log10(P_signal / P_noise);

% 显示结果
disp(['信噪比(SNR)为: ', num2str(SNR), ' dB']);

2、仿真结果

频域互相关计算结果如下：

信噪比结果如下：

从上面的结果看出，当前最大值点在第 14 个点，但是我们的信号在第 5 个点，因此频域互相关后并不能提取出 y 信号在 x 信号中的偏移位置，目前信噪比为 -7dB。

三、解决办法

我们将 x 信号中数值较大的改的小一些，也就是从提高信噪比方面来解决问题：

1、MATLAB 代码

% 假设 x 和 y 是你的两个输入信号
y =[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16];    % 举例的信号 y[n] 
x =[ 9 6 14 0 y 7 6 7 12 5 12 10 10 10 10 10 10];          % 举例的信号 x[n]

x_length = length(x);
y_length = length(y);

lags = 1 : x_length;

% 进行频域计算
X = fft(x, 32);
Y = fft(y, 32);
R_freq = ifft(X .* conj(Y));

figure(1);
plot(lags,R_freq,'r')
xlabel('时延')
ylabel('相关值幅度')
legend('频域计算结果')

% 计算噪声（x中去除y）
noise = x - [zeros(1, 4), y, zeros(1, length(x) - length(y) - 4)];

% 计算信号和噪声的功率
P_signal = sum(y .^ 2) / length(y);
P_noise = sum(noise .^ 2) / length(noise);

% 计算信噪比
SNR = 10 * log10(P_signal / P_noise);

% 显示结果
disp(['信噪比(SNR)为: ', num2str(SNR), ' dB']);

2、仿真结果

频域互相关计算结果如下：

信噪比结果如下：

从上面的结果看出，频域互相关后并能够提取出 y 信号在 x 信号中的偏移位置，也就是偏移 5 个位置，目前信噪比为 3.4242 dB。

四、结论

从上面的问题可以得出结论，我们在做信号互相关提取有用信息的时候一定要考虑信噪比这个问题，如果信噪比过小，我们将无法提取出我们的有用信息。

---

我的qq：2442391036，欢迎交流！

---