g 3uLEPA (3 mog.) ********** n"; << ; "BBEguTE E0 (B): n";
cin>> E0; <<"BBEguTE L (MuLurEHPu): n";>> LA; = LA/1000; <<"BBEguTE C ( MuKPOqpAPAgbI): n ";>> CA; = CA/1000000; <<" BBEguTE t0 (c): n ";>> t0; <<" BBEguTE t1 (c ); n ";>> t1; <<" BBEguTE t2 (c): n ";>> t2; <<" BBEguTE n: n ";
cin>> n; = (t2 - t0)/n; <<"War h =" <
{<<"BBEguTE R (" <> R [i]; p>
} <<"BBEguTE f (Gertz): n";>> f; <<"BBEguTE Phi (* PI PAguAH): n",> ;> PhiA; = PhiA * PI;
// Визначення початкових умов <<"BBEguTE I (0), (A): n";
cin>> I1; <<"BBEguTE U (0), (B): n";
cin>> U1; = t0; <<" n";
// Запис введених даних у файл (FOUT, "********** Метод Ейлера (3 мод.) ********** n" ); (FOUT, "Параметри ланцюга: n");
fprintf (FOUT, "E0 =% g (B) n", E0);
fprintf (FOUT, "L =% g (мГн) n", LA);
fprintf (FOUT, "C =% g (мкФ) n", CA);
fprintf (FOUT, "t0 =% g (c) nt1 =% g (c) nt2 =% g (c) n", t0, t1, t2); p>
fprintf (FOUT, "n =% d n", n); (FOUT, "h =% g n", h); (i = 1; i <= 6 ; i + +)
{(FOUT, "R (% d) =% g (Ом) n", i, R [i]);
} (FOUT, "f =% g (Гц) n", f); (FOUT, "ф =% g * п =% g (рад.) n", PhiA, Phi); (FOUT, " nТабліца значень: n"); (FOUT, "i tt tI tU n"); (FOUT, "0 t% g t% g t% g n ", xt1, I1, U1);
// Цикл розрахунку значень функцій 3 модифікацією методу Ейлера (i = 1; i <= n; i + +)
{
// Вичісленіезначеній = xt1 + h; = I1 + h * FDI (xt1, I1, U1);
pU = U1 + h * FDU (xt1, I1, U1);
I2 = I1 + (h/2) * (FDI (xt1, I1, U1) + FDI (xt2, pI, pU));
U2 = U1 + (h/2) * (FDU (xt1, I1, U1) + FDU (xt2, pI, pU));
// Висновок отриманих значень на екран <<"i =" <
// Висновок отриманих значень у файл (FOUT, "% d t% g t% g t% g n", i, xt2, I2, U2);
// Зсув значень для наступної ітерації = xt2;
I1 = I2; = U2;
} <<" n n n"; (FOUT, " nРасчетокончен n n n");
fclose (FOUT);
}
Аналіз результатів.
Чисельне рішення системи диференціальних рівнянь у програмах MathCad і C + + було реалізовано за допомогою методу Ейлера 1-й модифікації. p align="justify"> Чисельні результати двох програм збігаються. p> Отримано графіки залежності сили струму і напруги від часу відповідно до системи рівнянь і...