2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки


Правила форума


В этом разделе нельзя создавать новые темы.

Если Вы хотите задать новый вопрос, то не дописывайте его в существующую тему, а создайте новую в корневом разделе "Помогите решить/разобраться (М)".

Если Вы зададите новый вопрос в существующей теме, то в случае нарушения оформления или других правил форума Ваше сообщение и все ответы на него могут быть удалены без предупреждения.

Не ищите на этом форуме халяву, правила запрещают участникам публиковать готовые решения стандартных учебных задач. Автор вопроса обязан привести свои попытки решения и указать конкретные затруднения.

Обязательно просмотрите тему Правила данного раздела, иначе Ваша тема может быть удалена или перемещена в Карантин, а Вы так и не узнаете, почему.



Начать новую тему Ответить на тему На страницу 1, 2  След.
 
 метод наименьших квадратов, логарифм
Сообщение16.01.2009, 14:48 


16/01/09
14
есть набор точек $(x_i,y_i)$ измеренных с некоторой погрешностью
надо аппроксимировать их логарифмом $y=a\log(x)+b$ по методу наименьших квадратов
для этого нужна формула по которой следует оценивать расстояние от точки до логарифма (до ближайшей касательной к логарифму) но вывести ее мне что-то не удалось, найти тоже
случай вроде не из сложных, наверное где-то есть готовое решение

 Профиль  
                  
 
 
Сообщение16.01.2009, 14:54 
Супермодератор
Аватара пользователя


29/07/05
8248
Москва
 !  PAV:
Вы неправильно набираете формулы. Из-за этого неправильные шрифты. Каждую формулу нужно окружить знаками долларов, а тег math можно самому и не добавлять, он будет добавлен автоматически. Подробнее об этом можно прочитать во втором сообщении темы Что такое карантин и что нужно делать, чтобы там оказаться, раздел "Чем окружать формулы". Отредактируйте, пожалуйста, свое сообщение.


По сути вопроса: обычно в методе наименьших квадратов используют не "честное" расстояние от точки до графика, а лишь разность ординат при одинаковом значении аргумента, т.е. расстояние "по вертикали".

В Вашем случае, если идти по этому пути, то достаточно перейти к новым переменным $t_i = \log x_i$ и использовать обычный метод наименьших квадратов для пар $(t_i,y_i)$, т.е. аппроксимировать $y\approx at+b$.

 Профиль  
                  
 
 
Сообщение16.01.2009, 14:56 
Заслуженный участник


11/05/08
32166
Случай -- сложный. И обычно так задача не ставится -- минимизировать, дескать, кратчайшие расстояния до кривой.

Наверняка имелась в виду стандартная постановка задачи, когда минимизируются расстояния по вертикали.

 Профиль  
                  
 
 
Сообщение16.01.2009, 14:59 
Заслуженный участник
Аватара пользователя


13/08/08
14494
только не до ближайшей касательной, ибо найдется касательная даже проходящая через точку(лежащую выше кривой). Тогда уж по нормали.
Хотя если рассматривается достаточно пологий участок логарифмической фукнкции, может быть подойдут квадраты обычных отклонений?

 Профиль  
                  
 
 
Сообщение16.01.2009, 16:30 


16/01/09
14
Мне кажется что обычный метод, когда расстояние считается вдоль оси ординат, не совсем точен, так как функция скорей всего не пологая, поэтому я решил использовать более точное расстояние.

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

 Профиль  
                  
 
 
Сообщение16.01.2009, 16:35 


29/09/06
4552
PAV в сообщении #177957 писал(а):
По сути вопроса: обычно в методе наименьших квадратов используют не "честное" расстояние от точки до графика, а лишь разность ординат при одинаковом значении аргумента, т.е. расстояние "по вертикали".
Мне кажется, акценты не таковы. Не то, что бы это было "обычно", а так надо, и только так! Использовать "честное" расстояние было бы неверно. Так, если у вас график ток---напряжение или путь---время, то переходя от метров к километрам или от минут к милиамперам вы получите совсем другой результат! Иными словами "расстояние от точки до графика" есть в этом случае величина искусственно-неразумная.

Другое дело и другая задача, когда мы фитируем прямую (окружность, любой геом. объект) по её координатным измерениям и обязаны минимизировать сумму квадратов расстояний до прямой.
Собственно, оно не другое, а то же самое: и там мы смотрим расстояния "по вертикали" от точки $(x_i,y_i)$ до "графика" $z(x,y;\varphi)=x\cos\varphi-y\sin\varphi+L$ или $z(x,y;a,b,R)=(x-a)^2+(y-b)^2-R^2$ (последнее не есть расстояние, но при малых отклонениях пророрционально ему).

 Профиль  
                  
 
 
Сообщение16.01.2009, 16:42 
Заслуженный участник


11/05/08
32166
Я так понял, что задача чисто эмпирическая. Ну попробуйте поиграться с весовой функцией в стандартном МНК. В качестве неё можно взять

$${1\over\sqrt{1+{y'}^2}}={1\over\sqrt{1+{a^2\over x^2}}}\;.$

В качестве начального приближения для $a$ возьмите значение, полученное обычным МНК, а потом сделайте пару итераций с весовым МНК.

 Профиль  
                  
 
 
Сообщение16.01.2009, 16:43 


29/09/06
4552
GoldFinch в сообщении #177998 писал(а):
Есть ли какая-то методика для определения хорошего масштаба по осям?
Нет. Но нередко она легко возникает при конкретизации задачи --- что за точки, откуда берутся, хорошо ли на прямую ложатся, итп. Наверное, надо сначала исследовать чувствительность результатов к этому самому масштабу.

 Профиль  
                  
 
 
Сообщение16.01.2009, 16:44 
Супермодератор
Аватара пользователя


29/07/05
8248
Москва
Алексей К. в сообщении #178001 писал(а):
Мне кажется, акценты не таковы. Не то, что бы это было "обычно", а так надо, и только так!


Согласен

 Профиль  
                  
 
 
Сообщение16.01.2009, 16:45 
Заслуженный участник


11/05/08
32166
Алексей К. в сообщении #178001 писал(а):
Другое дело и другая задача, когда мы фитируем прямую (окружность, любой геом. объект) по её координатным измерениям и обязаны минимизировать сумму квадратов расстояний до прямой.
Собственно, оно не другое, а то же самое:

Оно не то же самое, а совсем другое, и алгоритмы решения, грубо говоря, не имеют ничего общего!

 Профиль  
                  
 
 
Сообщение16.01.2009, 17:06 


29/09/06
4552
GoldFinch, может, что-то из этой беседы Вам сгодится.Там расписана задачка про фитирование прямой по нормали.

Добавлено спустя 15 минут 33 секунды:

ewert в сообщении #178012 писал(а):
Оно не то же самое, а совсем другое
Оно всего лишь нелинейное. Найдите хорошее первое приближение, линеаризуйте формулки ($\varphi=\varphi_0+\Delta\varphi)$, и всё будет как раньше...

 Профиль  
                  
 
 
Сообщение16.01.2009, 17:09 
Заслуженный участник


11/05/08
32166
Если точки аппроксимируются прямой (вообще линейным многообразием в любых размерностях), то никакой линеаризации не надо -- задача допускает точное решение, причём довольно простое.

Ну или не допускает (при определённых условиях количество решений будет бесконечным).

 Профиль  
                  
 
 
Сообщение16.01.2009, 17:22 


29/09/06
4552
Это-то ясно, но тогда я просто не понял Вашего возражения про нечто "совсем другое". Дело, возможно, в неформализованном совсем. В любом случае факт пятницы радует.

 Профиль  
                  
 
 
Сообщение16.01.2009, 17:24 
Заслуженный участник


11/05/08
32166
Какая может быть пятница в сессию?...

 Профиль  
                  
 
 
Сообщение16.01.2009, 18:49 
Заслуженный участник
Аватара пользователя


01/03/06
13626
Москва
ewert в сообщении #178030 писал(а):
Какая может быть пятница в сессию?...
Та, которая обычно перед субботой в сессию бывает. :D

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 26 ]  На страницу 1, 2  След.

Модераторы: Модераторы Математики, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group