U=Ea*Rf/(Ra+Rf); % мкВ – амплитуда входного сигнала
P=U^2/Rf; % мкВт – мощность входного сигнала
% тепловой шум
Pha=4e6*k*Ta*df; % мкВт – мощность шума антенны
Php=4e6*k*Tp*df; % мкВт – мощность шума входных цепей приемника
Ph=Pha+Php; % мкВт – мощность теплового шума на входе приемника
% эфирный шум (атмосферные помехи)
Faa=Fama+dFama; % дБ – коэффициент атмосферного шума
Epd=Faa+20*log10(f)+10*log10(df)-99; % дБ/мкВ/м – напряженность поля помехи в точке приема для заданного df
Epm=10^(Epd/20); % мкВ/м – медианное изначение напряженности поля помехи в точке приема
ka=10^(A0/20); % коэффициент
Ep=Epm*(1+ka); % действующее изначение напряженности поля помехи в точке приема
Epa=Ep*hde; % мкВ – эдс помехи, наводимая в эквивалентной антенне
Upa=Kp^(1/2)*Epa*Rf/(Ra+Rf); % мкВ – амплитуда атмосферной помехи на выходе реальной антенны
Upam=Upa/(1+ka); % медианное изначение напряжения атмосферной помехи на выходе реальной антенны
siga=Upam*ka; % стандартное отклонение напряжения атмосферной помехи на выходе реальной антенны от медианного значения
sigax=2*siga; % стандартное отклонение моделируемого процесса с нормальным распределением
moax=4*(Upam-siga/(2*pi)^(1/2)); % мат. ожидание моделируемого процесса
Pap=Upa^2/Rf; % мкВт – мощность атмосферной помехи
Pp=(Ph+Pap); % мкВт – мощность атмосферной и тепловой помех на входе приемника
h02a=P/Pp; % раз отношение с/ш на входе приемника
Pad=10*log10(Pp/P); % дБ требуемый уровень помехи относительно сигнала для имитации канала связи
L=D/cos(a1); % км – длина пути луча
t=L/300; % сек. – время прохождения пути
disp([' Входной сигнал', Ap]);
disp(['U=', num2str(U),' мкВ (амплитуда входного сигнала)']);
disp(['P=', num2str(P),' мкВт (мощность входного сигнала)']);
disp(['Ph=', num2str(Ph),' мкВт (мощность тепловой помехи)']);
disp(['L=', num2str(L),' км (длина пути, пройденного лучем)']);
disp(['t=', num2str(t),' мсек. (время прохождения пути)']);
disp(' Атмосферная помеха. Логнормальное распределение.');
disp(['Pap=', num2str(Pap),' мкВт (мощность атмосферной помехи)']);
disp(['h02a=', num2str(h02a),' раз (отношение с/ш)']);
if Udisp([' Амплитуда входного сигнала меньше чувствительности приемника']);
else
disp(['sigax=', num2str(sigax),' (требуемое стандартное отклонение моделируемого процесса)']);
disp(['moax=', num2str(moax),' (требуемое мат. ожидание моделируемого процесса)']);
disp(['Pad=', num2str(Pad),' дБ (требуемая мощность атмосферной помехи относительно мощности входного сигнала для имитации канала связи)']);
