在高动态场景下,接收端双方的相对运动过大,不仅会有多普勒频偏,而且还会有多普勒变化率的频偏,这一节主要从图的角度观察多普勒和多普勒变化率的区别。
我前面只讲了如何加多普勒,其实加多普勒变化率是一样的道理,从单位就能看出来,多普勒的单位是Hz,多普勒变化率的单位是Hz/s,所以多普勒变化率肯定要乘以时间,代码变化如下:
ts1 = 0 / Fs : 1 / Fs : (ceil(Fs / Rc * (length(adata)) * Gp) - 1) / Fs;
ts2 = phase * (Fs / Rc) /Fs;
ts = ts1 + ts2;
Rc_rec = Rc * (1 - (f_doppler + 0.5 * f_accleration * ts1) / fc);
signal = data(floor(mod(Rc_rec .* ts,length(data))) + 1) .* exp(-1j * (2 * pi * f_doppler * ts + 0.5 * 2 * pi * f_accleration .* ts1 .* ts + theta));
黄色区域部分就是相对于只加上多普勒时代码的变化,至于为啥要乘以一个0.5,我之前问过老师,老师说如果加多普勒相当于vt的话,那么加多普勒变化率相当于0.5at^2。
下面两幅图分别是多普勒和多普勒加速度为100Hz和20kHz/s时signal的图
f_doppler = 100;
f_accleration = 0;
f_doppler = 0;
f_accleration = 2e4;
文章评论(0条评论)
登录后参与讨论