Добавил:
tinkapunka
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:kursach_RRL
.m clc
clear
close all
%% ИСХОДНЫЕ ДАННЫЕ
k = 0:0.1:1; % относительные координаты
R_zemli = 6370*10^3; % м, радиус Земли
R1_H = 17;
R2_H = 20;
R1 = R1_H * 10^3; % м, протяжнность пролета 1
R2 = R2_H * 10^3; % м, протяжнность пролета 2
h1 = [25 20 10 0 10 20 40 40 30 40 35]; % oтметки высот профиля 1
h2 = [35 20 30 20 30 20 40 40 30 40 35]; % oтметки высот профиля 2
mp1 = [0 9 9 9 20 20 20 20 0 0 0]; % местные предметы рельефа1
mp2 = [0 0 0 0 6 6 25 25 25 0 0]; % местные предметы рельефа2
f = 13; % рабочая частота в ГГц
H1_1 = 17.3 * sqrt((R1_H * 0.6 * (1 - 0.6)) / f)
H1_2 = 17.3 * sqrt((R2_H * 0.6 * (1 - 0.6)) / f)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% ЧАСТЬ 1
% ВЫБОР ОРИЕНТИРОВЧОНОЙ ВЫСОТЫ ПОДВЕСА АНТЕНН
% ПРОФИЛЬ 1
% Форма Земли
Y_1 = R1^2.*k.*(1-k)/(2*R_zemli);
f1 = figure;
f1.Position = [500 250 750 400];
plot(k, Y_1, 'k--', 'LineWidth', 1.5)
hold on
% С рельефом
Y_1_h = Y_1 + h1;
plot(k, Y_1_h, 'k', 'LineWidth', 1.2)
hold on
% с местными предметами рельефа (лес...)
Y_1_m = Y_1_h + mp1;
% plot(k, Y_1_m, 'k.')
X = [k(2) k(3) k(4) k(4) k(3) k(2)];
Y = [Y_1_m(2) Y_1_m(3) Y_1_m(4) Y_1_h(4) Y_1_h(3) Y_1_h(2)];
p1 = patch(X,Y,[0.7 0.7 0.7]);
p1.LineStyle = "-";
mp1 = [0 9 9 20 20 20 20 20 0 0 0]; % местные предметы рельефа1
Y_1_m = Y_1_h + mp1;
X = [k(4) k(5) k(6) k(7) k(8) k(8) k(7) k(6) k(5) k(4)];
Y = [Y_1_m(4) Y_1_m(5) Y_1_m(6) Y_1_m(7) Y_1_m(8) Y_1_h(8) Y_1_h(7) Y_1_h(6) Y_1_h(5) Y_1_h(4)];
p2 = patch(X,Y,[0.5 0.5 0.5]);
p2.LineStyle = "-";
H_1_y = Y_1_m(7) + H1_1
H_1_y = [67.5 75 80];
H_1_x = [0 0.6 1];
plot(H_1_x, H_1_y, 'k--', 'LineWidth', 1.2);
h_1 = H_1_y(1) - Y_1_h(1)
h_2 = H_1_y(3) - Y_1_h(11)
% !! отмечено точками - дорисовать отдельно
grid on
grid minor
xlabel('k')
ylabel('h,м')
ylim([0 120]) % ваша высота
legend('нулевой уровень','рельеф')
% ПРОФИЛЬ 2
% Форма Земли
Y_2 = R2^2.*k.*(1-k)/(2*R_zemli);
f2 = figure;
f2.Position = [500 250 750 400];
plot(k, Y_2, 'k--', 'LineWidth', 1.5)
hold on
% С рельефом
Y_2_h = Y_2 + h2;
plot(k, Y_2_h, 'k', 'LineWidth', 1.2)
hold on
% с лесом
Y_2_m = Y_2_h + mp2;
% plot(k,Y_2,'k.')
X = [k(5) k(6) k(6) k(5)];
Y = [Y_2_m(5) Y_2_m(6) Y_2_h(6) Y_2_h(5)];
patch(X,Y,[0.7 0.7 0.7])
mp2 = [0 0 0 0 6 25 25 25 25 0 0]; % местные предметы рельефа2
Y_2_m = Y_2_h + mp2;
X = [k(6) k(7) k(8) k(9) k(9) k(8) k(7) k(6)];
Y = [Y_2_m(6) Y_2_m(7) Y_2_m(8) Y_2_m(9) Y_2_h(9) Y_2_h(8) Y_2_h(7) Y_2_h(6)];
patch(X,Y,[0.5 0.5 0.5])
H_2_y = Y_2_m(7) + H1_2
H_2_y = [83 83 83];
H_2_x = [0 0.6 1];
plot(H_2_x, H_2_y, 'k--', 'LineWidth', 1.2);
h_1 = H_2_y(1) - Y_2_h(1)
h_2 = H_2_y(3) - Y_2_h(11)
% !! отмечено точками - дорисовать отдельно
grid on
grid minor
xlabel('k')
ylabel('h,м')
ylim([0 120]) % ваша высота
legend('нулевой уровень','рельеф')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% ЧАСТЬ 2
% ЛИНИЯ ПРЯМОЙ ВИДИМОСТИ
% k_max_1 = 0.6; %относительная координата наивысшей точки на трассе
% k_max_2 = 0.6;
%
% point_1 = H1_1 + max(Y_1); % точка через которую будет проходить линия прямой видимости
% point_2 = H1_2 + max(Y_2); % нужна только теоритически для определения высот антенн
%
% % !!!! линии прямой видимости для соответствующих профилей
% % по прошлому графику определяете с одной стороны точку где будет антенна
% % зная точку point через которую проходит прямая составляем уравнение
% % прямой
%
% line1 = 3*k+90;
% line2 = 3.4*k+85;
%
% Y_1 = R1^2.*k.*(1-k)/(2*R_zemli);
% Y_2 = R2^2.*k.*(1-k)/(2*R_zemli);
%
% % 1 профиль
% figure
% plot(k, Y_1, 'k--', 'LineWidth', 1.5)
% hold on
% Y_1 = Y_1 + h1;
% plot(k, Y_1, 'k', 'LineWidth', 1.5)
% hold on
% Y_1 = Y_1 + mp1;
% plot(k,Y_1,'k.')
% hold on
% plot(k, line1, 'k', 'LineWidth', 1)
% grid on
% grid minor
% xlabel('k')
% ylabel('h,м')
% ylim([0 120]) % ваша высота
%
% % 2 профиль
% figure
% plot(k, Y_2, 'k--', 'LineWidth', 1.5)
% hold on
% Y_2 = Y_2 + h2;
% plot(k, Y_2, 'k', 'LineWidth', 1.5)
% hold on
% Y_2 = Y_2 + mp2;
% plot(k,Y_2,'k.')
% hold on
% plot(k, line2, 'k', 'LineWidth', 1)
% grid on
% grid minor
% xlabel('k')
% ylabel('h,м')
% ylim([0 120]) % ваша высота
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% ЧАСТЬ 3
% С УЧЕТОМ РЕФРАКЦИИ
% 1 профиль
Y_1 = R1^2.*k.*(1-k)/(2*R_zemli); % без рефр
Y_norm = R1^2.*k.*(1-k)/(2*R_zemli*1.333); % нормальная рефракция
Y_sub = R1^2.*k.*(1-k)/(2*R_zemli*0.7); % субрефракция
% просто профиль
f3 = figure;
f3.Position = [500 250 750 400];
Y_1_plot = plot(k, Y_1, 'k', 'LineWidth', 1.5)
hold on
% С рельефом
Y_1_h = Y_1 + h1;
plot(k, Y_1_h, 'k', 'LineWidth', 1.2)
hold on
mp1 = [0 9 9 9 20 20 20 20 0 0 0];
% с местными предметами рельефа (лес...)
Y_1_m = Y_1_h + mp1;
clear p1_1
p1_1.Vertices = [k(2) Y_1_m(2); k(3) Y_1_m(3); k(4) Y_1_m(4); k(4) Y_1_h(4); k(3) Y_1_h(3); k(2) Y_1_h(2)];
p1_1.Faces = [1 2 3 4 5 6];
p1_1.LineStyle = "-";
p1_1.LineWidth = 1.2;
p1_1.FaceColor = [0.7 0.7 0.7];
p1_1.FaceAlpha = .3;
patch(p1_1)
mp1 = [0 9 9 20 20 20 20 20 0 0 0]; % местные предметы рельефа1
Y_1_m = Y_1_h + mp1;
clear p1_2
p1_2.Vertices = [k(4) Y_1_m(4); k(5) Y_1_m(5); k(6) Y_1_m(6); k(7) Y_1_m(7); k(8) Y_1_m(8); k(8) Y_1_h(8); k(7) Y_1_h(7); k(6) Y_1_h(6); k(5) Y_1_h(5); k(4) Y_1_h(4)];
p1_2.Faces = [1 2 3 4 5 6 7 8 9 10];
p1_2.LineStyle = "-";
p1_2.LineWidth = 1.2;
p1_2.FaceColor = [0.5 0.5 0.5];
p1_2.FaceAlpha = .3;
patch(p1_2)
% профиль с учётом рефракции
Y_1_norm_plot = plot(k, Y_norm, 'k--', 'LineWidth', 1.5);
hold on
% С рельефом
Y_1_h_norm = Y_norm + h1;
plot(k, Y_1_h_norm, 'k--', 'LineWidth', 1.2)
hold on
mp1 = [0 9 9 9 20 20 20 20 0 0 0];
% с местными предметами рельефа (лес...)
Y_1_m_norm = Y_1_h_norm + mp1;
clear p1_3
p1_3.Vertices = [k(2) Y_1_m_norm(2); k(3) Y_1_m_norm(3); k(4) Y_1_m_norm(4); k(4) Y_1_h_norm(4); k(3) Y_1_h_norm(3); k(2) Y_1_h_norm(2)];
p1_3.Faces = [1 2 3 4 5 6];
p1_3.LineStyle = "--";
p1_3.LineWidth = 1.2;
p1_3.FaceColor = [0.7 0.7 0.7];
p1_3.FaceAlpha = .3;
patch(p1_3)
mp1 = [0 9 9 20 20 20 20 20 0 0 0]; % местные предметы рельефа1
Y_1_m_norm = Y_1_h_norm + mp1;
clear p1_4
p1_4.Vertices = [k(4) Y_1_m_norm(4); k(5) Y_1_m_norm(5); k(6) Y_1_m_norm(6); k(7) Y_1_m_norm(7); k(8) Y_1_m_norm(8); k(8) Y_1_h_norm(8); k(7) Y_1_h_norm(7); k(6) Y_1_h_norm(6); k(5) Y_1_h_norm(5); k(4) Y_1_h_norm(4)];
p1_4.Faces = [1 2 3 4 5 6 7 8 9 10];
p1_4.LineStyle = "--";
p1_4.LineWidth = 1.2;
p1_4.FaceColor = [0.5 0.5 0.5];
p1_4.FaceAlpha = .3;
patch(p1_4)
% профиль с учётом субрефракции
Y_1_sub_plot = plot(k, Y_sub, 'k:', 'LineWidth', 1.5)
hold on
% С рельефом
Y_1_h_sub = Y_sub + h1;
plot(k, Y_1_h_sub, 'k:', 'LineWidth', 1.2)
hold on
mp1 = [0 9 9 9 20 20 20 20 0 0 0];
% с местными предметами рельефа (лес...)
Y_1_m_sub = Y_1_h_sub + mp1;
clear p1_5
p1_5.Vertices = [k(2) Y_1_m_sub(2); k(3) Y_1_m_sub(3); k(4) Y_1_m_sub(4); k(4) Y_1_h_sub(4); k(3) Y_1_h_sub(3); k(2) Y_1_h_sub(2)];
p1_5.Faces = [1 2 3 4 5 6];
p1_5.LineStyle = ":";
p1_5.LineWidth = 1.2;
p1_5.FaceColor = [0.7 0.7 0.7];
p1_5.FaceAlpha = .3;
patch(p1_5)
mp1 = [0 9 9 20 20 20 20 20 0 0 0]; % местные предметы рельефа1
Y_1_m_sub = Y_1_h_sub + mp1;
clear p1_6
p1_6.Vertices = [k(4) Y_1_m_sub(4); k(5) Y_1_m_sub(5); k(6) Y_1_m_sub(6); k(7) Y_1_m_sub(7); k(8) Y_1_m_sub(8); k(8) Y_1_h_sub(8); k(7) Y_1_h_sub(7); k(6) Y_1_h_sub(6); k(5) Y_1_h_sub(5); k(4) Y_1_h_sub(4)];
p1_6.Faces = [1 2 3 4 5 6 7 8 9 10];
p1_6.LineStyle = ":";
p1_6.LineWidth = 1.2;
p1_6.FaceColor = [0.5 0.5 0.5];
p1_6.FaceAlpha = .3;
patch(p1_6)
ylim([0 120])
grid on
grid minor
xlabel('k')
ylabel('h,м')
H_1_y = Y_1_m(7) + H1_1
H_1_y_norm = Y_1_m_norm(7) + H1_1
H_1_y_sub = Y_1_m_sub(7) + H1_1
H_1_y = [67.5 75 80];
H_1_y_norm = [66 73.5 78.5];
H_1_y_sub = [70 77.5 82.5];
H_1_x = [0 0.6 1];
plot(H_1_x, H_1_y, 'k-', 'LineWidth', 1.2);
plot(H_1_x, H_1_y_norm, 'k--', 'LineWidth', 1.2);
plot(H_1_x, H_1_y_sub, 'k:', 'LineWidth', 1.2);
h_1 = H_1_y(1) - Y_1_h(1)
h_2 = H_1_y(3) - Y_1_h(11)
legend([Y_1_plot Y_1_norm_plot Y_1_sub_plot], {'без рефракции K_а_т_м = 1', 'нормальная рефракция K_а_т_м = 1.333', 'субрефракция K_а_т_м = 0.7'})
% 2 профиль
% Форма Земли
Y_2 = R2^2.*k.*(1-k)/(2*R_zemli); % без рефр
Y_norm = R2^2.*k.*(1-k)/(2*R_zemli*1.333); % нормальная рефракция субрефракция
Y_sub = R2^2.*k.*(1-k)/(2*R_zemli*0.7); % нормальная рефракция субрефракция
% просто профиль
f4 = figure;
f4.Position = [500 250 750 400];
Y_2_plot = plot(k, Y_2, 'k', 'LineWidth', 1.5)
hold on
% С рельефом
Y_2_h = Y_2 + h2;
plot(k, Y_2_h, 'k', 'LineWidth', 1.2)
hold on
mp2 = [0 0 0 0 6 6 25 25 25 0 0]; % местные предметы рельефа2
% с местными предметами рельефа (лес...)
Y_2_m = Y_2_h + mp2;
clear p1_1
p1_1.Vertices = [k(5) Y_2_m(5); k(6) Y_2_m(6); k(6) Y_2_h(6); k(5) Y_2_h(5)];
p1_1.Faces = [1 2 3 4];
p1_1.LineStyle = "-";
p1_1.LineWidth = 1.2;
p1_1.FaceColor = [0.7 0.7 0.7];
p1_1.FaceAlpha = .3;
patch(p1_1)
mp2 = [0 0 0 0 6 25 25 25 25 0 0]; % местные предметы рельефа2
% с местными предметами рельефа (лес...)
Y_2_m = Y_2_h + mp2;
clear p1_2
p1_2.Vertices = [k(6) Y_2_m(6); k(7) Y_2_m(7); k(8) Y_2_m(8); k(9) Y_2_m(9); k(9) Y_2_h(9); k(8) Y_2_h(8); k(7) Y_2_h(7); k(6) Y_2_h(6)];
p1_2.Faces = [1 2 3 4 5 6 7 8];
p1_2.LineStyle = "-";
p1_2.LineWidth = 1.2;
p1_2.FaceColor = [0.5 0.5 0.5];
p1_2.FaceAlpha = .3;
patch(p1_2)
% профиль с учётом рефракции
Y_2_norm_plot = plot(k, Y_norm, 'k--', 'LineWidth', 1.5);
hold on
% С рельефом
Y_2_h_norm = Y_norm + h2;
plot(k, Y_2_h_norm, 'k--', 'LineWidth', 1.2)
hold on
mp2 = [0 0 0 0 6 6 25 25 25 0 0]; % местные предметы рельефа2
Y_2_m_norm = Y_2_h_norm + mp2;
clear p1_3
p1_3.Vertices = [k(5) Y_2_m_norm(5); k(6) Y_2_m_norm(6); k(6) Y_2_h_norm(6); k(5) Y_2_h_norm(5)];
p1_3.Faces = [1 2 3 4];
p1_3.LineStyle = "--";
p1_3.LineWidth = 1.2;
p1_3.FaceColor = [0.7 0.7 0.7];
p1_3.FaceAlpha = .3;
patch(p1_3)
mp2 = [0 0 0 0 6 25 25 25 25 0 0]; % местные предметы рельефа2
Y_2_m_norm = Y_2_h_norm + mp2;
clear p1_4
p1_4.Vertices = [k(6) Y_2_m_norm(6); k(7) Y_2_m_norm(7); k(8) Y_2_m_norm(8); k(9) Y_2_m_norm(9); k(9) Y_2_h_norm(9); k(8) Y_2_h_norm(8); k(7) Y_2_h_norm(7); k(6) Y_2_h_norm(6)];
p1_4.Faces = [1 2 3 4 5 6 7 8];
p1_4.LineStyle = "--";
p1_4.LineWidth = 1.2;
p1_4.FaceColor = [0.5 0.5 0.5];
p1_4.FaceAlpha = .3;
patch(p1_4)
% профиль с учётом субрефракции
Y_2_sub_plot = plot(k, Y_sub, 'k:', 'LineWidth', 1.5)
hold on
% С рельефом
Y_2_h_sub = Y_sub + h2;
plot(k, Y_2_h_sub, 'k:', 'LineWidth', 1.2)
hold on
mp2 = [0 0 0 0 6 6 25 25 25 0 0]; % местные предметы рельефа2
Y_2_m_sub = Y_2_h_sub + mp2;
clear p1_5
p1_5.Vertices = [k(5) Y_2_m_sub(5); k(6) Y_2_m_sub(6); k(6) Y_2_h_sub(6); k(5) Y_2_h_sub(5)];
p1_5.Faces = [1 2 3 4];
p1_5.LineStyle = ":";
p1_5.LineWidth = 1.2;
p1_5.FaceColor = [0.7 0.7 0.7];
p1_5.FaceAlpha = .3;
patch(p1_5)
mp2 = [0 0 0 0 6 25 25 25 25 0 0]; % местные предметы рельефа2
Y_2_m_sub = Y_2_h_sub + mp2;
clear p1_6
p1_6.Vertices = [k(6) Y_2_m_sub(6); k(7) Y_2_m_sub(7); k(8) Y_2_m_sub(8); k(9) Y_2_m_sub(9); k(9) Y_2_h_sub(9); k(8) Y_2_h_sub(8); k(7) Y_2_h_sub(7); k(6) Y_2_h_sub(6)];
p1_6.Faces = [1 2 3 4 5 6 7 8];
p1_6.LineStyle = ":";
p1_6.LineWidth = 1.2;
p1_6.FaceColor = [0.5 0.5 0.5];
p1_6.FaceAlpha = .3;
patch(p1_6)
ylim([0 120])
grid on
grid minor
xlabel('k')
ylabel('h,м')
H_2_y = Y_2_m(7) + H1_2
H_2_y_norm = Y_2_m_norm(7) + H1_2
H_2_y_sub = Y_2_m_sub(7) + H1_2
H_2_y = [83 83 83];
H_2_y_norm = [81 81 81];
H_2_y_sub = [86 86 86];
H_2_x = [0 0.6 1];
plot(H_2_x, H_2_y, 'k-', 'LineWidth', 1.2);
plot(H_2_x, H_2_y_norm, 'k--', 'LineWidth', 1.2);
plot(H_2_x, H_2_y_sub, 'k:', 'LineWidth', 1.2);
h_1 = H_2_y(1) - Y_2_h(1)
h_2 = H_2_y(3) - Y_2_h(11)
legend([Y_2_plot Y_2_norm_plot Y_2_sub_plot], {'без рефракции K_а_т_м = 1', 'нормальная рефракция K_а_т_м = 1.333', 'субрефракция K_а_т_м = 0.7'})
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% ЧАСТЬ 4
% построение диаграммы уровней для 1 пролета
f1 = 12.75;
f2 = 13.25;
f = (f2 + f1)/2;
D = 1.2;
L01 = 20*log10(4.189*10e4*R1_H*f)
L02 = 20*log10(4.189*10e4*R2_H*f)
Lr1 = (0.004+0.03)*17
Lr2 = (0.004+0.03)*20
Ldop = 1;
Lf1 = 0.5;
Lf2 = 0.5;
Lrf = 0;
Ppd = 28.5;
G1 = 20*log10(D) + 20*log10(f) + 17.5
G2 = G1;
Ppor = -81;
Ppr1 = Ppd + G1 + G2 - L01 - Lf1 - Lf2 - Lr1 - Lrf - Ldop
Ppr2 = Ppd + G1 + G2 - L02 - Lf1 - Lf2 - Lr2 - Lrf - Ldop
M1 = Ppr1 - Ppor
M2 = Ppr2 - Ppor
% Рпд -L0 -Lф1 +G1 -L0-Lг-Lдоп +G2 -Lф2
% y = [23 23 22.5 57.84 -85.942 -50.6 -51];
% y_M = [23 23 22.5 57.84 -127.442 -92.1 -92.5];
%
% x = [0 1 2 2.5 8.5 9 10];
% figure
% plot(x,y,':',LineWidth=1.3,Color="k")
% hold on
% plot(x,y_M,'--',LineWidth=1.2,Color="k")
%
% grid minor
% ylim([-150 150])
% ylabel('дБ')
% хуйня с затуханиями
a_g = 1.2;
a_v = 1.17;
b_g = 0.025;
b_v = 0.03;
J = 100;
yd_g = b_g * J^a_g
yd_v = b_v * J^a_v
R_ef_1 = R1_H / (1 + (R1_H / (35 * exp(-0.015*J))))
R_ef_2 = R2_H / (1 + (R2_H / (35 * exp(-0.015*J))))
A_g_1 = yd_g * R_ef_1
A_v_1 = yd_v * R_ef_1
A_g_2 = yd_g * R_ef_2
A_v_2 = yd_v * R_ef_2
T_v_1 = 10 ^ ((-0.546 + sqrt(0.295 - 0.172 * log10(M1/(0.12*A_v_1))))/0.086)
T_g_1 = 10 ^ ((-0.546 + sqrt(0.295 - 0.172 * log10(M1/(0.12*A_g_1))))/0.086)
T_v_2 = 10 ^ ((-0.546 + sqrt(0.295 - 0.172 * log10(M2/(0.12*A_v_2))))/0.086)
T_g_2 = 10 ^ ((-0.546 + sqrt(0.295 - 0.172 * log10(M2/(0.12*A_g_2))))/0.086)
K_kl = 2e-3;
b = 1.5;
d = 2;
c = 1;
Pint1 = K_kl * 1 * f^b * R1_H * c
Pint2 = K_kl * 1 * f^b * R2_H * c
CPC1 = Pint1 * 10^(-0.1 * M1)
CPC2 = Pint2 * 10^(-0.1 * M2)
f5 = figure;
f5.Position = [500 250 750 400];
x1 = [0 2 4 5 16 17 19 21];
y1 = [Ppd Ppd Ppd-Lf1 Ppd-Lf1+G1 Ppd-Lf1+G1-L02-Ldop-Lr1 Ppd-Lf1+G1-L02-Ldop-Lr1+G2 Ppd-Lf1+G1-L02-Ldop-Lr1+G2-Lf2 Ppr1];
y1_2 = [Ppd Ppd Ppd-Lf1 Ppd-Lf1+G1 Ppd-Lf1+G1-L02-Ldop-Lr1-M1 Ppd-Lf1+G1-L02-Ldop-Lr1+G2-M1 Ppd-Lf1+G1-L02-Ldop-Lr1+G2-Lf2-M1 Ppr1-M1];
line1_x = [0 21]
line1_y = [0 0];
plot(line1_x, line1_y, 'k-', 'LineWidth', 1.2)
hold on
for i = 2 : 7
line1_x = [x1(i) x1(i)]
line1_y = [y1_2(i) 100];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
end
line1_x = [5 11]
line1_y = [y1(4) y1(4)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
line1_x = [9 16]
line1_y = [y1(5) y1(5)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
line1_x = [9 16]
line1_y = [y1_2(5) y1_2(5)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
plot(x1, y1, 'k-', 'LineWidth', 1.2);
plot(x1, y1_2, 'k--', 'LineWidth', 1.2);
xlim([-4 25]);
ylim([-130 150]);
grid on
grid minor
f6 = figure;
f6.Position = [500 250 750 400];
x1 = [0 2 4 5 16 17 19 21];
y1 = [Ppd Ppd Ppd-Lf1 Ppd-Lf1+G1 Ppd-Lf1+G1-L01-Ldop-Lr2 Ppd-Lf1+G1-L01-Ldop-Lr2+G2 Ppd-Lf1+G1-L01-Ldop-Lr2+G2-Lf2 Ppr2];
y1_2 = [Ppd Ppd Ppd-Lf1 Ppd-Lf1+G1 Ppd-Lf1+G1-L01-Ldop-Lr2-M2 Ppd-Lf1+G1-L01-Ldop-Lr2+G2-M2 Ppd-Lf1+G1-L01-Ldop-Lr2+G2-Lf2-M2 Ppr2-M2];
line1_x = [0 21]
line1_y = [0 0];
plot(line1_x, line1_y, 'k-', 'LineWidth', 1.2)
hold on
for i = 2 : 7
line1_x = [x1(i) x1(i)]
line1_y = [y1_2(i) 100];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
end
line1_x = [5 11]
line1_y = [y1(4) y1(4)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
line1_x = [9 16]
line1_y = [y1(5) y1(5)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
line1_x = [9 16]
line1_y = [y1_2(5) y1_2(5)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
plot(x1, y1, 'k-', 'LineWidth', 1.2);
plot(x1, y1_2, 'k--', 'LineWidth', 1.2);
xlim([-4 25]);
ylim([-130 150]);
grid on
grid minor
Lr1 + L01 + Ldop
Lr2 + L02 + Ldop
f7 = figure;
f7.Position = [500 250 750 400];
line1_y = [h1(1) h1(11)];
line1_x = [0 10];
line2_y = [h2(1) h2(11)];
line2_x = [10 20];
H1 = [H_1_y(1) H_1_y(3)];
H2 = [H_2_y(1) H_2_y(3)];
line3_y = [0 H_1_y(1)];
line3_x = [0 0];
line4_y = [0 H_1_y(3)];
line4_x = [10 10];
line5_y = [0 H_2_y(1)];
line5_x = [10 10];
line6_y = [0 H_2_y(3)];
line6_x = [20 20];
line7_y = [Y_1_m(7) Y_1_m(7) + H1_1 + 1];
line7_x = [7 7];
line8_y = [Y_2_m(7) Y_2_m(7) + H1_2];
line8_x = [17 17];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1);
hold on
plot(line2_x, line2_y, 'k--', 'LineWidth', 1);
plot(line1_x, H1, 'k-', 'LineWidth', 0.5);
plot(line2_x, H2, 'k-', 'LineWidth', 0.5);
plot(line3_x, line3_y, 'k-', 'LineWidth', 2);
plot(line4_x, line4_y, 'k-', 'LineWidth', 2);
plot(line5_x, line5_y, 'k-', 'LineWidth', 2);
plot(line6_x, line6_y, 'k-', 'LineWidth', 2);
plot(line7_x, line7_y, 'k-', 'LineWidth', 0.5);
plot(line8_x, line8_y, 'k-', 'LineWidth', 0.5);
ylim([0 120]);
xlim([-5 25]);
grid on
grid minor
clear
close all
%% ИСХОДНЫЕ ДАННЫЕ
k = 0:0.1:1; % относительные координаты
R_zemli = 6370*10^3; % м, радиус Земли
R1_H = 17;
R2_H = 20;
R1 = R1_H * 10^3; % м, протяжнность пролета 1
R2 = R2_H * 10^3; % м, протяжнность пролета 2
h1 = [25 20 10 0 10 20 40 40 30 40 35]; % oтметки высот профиля 1
h2 = [35 20 30 20 30 20 40 40 30 40 35]; % oтметки высот профиля 2
mp1 = [0 9 9 9 20 20 20 20 0 0 0]; % местные предметы рельефа1
mp2 = [0 0 0 0 6 6 25 25 25 0 0]; % местные предметы рельефа2
f = 13; % рабочая частота в ГГц
H1_1 = 17.3 * sqrt((R1_H * 0.6 * (1 - 0.6)) / f)
H1_2 = 17.3 * sqrt((R2_H * 0.6 * (1 - 0.6)) / f)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% ЧАСТЬ 1
% ВЫБОР ОРИЕНТИРОВЧОНОЙ ВЫСОТЫ ПОДВЕСА АНТЕНН
% ПРОФИЛЬ 1
% Форма Земли
Y_1 = R1^2.*k.*(1-k)/(2*R_zemli);
f1 = figure;
f1.Position = [500 250 750 400];
plot(k, Y_1, 'k--', 'LineWidth', 1.5)
hold on
% С рельефом
Y_1_h = Y_1 + h1;
plot(k, Y_1_h, 'k', 'LineWidth', 1.2)
hold on
% с местными предметами рельефа (лес...)
Y_1_m = Y_1_h + mp1;
% plot(k, Y_1_m, 'k.')
X = [k(2) k(3) k(4) k(4) k(3) k(2)];
Y = [Y_1_m(2) Y_1_m(3) Y_1_m(4) Y_1_h(4) Y_1_h(3) Y_1_h(2)];
p1 = patch(X,Y,[0.7 0.7 0.7]);
p1.LineStyle = "-";
mp1 = [0 9 9 20 20 20 20 20 0 0 0]; % местные предметы рельефа1
Y_1_m = Y_1_h + mp1;
X = [k(4) k(5) k(6) k(7) k(8) k(8) k(7) k(6) k(5) k(4)];
Y = [Y_1_m(4) Y_1_m(5) Y_1_m(6) Y_1_m(7) Y_1_m(8) Y_1_h(8) Y_1_h(7) Y_1_h(6) Y_1_h(5) Y_1_h(4)];
p2 = patch(X,Y,[0.5 0.5 0.5]);
p2.LineStyle = "-";
H_1_y = Y_1_m(7) + H1_1
H_1_y = [67.5 75 80];
H_1_x = [0 0.6 1];
plot(H_1_x, H_1_y, 'k--', 'LineWidth', 1.2);
h_1 = H_1_y(1) - Y_1_h(1)
h_2 = H_1_y(3) - Y_1_h(11)
% !! отмечено точками - дорисовать отдельно
grid on
grid minor
xlabel('k')
ylabel('h,м')
ylim([0 120]) % ваша высота
legend('нулевой уровень','рельеф')
% ПРОФИЛЬ 2
% Форма Земли
Y_2 = R2^2.*k.*(1-k)/(2*R_zemli);
f2 = figure;
f2.Position = [500 250 750 400];
plot(k, Y_2, 'k--', 'LineWidth', 1.5)
hold on
% С рельефом
Y_2_h = Y_2 + h2;
plot(k, Y_2_h, 'k', 'LineWidth', 1.2)
hold on
% с лесом
Y_2_m = Y_2_h + mp2;
% plot(k,Y_2,'k.')
X = [k(5) k(6) k(6) k(5)];
Y = [Y_2_m(5) Y_2_m(6) Y_2_h(6) Y_2_h(5)];
patch(X,Y,[0.7 0.7 0.7])
mp2 = [0 0 0 0 6 25 25 25 25 0 0]; % местные предметы рельефа2
Y_2_m = Y_2_h + mp2;
X = [k(6) k(7) k(8) k(9) k(9) k(8) k(7) k(6)];
Y = [Y_2_m(6) Y_2_m(7) Y_2_m(8) Y_2_m(9) Y_2_h(9) Y_2_h(8) Y_2_h(7) Y_2_h(6)];
patch(X,Y,[0.5 0.5 0.5])
H_2_y = Y_2_m(7) + H1_2
H_2_y = [83 83 83];
H_2_x = [0 0.6 1];
plot(H_2_x, H_2_y, 'k--', 'LineWidth', 1.2);
h_1 = H_2_y(1) - Y_2_h(1)
h_2 = H_2_y(3) - Y_2_h(11)
% !! отмечено точками - дорисовать отдельно
grid on
grid minor
xlabel('k')
ylabel('h,м')
ylim([0 120]) % ваша высота
legend('нулевой уровень','рельеф')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% ЧАСТЬ 2
% ЛИНИЯ ПРЯМОЙ ВИДИМОСТИ
% k_max_1 = 0.6; %относительная координата наивысшей точки на трассе
% k_max_2 = 0.6;
%
% point_1 = H1_1 + max(Y_1); % точка через которую будет проходить линия прямой видимости
% point_2 = H1_2 + max(Y_2); % нужна только теоритически для определения высот антенн
%
% % !!!! линии прямой видимости для соответствующих профилей
% % по прошлому графику определяете с одной стороны точку где будет антенна
% % зная точку point через которую проходит прямая составляем уравнение
% % прямой
%
% line1 = 3*k+90;
% line2 = 3.4*k+85;
%
% Y_1 = R1^2.*k.*(1-k)/(2*R_zemli);
% Y_2 = R2^2.*k.*(1-k)/(2*R_zemli);
%
% % 1 профиль
% figure
% plot(k, Y_1, 'k--', 'LineWidth', 1.5)
% hold on
% Y_1 = Y_1 + h1;
% plot(k, Y_1, 'k', 'LineWidth', 1.5)
% hold on
% Y_1 = Y_1 + mp1;
% plot(k,Y_1,'k.')
% hold on
% plot(k, line1, 'k', 'LineWidth', 1)
% grid on
% grid minor
% xlabel('k')
% ylabel('h,м')
% ylim([0 120]) % ваша высота
%
% % 2 профиль
% figure
% plot(k, Y_2, 'k--', 'LineWidth', 1.5)
% hold on
% Y_2 = Y_2 + h2;
% plot(k, Y_2, 'k', 'LineWidth', 1.5)
% hold on
% Y_2 = Y_2 + mp2;
% plot(k,Y_2,'k.')
% hold on
% plot(k, line2, 'k', 'LineWidth', 1)
% grid on
% grid minor
% xlabel('k')
% ylabel('h,м')
% ylim([0 120]) % ваша высота
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% ЧАСТЬ 3
% С УЧЕТОМ РЕФРАКЦИИ
% 1 профиль
Y_1 = R1^2.*k.*(1-k)/(2*R_zemli); % без рефр
Y_norm = R1^2.*k.*(1-k)/(2*R_zemli*1.333); % нормальная рефракция
Y_sub = R1^2.*k.*(1-k)/(2*R_zemli*0.7); % субрефракция
% просто профиль
f3 = figure;
f3.Position = [500 250 750 400];
Y_1_plot = plot(k, Y_1, 'k', 'LineWidth', 1.5)
hold on
% С рельефом
Y_1_h = Y_1 + h1;
plot(k, Y_1_h, 'k', 'LineWidth', 1.2)
hold on
mp1 = [0 9 9 9 20 20 20 20 0 0 0];
% с местными предметами рельефа (лес...)
Y_1_m = Y_1_h + mp1;
clear p1_1
p1_1.Vertices = [k(2) Y_1_m(2); k(3) Y_1_m(3); k(4) Y_1_m(4); k(4) Y_1_h(4); k(3) Y_1_h(3); k(2) Y_1_h(2)];
p1_1.Faces = [1 2 3 4 5 6];
p1_1.LineStyle = "-";
p1_1.LineWidth = 1.2;
p1_1.FaceColor = [0.7 0.7 0.7];
p1_1.FaceAlpha = .3;
patch(p1_1)
mp1 = [0 9 9 20 20 20 20 20 0 0 0]; % местные предметы рельефа1
Y_1_m = Y_1_h + mp1;
clear p1_2
p1_2.Vertices = [k(4) Y_1_m(4); k(5) Y_1_m(5); k(6) Y_1_m(6); k(7) Y_1_m(7); k(8) Y_1_m(8); k(8) Y_1_h(8); k(7) Y_1_h(7); k(6) Y_1_h(6); k(5) Y_1_h(5); k(4) Y_1_h(4)];
p1_2.Faces = [1 2 3 4 5 6 7 8 9 10];
p1_2.LineStyle = "-";
p1_2.LineWidth = 1.2;
p1_2.FaceColor = [0.5 0.5 0.5];
p1_2.FaceAlpha = .3;
patch(p1_2)
% профиль с учётом рефракции
Y_1_norm_plot = plot(k, Y_norm, 'k--', 'LineWidth', 1.5);
hold on
% С рельефом
Y_1_h_norm = Y_norm + h1;
plot(k, Y_1_h_norm, 'k--', 'LineWidth', 1.2)
hold on
mp1 = [0 9 9 9 20 20 20 20 0 0 0];
% с местными предметами рельефа (лес...)
Y_1_m_norm = Y_1_h_norm + mp1;
clear p1_3
p1_3.Vertices = [k(2) Y_1_m_norm(2); k(3) Y_1_m_norm(3); k(4) Y_1_m_norm(4); k(4) Y_1_h_norm(4); k(3) Y_1_h_norm(3); k(2) Y_1_h_norm(2)];
p1_3.Faces = [1 2 3 4 5 6];
p1_3.LineStyle = "--";
p1_3.LineWidth = 1.2;
p1_3.FaceColor = [0.7 0.7 0.7];
p1_3.FaceAlpha = .3;
patch(p1_3)
mp1 = [0 9 9 20 20 20 20 20 0 0 0]; % местные предметы рельефа1
Y_1_m_norm = Y_1_h_norm + mp1;
clear p1_4
p1_4.Vertices = [k(4) Y_1_m_norm(4); k(5) Y_1_m_norm(5); k(6) Y_1_m_norm(6); k(7) Y_1_m_norm(7); k(8) Y_1_m_norm(8); k(8) Y_1_h_norm(8); k(7) Y_1_h_norm(7); k(6) Y_1_h_norm(6); k(5) Y_1_h_norm(5); k(4) Y_1_h_norm(4)];
p1_4.Faces = [1 2 3 4 5 6 7 8 9 10];
p1_4.LineStyle = "--";
p1_4.LineWidth = 1.2;
p1_4.FaceColor = [0.5 0.5 0.5];
p1_4.FaceAlpha = .3;
patch(p1_4)
% профиль с учётом субрефракции
Y_1_sub_plot = plot(k, Y_sub, 'k:', 'LineWidth', 1.5)
hold on
% С рельефом
Y_1_h_sub = Y_sub + h1;
plot(k, Y_1_h_sub, 'k:', 'LineWidth', 1.2)
hold on
mp1 = [0 9 9 9 20 20 20 20 0 0 0];
% с местными предметами рельефа (лес...)
Y_1_m_sub = Y_1_h_sub + mp1;
clear p1_5
p1_5.Vertices = [k(2) Y_1_m_sub(2); k(3) Y_1_m_sub(3); k(4) Y_1_m_sub(4); k(4) Y_1_h_sub(4); k(3) Y_1_h_sub(3); k(2) Y_1_h_sub(2)];
p1_5.Faces = [1 2 3 4 5 6];
p1_5.LineStyle = ":";
p1_5.LineWidth = 1.2;
p1_5.FaceColor = [0.7 0.7 0.7];
p1_5.FaceAlpha = .3;
patch(p1_5)
mp1 = [0 9 9 20 20 20 20 20 0 0 0]; % местные предметы рельефа1
Y_1_m_sub = Y_1_h_sub + mp1;
clear p1_6
p1_6.Vertices = [k(4) Y_1_m_sub(4); k(5) Y_1_m_sub(5); k(6) Y_1_m_sub(6); k(7) Y_1_m_sub(7); k(8) Y_1_m_sub(8); k(8) Y_1_h_sub(8); k(7) Y_1_h_sub(7); k(6) Y_1_h_sub(6); k(5) Y_1_h_sub(5); k(4) Y_1_h_sub(4)];
p1_6.Faces = [1 2 3 4 5 6 7 8 9 10];
p1_6.LineStyle = ":";
p1_6.LineWidth = 1.2;
p1_6.FaceColor = [0.5 0.5 0.5];
p1_6.FaceAlpha = .3;
patch(p1_6)
ylim([0 120])
grid on
grid minor
xlabel('k')
ylabel('h,м')
H_1_y = Y_1_m(7) + H1_1
H_1_y_norm = Y_1_m_norm(7) + H1_1
H_1_y_sub = Y_1_m_sub(7) + H1_1
H_1_y = [67.5 75 80];
H_1_y_norm = [66 73.5 78.5];
H_1_y_sub = [70 77.5 82.5];
H_1_x = [0 0.6 1];
plot(H_1_x, H_1_y, 'k-', 'LineWidth', 1.2);
plot(H_1_x, H_1_y_norm, 'k--', 'LineWidth', 1.2);
plot(H_1_x, H_1_y_sub, 'k:', 'LineWidth', 1.2);
h_1 = H_1_y(1) - Y_1_h(1)
h_2 = H_1_y(3) - Y_1_h(11)
legend([Y_1_plot Y_1_norm_plot Y_1_sub_plot], {'без рефракции K_а_т_м = 1', 'нормальная рефракция K_а_т_м = 1.333', 'субрефракция K_а_т_м = 0.7'})
% 2 профиль
% Форма Земли
Y_2 = R2^2.*k.*(1-k)/(2*R_zemli); % без рефр
Y_norm = R2^2.*k.*(1-k)/(2*R_zemli*1.333); % нормальная рефракция субрефракция
Y_sub = R2^2.*k.*(1-k)/(2*R_zemli*0.7); % нормальная рефракция субрефракция
% просто профиль
f4 = figure;
f4.Position = [500 250 750 400];
Y_2_plot = plot(k, Y_2, 'k', 'LineWidth', 1.5)
hold on
% С рельефом
Y_2_h = Y_2 + h2;
plot(k, Y_2_h, 'k', 'LineWidth', 1.2)
hold on
mp2 = [0 0 0 0 6 6 25 25 25 0 0]; % местные предметы рельефа2
% с местными предметами рельефа (лес...)
Y_2_m = Y_2_h + mp2;
clear p1_1
p1_1.Vertices = [k(5) Y_2_m(5); k(6) Y_2_m(6); k(6) Y_2_h(6); k(5) Y_2_h(5)];
p1_1.Faces = [1 2 3 4];
p1_1.LineStyle = "-";
p1_1.LineWidth = 1.2;
p1_1.FaceColor = [0.7 0.7 0.7];
p1_1.FaceAlpha = .3;
patch(p1_1)
mp2 = [0 0 0 0 6 25 25 25 25 0 0]; % местные предметы рельефа2
% с местными предметами рельефа (лес...)
Y_2_m = Y_2_h + mp2;
clear p1_2
p1_2.Vertices = [k(6) Y_2_m(6); k(7) Y_2_m(7); k(8) Y_2_m(8); k(9) Y_2_m(9); k(9) Y_2_h(9); k(8) Y_2_h(8); k(7) Y_2_h(7); k(6) Y_2_h(6)];
p1_2.Faces = [1 2 3 4 5 6 7 8];
p1_2.LineStyle = "-";
p1_2.LineWidth = 1.2;
p1_2.FaceColor = [0.5 0.5 0.5];
p1_2.FaceAlpha = .3;
patch(p1_2)
% профиль с учётом рефракции
Y_2_norm_plot = plot(k, Y_norm, 'k--', 'LineWidth', 1.5);
hold on
% С рельефом
Y_2_h_norm = Y_norm + h2;
plot(k, Y_2_h_norm, 'k--', 'LineWidth', 1.2)
hold on
mp2 = [0 0 0 0 6 6 25 25 25 0 0]; % местные предметы рельефа2
Y_2_m_norm = Y_2_h_norm + mp2;
clear p1_3
p1_3.Vertices = [k(5) Y_2_m_norm(5); k(6) Y_2_m_norm(6); k(6) Y_2_h_norm(6); k(5) Y_2_h_norm(5)];
p1_3.Faces = [1 2 3 4];
p1_3.LineStyle = "--";
p1_3.LineWidth = 1.2;
p1_3.FaceColor = [0.7 0.7 0.7];
p1_3.FaceAlpha = .3;
patch(p1_3)
mp2 = [0 0 0 0 6 25 25 25 25 0 0]; % местные предметы рельефа2
Y_2_m_norm = Y_2_h_norm + mp2;
clear p1_4
p1_4.Vertices = [k(6) Y_2_m_norm(6); k(7) Y_2_m_norm(7); k(8) Y_2_m_norm(8); k(9) Y_2_m_norm(9); k(9) Y_2_h_norm(9); k(8) Y_2_h_norm(8); k(7) Y_2_h_norm(7); k(6) Y_2_h_norm(6)];
p1_4.Faces = [1 2 3 4 5 6 7 8];
p1_4.LineStyle = "--";
p1_4.LineWidth = 1.2;
p1_4.FaceColor = [0.5 0.5 0.5];
p1_4.FaceAlpha = .3;
patch(p1_4)
% профиль с учётом субрефракции
Y_2_sub_plot = plot(k, Y_sub, 'k:', 'LineWidth', 1.5)
hold on
% С рельефом
Y_2_h_sub = Y_sub + h2;
plot(k, Y_2_h_sub, 'k:', 'LineWidth', 1.2)
hold on
mp2 = [0 0 0 0 6 6 25 25 25 0 0]; % местные предметы рельефа2
Y_2_m_sub = Y_2_h_sub + mp2;
clear p1_5
p1_5.Vertices = [k(5) Y_2_m_sub(5); k(6) Y_2_m_sub(6); k(6) Y_2_h_sub(6); k(5) Y_2_h_sub(5)];
p1_5.Faces = [1 2 3 4];
p1_5.LineStyle = ":";
p1_5.LineWidth = 1.2;
p1_5.FaceColor = [0.7 0.7 0.7];
p1_5.FaceAlpha = .3;
patch(p1_5)
mp2 = [0 0 0 0 6 25 25 25 25 0 0]; % местные предметы рельефа2
Y_2_m_sub = Y_2_h_sub + mp2;
clear p1_6
p1_6.Vertices = [k(6) Y_2_m_sub(6); k(7) Y_2_m_sub(7); k(8) Y_2_m_sub(8); k(9) Y_2_m_sub(9); k(9) Y_2_h_sub(9); k(8) Y_2_h_sub(8); k(7) Y_2_h_sub(7); k(6) Y_2_h_sub(6)];
p1_6.Faces = [1 2 3 4 5 6 7 8];
p1_6.LineStyle = ":";
p1_6.LineWidth = 1.2;
p1_6.FaceColor = [0.5 0.5 0.5];
p1_6.FaceAlpha = .3;
patch(p1_6)
ylim([0 120])
grid on
grid minor
xlabel('k')
ylabel('h,м')
H_2_y = Y_2_m(7) + H1_2
H_2_y_norm = Y_2_m_norm(7) + H1_2
H_2_y_sub = Y_2_m_sub(7) + H1_2
H_2_y = [83 83 83];
H_2_y_norm = [81 81 81];
H_2_y_sub = [86 86 86];
H_2_x = [0 0.6 1];
plot(H_2_x, H_2_y, 'k-', 'LineWidth', 1.2);
plot(H_2_x, H_2_y_norm, 'k--', 'LineWidth', 1.2);
plot(H_2_x, H_2_y_sub, 'k:', 'LineWidth', 1.2);
h_1 = H_2_y(1) - Y_2_h(1)
h_2 = H_2_y(3) - Y_2_h(11)
legend([Y_2_plot Y_2_norm_plot Y_2_sub_plot], {'без рефракции K_а_т_м = 1', 'нормальная рефракция K_а_т_м = 1.333', 'субрефракция K_а_т_м = 0.7'})
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% ЧАСТЬ 4
% построение диаграммы уровней для 1 пролета
f1 = 12.75;
f2 = 13.25;
f = (f2 + f1)/2;
D = 1.2;
L01 = 20*log10(4.189*10e4*R1_H*f)
L02 = 20*log10(4.189*10e4*R2_H*f)
Lr1 = (0.004+0.03)*17
Lr2 = (0.004+0.03)*20
Ldop = 1;
Lf1 = 0.5;
Lf2 = 0.5;
Lrf = 0;
Ppd = 28.5;
G1 = 20*log10(D) + 20*log10(f) + 17.5
G2 = G1;
Ppor = -81;
Ppr1 = Ppd + G1 + G2 - L01 - Lf1 - Lf2 - Lr1 - Lrf - Ldop
Ppr2 = Ppd + G1 + G2 - L02 - Lf1 - Lf2 - Lr2 - Lrf - Ldop
M1 = Ppr1 - Ppor
M2 = Ppr2 - Ppor
% Рпд -L0 -Lф1 +G1 -L0-Lг-Lдоп +G2 -Lф2
% y = [23 23 22.5 57.84 -85.942 -50.6 -51];
% y_M = [23 23 22.5 57.84 -127.442 -92.1 -92.5];
%
% x = [0 1 2 2.5 8.5 9 10];
% figure
% plot(x,y,':',LineWidth=1.3,Color="k")
% hold on
% plot(x,y_M,'--',LineWidth=1.2,Color="k")
%
% grid minor
% ylim([-150 150])
% ylabel('дБ')
% хуйня с затуханиями
a_g = 1.2;
a_v = 1.17;
b_g = 0.025;
b_v = 0.03;
J = 100;
yd_g = b_g * J^a_g
yd_v = b_v * J^a_v
R_ef_1 = R1_H / (1 + (R1_H / (35 * exp(-0.015*J))))
R_ef_2 = R2_H / (1 + (R2_H / (35 * exp(-0.015*J))))
A_g_1 = yd_g * R_ef_1
A_v_1 = yd_v * R_ef_1
A_g_2 = yd_g * R_ef_2
A_v_2 = yd_v * R_ef_2
T_v_1 = 10 ^ ((-0.546 + sqrt(0.295 - 0.172 * log10(M1/(0.12*A_v_1))))/0.086)
T_g_1 = 10 ^ ((-0.546 + sqrt(0.295 - 0.172 * log10(M1/(0.12*A_g_1))))/0.086)
T_v_2 = 10 ^ ((-0.546 + sqrt(0.295 - 0.172 * log10(M2/(0.12*A_v_2))))/0.086)
T_g_2 = 10 ^ ((-0.546 + sqrt(0.295 - 0.172 * log10(M2/(0.12*A_g_2))))/0.086)
K_kl = 2e-3;
b = 1.5;
d = 2;
c = 1;
Pint1 = K_kl * 1 * f^b * R1_H * c
Pint2 = K_kl * 1 * f^b * R2_H * c
CPC1 = Pint1 * 10^(-0.1 * M1)
CPC2 = Pint2 * 10^(-0.1 * M2)
f5 = figure;
f5.Position = [500 250 750 400];
x1 = [0 2 4 5 16 17 19 21];
y1 = [Ppd Ppd Ppd-Lf1 Ppd-Lf1+G1 Ppd-Lf1+G1-L02-Ldop-Lr1 Ppd-Lf1+G1-L02-Ldop-Lr1+G2 Ppd-Lf1+G1-L02-Ldop-Lr1+G2-Lf2 Ppr1];
y1_2 = [Ppd Ppd Ppd-Lf1 Ppd-Lf1+G1 Ppd-Lf1+G1-L02-Ldop-Lr1-M1 Ppd-Lf1+G1-L02-Ldop-Lr1+G2-M1 Ppd-Lf1+G1-L02-Ldop-Lr1+G2-Lf2-M1 Ppr1-M1];
line1_x = [0 21]
line1_y = [0 0];
plot(line1_x, line1_y, 'k-', 'LineWidth', 1.2)
hold on
for i = 2 : 7
line1_x = [x1(i) x1(i)]
line1_y = [y1_2(i) 100];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
end
line1_x = [5 11]
line1_y = [y1(4) y1(4)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
line1_x = [9 16]
line1_y = [y1(5) y1(5)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
line1_x = [9 16]
line1_y = [y1_2(5) y1_2(5)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
plot(x1, y1, 'k-', 'LineWidth', 1.2);
plot(x1, y1_2, 'k--', 'LineWidth', 1.2);
xlim([-4 25]);
ylim([-130 150]);
grid on
grid minor
f6 = figure;
f6.Position = [500 250 750 400];
x1 = [0 2 4 5 16 17 19 21];
y1 = [Ppd Ppd Ppd-Lf1 Ppd-Lf1+G1 Ppd-Lf1+G1-L01-Ldop-Lr2 Ppd-Lf1+G1-L01-Ldop-Lr2+G2 Ppd-Lf1+G1-L01-Ldop-Lr2+G2-Lf2 Ppr2];
y1_2 = [Ppd Ppd Ppd-Lf1 Ppd-Lf1+G1 Ppd-Lf1+G1-L01-Ldop-Lr2-M2 Ppd-Lf1+G1-L01-Ldop-Lr2+G2-M2 Ppd-Lf1+G1-L01-Ldop-Lr2+G2-Lf2-M2 Ppr2-M2];
line1_x = [0 21]
line1_y = [0 0];
plot(line1_x, line1_y, 'k-', 'LineWidth', 1.2)
hold on
for i = 2 : 7
line1_x = [x1(i) x1(i)]
line1_y = [y1_2(i) 100];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
end
line1_x = [5 11]
line1_y = [y1(4) y1(4)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
line1_x = [9 16]
line1_y = [y1(5) y1(5)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
line1_x = [9 16]
line1_y = [y1_2(5) y1_2(5)];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1.2)
plot(x1, y1, 'k-', 'LineWidth', 1.2);
plot(x1, y1_2, 'k--', 'LineWidth', 1.2);
xlim([-4 25]);
ylim([-130 150]);
grid on
grid minor
Lr1 + L01 + Ldop
Lr2 + L02 + Ldop
f7 = figure;
f7.Position = [500 250 750 400];
line1_y = [h1(1) h1(11)];
line1_x = [0 10];
line2_y = [h2(1) h2(11)];
line2_x = [10 20];
H1 = [H_1_y(1) H_1_y(3)];
H2 = [H_2_y(1) H_2_y(3)];
line3_y = [0 H_1_y(1)];
line3_x = [0 0];
line4_y = [0 H_1_y(3)];
line4_x = [10 10];
line5_y = [0 H_2_y(1)];
line5_x = [10 10];
line6_y = [0 H_2_y(3)];
line6_x = [20 20];
line7_y = [Y_1_m(7) Y_1_m(7) + H1_1 + 1];
line7_x = [7 7];
line8_y = [Y_2_m(7) Y_2_m(7) + H1_2];
line8_x = [17 17];
plot(line1_x, line1_y, 'k--', 'LineWidth', 1);
hold on
plot(line2_x, line2_y, 'k--', 'LineWidth', 1);
plot(line1_x, H1, 'k-', 'LineWidth', 0.5);
plot(line2_x, H2, 'k-', 'LineWidth', 0.5);
plot(line3_x, line3_y, 'k-', 'LineWidth', 2);
plot(line4_x, line4_y, 'k-', 'LineWidth', 2);
plot(line5_x, line5_y, 'k-', 'LineWidth', 2);
plot(line6_x, line6_y, 'k-', 'LineWidth', 2);
plot(line7_x, line7_y, 'k-', 'LineWidth', 0.5);
plot(line8_x, line8_y, 'k-', 'LineWidth', 0.5);
ylim([0 120]);
xlim([-5 25]);
grid on
grid minor
Соседние файлы в предмете Спутниковые и радиорелейные линии связи