Добро пожаловать!

Войдите или зарегистрируйтесь сейчас!

Войти

  1. С 13.10.2020 снова работает авторизация/регистрация через социальные сети: VK (Вконтакте) и Facebook.
    Скрыть объявление

Гауссова нормализация переменных в МНК

Тема в разделе "Геодезия как наука", создана пользователем zvezdochiot, 8 ноя 2020.

  1. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    При составлении матрицы нормальных уравнений:

    https://geodesist.ru/threads/iteracionnyj-algoritm-konformnogo-preobrazovanija.85528/#post-972309

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

    Для преодоления "плохости" матрицы нормальных уравнений "рекомендуют" использовать SVD. Но не всегда и не всем он доступен (в офисных пакетах его нет).

    Но существует такая методика, как "Гауссова нормализация переменных" перед составлением нормальных уравнений:

    https://geodesist.ru/threads/iteracionnyj-algoritm-konformnogo-preobrazovanija.85528/#post-972330

    - позволяющая получать "неплохие" матрицы нормальных уравнений на тех же данных. Что в свою очередь, позволяет применять стандартные методы обращения матриц (включая офисные пакеты).

    Слабая "освещённость" этой методики (не видел ни одной толковой статьи по этой теме) считаю диверсией. Если у кого-нибудь имеется соответствующий материал, будьте добры, окажите любезность, скиньте ссылочку.

    Оффтоп
    Модераторы! Поправьте название, пожалуйста. Опечатался. Cast @Qvinto
     
    #1
    Qvinto нравится это.
  2. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    493
    Симпатии:
    421
    Адрес:
    Любимый Иркутск - середина земли.
    Несколько ремарок.
    1. Идеальным средством ухода от плохой обусловленности матрицы нормальных уравнений параметрического способа уравнивания является использование коррелатноой формы МНК. Правда, надо понимать о каких измерениях и сетях идет речь. Вы ничего по этому поводу не сказали.
    2. То , что вы называете конформным преобразованием по сути является решением задачи определения параметров преобразования координат двух систем координат. Если системы координат обычные для отечественной геодезии и топографии, то они по определению конформные, т.к. построены от ГСК, т.е в логике плоских прямоугольных координат в проекции Гаусса-Крюгера. Для перехода необходимо знать параметры плоского разворота и относительного смещения центра систем плюс относительное изменение масштаба, обусловленное различием в положении осевого меридиана рассматриваемых СК. Поэтому искать афиноры для таких систем - излишнее и совершенно никчемное усложнение.
    3. Для решения задачи в контексте п.2 хорошим средством улучшения обусловленности матрицы нормальных уравнений является использование центрированных координат.
    Успехов!
     
    #2
  3. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    1. Гауссова нормализация не применима к коррелатноой форме МНК, чтоле?
    2. Гауссова нормализация не содержит "центрирования" по мат. ожиданию, чтоле?
    3. Только центрирование даёт на порядки худший результат, чем нормировка (отношение чисел обусловленности), при тех же осложнениях с дальнейшей обработкой (децентрализацией результата).
     
    #3
  4. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    493
    Симпатии:
    421
    Адрес:
    Любимый Иркутск - середина земли.
    Ну почему же не 1, 2 и 3? Я говорю о соразмерности усложнения решения в зависимости от задачи определения параметров перехода. Ваши навороты с прибамбасами могут быть не сообразными с решением этой задачи. Только и всего.
     
    #4
  5. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    Ну да! Это вот верно! Полностью соглашаюсь!

    Но речь не о конкретной задаче (и чего Вы привязались к линейному конформному преобразованию, там нормализация играет ту же роль, что и централизация, матрица один чорт единичная), а о всех видах задач. О ВСЕХ!

    В общем случае, нормализация даёт на порядки (НА ПОРЯДКИ!) менее плохой результат по соотношению чисел обусловленнности, чем простая центрировка, при тех же (ТЕХ ЖЕ!) сложностях в дальнейшей обработке (децентрировке).

    PS: Но можно конечно замахнуться и назвать SVD простым.
     
    #5
  6. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    493
    Симпатии:
    421
    Адрес:
    Любимый Иркутск - середина земли.
    Возможно, я не совсем понял смысл темы, на которую вы сослались. А говорить о ВСЕХ задачах - это вы круто замахнулись со всего плеча. Лично я воздержусь. У вас это диссертация, или ковырялки саморазвития?!
    Матричные вычисления сейчас, при современных вычислительных возможностях и при наличии свободно доступных программных средств, на мой взгляд, вовсе не проблема. В исследовательском плане применения в геодезии - весьма интересно и полезно.
    Кстати, говоря о возможном упрощении, или обходе проблемы обращения плохо обусловленных матриц смею напомнить рекуррентную схему решения задачи МНК. Первоначально формируется нормально обусловленная обратная матрица нормальных уравнений (ковариационная матрица), например по необходимым измерениям, а потом она обновляется по блокам, или даже отдельным измерениям. Задача обращения уходит в тень, но вычислений поболее будет. Но это, как я уже обмолвился, сейчас не проблема.
     
    #6
  7. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    Проблема есть. Усложняет по полной. Даже сведущий может запутаться. Не вижу никаких преимуществ перед нормализацией (ни в плане простоты, ни в плане вычислений).

    PS: Возникает также вопрос "плохого" начального приближения.
     
    #7
  8. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    493
    Симпатии:
    421
    Адрес:
    Любимый Иркутск - середина земли.
    Не буду спорить, поскольку вижу, что вы уже перебрали какой-то набор вариантов решения. Понаблюдаю за вашими успехами. Не забудьте поделиться ими. Один заинтересованный наблюдатель у вас уже есть.
    Мне довелось использовать отмеченные мною варианты регуляризации, а также технику выравнивания значений коэффициентов матрицы параметрических уравнений и стохастическую регуляризацию. Но решение было для конкретных случаев и не ВСЕохватными средствами. Вариантность подхода давала локальный успех при разумных затратах. Универсального метода я не нашел, да и не ставил задачи такого поиска.
     
    #8
  9. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    #9
  10. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    493
    Симпатии:
    421
    Адрес:
    Любимый Иркутск - середина земли.
    Спасибо. Поделился, так поделился гнутыми разработками под линухами. Как-то не очень мне это помогло, интерес отбило и не убедило. Я то настроился на науку, публикации. А тут куча специфического софта. К этой куче комментарии нужны. Но лучше не сейчас.
    Успехов!
     
    #10
    мирось и zvezdochiot нравится это.
  11. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    #11

Поделиться этой страницей

  1. Этот сайт использует файлы cookie. Продолжая пользоваться данным сайтом, Вы соглашаетесь на использование нами Ваших файлов cookie.
    Скрыть объявление