Моделирование канала коротковолновой радиосвязи | страница 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); % дБ – медианное значение коэффициента промышленного шума для сельской местности