On 15 Wrz, 21:42, Mariusz Marszałkowski <m...@g...com> wrote:
> On 15 Wrz, 03:52, bartekltg <b...@g...com> wrote:> On 14 Wrz, 20:31, Mariusz
Marszałkowski <m...@g...com> wrote:
>
> > > Link do przykładu:http://www.przeklej.pl/plik/eq-xls-0020nk3
8s5a5
>
> > Ześ sobie format znalzal;) Ściagnalem dane, makr nie ruszam:)
>
> Makra to tylko copy-paste z komórki do komórki i przeliczenie arkusza.
> Nie sądzę abym miał jakiegoś wirusa który dokleił coś samemu.
>
> > Zanim zaczniesz szukaszukać błędów w stabilnosci  numerycznej
> > i prezycji obliczeń, zerknij na niezalesnosc wektorów;)
>
> > Twoje wektroki(pionowe z macierzy data) _nie_ są liniowo niezalezne!
> > Rank (data) = 13. A wymiar 16. Masz trzy stopnie swobody w zapisaniu
> > wyniku:)
>
> Danych jest (było gdy wysyłałem) 10tys.

Nie zrozumiales. Macierz jest 10 000 na 16. Ale jej
rząd to 13. Maceirz ejst osobliwa, liniowo zalezna.
To stad w pierwszej kolejnosci biorą sie zypelnie
rozne wektroy dajace ten sam wynkik.

BTW, kopnales sie w wyliczaniach:

> Co z czasem obliczeń?
> Załóżmy że mamy 10^9 danych i trzeba znaleźć 10^6
> parametrów. Daje to 10^9 * 10^6 * 3 ~=~ 10^16 operacji.
> Czyli tak czy siak odpada.

Mamy M(10^9) wektorkow dlugośći N=(M/1000), kazdy ma (k=4)
niezerowych skladowych. Daje to S=4*10^9 liczb (satrczy dwa bajty
całkowitoliczbowej;) w data i 10^9 na wektor wyników(tu juz float/
double).

Co prawda Wykonanie mnozenia A*x zajmuje rzędu S operacji,
ale nie trzeba tego robić co iteracje! A interesujace nas (A*e_i, b)
da sie policzyć w zlozonosci k*(M/N) (*powiedzmy 8 operacji).

Moze poznym wieczorkeim wrzuce kod (od razu ostrzegam,
z lenistwa matlabowski;)

pozdrawiam
bartekltg