Исследование диаграмм направленности и коэффициента направленного действия апертурных антенн СВЧ диапазона. Струков И.Ф - 35 стр.

                                               35
begin Line (GetMaxX div 2, 0, Xm div 2, Ym);
        Line (0, Round(0.9*Ym), Xm, Round(0.9*Ym));
        OutTextXY (Xm div 2+8, Round(Ym*0.45), '0.5');
        OutTextXY (Xm div 2+8, Round(Ym*0.26), '0.7');
        OutTextXY (Xm-35, Round(0.9*Ym)-10, 'Teta');
        OutTextXY (0, Round(0.9*Ym)+5, '-pi');
        OutTextXY (Xm div 2 +5, Round(0.9*Ym)+5, '0');
        OutTextXY (Xm-15, Round(0.9*Ym)+5, 'pi');
end;
BEGIN Assign(Q, 'e:\FileOfQ.txt');
        Rewrite(Q);
        Writeln('Д Н волноводной илир у пор ной а нтенны пр и лю бы х ква др а тичны х
ф а зовы х иска ж ения х');
        WriteLn('Введите р а змер ы а нтенны идлину волны в одних единиц а х:');
        Write('D1='); ReadLn(D1);
        Write('D2='); ReadLn(D2);
        Write('R1='); Readln(R1);
        Write('R2='); Readln(R2);
        Write('lambda='); ReadLn(lambda);
        D1:=D1/lambda; {э лектр ические р а змер ы а нтенны }
        D2:=D2/lambda;
        R1:=R1/lambda;
        R2:=R2/lambda;
          {==============П остр оение за висимости Fe(Teta)=============}
          i:=-1;
          di:=0.005;
          Max:=0;
          while i<=1 do
            begin v1:=sqrt(2/R2)*(-D2/2-R2*sin(pi*i));
                     v2:=sqrt(2/R2)*(D2/2-R2*sin(pi*i));
                     y:=sqr(cos(pi*i/2))*sqrt(sqr(C(v2)-C(v1))+sqr(S(v2)-S(v1)));
                     Write(Q, y);
                     if y>Max then Max:=y;
                     i:=i+di;
            end;
          {-------------------------------------------}
          init;
          Xm:=GetMaxX; Ym:=GetMaxY;
          SetColor (White);
          XYplot;
          OutTextXY (Xm div 2+10, 5, 'Fe(Teta)');
          i:=-1;
          Reset(Q);
          Read(Q, y1);
          While i<=1 do begin Read(Q, y);
                     Line(Round(Xm*(1+i)/2), Round(Ym*0.9*(1-y1/Max)), Round(Xm*(1+i+di)/2),
                     Round(Ym*0.9*(1-y/Max))); y1:=y;
                     i:=i+di;
            end;
          OutTextXY (10, Ym-15, 'Press any key');
          ReadLn;
      {===============П остр оение за висимости Fh(Teta)==============}
          Rewrite(Q); i:=-1; di:=0.005; Max:=0;
          while i<=1 do
            begin v3:=sqrt(2/R1)*(D1/2+R1/2*(1/D1-2*sin(pi*i)));
                   v4:=sqrt(2/R1)*(-D1/2+R1/2*(1/D1-2*sin(pi*i)));
                   v5:=sqrt(2/R1)*(D1/2-R1/2*(1/D1+2*sin(pi*i)));
                   v6:=sqrt(2/R1)*(-D1/2-R1/2*(1/D1+2*sin(pi*i)));
                   v1:=cos(pi*R1/4*sqr(1/D1-2*sin(pi*i)))*(C(v3)-
                     C(v4))+cos(pi*R1/4*sqr(1/D1+2*sin(pi*i)))*(C(v5)-
                     C(v6))+sin(pi*R1/4*sqr(1/D1-2*sin(pi*i)))*(S(v3)-
                     S(v4))+sin(pi*R1/4*sqr(1/D1+2*sin(pi*i)))*(S(v5)-S(v6));
                     v2:=-sin(pi*R1/4*sqr(1/D1-2*sin(pi*i)))*(C(v3)-C(v4))-
                     sin(pi*R1/4*sqr(1/D1+2*sin(pi*i)))*(C(v5)-C(v6))+cos(pi*R1/4*sqr(1/D1-
                     2*sin(pi*i)))*(S(v3)-S(v4))+cos(pi*R1/4*sqr(1/D1+2*sin(pi*i)))*(S(v5)-(v6));
                   y:=sqr(cos(pi*i/2))*sqrt(sqr(v1)+sqr(v2));