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

Чтение книги онлайн.

СКАЧАТЬ abs(R-200)<abs(R-150)

      R1=200;

      elseif abs(R-150)<abs(R-100)

      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)<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=', СКАЧАТЬ