Текст програми.
Program mzkt;
uses crt;
label 1, 2;
const E=7; tau1=1; tau2=4; T=10.0; n=20; Ut0=0;
var fajl : Text;
omega, dt, te : real;
C, B, A, fi, Ut, f : array[1..n] of real;
z1, z2, z3, Wser, C0 : real;
i, k : integer;
Begin
clrscr;
Assign(fajl, 'D:\Docs\mzoleg.txt');
rewrite(fajl);
omega:=2*pi/T;
for i:=1 to n do f[i]:=i/T;
for k:=1 to n do
begin
z1:=(1.0/(k*omega*tau1))*(cos(k*omega*tau1-1));
z2:=((sqr(tau1)-tau1+tau2)/tau2)*sin(k*omega*tau1);
z3:=(1.0/(k*omega*tau2))*(cos(k*omega*tau1)-cos(k*omega*(tau1+tau2)));
C[k]:=(2*E/(k*omega*T))*(z1+z2+z3)
end;
for k:=1 to n do
begin
z1:=(1.0/tau1)*sin(k*omega*tau1);
z2:=(1.0/tau2)*(sin(k*omega*(tau1+tau2))-sin(k*omega*tau1));
B[k]:=((-1.0)*2*E)/(sqr(k*omega)*T)*(z1+z2)
end;
for k:=1 to n do
begin
A[k]:=sqrt(sqr(C[k])+sqr(B[k]));
fi[k]:=arctan(B[k]/C[k])
end;
Wser:=0;
for k:=1 to n do Wser:=Wser+sqr(A[k]);
Wser:=1.0/2.0*Wser;
for i:=1 to n do writeln(fajl, i, 'f=',f[i]);
writeln(fajl,'');
for k:=1 to n do writeln(fajl,'C',k,'=',C[k]);
writeln(fajl,'');
for k:=1 to n do writeln(fajl,'B',k,'=',B[k]);
writeln(fajl,'');
for k:=1 to n do writeln(fajl, 'A',k,'=',A[k]);
writeln(fajl,'');
for k:=1 to n do writeln(fajl, 'fi',k,'=',fi[k]);
writeln(fajl,'');
writeln(fajl,' Wser=',Wser);
dt:=1.0*T/n;
te:=0;
k:=1;
C0:=(2*E/T)*(3.0/2.0*tau1+tau2-tau1/2.0*(2.0*tau1+tau2));
C0:=C0/2.0;
1: te:=te+dt;
if te>T then goto 2;
Ut[k]:=0;
Ut[k]:=Ut[k]+C0;
for i:=1 to n do Ut[k]:=Ut[k]+C[i]*cos(i*omega*te)+B[i]*sin(i*omega*te);
k:=k+1;
goto 1;
2: writeln(fajl,'');
te:=0;
for i:=1 to n do
begin
te:=te+dt;
writeln(fajl,'U(t=',te:4:2,')= ',Ut[i])
end;
writeln(fajl,'');
writeln(fajl,' END');
close(fajl)
end.
Результати роботи програми.
1f= 1.0000000000E-01
2f= 2.0000000000E-01
3f= 3.0000000000E-01
4f= 4.0000000000E-01
5f= 5.0000000000E-01
6f= 6.0000000000E-01
7f= 7.0000000000E-01
8f= 8.0000000000E-01
9f= 9.0000000000E-01
10f= 1.0000000000E+00
11f= 1.1000000000E+00
12f= 1.2000000000E+00
13f= 1.3000000000E+00
14f= 1.4000000000E+00
15f= 1.5000000000E+00
16f= 1.6000000000E+00
17f= 1.7000000000E+00
18f= 1.8000000000E+00
19f= 1.9000000000E+00
20f= 2.0000000000E+00
C1= 6.2175846582E+00
C2= 1.7639325551E+00
C3= 1.0239846668E+00
C4= 2.3992577411E-01
C5=-7.6641696819E-02
C6=-3.5461137902E-01
C7=-3.6023070248E-01
C8=-3.0955369066E-01
C9=-1.2823749842E-01
C10= 1.9160424207E-02
C11= 1.5962342671E-01
C12= 1.9615887932E-01
C13= 1.7992313848E-01
C14= 8.6406440604E-02
C15=-8.5157440943E-03
C16=-1.0102679538E-01
C17=-1.3440295095E-01
C18=-1.2655118705E-01
C19=-6.5052873531E-02
C20= 4.7901060538E-03
B1=-1.5633213091E+00
B2=-6.3237675337E-01
B3=-2.8105633483E-01
B4=-9.7707581820E-02
B5= 1.1611036042E-12
B6= 4.3425591919E-02
B7= 5.1622592112E-02
B8= 3.9523547085E-02
B9= 1.9300263076E-02
B10=-7.0927782976E-13
B11=-1.2920010819E-02
B12=-1.7566020927E-02
B13=-1.4967497120E-02
B14=-7.9761291282E-03
B15= 8.0273829427E-13
B16= 6.1067238635E-03
B17= 8.7526194241E-03
B18= 7.8071204116E-03
B19= 4.3305299422E-03
B20=-3.5463891488E-13
A1= 6.4111100831E+00
A2= 1.8738618992E+00
A3= 1.0618555745E+00
A4= 2.5905819544E-01
A5= 7.6641696819E-02
A6= 3.5726042625E-01
A7= 3.6391077344E-01
A8= 3.1206665662E-01
A9= 1.2968174951E-01
A10= 1.9160424207E-02
A11= 1.6014544962E-01
A12= 1.9694382709E-01
A13= 1.8054462531E-01
A14= 8.6773795663E-02
A15= 8.5157440943E-03
A16= 1.0121119237E-01
A17= 1.3468764446E-01
A18= 1.2679177447E-01
A19= 6.5196854557E-02
A20= 4.7901060538E-03
fi1=-2.4632923496E-01
fi2=-3.4423052506E-01
fi3=-2.6787637102E-01
fi4=-3.8673288209E-01
fi5=-1.5149763802E-11
fi6=-1.2185296801E-01
fi7=-1.4233519614E-01
fi8=-1.2699203524E-01
fi9=-1.4938287109E-01
fi10=-3.7017856291E-11
fi11=-8.0764502187E-02
fi12=-8.9311738034E-02
fi13=-8.2997181049E-02
fi14=-9.2048570126E-02
fi15=-9.4265196956E-11
fi16=-6.0373117289E-02
fi17=-6.5030401072E-02
fi18=-6.1613319947E-02
fi19=-6.6471314364E-02
fi20=-7.4035712550E-11
Wser= 2.3171107808E+01
U(t=0.50)= 8.6806163210E+00
U(t=1.00)= 5.1846561821E+00
U(t=1.50)= 1.6107332443E+00
U(t=2.00)= 1.7871462054E-01
U(t=2.50)= -9.2025528054E-01
U(t=3.00)= -1.8007564784E+00
U(t=3.50)= -2.5025583464E+00
U(t=4.00)= -3.0443319260E+00
U(t=4.50)= -3.4350621398E+00
U(t=5.00)= -3.6394032966E+00
U(t=5.50)= -2.9101006622E+00
U(t=6.00)= -1.9944184772E+00
U(t=6.50)= -9.2771763071E-01
U(t=7.00)= 2.9895325028E-01
U(t=7.50)= 1.7041701197E+00
U(t=8.00)= 3.3273356782E+00
U(t=8.50)= 5.2805671394E+00
U(t=9.00)= 9.2520901986E+00
U(t=9.50)= 1.0775904054E+01
U(t=10.00)= 9.9766655511E+00
END
Графіки