dFamp=7.27; % дБ – стандартное отклонение коэффициента промышленного шума для сельской местности
Fap=Famp+7.27; % дБ – коэффициент промышленного шума для сельской местности
Mp=' (сельская местность) ';
end
Eppd=Fap+20*log10(f)+10*log10(df)-99; % дБ/мкВ/м – напряженность поля помехи в точке приема для заданного df
Eppm=Famp+20*log10(f)+10*log10(df)-99; % дБ/мкВ/м – медианное значение напряженности поля помехи в точке приема для заданного df
dEpp=dFamp+20*log10(f)+10*log10(df)-99; % дБ/мкВ/м – стандортное отклонение от медианного значения напряженности поля помехи в точке приема для заданного df
kp=10^(dEpp/20)/10^(Eppm/20); % коэффициент
Epp=10^(Eppd/20); % мкВ/м – напряженность поля помехи в точке приема
Eppa=Epp*hde; % мкВ – эдс помехи, наводимая в эквивалентной антенне
Upp=Kp^(1/2)*Eppa*Rf/(Ra+Rf); % мкВ – амплитуда промышленной помехи на выходе реальной антенны
Ppp=Upp^2/Rf; % мкВт – мощность промышленной помехи
Uppm=Upp/(1+kp); % медианное изначение напряжения промышленной помехи на выходе реальной антенны
sigp=Uppm*kp; % стандартное отклонение напряжения промышленной помехи на выходе реальной антенны от медианного значения
sigpx=2*sigp; % стандартное отклонение моделируемого процесса с нормальным распределением
mopx=4*(Uppm-sigp/(2*pi)^(1/2)); % мат. ожидание моделируемого процесса
h02p=P/Ppp; % раз отношение с/ш на входе приемника
Ppd=10*log10(1/h02p); % дБ требуемый уровень помехи относительно сигнала для имитации канала связи
disp([' Промышленная помеха. Нормальное распределение.', Mp]);
disp(['Ppp=', num2str(Ppp),' мкВт (мощность промышленной помехи)']);
disp(['h02p=', num2str(h02p),' раз (отношение с/ш)']);
disp(['sigpx=', num2str(sigpx),' (требуемое стандартное отклонение моделируемого процесса)']);
disp(['mopx=', num2str(mopx),' (требуемое мат. ожидание моделируемого процесса)']);
disp(['Ppd=', num2str(Ppd),' дБ (требуемая мощность промышленной помехи относительно мощности входного сигнала для имитации канала связи)']);
% суммарная помеха
h02s=P/(Pp+Ppg+Ppp); % раз отношение с/ш на входе приемника по упрощенной формуле
c=4.343;
a1=exp(Fama/c+dFama^2/(2*c^2)); % частные промежуточные коэффициенты
a2=exp(Famg/c+dFamg^2/(2*c^2));
a3=exp(Famp/c+dFamp^2/(2*c^2));
at=a1+a2+a3; % промежуточный коэффициент
bt=a1^2*(exp(dFama^2/(2*c^2))-1)+a2^2*(exp(dFamg^2/(2*c^2))-1)+a3^2*(exp(dFamp^2/(2*c^2))-1); % промежуточный коэффициент
sigmt=c*(log(1+bt/at^2))^1/2; % стандартное отклонение коэффициента шума от медианного значения
