Famt=c*(log(at)-sigmt^2/(2*c^2)); % медианное значение коэффициента шума
Fat=Famt+sigmt;% действующее значение коэффициента шума
Eptd=Fat+20*log10(f)+10*log10(df)-99; % дБ/мкВ/м – действующее значение напряженности поля помехи в точке приема для заданного df
Eptm=Famt+20*log10(f)+10*log10(df)-99; % дБ/мкВ/м – медианное значение напряженности поля помехи в точке приема для заданного df
dEpt=sigmt+20*log10(f)+10*log10(df)-99; % дБ/мкВ/м – стандортное отклонение от медианного значения напряженности поля помехи в точке приема для заданного df
kt=10^(dEpt/20)/10^(Eptm/20); % коэффициент
Ept=10^(Eptd/20); % мкВ/м – действующее значение напряженности поля помехи в точке приема для заданного df
Epta=Ept*hde; % мкВ – эдс суммарной помехи, наводимая в эквивалентной антенне
Upt=Kp^(1/2)*Epta*Rf/(Ra+Rf); % мкВ – амплитуда суммарной помехи на выходе реальной антенны
Ppt=Upt^2/Rf; % мкВт – мощность суммарной помехи
Uptm=Upt/(1+kt); % медианное изначение напряжения суммарной помехи на выходе реальной антенны
sigt=Uptm*kt; % стандартное отклонение напряжения суммарной помехи на выходе реальной антенны от медианного значения
sigtx=sigt; % стандартное отклонение моделируемого процесса с нормальным распределением
motx=2*(Uptm-sigt/(2*pi)^(1/2)); % мат. ожидание моделируемого процесса
h02t=P/Ppt;
Ptd=10*log10(1/h02t); % дБ требуемый уровень помехи относительно сигнала для имитации канала связи
disp(' Суммарная помеха. Полунормальное распределение.');
disp(['h02s=', num2str(h02s),' раз (отношение с/ш по упрощенной формуле)']);
disp(['Ppt=', num2str(Ppt),' мкВт (мощность суммарной помехи)']);
disp(['h02t=', num2str(h02t),' раз (отношение с/ш по формуле для суммарной помехи)']);
disp(['sigtx=', num2str(sigtx),' (требуемое стандартное отклонение моделируемого процесса)']);
disp(['motx=', num2str(motx),' (требуемое мат. ожидание моделируемого процесса)']);
disp(['Ptd=', num2str(Ptd),' дБ (требуемая мощность суммарной помехи относительно мощности входного сигнала для имитации канала связи)']);
Приложение Б
Инструкция по работе с программой
KVkanal
Программа KVkanal работает в среде MATLAB 2011. Для подготовки исходных данных нужно иметь распечатку оперативного прогноза по данным с сети ионосферных станций ИЗМИРАН и рекомендации МСЭ-R P.372-9 (Радиошум).
Подготовка исходных данных заключается в определении по распечатке оперативного прогноза даты и времени проведения сеанса связи, координат размещения радиоприемного устройства, дальности связи, рабочей частоты, напряженности поля сигнала в точке приема, угла прихода сигнала, и определении по рекомендациям МСЭ-R P.372-9 медианного значения коэффициента атмосферного шума, стандартного отклонения коэффициента атмосферного шума, и превышение напряжения помехи над медианным значением. Пример распечатки оперативного прогноза приведен на рис. 1.
