如何计算2信号之间的多径信号,延迟,角度?

时间:2016-12-13 16:27:30

标签: gps signal-processing radio broadcast telecommunication

我试图了解如何计算相位角,延迟和多径信号。我希望每个人都可以帮助我提出想法。

我的场景:我有一个带rx1和rx2的无线电硬件(SDR)和电视台(中心频率= 533MHZ,带宽= 6MHz),离我4.8公里(从谷歌地图测量距离) 。从FCC数据库中我可以获得TX的信息(纬度,经度,天线高度和型号......)。

首先,我计算两个天线之间的相位:

%%%%%%

x = .... ; % complex data from rx1

y = .... ; % complex data from rx2

x = x - mean(x);  X = fft(x);
y = y - mean(y);  Y = fft(y);

[mag_x idx_x] = max(abs(X));
[mag_y idx_y] = max(abs(Y));

px = angle(X(idx_x));
py = angle(Y(idx_y));
phase_lag = py - px;

%%%%%

是对的吗?然后尝试计算多径信号:

speed_of_light = 3e8 ; % m/s
distance = 3 * 1.6;
freq = 533e6;
lamda = speed_of_light / freq;
k = 2 * pi / lamda;
x = distance / 2;
y = 0;
w = 2 * pi * freq;
phase = 0;
delta = 1e-6; % delay
i = sqrt(-1);
t = 0:0.01:1;
degree = 0;
mutilpath = zeros(1,length(t));

index = 10;

for j = 1: index

   mutilpath  = mutilpath + exp(i * (k * (x * cos(phase+degree) + y * sin(phase+degree)) - w * (t - delta)));
    degree    = degree + pi/6;

end

figure(1), plot(t,abs(mutilpath));
figure(2), plot(t,angle(mutilpath));

%figure(3), plot(t,abs(mutilpath));
%figure(4), plot(t,angle(mutilpath));

%%%%%%%%%

当我设置index = 1时;我有一个情节abs(mutilpath)是一条直线(图1)

当我设置index =! 1;我有一个情节abs(mutilpath)是一个不同的线(图3)

我有一种不好的感觉,我做得不对。

如果每个人都有想法,请推荐。

非常感谢你。

0 个答案:

没有答案