i + k-2)); [i, k]: = s;; i: = 1 to m +1 do: = 0; j: = 1 to n do: = s + ((power (x [j], i-1)) * y [j]); [i]: = s;;
{Метод Гаусса} i: = 1 to (m +1) do [i, i]: = 1; [1, i]: = a [1, i]/a [1,1]; [i , 1]: = a [i, 1];; [1]: = d [1]/a [1,1]; i: = 2 to (m +1) dok: = i to (m +1) do: = 0; j: = 1 to i-1 do: = s + c [k, j] * b [j, i]; [k, i]: = a [k, i]-s;; k : = i +1 to (m +1) do: = 0; j: = 1 to i-1 do: = s + c [i, j] * b [j, k]; [i, k]: = (a [i, k] - s)/c [i, i];;; i: = 2 to (m +1) doj: = 1 to i-1 do: = s + c [i, j] * YY [j]; [i]: = (d [i]-s)/c [i, i];; [n]: = YY [n]; i: = 1 to n-1 do: = 0; j: = n-i +1 to n do: = s + b [ni, j] * XX [j]; [ni]: = YY [ni]-s;; i: = 1 to m +1 do_StringGrid1. Cells [Kovalenko_StringGrid1.colCount-
q: = x [1]; _chart1.SeriesList.Series [0]. Clear; k: = 1 to 1500 doi: = 1 to (m +1) do: = s + XX [i] * (power (q, i-1)); {y}: = s;
1,0]: = FloatToStrF (q, ffGeneral, 6,6); _StringGrid2.Cells [Kovalenko_StringGrid2.colCount-
1,1]: = FloatToStrF (f, ffGeneral, 6,6); _StringGrid2.ColCount: = Kovalenko_StringGrid2.ColCount +1;;;.
3.3 Блок-схема
3.4 Результати
В
Малюнок 3
В
В
Стартові точки апроксимуючої функції:
, 5; 5.5, 1; 0.55, 1.5; -1, 2; 11.5, 2.5; 9.5, 3; -10.2, 3.5; 1.5.
апроксимується значення функції:
.15, 2281.1, -3643.52, 2768.05, -1075.68, 206.173, -15.4311;
Висновок: при розрахунку даним методом досягається високий ступінь точності для апроксимації різних функцій. Метод може застосовуватися для вирішення інженерних завдань всіляких видів. p align="justify"> апроксимація золотий перетин градієнтний
4. Одновимірна оптимізація багатоекстремального функцій. Метод золотого перерізу
.1 Короткі теоретичні відомості
Якщо на заданому інтервалі функції не унімодальних, то при використанні цього методу спочатку вихідний інтервал існування функції звужують до інтервалу унімодальних (існування одного екстремуму). Золотим перетином відрізка називається розподіл відрізка на дві нерівні частини так, відношення довжини всього відрізка до довжини більшої його частини дорівнює відношенню довжини більшої частини до меншої. p> Очевидно, що існує два варіанти розбиття:
В
З урахуванням використання лічильника ітерацій в точці
,.
Алгоритм реалізації методу:
. Приймають базові значення інтервалу унімодальних,, лічильник.
. Задають значення точності.
. Обчислюють допоміжні значення
,
4.Якщо, то,,. Інакше,,. p>. Обчислюють довжину звужуватися інтервалу.
. Якщо, то оптимальною точкою буде. Інакше обчислення продовжують, починаючи з пункту 3. p ali...