Re: Алгоритм Левенберга-Маркуарта

From

Victor Kremin ()

To

Alex Astafiev

Date

2002-12-21T22:02:44Z

Area

RU.ALGORITHMS

From: "Victor Kremin" <vkremin@email.lviv.ua>

Есть у меня Numerical Recipes in C (С, PASCAL, FORTRAN)
Могу кинуть архив, 300 кб.
А что интересует по оптимизиции?
Когда-то в дисере активно использовал.

--Пишите на почту.

"Alex Astafiev" <Alex.Astafiev@p16.f228.n5000.z2.fidonet.org> wrote in
message news:1039882119@p16.f228.n5000.z2.fidonet.ftn...
> Народ! Я так ничего толком и не нашел на тему сабжа. В интернете он
попадается либо на фортране, либо в составе платных библиотек
> (Numerical Receipes)
> Пошукайте в своих библиотеках, может быть где-нибудь среди файлов или в
книгах
> найдется?  Ф-ция минимизации LinearFit (Levenberg-Marquardt)
> либо может у кого есть архив исходников от Numerical Receipes?
>
> интересуют также остальные алгоритмы минимизации.
>
>
>
> У меня есть несколько точек .
> Мне необходимо интерполировать эти точки на 4096 других точек,
составляющих
> таблицу.
> Попытаюсь изобразить графически:
>
> |*
> | *
> |   *
> |      * x*   *
> |                *   *x
> |                       *
> |                           * x
> |                               *
> |                                  * x
> |                                     *
> |
> +----------------------------------------
> a                                        b
>
>
> (нарисовалось очень коряво, но суть в том что функция либо линейная либо
> степенная)
> где x - это известные мне пары "аргумент"-"значение" функции.
> a- начало таблицы, b - ее конец.
>
> Сабжевый алгоритм итеративно находит по нескольким парам коэффициенты
> произвольной функции (samplefunction), образец которой ему передается.
Заодно
> он определяет среднеквадратичное отколнение реальных значений от функции.
> (Алгоритм может и не решить предложеную ему задачу).
>
> Отыскав коэффициенты интерполяционной функции, я могу сделать не
> кусочно-линейную интерполяцию функции,  а вызовами samplefunction() с
найдеными
> коэффициентами заполню таблицу. Эта самая samplefunction() может быть и
> кубической сплайн-функцией. В этом случае получится интерполяция сплайном.
>
> В случае более простого кусочно-линейного метода, мне придется еще и
делать
> экстраполяцию ф-ции, ибо заполнить таблицу нужно от начала до конца, а
натурные
> коэффициенты могут быть где-то внутри таблицы, а не на концах, так вот
функцию
> до концов придется еще и экстраполировать.
>
> Я пока что сделал кусочно-линейным методом с экстраполяцией. Но хочу
сделать
> поиск коэффициентов интерполяционной функции.
>
> help.
>
>


--- ifmail v.2.15dev5
 * Origin: Demos online service (2:5020/400)