互相关法求时延

bqgup

互相关法求时延 [复制链接]

# 互相关法求时延

#### 假设已知信号为起始频率为250Hz，终止频率为2000Hz的Chirp信号，接收信号为麦克风利用ADC采集到的在幅度上有一定衰减的Chirp信号。

# 互相关法求时延
```
%%% 互相关法求时延 %%%
N=2048;     % 长度
Fs=10000;   % 采样频率
n=0:N-1;
t=n/Fs;     % 时间序列
%% Chirp信号
x = chirp(t,250,0.2047,2000);
x = [x,zeros(1,2*N)];
      % 要保证时延在两倍信号长度之内，不然就调整2*N的2
%% 接收到的Chirp信号
ADC1 = 0.78 * chirp(t,250,0.2047,2000);
      % 接收后的信息假设存在幅度衰减
tao = 0.15;   % 延时
Ndelay = fix(tao*Fs);
      % 换算成延时点数
ADC1 = [zeros(1,Ndelay),ADC1,zeros(1,length(x)-Ndelay-length(ADC1))];%,zeros(1,2048)];
%% 互相关函数
[c,lags]=xcorr(x,ADC1);
       % 互相关
subplot(211);
xaxis = 1:1:length(x);
plot(xaxis,x,'r');
hold on;
plot(xaxis,ADC1,':');
legend('x信号', 'ADC1信号');
xlabel('时间/s');ylabel('x(t) ADC1(t)');
title('原始信号');grid on;
hold off
subplot(212);
plot(lags/Fs,c,'r');
title('相关函数');
xlabel('时间(s)');ylabel('Rxadc1');
grid on
[Am,Lm]=max(c);
d = Lm - (length(c)+1)/2;
phy=(2*10*d*180/Fs);
phy = rem(phy,360)% 取余360

phy =

-180

Delay=d/Fs

Delay =

-0.1500
```
#### 运行结果正确，图形如下:

<ignore_js_op>

<a href="javascript:;" onclick="showWindow(this.id, this.getAttribute('url'), 'get', 0);" id="savephoto_466660" url="home.php?mod=spacecp&ac=album&op=saveforumphoto&aid=466660&handlekey=savephoto_466660">保存到相册</a>

</p>

</div>
<div class="tip_horn"></div>
</div>

</ignore_js_op>

littleshrimp

谢谢分享，在玩声源定位吗?

bqgup

littleshrimp 发表于 2020-3-26 11:52 谢谢分享，在玩声源定位吗?

是的，感觉这个挺有意思的，还望大神指教

littleshrimp

bqgup 发表于 2020-3-26 18:17 是的，感觉这个挺有意思的，还望大神指教

都研究一这步了，应该向你学习。

Smile丶

大佬，请问最后phy的计算结果代表的是什么意思

bqgup

Smile丶发表于 2020-5-18 14:10 大佬，请问最后phy的计算结果代表的是什么意思

相位差

互相关法求时延 [复制链接]

最新回复

赞赏

点评

点评

点评