On 8 Wrz, 23:15, bartekltg <b...@g...com> wrote:
> On 8 Wrz, 19:35, Mariusz Marszałkowski <m...@g...com> wrote:
> > Obawiam się że nie można tego faktu z istotnym zyskiem
> > wykorzystać, przynajmniej mi nic istotnego nie przychodzi
> > do głowy. Macierz układu równań już nie będzie miała zer.
>
> Jak to nie.
To mówimy o jakimś innym układzie :) W takim razie z
Twoich słów wynika że jest szybsza metoda od tej którą
to zadanie zawsze liczyłem. Na razie nic o niej nie wiem.

Zajrzałem do Cormena, strona 864, widzę jakąś metodę z
macierzą pseudo-odwrotną. Wzór na policzenie macierzy
pseudoodwrotnej: ( ( A^t x A ) ^ -1 ) x A^t - hmmm nie
widzę tutaj przyspieszenia.

Na stronie 864 mam " W praktyce równanie normalne
rozwiązujemy mnożąc y przez A^t, a następnie znajdując
rozkład LU macierzy A^t x A" - wydaje mi się że zawsze
tak to robiłem. Macierz A^t x A nie będzie miała zer, tzn
będzie je miała z bardzo małym prawdopodobieństwem.

Nie rozumiem o czym mówicie.
Pozdrawiam








> Macierz układu rownań
> zminimalizuj norm(Ax-b), czyli macierz A sklada
> sie z tych wektorkow, czyli glownie z jedynek.
Ja to widzę tak, że elementy macierzy A równania Ax=b
w postaci normalnej będą się składały z odpowiednich
sum iloczynów, zsumują się do dużych wartości.

> Dopiero macierz 'rownania normalnego' A'A bedzie pełna.


>
> Poszperaj za metodami dla pierwotnego zagadnienia,
> czyli najmniejszych kwadratow dla rzedkiej macierzy.
> Cos takiego widzialem (ale nie wiem czy dokladnie to
> i czy sie przyda).wiecej popisze jak bede
> miec troche czasu.
>
> pozdrawiam
> bartekltg