програмний код чисельне інтегрування В
Третій ділянку від t3 = 0.0006 до t4 = 0.0099
В
Поєднання отриманих апроксимуючих функцій
В
Аналіз результатів.
Рішення завдання апроксимації було проведено в програмі MathCad методом найменших квадратів і в пакеті Excel з використанням майстра діаграм з висновком рівняння лінії тренду. Результати двох програм збігаються. Отримано аналітичну формула залежності сили струму від часу на інтервалі 0 Чисельне інтегрування
.1 Реалізація в програмі С + +. Метод лівих прямокутників. p> 3.1.1 Блок-схема
В
Малюнок 16 - Блок-схема. Метод лівих прямокутників
3.1.2 Код програми
# include "stdafx. h "
# include
# include
# include ;
// Квадрат аналітичної функції I (t), отриманої в пакеті Excel
double I (double t)
{doubleIt, Res; (t <0.0004)
{It = 4 * (1E +9) * (pow (t, 3)) - 5 * (1E +6) * pow (t, 2) +1873.5 * t +0.0002;}
else (t <0.001)
{It = 2 * (1E +8) * pow (t, 3) - 513 156 * pow (t, 2) - 346.61 * t +0.1905;}
{It = -2 * (1E +6) * pow (t, 3) +19110 * pow (t, 2) - 99.733 * t +0.3122;}
Res = pow (It, 2); Res;}
void main ()
{double R4, T1, T2, h, S, Integ, Q;
inti, n; <
cout <<"Vveditesoprotivltnie R4 (OM)" <> R4; <<"Vvedite T1 (c):" << ; endl;>> T1; <<"Vvedite T2 (c):" <> T2; <<"Vveditekolichectvorasbienii:" <> n; = (T2 - T1)/n; <<"h =" <
S = 0; (i = 1; i <= (n-1); i + +)
{S = S + I (T1 + h * i);} = h * S; = R4 * Integ; <<"Integ =" <
В
Метод правих прямокутників.
# include "stdafx. h "
# include
# include
# include ;
// Квадрат аналітичної функції I (t), отриманої в пакеті ExcelI (double t)
{double It, Res; (t <0.0004)
{It = 4 * (1E +9) * (pow (t, 3)) - 5 * (1E +6) * pow (t, 2) +1873.5 * t +0.0002;} (t <0.001)
{It = 2 * (1E +8) * pow (t, 3) - 513 156 * pow (t, 2) - 346.61 * t +0.1905;}
{It = -2 * (1E +6) * pow (t, 3) +19110 * pow (t, 2) - 99.733 * t +0.3122;}
Res = pow (It, 2); Res;}