a1=10.09;
b1=0.2266;
c1=-10.33;
d1=-1.994;
elseif R1==150
a1=9.399;
b1=0.1913;
c1=-9.604;
d1=-1.576;
elseif R1==200
a1=7.04;
b1=0.2457;
c1=-7.011;
d1=-1.612;
elseif R1==250
a1=7.807;
b1=0.1683;
c1=-7.757;
d1=-1.013;
elseif R1==300
a1=7.073;
b1=0.2235;
c1=-7.159;
d1=-1.404;
end
Lp=a1*exp(b1*delta)+c1*exp(d1*delta);
X=h*100/l;
if X<12
if abs(R-400)<abs(R-300)
R2=400;
elseif abs(R-300)<abs(R-200)
R2=300;
elseif abs(R-200)<abs(R-150)
R2=200;
elseif abs(R-150)<abs(R-100)
R2=150;
else
R2=100;
end
if R2==100
a1=1.983e15;
b1=-54.05;
c1=9.514;
a2=18.74;
b2=-6.452;
c2=11.95;
elseif R2==150
a1=1.661e15;
b1=-45.48;
c1=8.023;
a2=14.88;
b2=-3.284;
c2=8.104;
elseif R2==200
a1=1.439e9;
b1=-52.15;
c1=12.28;
a2=3.78;
b2=3.267;
c2=3.428;
elseif R2==300
a1=7835;
b1=-4.066;
c1=3.974;
a2=-7425;
b2=-3.92;
c2=3.917;
elseif R2==400
a1=8.829e13;
b1=-47.8;
c1=9.207;
a2=-4457;
b2=-6.454;
c2=3.497;
end
Lh=a1*exp(-((X-b1)/c1)^2)+a2*exp(-((X-b2)/c2)^2); else
Lh=0;
end
La=0.07*exp(0.055*(Gr+Gp));
Ldop=Lmed+Lz+Lp+Lh+La+Lk;
L1m=L0r+Ldop;
Lm1m=L1m-Lz;
Lm1c=Lm1m+4.676*exp(-((R-124.9)/332.4)^2)+3.165*exp(-((R-3828)/508.3)^2);
Lm1g=Lm1c-14.96*exp(-((R-74.66)/492.7)^2
if Tp==4
Lm1z=Lm1g;
elseif Tp==3 ||Tp==1
Lm1z=Lm1m;
else
Lm1z=Lm1c;
end
L1=Lm1z+Lz;
disp(['L1=', num2str(L1),' дБ (полное затухание)']);
Pvh1=10*log10(P)-L1;
Pvh1vt=10^(Pvh1/10); \
disp(['Pvh1=', num2str(Pvh1),' дБ (мощность сигнала на входе приемника)']);
disp(['Pvh1=', num2str(Pvh1vt),' вт (мощность сигнала на входе приемника)']);
if kan==1|| tipcan==1 \
disp(['L1med=', num2str(Lm1z),' дБ (медианные потери для заданного периода оценки)']);
Pvh1med=10^((10*log10(P)-Lm1z)/10);
disp(['Pvh1med=', num2str(Pvh1med),' вт (медианная мощность сигнала на входе приемника)']);
Pch1=4e-21*10^(Fe/10)*Nk*3100;
Pcht1=6975e-24*10^(Fe/10)*(Fk/dfk)^2/Pvh1vt;
h0t1=10*log10(1e-3/Pcht1);
h01=1e-3/Pcht1;
disp(['h01=', num2str(h01),'(отношение с/ш на входе детектора)']);
disp(['h0t1=', num2str(h0t1),' дБ (отношение с/ш на выходе тлф канала)']);
disp(['Pcht1=', num2str(Pcht1),' Вт (мощность шума в тлф канале, в точке, где Рс=1мВт)']);
disp(['Pch1=', num2str(Pch1),' Вт (мощность шума на входе приемника при частотном уплотнении)']);
elseif kan==2 && tipcan==2
Pch1=4e-21*10^(Fe/10)*V*1000;
h01=(10^(Pvh1/10))/Pch1;
disp(['h01=', num2str(h01),'(отношение с/ш на входе приемника и на входе детектора)']);
L1m=Lm1z+Lmz;
disp(['L1med=', num2str(L1m),' дБ (затухание тлг канала без быстрых замираний)']);
Pvh1med=10^((10*log10(P)-L1m)/10);
disp(['Pvh1med=', num2str(Pvh1med),' Вт (медианная мощность сигнала на входе приемника тлг канала)']);
disp(['Pch1=', num2str(Pch1),' Вт (мощность шума на входе приемника в Вт при непосредственной модуляции)']);
end
Lpred=10*log10(P)-10*log10(Pch1)-B;
disp(['Lpred=', num2str(Lpred),' дБ (допустимое значение полных потерь)']);
Zap1=Lpred-L1;
disp(['Zap1=', num2str(Zap1),' дБ (энергетический запас трассы)']);
if Zap1<1 && Zap1>0.5
Rpred=R;
else
Lp=Ldop-Lmed;
if abs(Zap1)<3
Kd=L0r/((L1-Lp)*2);
else
Kd=L0r/(L1-Lp);
end
Rpred=l*1e-5*10^((L0+Kd*(Zap1-0.7))/20)/(4*pi);
end
disp(['Rpred=',