Моделирование канала коротковолновой радиосвязи | страница 16
end
% галактический шум
if f>10 || f==10
Famg=52.25-10.296*log(f); % дБ – медианное значение коэффициента галактического шума
dFamg=1.56;
Fag=Famg+1.56; % дБ – коэффициент галактического шума
Epgd=Fag+20*log10(f)+10*log10(df)-99; % дБ/мкВ/м – напряженность поля помехи в точке приема для заданного df
Epgm=Famg+20*log10(f)+10*log10(df)-99; % дБ/мкВ/м – медианное значение напряженности поля помехи в точке приема для заданного df
dEpg=dFamg+20*log10(f)+10*log10(df)-99; % дБ/мкВ/м – стандортное отклонение от медианного значения напряженности поля помехи в точке приема для заданного df
kg=10^(dEpg/20)/10^(Epgm/20); % коэффициент
Epg=10^(Epgd/20); % мкВ/м – напряженность поля помехи в точке приема
Epga=Epg*hde; % мкВ – эдс помехи, наводимая в эквивалентной антенне
Upg=Kp^(1/2)*Epga*Rf/(Ra+Rf); % мкВ – амплитуда галактической помехи на выходе реальной антенны
Ppg=Upg^2/Rf; % мкВт – мощность галактической помехи
Upgm=Upg/(1+kg); % медианное изначение напряжения галактической помехи на выходе реальной антенны
sigg=Upgm*kg; % стандартное отклонение напряжения галактической помехи на выходе реальной антенны от медианного значения
siggx=2*sigg; % стандартное отклонение моделируемого процесса с нормальным распределением
mogx=4*(Upgm-sigg/(2*pi)^(1/2)); % мат. ожидание моделируемого процесса
h02g=P/Ppg; % раз отношение с/ш на входе приемника
Pgd=10*log(1/h02g); % дБ требуемый уровень помехи относительно сигнала для имитации канала связи
disp(' Галактический шум. Нормальное распределение.');
disp(['Ppg=', num2str(Ppg),' мкВт (мощность галактической помехи)']);
disp(['h02g=', num2str(h02g),' раз (отношение с/ш)']);
disp(['siggx=', num2str(siggx),' (требуемое стандартное отклонение моделируемого процесса)']);
disp(['mogx=', num2str(mogx),' (требуемое мат. ожидание моделируемого процесса)']);
disp(['Pgd=', num2str(Pgd),' дБ (требуемая мощность галактической помехи относительно мощности входного сигнала для имитации канала связи)']);
else
Ppg=0;
Famg=0;
dFamg=0;
disp(['Ppg=', num2str(Ppg),' ( галактическая помеха отсутствует)']);
end
% промышленный шум
if M==1
Famp=71.52-11.987*log(f); % дБ – медианное значение коэффициента промышленного шума для жилого района
dFamp=8.18; % дБ – стандартное отклонение коэффициента промышленного шума для жилого района
Fap=Famp+8.18; % дБ – коэффициент промышленного шума для жилого района
Mp=' (жилой район) ';
elseif M==2
Famp=66.19-11.987*log(f); % дБ – медианное значение коэффициента промышленного шума для сельской местности