[Octave] Iloczyn macierzy trójkątnej i mającej 1 kolumnę

3mts3 · Post autor: **3mts3** » 5 gru 2011, o 02:22

Witam potrzebuje pomocy z pewnym zadaniem, z góry dziękuje za pomoc.

Napisać funkcję obliczającą iloczyn macierzy trójkątnej i 1 kolumnowej ub=c
Do obliczeń przyjąć:

Kod: Zaznacz cały

u=[5,2,3,1,3;
    0,2,1,4,5;
    0,0,1,3,0;
    0,0,0,6,7;
    0,0,0,0,10]

b=[5;
    2;
    3;
    1;
    5]

Algorytm dostosować do postaci macierzy.
Chodzi o to żeby nie mnożyć przez elementy zerowe.

Xitami · Post autor: **Xitami** » 5 gru 2011, o 07:49

to mój pierwszy kawałeczek kodu Octavowego więc może nie wolny od błędów
posłużyłem się metodą J. F. Champollion'a

Kod: Zaznacz cały

http://rosettacode.org/wiki/Matrix_multiplication#Octave

Kod: Zaznacz cały

u=[5,2,3,1,3;0,2,1,4,5;0,0,1,3,0;0,0,0,6,7;0,0,0,0,10];
b=[5;2;3;1;5];
c=[0;0;0;0;0];
for i=0:4
	for k=i:4
		c(i) = c(i) + u(i,k)*b(k);
	endfor
endfor

3mts3 · Post autor: **3mts3** » 5 gru 2011, o 08:27

trochę nie tak to miało być, wynik wychodzi raczej inny.

Dzięki że spróbowałeś

Xitami · Post autor: **Xitami** » 5 gru 2011, o 08:36

mam pomysła, bo Octavy nie mam

Kod: Zaznacz cały

for i=0:4
   s=0;    
   for k=i:4
      s = s + u(i,k)*b(k);
   endfor
   c(i)=s;
endfor

-----edit---------------------------------
e... guzik, przecie to to samo
a wynik poprawny[54;36;6;41;50]
wklej jakie jest "c" z tego mojego programiku

3mts3 · Post autor: **3mts3** » 5 gru 2011, o 10:29

tak wynik ma być taki jak napisałeś, a ten program który napisałeś niestety nie chce działać.

Xitami · Post autor: **Xitami** » 5 gru 2011, o 10:47

dodałem kilka średników, sprawdź
275119.htm#p4831382