我试图了解如何计算相位角,延迟和多径信号。我希望每个人都可以帮助我提出想法。
我的场景:我有一个带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)
我有一种不好的感觉,我做得不对。
如果每个人都有想法,请推荐。
非常感谢你。