无线电系统的主要特性和参数的计算
该代码演示了无线电系统的主要特性和参数的计算，包括无线电通信范围、信号衰减水平、信噪比和检测无线电信号的概率。 考虑了天线参数、大气衰减、频谱噪声密度等各种因素的影响。
计算的关键步骤和原则包括以下步骤。
- 计算天线高度、直接无线电能见度范围和发射机和接收机之间的倾斜范围.
- 无线电波在自由空间和大气中的衰减计算。
- 基于发射机功率、天线增益和噪声电平确定SNR。
In [ ]:
# Определение исходных данных
Rз = 6371.0 # Радиус Земли, км
RЭ = 8500.0 # Эквивалентный радиус Земли, км
f = 49.0 # Частота, МГц
P = 100.0 # Мощность, Вт
Δfc = 0.025 # Ширина полосы частот, МГц
φN = 0.965167 # Широта N, рад
λN = 1.448623 # Долгота N, рад
ΔhN0 = 0.162 # Высота над уровнем моря N, км
hN = 0.016 # Высота антенны передатчика, км
Gφ = -2.0 # Коэффициент усиления антенны передатчика, дБ
ha = 42000.0 # Высота апогея, км
hp = 32000.0 # Высота перигея, км
i = 0.17453 # Наклон орбиты, рад
φM = 0.0 # Широта M, рад
λM = 1.666789 # Долгота M, рад
ΔhM = 0.0 # Высота над уровнем моря M, км
N0 = -201.0 # Спектральная плотность мощности шума приемного устройства, дБ
Gp1 = 22.0 # Коэффициент усиления антенны на частоте f1
Gp2 = 35.0 # Коэффициент усиления антенны на частоте f2
pfa = 1e-5 # Вероятность ложной тревоги
ν = 1 # Число частотных параметров
println("Инициализация выполнена")
地球表面的PRD和PRM之间的距离
In [ ]:
θ0 = 0.981489 # Геоцентрический угол, рад
D = Rз * θ0 # Дальность между ПРД и ПРМ по поверхности Земли, км
println("Дальность между ПРД и ПРМ: $D км")
探测器在海平面以上的高度
In [ ]:
hM = 36423.578 # Высота M над Землей, км
hM0 = ΔhM + hM # Высота над уровнем моря
println("Высота обнаружителя над уровнем моря: $hM0 км")
发射器的天线高度高于海平面
In [ ]:
hN0 = ΔhN0 + hN # Высота антенны передатчика над уровнем моря
println("Высота антенны передатчика над уровнем моря: $hN0 км")
直接无线电能见度总范围
In [ ]:
DM0 = 9187.426 # Дальность прямой радиовидимости для обнаружителя, км
DN0 = 55.009 # Дальность прямой радиовидимости для ПРД, км
D0 = DM0 + DN0 # Общая дальность прямой радиовидимости, км
println("Общая дальность прямой радиовидимости: $D0 км")
无线电波在自由空间中的衰减
In [ ]:
Dн = 39609.292 # Наклонная дальность, км
Z0 = -32.45 - 20 * log10(Dн * f) # Ослабление в свободном пространстве, дБ
println("Ослабление радиоволн в свободном пространстве: $Z0 дБ")
大气和气象现象期间无线电信号的衰减
In [ ]:
Va = -0.02 # Затухание в атмосфере, дБ
Vг = 0.0 # Затухание в гидрометеообразованиях, дБ
V2 = Va + Vг # Общее затухание
println("Общее затухание в атмосфере и гидрометеообразованиях: $V2 дБ")
空间轨道上的信号衰减
In [ ]:
Z = Z0 + V2 # Полное ослабление сигнала
println("Полное ослабление сигнала: $Z дБ")
# 1.24.2 Уровень атмосферных шумов
Nш = -177.2 - 20 * log10(f) + 47.2 * (2.34 + 0.78 * log10(f))^(-2/3) + Gp1 # Атмосферные шумы, дБ
println("Уровень атмосферных шумов: $Nш дБ")
检测器接收器输出端的信噪比
In [ ]:
PВт = 100.0 # Мощность передатчика в ваттах
Gp = 22.0 # Коэффициент усиления антенны по главному лепестку
Z = Z0 + V2 # Полное ослабление сигнала
# Расчет отношения сигнал/шум
SNR = 10 * log10(PВт) + Gp - Z - N0 # Отношение сигнал/шум, дБ
println("Отношение сигнал/шум на выходе приемного устройства: $SNR дБ")
归一化阈值水平
In [ ]:
Qф = sqrt(2 * log(1 / pfa)) # Вероятность ложной тревоги pfa
println("Нормированный пороговый уровень: $Qф")
函数参数X
In [ ]:
X = SNR - Qф # Функциональный параметр
println("Функциональный параметр X: $X")
基于所述计数登记无线电信号的概率。 如果X<0，我们使用概率函数的近似值。
In [ ]:
if X < 0
Y = (0.707 * abs(X))^2
Pн = exp(-Y) # Вероятность регистрации
println("Вероятность регистрации радиосигнала по отсчетам: $Pн")
else
# Если X >= 0, необходимо использовать другую формулу (можно добавить при необходимости)
println("Функциональный параметр X >= 0, требуется другая аппроксимация")
end
# Если выполняем расчет для частоты f = 50 МГц
f2 = 50.0 # Новая частота
Nш2 = -177.2 - 20 * log10(f2) + 47.2 * (2.34 + 0.78 * log10(f2))^(-2/3) + Gp1 # Уровень шумов на 50 МГц
SNR2 = 10 * log10(PВт) + Gp - Z - N0 # Обновленное отношение сигнал/шум для f = 50 МГц
X2 = SNR2 - Qф
if X2 >= 0
println("На частоте 50 МГц вероятность регистрации = 1.0")
else
Y2 = (0.707 * abs(X2))^2
Pн2 = exp(-Y2)
println("На частоте 50 МГц вероятность регистрации радиосигнала: $Pн2")
end
结论
计算结果表明，:
-
系统中的信号稳定。 在衰减Z=-211.4dB时信噪比SNR=111dB提供条件X>0。
-
衰减和噪声最小。 大气衰减（-0.02dB）和噪声水平为信号传输提供了有利条件。
实际意义:
该系统的高水平可靠性证实了其在困难条件下运行的适用性，这对于远程探测系统和卫星通信尤其重要。