Проектирование электрической цепи
Дипломная работа - Компьютеры, программирование
Другие дипломы по предмету Компьютеры, программирование
);i:=1 to M do:=i*h;_An_t(t);_B_t(t);(n,A,B,h,X);:=X0+round(kx*t);:=Y01-round(ky1*x[1]);:=Y02-round(ky2*x[2]);:=Y01-round(ky1*d[1]);:=Y02-round(ky2*d[2]);:=Y03-round(ky3*U5(X));:=Y03-round(ky3*U5(D));(Red);(xx1,y11,xx,y1);(xx1,y21,xx,y2);(xx1,y51,xx,y5);(Green);(xx1,y31,xx,y3);(xx1,y41,xx,y4);(xx1,y61,xx,y6); xx1:=xx;y11:=y1;y21:=y2;y31:=y3;y41:=y4;y51:=y5;y61:=y6;;;;;
{****************************Метод установления****************}Ust(P:Vect);_A_B(P);[1]:=0;x[2]:=0;:=0;:=j+1;:=X;i:=1 to M do:=i*h;_B_t(t);(n,A,B,h,X);;:=true;i:=1 to n doabs(x[i]-x1[i])-delt*abs(x[i]+x1[i])/2>=0 then def:=false; until def;;:=(GetMaxX-X0-10)/TT/j;:=120/50;:=120/200;:=-100/2;[1]:=0;x[2]:=0;;:=X0;y11:=Y01;y21:=Y02;y31:=Y03;k:=0 to j-1 doi:=1 to M do:=i*h+k*TT;_B_t(t);(n,A,B,h,X);:=X0+round(kx*t);:=Y01-round(ky1*x[1]);:=Y02-round(ky2*x[2]);:=Y03-round(ky3*U5(X));(Red);(xx1,y11,xx,y1);(xx1,y21,xx,y2);(xx1,y31,xx,y3);:=xx;y11:=y1;y21:=y2;y31:=y3;;;;;
{*Определение коэффициента пульсаций по решению непосредственным методом} function Kp_nep(P:Vect):real;Usr,Umax,Umin:real;_A_B(P);_nep(X);:=0;Umin:=U5(X);Umax:=Umin;i:=1 to M do:=i*h;_B_t(t);(n,A,B,h,X);:=Usr+U5(X);U5(X)>Umax then Umax:=U5(X);U5(X)<Umin then Umin:=U5(X); end;_nep:=(Umax-Umin)*M/Usr/2;;
{***************Метод золотого сечения }Gold(P:Vect;i:byte):real;
const g=0.618034;r,xl,xp:real;,X2:Vect;:=pp[i];xl:=pl[i];X1:=P;X2:=P;r:=g*(xp-xl);x1[i]:=xl+r;x2[i]:=xp-r;kp_nep(X1)<kp_nep(X2) then xp:=x1[i] else xl:=x2[i]abs(xp-xl)-delt*abs(xp+xl)/2<0;:=(xp+xl)/2;;
{*************** Метод координатного спуска } Procedure Koord(var Pmin:Vect; var Kpmin:real);i:byte;,PL:Vect;:boolean;i:=1 to s do[i]:=(pl[i]+pp[i])/2;:=P;i:=1 to s do[i]:=Gold(P,i);;:=true;i:=1 to s doabs(p[i]-pl[i])-delt*abs(p[i]+pl[i])/2>=0 then def:=false; until def;i:=1 to s do[i]:=(p[i]+pl[i])/2;:=Kp_nep(Pmin);;
{************************ Управляющая программа *********** };(P,kpmax);{оптимизация методом координатного спуска} Def_A_B(P);{формирование матриц А,В математической модели} Det_L(A,Lam,def);{определение собственных чисел}('Оптимальная точка:'); writeln('p[1]=',p[1],'p[2]=',p[2],'kpmax=',kpmax); writeln('Собственные числа:'); {вывод собственных чисел на экран} if not def then
begin('l[1]=',lam[1]);writeln('l[2]=',lam[2]); writeln('Постоянные времени:'); writeln('tau1=',1/abs(lam[1]));writeln('tau2',1/abs(lam[2])); end
else('l[1]=',lam[1],'+j',lam[2]);writeln('l[1]=',lam[1],'-j',lam[2]); writeln('Постоянная времени: tau=',1/abs(lam[1]));('Собственная частота: fc=',lam[2]/2/pi);;;(P);{определение периодического решения методом установления} Per_Solve(P);{определение периодического решения непосредственным методом}.