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

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

СКАЧАТЬ style="font-size:15px;">      c2=14.05;

      a3=7302;

      b3=231.9;

      c3=50.68;

      elseif sigma==4

      a1=234;

      b1=100.8;

      c1=0.4421;

      a2=6.382e16;

      b2=179.1;

      c2=12.98;

      a3=1615;

      b3=219;

      c3=52.85;

      elseif sigma==5

      a1=4.173e15;

      b1=130.9;

      c1=5.308;

      a2=69.67;

      b2=140.1;

      c2=30.56;

      a3=0;

      b3=1;

      c3=1;

      elseif sigma==6

      a1=5.952e14;

      b1=113.5;

      c1=2.384;

      a2=4.38;

      b2=100.4;

      c2=2.216;

      a3=15.26;

      b3=112.7;

      c3=24.18;

      elseif sigma==7

      a1=2.977;

      b1=99.89;

      c1=0.0954;

      a2=1.88e15;

      b2=182.8;

      c2=14.36;

      a3=12400;

      b3=568;

      c3=126.3;

      end

      elseif n==4

      if sigma==3

      a1=1.386;

      b1=99.91;

      c1=0.1078;

      a2=-24.82;

      b2=103.8;

      c2=5.008;

      a3=67530;

      b3=143.8;

      c3=15.46;

      elseif sigma==4

      a1=2271;

      b1=102.4;

      c1=0.9524;

      a2=179.7;

      b2=114.6;

      c2=13.67;

      a3=-86.2;

      b3=108.2;

      c3=10.75;

      elseif sigma==5

      a1=9.461e14;

      b1=124.6;

      c1=4.296;

      a2=0.4379;

      b2=91.79;

      c2=3.112;

      a3=1429;

      b3=189.1;

      c3=40.17;

      elseif sigma==6

      a1=6.402e15;

      b1=133.6;

      c1=5.704;

      a2=0;

      b2=100.5;

      c2=2.282;

      a3=122.7;

      b3=160.2;

      c3=39.06;

      elseif sigma==7

      a1=6.102e15;

      b1=129.6;

      c1=5.049;

      a2=0;

      b2=59.72;

      c2=25.03;

      a3=129.3;

      b3=161.5;

      c3=40.62;

      end

      end

      Lz=a1*exp(-((x-b1)/c1)^2)+a2*exp(-((x-b2)/c2)^2)+a3*exp(-((x-b3)/c3)^2);

      end

      kramp.m

      function [Kr] = kramp(x)

      K=2/sqrt(2*pi)*int(sym('exp(-t^2/2)'),'t',0,x);

      Kr=double(vpa(K,5));

      end

      poly_approx.m

      clc

      clear

      x = [100 150 200 250 300 350 400];

      y = [15 14.5 14 13.2 12.2 11 9.6];

      p = polyfit(x, y, 6)

      f = polyval(p, x);

      f = polyval(p, x);

      plot(x, y, 'ob', x, f, '-g')

      cftool

      LK_data_AM.m

      function varargout = LK_data_AM(varargin)

      gui_Singleton = 1;

      gui_State = struct('gui_Name', mfilename, …

      'gui_Singleton', gui_Singleton, …

      'gui_OpeningFcn', @LK_data_AM_OpeningFcn, …

      'gui_OutputFcn', @LK_data_AM_OutputFcn, …

      'gui_LayoutFcn', СКАЧАТЬ