?????????????? ?????? ? ????????? ?????????? ??????? ?????

Lj_1=Lj; //новая индуктивность становится начальной внутри итерации
end; //окончание цикла по итерациям
uL_0=zL*iL-eL; // расчет напряжения на индуктивности
iLg(m)=iL_0; // вывод тока во внешний массив
iL_0=iL; //новый ток в индуктивности
L_0=Lj; //индуктивность для следущего шага
ietalon(m)=tok(t(m)); //запись эталонного тока во внешний массив
end; //окончание внешнего цикла по времени
plot2d(t',[iLg, ietalon]) //построение графиков
xgrid(1); // построение сетки
Результаты расчетов для эталонного решения, а также для метода Эйлера
1-го порядка и метода трапеций при dt=0.002 и eps=0.000002 приведены
ниже:
1.0 L, Гн
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
i, A
0.0
0 1 2 3 4 5 6 7 8 9 10
зависимость индуктивности от тока
i , A
L 6
5
4
3
2
1
0
t, c
0.000 0.005 0.010 0.015 0.020 0.025 0.030
Результаты расчета переходных процессов:
Сплошная линия − аналитическое решение
• − метод трапеций;
▼ − неявный метод 1-го порядка.
53