* [room] Библиотеку математическую не подскажите?
@ 2008-08-14 18:07 Alexei V. Mezin
2008-08-15 7:32 ` Vyatcheslav Perevalov
` (2 more replies)
0 siblings, 3 replies; 17+ messages in thread
From: Alexei V. Mezin @ 2008-08-14 18:07 UTC (permalink / raw)
To: smoke-room
Весь Гугль перекопал, ничего толком найти не могу :(
Задача: есть набор точек (x,y,z), то есть поверхность в 3Д. Хочется
скормить эти точки какой-нить функции, которая вычислит оптимальную
плоскость или сферическую поверхность, проходящую через эти точки
(аналог "наименьших квадратов", но для 3Д).
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-14 18:07 [room] Библиотеку математическую не подскажите? Alexei V. Mezin
@ 2008-08-15 7:32 ` Vyatcheslav Perevalov
2008-08-15 10:19 ` Alexei V. Mezin
2008-08-15 12:07 ` Michael Shigorin
2008-08-25 13:44 ` Andrei Lomov
2 siblings, 1 reply; 17+ messages in thread
From: Vyatcheslav Perevalov @ 2008-08-15 7:32 UTC (permalink / raw)
To: Культурный
офтопик
В сообщении от 15 августа 2008 Alexei V. Mezin написал(a):
> Хочется
> скормить эти точки какой-нить функции, которая вычислит оптимальную
> плоскость или сферическую поверхность, проходящую через эти точки
> (аналог "наименьших квадратов", но для 3Д).
Таки наименьшие квадраты и есть, только вычислять сумму квадратов расстояния
от точки до поверхности.
--
Всего хорошего
/vip
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-15 7:32 ` Vyatcheslav Perevalov
@ 2008-08-15 10:19 ` Alexei V. Mezin
2008-08-15 16:18 ` Vyatcheslav Perevalov
0 siblings, 1 reply; 17+ messages in thread
From: Alexei V. Mezin @ 2008-08-15 10:19 UTC (permalink / raw)
To: Культурный
офтопик
Vyatcheslav Perevalov пишет:
> В сообщении от 15 августа 2008 Alexei V. Mezin написал(a):
>> Хочется
>> скормить эти точки какой-нить функции, которая вычислит оптимальную
>> плоскость или сферическую поверхность, проходящую через эти точки
>> (аналог "наименьших квадратов", но для 3Д).
>
> Таки наименьшие квадраты и есть, только вычислять сумму квадратов расстояния
> от точки до поверхности.
>
Ну вот к примеру GNU Scientific Library умеет делать "квадраты" для
набора (x, y, ошибка y), то есть двумерных данных. А надо бы (x, y, z,
ошибка z).
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-15 10:19 ` Alexei V. Mezin
@ 2008-08-15 16:18 ` Vyatcheslav Perevalov
2008-08-15 17:23 ` v.n.belyaev
2008-08-15 17:31 ` Alexei V. Mezin
0 siblings, 2 replies; 17+ messages in thread
From: Vyatcheslav Perevalov @ 2008-08-15 16:18 UTC (permalink / raw)
To: Культурный
офтопик
В сообщении от 15 августа 2008 Alexei V. Mezin написал(a):
> Ну вот к примеру GNU Scientific Library умеет делать "квадраты" для
> набора (x, y, ошибка y), то есть двумерных данных. А надо бы (x, y, z,
> ошибка z).
Я подобные задачи решаю по старинке - беру Справочник по высшей математике
Выгодского в руки, трачу на пролистывание и поиск решения максимум час - и
вуаля.
--
Всего хорошего
/vip
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-15 16:18 ` Vyatcheslav Perevalov
@ 2008-08-15 17:23 ` v.n.belyaev
2008-08-15 18:56 ` Vyatcheslav Perevalov
2008-08-15 17:31 ` Alexei V. Mezin
1 sibling, 1 reply; 17+ messages in thread
From: v.n.belyaev @ 2008-08-15 17:23 UTC (permalink / raw)
To: Культурный
офтопик
Здравствуйте, Vyatcheslav.
Вы писали 15 августа 2008 г., 22:18:27:
VP> В сообщении от 15 августа 2008 Alexei V. Mezin написал(a):
>> Ну вот к примеру GNU Scientific Library умеет делать "квадраты" для
>> набора (x, y, ошибка y), то есть двумерных данных. А надо бы (x, y, z,
>> ошибка z).
VP> Я подобные задачи решаю по старинке - беру Справочник по высшей математике
VP> Выгодского в руки, трачу на пролистывание и поиск решения максимум час - и
VP> вуаля.
Угу, и интегралы по сложной поверхности в уме
берете?
--
С уважением,
v mailto:v.n.belyaev@gmail.com
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-15 16:18 ` Vyatcheslav Perevalov
2008-08-15 17:23 ` v.n.belyaev
@ 2008-08-15 17:31 ` Alexei V. Mezin
2008-08-15 19:04 ` Vyatcheslav Perevalov
1 sibling, 1 reply; 17+ messages in thread
From: Alexei V. Mezin @ 2008-08-15 17:31 UTC (permalink / raw)
To: smoke-room
Vyatcheslav Perevalov пишет:
> В сообщении от 15 августа 2008 Alexei V. Mezin написал(a):
>> Ну вот к примеру GNU Scientific Library умеет делать "квадраты" для
>> набора (x, y, ошибка y), то есть двумерных данных. А надо бы (x, y, z,
>> ошибка z).
>
> Я подобные задачи решаю по старинке - беру Справочник по высшей математике
> Выгодского в руки, трачу на пролистывание и поиск решения максимум час - и
> вуаля.
>
Я на решение такой задачи трачу 10 минут: беру Гугль и вуаля -- вот она
формула в общем виде, и даже с объяснениями. Че там, делов то всего
ничего: просуммировать квадраты разностей для каждой точки, да потом
взять три частных производных, найти их минимумы, и получить коэффициенты.
Только вот надо все это в виде кода (и суммирование нескольких сотен
"квадратов", и численный поиск минимума частной производной, и
желательно чтоб с оптимизацией). В вашем справочнике, случайно, этого
всего не написано?
Причем код этот наверняка написан еще во времена ENIAC, и нету никакого
смысла писать его еще раз.
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-15 17:31 ` Alexei V. Mezin
@ 2008-08-15 19:04 ` Vyatcheslav Perevalov
0 siblings, 0 replies; 17+ messages in thread
From: Vyatcheslav Perevalov @ 2008-08-15 19:04 UTC (permalink / raw)
To: alexei-mezin,
Культурный
офтопик
В сообщении от 16 августа 2008 Alexei V. Mezin написал(a):
> Я на решение такой задачи трачу 10 минут: беру Гугль и вуаля -- вот она
> формула в общем виде, и даже с объяснениями.
Ни разу не встречал. Да и не надо было.
> Че там, делов то всего
> ничего: просуммировать квадраты разностей для каждой точки,
Не квадраты разностей, а квадраты расстояний.
> да потом
> взять три частных производных, найти их минимумы, и получить
> коэффициенты.
>
Минимумы не производных, а самих функций. А если вспомнить, что производная
квадратичной функции есть линейная функция, да и Ваши x,y и z ортогональны
(строго), то всё сводится к системе из четырёх линейных уравнений. Решение
такой системы самостоятельно накодить смогёте? (четыре уравнения - для
случая сферы - x, y, z центра сферы и её радиус)
> Только вот надо все это в виде кода (и суммирование нескольких сотен
> "квадратов", и численный поиск минимума частной производной, и
> желательно чтоб с оптимизацией). В вашем справочнике, случайно, этого
> всего не написано?
>
> Причем код этот наверняка написан еще во времена ENIAC, и нету никакого
> смысла писать его еще раз.
--
Всего хорошего
/vip
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-14 18:07 [room] Библиотеку математическую не подскажите? Alexei V. Mezin
2008-08-15 7:32 ` Vyatcheslav Perevalov
@ 2008-08-15 12:07 ` Michael Shigorin
2008-08-15 14:12 ` Alexei V. Mezin
2008-08-25 13:44 ` Andrei Lomov
2 siblings, 1 reply; 17+ messages in thread
From: Michael Shigorin @ 2008-08-15 12:07 UTC (permalink / raw)
To: smoke-room
On Thu, Aug 14, 2008 at 10:07:00PM +0400, Alexei V. Mezin wrote:
> Весь Гугль перекопал, ничего толком найти не могу :(
Когда-то была хорошая подборка ссылок на sal.kachinatech.com
-- сейчас только через web.archive.org.
--
---- WBR, Michael Shigorin <mike@altlinux.ru>
------ Linux.Kiev http://www.linux.kiev.ua/
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-15 12:07 ` Michael Shigorin
@ 2008-08-15 14:12 ` Alexei V. Mezin
0 siblings, 0 replies; 17+ messages in thread
From: Alexei V. Mezin @ 2008-08-15 14:12 UTC (permalink / raw)
To: Культурный
офтопик
Michael Shigorin пишет:
> On Thu, Aug 14, 2008 at 10:07:00PM +0400, Alexei V. Mezin wrote:
>> Весь Гугль перекопал, ничего толком найти не могу :(
>
> Когда-то была хорошая подборка ссылок на sal.kachinatech.com
> -- сейчас только через web.archive.org.
>
SAL (scientific applications on linux) вполне и без веб-архива работает.
Копался там, но ничего не нашел :( Видимо плохо или не так ищу, не
математик совсем.
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-14 18:07 [room] Библиотеку математическую не подскажите? Alexei V. Mezin
2008-08-15 7:32 ` Vyatcheslav Perevalov
2008-08-15 12:07 ` Michael Shigorin
@ 2008-08-25 13:44 ` Andrei Lomov
2008-08-25 13:45 ` Alexei V. Mezin
2 siblings, 1 reply; 17+ messages in thread
From: Andrei Lomov @ 2008-08-25 13:44 UTC (permalink / raw)
To: smoke-room
Alexei V. Mezin wrote:
> Весь Гугль перекопал, ничего толком найти не могу :(
>
> Задача: есть набор точек (x,y,z), то есть поверхность в 3Д. Хочется
> скормить эти точки какой-нить функции, которая вычислит оптимальную
> плоскость или сферическую поверхность, проходящую через эти точки
> (аналог "наименьших квадратов", но для 3Д).
Я бы сделал в GNU Octave.
Посмотрите в wiki что-нибудь про МНК
(least squares method)
Там размерность -- 2, 3, ... -- не играет роли,
формулки матричные несложные (одно обращение).
Если чем-то конкретно надо помочь,
напишите в личку.
lomov собачка math.nsc.ru
--
Всего доброго,
А.Л.
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-25 13:44 ` Andrei Lomov
@ 2008-08-25 13:45 ` Alexei V. Mezin
2008-08-25 15:10 ` Andrei Lomov
0 siblings, 1 reply; 17+ messages in thread
From: Alexei V. Mezin @ 2008-08-25 13:45 UTC (permalink / raw)
To: Культурный
офтопик
Andrei Lomov пишет:
> Я бы сделал в GNU Octave.
Мне бы интегрировать математику в самописный софт. Octave можно
подключать в виде библиотеки?
Пока позаимствовал часть кода из Gwyddion и смотрю в сторону GNU
Scientific Library (плагин ndlinear).
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] Библиотеку математическую не подскажите?
2008-08-25 13:45 ` Alexei V. Mezin
@ 2008-08-25 15:10 ` Andrei Lomov
2008-08-26 11:46 ` [room] SL Michael Shigorin
0 siblings, 1 reply; 17+ messages in thread
From: Andrei Lomov @ 2008-08-25 15:10 UTC (permalink / raw)
To: smoke-room
Alexei V. Mezin wrote:
> Andrei Lomov пишет:
>
>> Я бы сделал в GNU Octave.
>
> Мне бы интегрировать математику в самописный софт. Octave можно
> подключать в виде библиотеки?
Это скорее вещь в себе (imho).
Насчет библиотек -- написано на разных языках много всего,
в частности, по матричным вычислениям,
но где и как это можно найти в удобной упаковке ...
На си и фортране есть разные "линпаки",
но эти гусли сами не поют.
Может, имеет смысл посмотреть на scientific linux
на основе red hat
--
Всего доброго,
А.Л.
^ permalink raw reply [flat|nested] 17+ messages in thread
* [room] SL
2008-08-25 15:10 ` Andrei Lomov
@ 2008-08-26 11:46 ` Michael Shigorin
2008-08-26 18:02 ` Alexei V. Mezin
2008-08-28 11:47 ` Andrei Lomov
0 siblings, 2 replies; 17+ messages in thread
From: Michael Shigorin @ 2008-08-26 11:46 UTC (permalink / raw)
To: smoke-room
On Mon, Aug 25, 2008 at 10:10:47PM +0700, Andrei Lomov wrote:
> Может, имеет смысл посмотреть на scientific linux
> на основе red hat
Гадость редкая. Впрочем, это было бы несложно предсказать
любому, кто делал дистрибутивы и плотно работал с Настоящими
Учёными(TM).
--
---- WBR, Michael Shigorin <mike@altlinux.ru>
------ Linux.Kiev http://www.linux.kiev.ua/
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] SL
2008-08-26 11:46 ` [room] SL Michael Shigorin
@ 2008-08-26 18:02 ` Alexei V. Mezin
2008-09-02 21:28 ` Michael Shigorin
2008-08-28 11:47 ` Andrei Lomov
1 sibling, 1 reply; 17+ messages in thread
From: Alexei V. Mezin @ 2008-08-26 18:02 UTC (permalink / raw)
To: Культурный
офтопик
Michael Shigorin пишет:
> On Mon, Aug 25, 2008 at 10:10:47PM +0700, Andrei Lomov wrote:
>> Может, имеет смысл посмотреть на scientific linux
>> на основе red hat
>
> Гадость редкая. Впрочем, это было бы несложно предсказать
> любому, кто делал дистрибутивы и плотно работал с Настоящими
> Учёными(TM).
>
А можно чуть развить тему? Интересно для общего развития. Чем так плохи
ученые? :) И зачем им какой-то отдельный дистрибутив? Почему нельзя
просто собрать под любой мэйнстрим-дистриб нужные пакеты?
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] SL
2008-08-26 18:02 ` Alexei V. Mezin
@ 2008-09-02 21:28 ` Michael Shigorin
0 siblings, 0 replies; 17+ messages in thread
From: Michael Shigorin @ 2008-09-02 21:28 UTC (permalink / raw)
To: Культурный
офтопик
On Tue, Aug 26, 2008 at 10:02:46PM +0400, Alexei V. Mezin wrote:
> >>Может, имеет смысл посмотреть на scientific linux
> >>на основе red hat
> >Гадость редкая. Впрочем, это было бы несложно предсказать
> >любому, кто делал дистрибутивы и плотно работал с Настоящими
> >Учёными(TM).
> А можно чуть развить тему? Интересно для общего развития.
> Чем так плохи ученые? :)
Не плохи, а специфичны. Как, скажем, настоящие военные(ТМ).
> И зачем им какой-то отдельный дистрибутив?
Кому для чего... линейной алгебры из коробки надо, например.
> Почему нельзя просто собрать под любой мэйнстрим-дистриб нужные
> пакеты?
Вот из таких addons обычно и возникают производные дистры.
--
---- WBR, Michael Shigorin <mike@altlinux.ru>
------ Linux.Kiev http://www.linux.kiev.ua/
^ permalink raw reply [flat|nested] 17+ messages in thread
* Re: [room] SL
2008-08-26 11:46 ` [room] SL Michael Shigorin
2008-08-26 18:02 ` Alexei V. Mezin
@ 2008-08-28 11:47 ` Andrei Lomov
1 sibling, 0 replies; 17+ messages in thread
From: Andrei Lomov @ 2008-08-28 11:47 UTC (permalink / raw)
To: smoke-room
Michael Shigorin wrote:
> On Mon, Aug 25, 2008 at 10:10:47PM +0700, Andrei Lomov wrote:
>> Может, имеет смысл посмотреть на scientific linux
>> на основе red hat
>
> Гадость редкая. Впрочем, это было бы несложно предсказать
> любому, кто делал дистрибутивы и плотно работал с Настоящими
> Учёными(TM).
Вопрос именно о списке пакетов,
а не о работе с "Н. У." ;)
--
Всего доброго,
А.Л.
^ permalink raw reply [flat|nested] 17+ messages in thread
end of thread, other threads:[~2008-09-02 21:28 UTC | newest]
Thread overview: 17+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2008-08-14 18:07 [room] Библиотеку математическую не подскажите? Alexei V. Mezin
2008-08-15 7:32 ` Vyatcheslav Perevalov
2008-08-15 10:19 ` Alexei V. Mezin
2008-08-15 16:18 ` Vyatcheslav Perevalov
2008-08-15 17:23 ` v.n.belyaev
2008-08-15 18:56 ` Vyatcheslav Perevalov
2008-08-15 17:31 ` Alexei V. Mezin
2008-08-15 19:04 ` Vyatcheslav Perevalov
2008-08-15 12:07 ` Michael Shigorin
2008-08-15 14:12 ` Alexei V. Mezin
2008-08-25 13:44 ` Andrei Lomov
2008-08-25 13:45 ` Alexei V. Mezin
2008-08-25 15:10 ` Andrei Lomov
2008-08-26 11:46 ` [room] SL Michael Shigorin
2008-08-26 18:02 ` Alexei V. Mezin
2008-09-02 21:28 ` Michael Shigorin
2008-08-28 11:47 ` Andrei Lomov
Культурный офтопик
This inbox may be cloned and mirrored by anyone:
git clone --mirror http://lore.altlinux.org/smoke-room/0 smoke-room/git/0.git
# If you have public-inbox 1.1+ installed, you may
# initialize and index your mirror using the following commands:
public-inbox-init -V2 smoke-room smoke-room/ http://lore.altlinux.org/smoke-room \
smoke-room@lists.altlinux.org smoke-room@lists.altlinux.ru smoke-room@lists.altlinux.com smoke-room@altlinux.ru smoke-room@altlinux.org smoke-room@altlinux.com
public-inbox-index smoke-room
Example config snippet for mirrors.
Newsgroup available over NNTP:
nntp://lore.altlinux.org/org.altlinux.lists.smoke-room
AGPL code for this site: git clone https://public-inbox.org/public-inbox.git