n> i + h Г— k 2 /2), 4 = f (x i + h, y i + h Г— k 3 ).
9. Практична частина
1 Чисельна реалізація рішення системи диференціальних рівнянь
.1 Реалізація в програмі MathCad. Методу Рунге-Кутта. br/>В
Метод Рунге-Кутта
В
Таблиця 1. Результати роботи програми MathCad. <В
1.2 Реалізація в програмі MathCad. Метод Ейлера 1-а модифікація
В
Таблиця 2. Результати роботи програми MathCad. <В
В
Рисунок 2 - Графік залежності сили струму від часу
В
Рисунок 3 - Графік Залежності напруги від часу
1.3 Реалізація рішення на мові програмування високого рівня (C + +) методом Ейлера (1 модифікація).
Блок-схема
В
Програмний код
# include
# include
# include
# include double PI = 3.1415926;
double f, PhiA, Phi, E0, LA, L, CA, C , T0, t1, t2, h;
double R [7]; n, i;
// Функція гармонійного джерела струму, що враховує перемикання ключа
double E (double t)
{; (t
{= E0 + E0 * sin (2 * PI * f * t + Phi);
}
{= 0;
};
}
// Похідна від струму I за часом t
double FDI (double t, double I, double U)
{A, B, D, ResultFDI; = R [2]/(R [1] + R [2]); = (R [1] * R [2])/(R [1] + R [2]); = (R [5] + R [6])/(R [3] + R [5] + R [6]); = (1/L) * ((A * E (t)) - (I * (R [4] + (R [3] * D) + B)) - (U * D));
returnResultFDI;
}
// Похідна від напруги U за часом t
double FDU (double t, double I, double U)
{A, B, ResultFDU; = (R [5] + R [6])/(R [3] + R [5] + R [6]); = 1/(R [3] + R [5] + R [6]); = (1/C) * ((I * A) - (U * B));;
} main ()
{* FOUT; xt1, xt2, pI, I1, I2, pU, U1, U2; = fopen ("C: student Result08. txt", "a");
// Вводпараметровцепі <<"********** METO...
