Уравнений на основе итерационной процедуры метода Ньютона
Поможем в ✍️ написании учебной работы
Поможем с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой
0 1,681 -0,901 -2,211 2,539 -35,61 68,130 39,01 -90,76    
1 1,633 -0,894 -0,034 0,045 -35,05 62,53 38,26 -83,31 0,0481 0,0073
2 1,633 -0,894 0,000 0,000 -35,01 62,40 38,20 -83,12 0,0000 0,0005
3 1,633 -0,894 0,000 0,000 -35,01 62,40 38,20 -83,12 0,0000 0,0000

Таким образом, минимум суммы квадратов отклонения значений  экспоненциальной зависимости от табличных данных  достигается при значениях её параметров  и  (с погрешностью, не превышающей 0,001). При этом остаточная сумма квадратов составляет

,

что в относительных единицах равно 14,0 %.

Полученная величина  может служить критерием для сравнения с остаточной суммой квадратов, соответствующей какой-либо другой оценке параметров экспоненциальной зависимости, вычисленной тем или иным способом.

Рассматриваемая задача однофакторной нелинейной регрессии была сведена к решению нормальной системы нелинейных уравнений методом Ньютона. Другим известным подходом к решению задачи нелинейного оценивания является линейная параметризация математической модели [6].

Разложим функцию  в ряд Тейлора в окрестности точки :

Ограничиваясь линейной частью разложения, получаем аппроксимацию:

,

где , .

В этом случае сумма квадратов остатков (1.2) будет описываться выражением вида

. (1.6)

Очевидно, что функционал (1.6) представляет собой квадратичную форму относительно параметров  и  и, следовательно, система нормальных уравнений  и  будет линейной относительно этих параметров:

;

.

Отсюда имеем:

;

.

Вводя обозначения , , , , , , получаем систему уравнений

                   (1.7)

из решения которой находим

, .

Вычисляя новые приближения по формулам  и , получаем итерационный процесс уточнения параметров заданной нелинейной зависимости,  – номер итерации. Условием завершения итерационной процедуры может служить выполнение неравенств: , , где  – заданная величина погрешности вычислений.

В качестве начального приближения  можно взять значения, полученные выше по формулам (1.5): , .

Данные для вычисления коэффициентов нормальной системы уравнений (1.7) на первой итерации  представлены в табл. 1.4.

Таблица 1.4

Данные для вычисления коэффициентов нормальной системы
уравнений (1.7) на первом шаге итерационной процедуры

1 -1,5 6,50 14,924 -37,632 -22,387 56,448 0,023 -0,034
2 -1,2 4,30 8,692 -17,533 -10,430 21,040 -1,934 2,321
3 -0,9 3,53 5,062 -7,659 -4,556 6,893 -0,567 0,511
4 -0,6 2,93 2,948 -2,974 -1,769 1,784 0,075 -0,045
5 -0,3 2,42 1,717 -0,866 -0,515 0,260 0,285 -0,085
6 0,0 2,01 1,000 0,000 0,000 0,000 0,329 0,000
7 0,3 0,68 0,582 0,294 0,175 0,088 -0,460 -0,138
8 0,6 1,67 0,339 0,342 0,204 0,205 0,402 0,241
9 0,9 0,03 0,198 0,299 0,178 0,269 -0,319 -0,287
10 1,2 0,57 0,115 0,232 0,138 0,278 0,000 0,000
  35,578 -65,497 -38,963 87,266 -2,166 2,483

 

Система нормальных уравнений

имеет решение:  и . Отсюда новые уточненные оценки параметров:  и .

Данные для вычисления коэффициентов нормальной системы уравнений (1.7) на второй итерации  представлены в табл. 1.5.

Таблица 1.5

Данные для вычисления коэффициентов нормальной системы
уравнений (1.7) на втором шаге итерационной процедуры

1 -1,5 6,50 14,608 -35,789 -21,912 53,683 0,984 -1,477
2 -1,2 4,30 8,544 -16,746 -10,253 20,095 -1,386 1,663
3 -0,9 3,53 4,998 -7,346 -4,498 6,611 -0,271 0,244
4 -0,6 2,93 2,923 -2,864 -1,754 1,719 0,235 -0,141
5 -0,3 2,42 1,710 -0,838 -0,513 0,251 0,372 -0,112

 

Окончание табл. 1.5

 

6 0,0 2,01 1,000 0,000 0,000 0,000 0,377 0,000
7 0,3 0,68 0,585 0,287 0,175 0,086 -0,435 -0,131
8 0,6 1,67 0,342 0,335 0,205 0,201 0,418 0,251
9 0,9 0,03 0,200 0,294 0,180 0,265 -0,313 -0,282
10 1,2 0,57 0,117 0,229 0,140 0,275 0,004 0,005
  35,027 -62,438 -38,229 83,187 -0,015 0,020

 

Система нормальных уравнений ;

имеет решение:  и . Отсюда новые, уточненные оценки параметров:  и . Так как  и , эти оценки с погрешностью, не превышающей , обеспечивают минимум остаточной суммы квадратов (1.6). Очевидно, что полученный результат  и  совпадает с результатом вычислений методом Ньютона, причем в обоих случаях потребовалось по две итерации. Следовательно, и остаточные суммы квадратов будут одинаковыми: .

Рассмотрим решение задачи из примера 1.1 на основе линеаризующих преобразований нелинейной зависимости посредством логарифмирования [2, 6]. Целью таких преобразований является сведение системы нормальных уравнений к линейному виду, что существенно упрощает её решение.

Представим табличные данные  в виде

, .      (1.8)

Прологарифмируем обе части выражения (1.8):

Обозначая , ,   и , получаем линейную по параметрам  и  модель, описывающую табличные данные :

,   .          (1.9)

Минимизируя квадрат “невязки”: , получаем линейную систему нормальных уравнений

 

                       (1.10)

 

коэффициенты которой вычисляются по формулам:

 

; . (1.11)

 

С учетом введенных обозначений оценки параметров экспоненциальной зависимости находятся по формулам  и .

Для вычисления коэффициентов в системе уравнений (1.10) используем данные из табл. 1.2 и формируем табл. 1.6.


Таблица 1.6




Дата: 2019-07-30, просмотров: 167.