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 - дисперсия экспериментальных данных, вычисляемая по формуле