Page 129 - 480

Поточний час знаходиться за формулою: t  t  t  ,
i i1
i  n , 1 ; t  0 .
0
У випадку, якщо відлік значень F i відбувається у довільні
моменти часу, то повинен бути сформований також масив значень
моментів t i відліку значень F i .
Алгоритм розв’язку задачі можна представити в такій
послідовності.
Крок 1. Розраховуються початкові стартові значення
параметрів рівняння (2.9). З цією метою розв’язується нелінійне
рівняння (2.24) відносно , а потім обраховується К з (2.23). Для
числового рішення рівняння (2.24) використовується метод
дихотомії.
Величина Е забезпечує задану точність розв’язку рівняння
(2.24). Рекомендується наступне значення величини Е=0,001.
1 N 2
Крок 2. Знаходимо значення Q   F  F  tc,  .
1 i i
2  i 1
Крок 3. Вираховується норма градієнта функції нев’язки.
Якщо  J  ,tc   E , то кінець обрахунків. В іншому
i 2
випадку проводиться перехід до наступного кроку. Величина Е
визначає точність розв’язку задачі. Рекомендується приймати
E  1 . 0 .
2
Крок 4. Визначає оптимальну довжину кроку шляхом
розв’язку задачі min J c  p, i  t .

Для розв’язку такої задачі використовується метод “золотого
перерізу”.
Підпрограма складається з двох блоків. Перший блок
здійснює пошук відрізка  ,0 , що містить точку мінімуму. З цією
метою послідовним збільшенням змінної  досягають виконання
співвідношення   JJ   . Другий блок здійснює пошук
i  i 1
мінімуму функції cJ   p,  t по змінній  . Процес обрахунку
i
закінчується після того, як точка мінімуму локалізована на відрізку,
довжина якого не перевищує E . Рекомендовано приймати
1
E 1=0.001.
Крок 5. Визначаємо значення К та  на наступній ітерації.
K  K   p ,      p , n=0,1,2,...
i 1 i i i 1 i 1 i i 2i
де складові p та p вектора p визначаються з рівняння (11).
1 2
126