[P&AM Lab] Linear Least Square -- matrix inverse LU пример программы

[Рябов Евгений]

06.11.2011, 04:33, "Grigoriy A. Sitkarev" <sitkarev на komitex.ru>:
> Приветствую всех!
>
> В эту субботу было много информации поэтому решил сделать пример для
> нахождения коэффициентов методом нахождения обратной матрицы. Это то что
> мы рассмотрели в субботу детально. Я думаю что пример на Си будет очень
> кстати.
>
> Чтобы с ним поиграться, нужно установить библиотеку GSL (GNU Scientific
> Library). В ОС Linux GNU/Debian это пакеты libsgl0-dev и libgsl0ldbl.
>
> Исходник находится на репозитории в каталоге /dsp/lsfit.c.
>
> Т.к. у нас нет настоящих измерений, я заранее придумал коэффиценты (от
> балды) которые нам как раз предстоит отыскать. С помощью этих
> коэффициентов для пар значений [x,y], расположенных на сетке с шагом
> [0.25,0.25] генерируются значения u = U(x,y) и v = V(x,y). Дальше, наша
> задача найти коэффициенты [a0,...,a5] и [b0,...,b5], пользуясь только
> известными значениями координат точек [u,v] и [x,y].
>
> Программа находит коэффициенты и печатает их.
>
> Удачных экспериментов.
>
> --
> Г.А.
>
> _______________________________________________
> Lab mailing list
> Lab на wiki.syktsu.ru
> http://wiki.syktsu.ru/cgi-bin/mailman/listinfo/lab

Есть вопрос по поводу LU факторизации матрицы
gsl_linalg_LU_decomp(&lu.matrix, p, &s); 

это когда матрица распиливается в произведение треугольных матриц. Это lu и s ? 
и  gsl_permutation *p что это? 

typedef struct { size_t size; size_t * data; } gsl_permutation; всё что нашел.

Заранее спасибо

-- 
Заходим на melinux.ru - всё про Linux и не только