105
, (72)
где
- вектор столбец неизвестных параметров;
- экспериментальное относительное давление в момент времени tk;
- расчетное относительное давление в момент времени tk, получаемое численным интегрированием уравнения (65);
N - число точек эксперимента, пo которым производят оптимизацию.
1.3.3.1. Метод нелинейных оценок.
В данном методе минимизацию функционала (72) производят итеративно, причем
, (73)
где L - номер итерации. Значение поправки
определяют при каждой итерации из выражения
, (74)
где
- матрица размерностью 2хN;
- транспонированная матрица
;
- вектор-столбец размерностью N;
- матрица, обратная матрице
.
Транспонирование, перемножение и обращение матриц осуществляют по стандартным процедурам, имеющимся в библиотеках подпрограмм ЭВМ.
Итерационный процесс прекращается при выполнении условия
, (75)
где
.- Элементы матрицы
определяются численным интегрированием уравнения
, (76)
с. начальными условиями
. В уравнении (76) символ f использован для обозначения правой части уравнения (65).
В процессе вычислений для нахождения доверительных интервалов (см. формулу (81) найденных параметров получают информационную матрицу Фишера
, (77)
где S2 - дисперсия экспериментальных данных, вычисляемая по формуле