Определение ускорения свободного падения с помощью физического маятника. Период колебания физического маятника. Нахождение ускорения свободного падения методом наименьших квадратов. Решение задач методами Гаусса-Ньютона и квазиньютоновскими методами.
Аннотация к работе
Расчет методом Левенберга-Марквардта Условия останова: Достигнуто значение нормы градиента: 1,000E-3 Начальное значение параметра мю: 10000,000 Количество параметров: 1 Количество уравнений: 9 Начальное значение параметров: A = (6,000)------------------------------------------------Итерация 0 Текущее приближение: A = (6,000) Вектор невязок: R = (7,963 7,026 5,771 5,898 6,208 8,870 5,652 5,438 9,469) Значение невязки: 443,532 Якобиан:-3,031-2,663-2,159-2,213-2,336-3,386-2,109-2,020-3,618 Норма градиента: 170,058 Значение параметра мю: 5000,000 Итерация 1 Текущее приближение: A = (6,017) Вектор невязок: R = (7,912 6,981 5,734 5,861 6,168 8,813 5,616 5,404 9,408) Значение невязки: 432,467 Якобиан:-3,018-2,652-2,150-2,204-2,326-3,372-2,100-2,011-3,603 Норма градиента: 168,261 Значение параметра мю: 2500,000 Итерация 2 Текущее приближение: A = (6,050) Вектор невязок: R = (7,812 6,893 5,663 5,788 6,091 8,701 5,547 5,338 9,289) Значение невязки: 411,699 Якобиан:-2,993-2,630-2,132-2,186-2,307-3,344-2,083-1,995-3,573 Норма градиента: 164,782 Значение параметра мю: 1250,000 Итерация 3 Текущее приближение: A = (6,114) Вектор невязок: R = (7,621 6,726 5,527 5,649 5,944 8,488 5,414 5,210 9,061) Значение невязки: 374,909 Якобиан:-2,946-2,589-2,099-2,151-2,270-3,291-2,050-1,963-3,517 Норма градиента: 158,247 Значение параметра мю: 625,000 Итерация 4 Текущее приближение: A = (6,235) Вектор невязок: R = (7,271 6,418 5,277 5,393 5,674 8,096 5,170 4,977 8,642) Значение невязки: 315,981 Якобиан:-2,861-2,514-2,038-2,089-2,205-3,196-1,991-1,907-3,416 Норма градиента: 146,646 Значение параметра мю: 312,500 Итерация 5 Текущее приближение: A = (6,450) Вектор невязок: R = (6,671 5,891 4,850 4,955 5,212 7,427 4,753 4,577 7,927) Значение невязки: 235,516 Якобиан:-2,719-2,389-1,937-1,986-2,095-3,038-1,892-1,812-3,246 Норма градиента: 127,955 Значение параметра мю: 156,250 Итерация 6 Текущее приближение: A = (6,801) Вектор невязок: R = (5,753 5,085 4,197 4,285 4,505 6,402 4,114 3,966 6,831) Значение невязки: 147,427 Якобиан:-2,512-2,207-1,789-1,834-1,935-2,805-1,748-1,674-2,998 Норма градиента: 102,020 Значение параметра мю: 78,125 Итерация 7 Текущее приближение: A = (7,310) Вектор невязок: R = (4,543 4,021 3,334 3,401 3,572 5,049 3,271 3,159 5,386) Значение невязки: 73,516 Якобиан:-2,254-1,980-1,606-1,646-1,737-2,518-1,568-1,502-2,690 Норма градиента: 72,432 Значение параметра мю: 39,063 Итерация 8 Текущее приближение: A = (7,947) Вектор невязок: R = (3,195 2,837 2,374 2,417 2,533 3,543 2,333 2,261 3,777) Значение невязки: 27,360 Якобиан:-1,988-1,747-1,416-1,452-1,532-2,221-1,384-1,325-2,374 Норма градиента: 45,118 Значение параметра мю: 19,531 Итерация 9 Текущее приближение: A = (8,623) Вектор невязок: R = (1,931 1,726 1,474 1,494 1,559 2,132 1,454 1,418 2,268) Значение невязки: 6,854 Якобиан:-1,759-1,546-1,253-1,285-1,356-1,965-1,224-1,172-2,100 Норма градиента: 24,342 Значение параметра мю: 9,766 Итерация 10 Текущее приближение: A = (9,213) Вектор невязок: R = (0,942 0,858 0,770 0,772 0,798 1,028 0,766 0,760 1,088) Значение невязки: 0,998 Якобиан:-1,593-1,399-1,135-1,163-1,227-1,779-1,108-1,061-1,901 Норма градиента: 11,006 Значение параметра мю: 4,883 Итерация 11 Текущее приближение: A = (9,613) Вектор невязок: R = (0,326 0,316 0,330 0,322 0,323 0,339 0,337 0,349 0,352) Значение невязки: 0,100 Якобиан:-1,495-1,313-1,065-1,091-1,152-1,670-1,040-0,996-1,784 Норма градиента: 3,871 Значение параметра мю: 2,441 Итерация 12 Текущее приближение: A = (9,801) Вектор невязок: R = (0,048 0,072 0,133 0,119 0,109 0,029 0,144 0,164 0,021) Значение невязки: 0,042 Якобиан:-1,452-1,275-1,034-1,060-1,119-1,622-1,010-0,967-1,733 Норма градиента: 0,935 Значение параметра мю: 1,221 Итерация 13 Текущее приближение: A = (9,856) Вектор невязок: R = (-0,030 0,003 0,077 0,062 0,048-0,059 0,089 0,112-0,073) Значение невязки: 0,041 Якобиан:-1,440-1,265-1,026-1,051-1,109-1,608-1,002-0,959-1,719 Норма градиента: 0,135 Значение параметра мю: 0,610 Итерация 14 Текущее приближение: A = (9,864) Вектор невязок: R = (-0,043-0,007 0,068 0,053 0,039-0,073 0,081 0,103-0,087) Значение невязки: 0,041 Якобиан:-1,438-1,263-1,024-1,050-1,108-1,606-1,001-0,958-1,716 Норма градиента: 0,011 Значение параметра мю: 0,305------------------------------------------------Расчет окончен Достигнута заданная норма градиента.