Исследование эффекта автодинного детектирования в многоконтурном генераторе на диоде Ганна

Страница: 10/13

outtextxy(540,75,s);

end;

end

else begin

if (visir_f<>0) then

begin

setcolor(0);

outtextxy(540,75,' _');

setcolor(13);

str((visir_f1/visir_f):5:3,s);

outtextxy(540,75,s);

end;

end;

end;

procedure v12; { вывод информации физиров 1 и 2 }

begin

d_visir:=1e-9*abs(visir_2-visir_1)*(xmax-xmin)/(xgmax- xgmin);

setcolor(0);

outtextxy(540,255,' _');

outtextxy(540,35,' _');

setcolor(15);

if(d_visir<>0) then begin

an2;

line(trunc(visir_s),ygmin,trunc(visir_s),ygmax);

visir_s:=xgmax-trunc((xmax-1/(d_visir*1e9))*(xgmax- xgmin)/(xmax-xmin));

line(trunc(visir_s),ygmin,trunc(visir_s),ygmax);

str((1e-9/d_visir):5:3,s);

outtextxy(540,35,s+' GHz');

end;

str(d_visir*1e9:5:4,s);

outtextxy(540,255,s+' ns');

end;

BEGIN

oldc1:=0;

oldc2:=0;

gd:=0;

InitGraph(gd,gm,'E: p-7gi');

an2; scal;

an4; scal;

an3; scal;

setcolor(11);

current;

kzp;

{ Начальные условия }

dh:=4;

dj:=2;

x:=0;

h:=8e-13;

y[1]:=eds;

w[1]:=eds;

y[3]:=eds; y[6]:=iv1;

w[3]:=eds; w[6]:=iv1;

y[2]:=eds; y[7]:=iv1;

w[2]:=eds; w[7]:=iv1;

y[5]:=eds; y[8]:=iv1;

w[5]:=eds; w[8]:=iv1;

y[4]:=eds; y[6]:=iv1;

w[4]:=eds; w[6]:=iv1;

y[11]:=eds; y[10]:=0;

y[9]:=iv1; w[9]:=iv1;

w[11]:=eds; w[10]:=0;

y[12]:=0; w[12]:=y[12];

y[13]:=eds; w[13]:=y[13];

y[14]:=0; w[14]:=y[14];

y[15]:=0; w[15]:=y[15];

loop:=1; { номеp pазвеpтки тока }

phas_x:=0; phas_y:=0; { сдвиг фазового поpтpета }

size_x:=1;size_y:=1; { масштаб фазового портрета }

an2;

visir_s:=800;

visir_3:=xgmin;

visir_f:=0;

visir_4:=xgmin;

visir_f1:=0;

an3;

visir_1:=xgmin;

visir_2:=xgmin; { визиры }

count:=1;

mark:=0;

round:=0;

old_cur:=iv1;

fcount:=0;

fsign:=1;

fpoint:=1;

frequency:=1e10;

old_f:=1e10;

Smax:=0;

power:=0;

oldx:=xgmax-trunc((xmax-0)*(xgmax-xgmin)/(xmax-xmin));

for aaa:=1 to 10 do

oldy[aaa]:=ygmin-trunc((ymax-y[8]*10)*(ygmin- ygmax)/(ymax-ymin));

{ Рунге-Кутт }

for iii1:=-249 to maxpoint do begin

for iii:=0 to 4 do begin

anna(y,f);

for k:=1 to n do begin

K1[k]:=f[k]*h;

y[k]:=w[k]+h*f[k]/2;

end;

x:=x+h/2;

anna(y,f);

for k:=1 to n do begin

K1[k]:=K1[k]+2*f[k]*h;

y[k]:=w[k]+f[k]*h/2;

end;

anna(y,f);

for k:=1 to n do begin

K1[k]:=K1[k]+2*f[k]*h;

y[k]:=w[k]+f[k]*h;

end;

x:=x+h/2;

anna(y,f);

for k:=1 to n do begin

y[k]:=w[k]+(K1[k]+f[k]*h)/6;

w[k]:=y[k];

end;

end;

{ вычисление мощности }

power:=power+y[8]*y[2];

{ вычисление частоты по изменению знака производной }

if fsign > 0 then begin

if y[8]-old_cur <= 0 then begin

if fcount = 0 then fpoint:=iii1;

fcount:=fcount+1;

fsign:=-1;

end;

end

else begin

if y[8]-old_cur >= 0 then begin

if fcount = 0 then fpoint:=iii1;

fcount:=fcount+1;

fsign:=1;

end;

end;

old_cur:=y[8];

if fcount = 15 then begin { Частота сигнала }

fcount:=1;

mark:=1;

old_f:=frequency;

frequency:=(iii1-fpoint)/(h*4.2e3 * 5);

fpoint:=iii1;

power:=power *h*frequency/5;

str(power:5:4,s);

power:=0;

setcolor(0); outtextxy(250,460,' ');

setcolor(11);

outtextxy(250,460,'Puhf = '+s+' W');

end;

{ вывод графиков токов и напряжений }

if(iii1>0) then begin

an3;

if(iii1=loop*1000) then begin

loop:=loop+1;

setfillstyle(0,0);

bar(xgmin+1,ygmin+1,xgmax-1,ygmax-1);

scal;

setwritemode(XORput);

setcolor(15);

line(visir_1,ygmin,visir_1,ygmax);

line(visir_2,ygmin,visir_2,ygmax);

setwritemode(COPYput);

str(d_visir*1e9:5:4,s);

outtextxy(540,255,s+' ns');

round:=round+1;

setcolor(0);

outtextxy(50,460,' ');

Реферат опубликован: 30/10/2008