Энергетический расчет канала дальней тропосферной радиосвязи (Программа для ЭВМ) - стр. 5

R1=250;

elseif abs(R-200)

R1=200;

elseif abs(R-150)

R1=150;

else

R1=100;

end

if R1==100

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)

R2=400;

elseif abs(R-300)

R2=300;

elseif abs(R-200)

R2=200;

elseif abs(R-150)

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),' Вт (медианная мощность сигнала на входе приемника тлг канала)']);