Отправлено Max_Fly 30 июня 2004 г. 11:31

если надо - напишу полные выкладки - в кратце суть следующая. имеем рассеяние точек. ТОЧНО знаю что это рассеяние суть окружность. надо определить параметры этой окружности. по точкам делаю окружность методом наименьших квадратов. проверяю на данных с ошибками и чётком математическом круге (нихрена нет там ошибки). на некоторых радиусах получается всё как в копейку - на некоторых радиус нуль :(((. я лох или где-то тут грабли? При чём сдвигаешь радиус от критического хоть на единичку - всё в норме. грубо говоря радиус 20 точек - нихрена не берёт. 19 и 21 - нет вопросов. я просто в замешательстве.

• Задача мне тоже актуальная. Есть такие мысли. — AlexandrY   (30.06.2004 14:15, пустое, картинка)
  • гм... (+) — Max_Fly   (30.06.2004 14:28, 351 байт)
    • x,y,r здесь неизвестные параметры окружности, а не переменные в уравнении окружности — AlexandrY   (30.06.2004 14:36, 151 байт)
      • понял :). сорри. — Max_Fly   (30.06.2004 14:38, пустое)
  • Поясню. — AlexandrY   (30.06.2004 14:24, 430 байт)
    • в чистом виде нельзя (+) — Max_Fly   (30.06.2004 14:31, 352 байт)
      • Хм.... (+) — SM   (30.06.2004 14:37, 315 байт)
        • вот и думаю я о том что я таки попробую построить (+) — Max_Fly   (30.06.2004 14:42, 80 байт)
          • Врядли вы откроете америку. Вот еще можете решать систему нелинейных уравнений для точного ответа. — AlexandrY   (30.06.2004 17:37, пустое, картинка)
            • попробуйте прогнать такие точки — Max_Fly   (30.06.2004 18:42, 330 байт)
              • Нет проблем. Минимизатор решил, а вот система нелинейных не решилась. — AlexandrY   (30.06.2004 19:22, пустое, картинка)
                • То есть нет, извиняюсь. Нелинейные тоже решились, но при хорошем начальном приближении, они дали ответ 250,250,20 — AlexandrY   (30.06.2004 19:27, пустое)
                  • Причем с радиусом можете ошибаться сильно, но с координатами не более 20% — AlexandrY   (30.06.2004 19:31, пустое)
• А что является нормой ошибки? — GF(6)   (30.06.2004 13:41, 144 байт)
  • работаешь по 2-м изображениям? — Max_Fly   (30.06.2004 14:11, 80 байт)
    • Нет речь в моем случае идет не об изображении — GF(6)   (30.06.2004 14:17, 54 байт)
      • а.. в этом смысле :) — Max_Fly   (30.06.2004 14:33, 72 байт)
  • не совсем понял слово норма. — Max_Fly   (30.06.2004 13:51, 191 байт)
    • т. е. критерий ошибки какой? сумма квадратов расстояний от точек до искомой окружности? — GF(6)   (30.06.2004 13:57, пустое)
      • да. — Max_Fly   (30.06.2004 14:01, пустое)
        • :) или как в старых учебниках их называют "невязок" — Max_Fly   (30.06.2004 14:02, пустое)
          • Просто эта проблема меня очень интересует, как я выше сказал. Скопировал все сообщения - хочу на досуге подумать. А щас, увы, работать надо :( — GF(6)   (30.06.2004 14:12, пустое)
    • Вроде как норма ошибки это (+) — SM   (30.06.2004 13:55, 207 байт)
• Ну математик я в общем-то никакой:)... — POV   (30.06.2004 11:37, 170 байт)
  • может... отошлите меня в какой-нибудь тривиальный учебник по арифметике - должен жу быть какой-то признак что решение тривиально! — Max_Fly   (30.06.2004 11:59, пустое)
    • Вряд ли это тривиальная арифметика... — POV   (30.06.2004 12:47, 691 байт)
      • Ответ: Именно - элементарная арифметика... — Malya   (30.06.2004 13:28, 546 байт)
        • :) ну а если должно получиться круг и апроксимируем мы его кругом? ошибки не будет нигде и никогда? — Max_Fly   (30.06.2004 13:53, пустое)
      • :) радиолокацией меня тоже кормили как студента. (равно как и радионавигацией и т.п. :) ) — Max_Fly   (30.06.2004 12:51, 30 байт)
• Дык. Глюк наверное в программе. — SM   (30.06.2004 11:35, пустое)
  • этого я и ожидал :). сам так думал. привожу полные выкладки (+) — Max_Fly   (30.06.2004 11:50, 1193 байт)
    • Неправильные у тебя уравнения дядя Фё... Max_Fly — GF(6)   (01.07.2004 08:55, 475 байт)
    • По моему не совсем правильно. — -=ВН=-   (30.06.2004 15:30, 579 байт)
    • Если я чего-то понимаю, то немного не так все (+) — SM   (30.06.2004 12:13, 683 байт)
      • я так пытался. (+) — Max_Fly   (30.06.2004 12:28, 458 байт)
        • Ответ: допустим три параметра — -Tумблер-   (30.06.2004 14:35, 498 байт)
          • что-то подобное я делал (+) — Max_Fly   (30.06.2004 15:33, 430 байт)
            • Ответ: это означает, что — -Tумблер-   (30.06.2004 15:44, 461 байт)
              • ваши слова тоже содержат истину... действительно когда-то я вычислял геометрический центр изображения... надо вспомнить буйную юность. — Max_Fly   (30.06.2004 18:57, пустое)
        • ну и что, что три параметра (+) — SM   (30.06.2004 12:55, 358 байт)
          • Немного не так (+) — SM   (30.06.2004 13:19, 345 байт)
            • мдя. :) а задач для DSP становится всё больше и больше. чует моё сердце - скоро надо будет 67-ю серию брать :) — Max_Fly   (30.06.2004 13:22, пустое)
              • Это скорее не для DSP задача, а для mapple. Ну в крайнем случае для ручки с бумажкой (и мозгов) — SM   (30.06.2004 13:27, пустое)
                • ну первый раз - да. а потом мне надо раз 5-7 в секунду считать. :). И это только одна из задач. — Max_Fly   (30.06.2004 13:29, пустое)
                  • Так когда оно решено аналитически (+) — SM   (30.06.2004 13:33, 230 байт)
                    • ситуация такая (+) — Max_Fly   (30.06.2004 13:43, 735 байт)
                      • Кстати - а вот если (+) — SM   (30.06.2004 14:29, 398 байт)
                        • завтра я на выезд в цех - будет время подумать. — Max_Fly   (30.06.2004 14:36, 1 байт)
  • А что собственно происходит-то? (+) — SM   (30.06.2004 11:37, 75 байт)
    • да всё проще — Max_Fly   (30.06.2004 12:06, 346 байт)
    • Или Вы сами аналитически (+) — SM   (30.06.2004 11:44, 90 байт)
      • сам аналитически. :(. — Max_Fly   (30.06.2004 11:53, 171 байт)
        • Бывают ситуации, когда решение неустойчиво - сильно зависит от точности решения. — Крок   (30.06.2004 12:33, 68 байт)