Курсовая работа: Программирование и основы алгоритмизации
h – шаг сетки;
После задания всем переменным определенных значений происходит последовательное вычисление значений функции
while t0<=tk do
begin
Eu(t0,t,h,g,y0,y);
writeln(t,y[1],y[2]);
for i:=1 to g do
y0[i]:=y[i];
end;
readln;
end.
Для приведённой выше схемы составим систему дифференциальных уравнений:
Выразим неизвестные величины через данные начальных условий
1. Но
Итак,
2. Аналогично
3.
4.
5. Пусть
тогда
Таким образом имеем систему из пяти дифференциальных уравнений:
Итак, составим программу для решения дифференциальных уравнений и построения графиков переходных процессов.
Сначала сделаем некоторые переобозначения:
F — массив значений правых частей
F[1]=dZ1 /dt
F[2]=X6
F[3]=dα/dt