Witaj Gościu! ( Zaloguj | Rejestruj )

Forum PHP.pl

 
 Reply to this topicStart new topic
> Rozkład Macierzy LU - Metoda Dollitle'a
kudlaty12
post
Post #1





Grupa: Zarejestrowani
Postów: 2
Pomógł: 0
Dołączył: 24.03.2010

Ostrzeżenie: (0%)
-----

Witajcie.
Mam pewien problem, ponieważ napisałem skrypt do liczenia macierzy (rozkład LU) metodą Dollitle'a. Jest on oparty na przykładzie zawartym na wikipedia.pl , jednak nie dziala do końca tak jak powinien. Wartości bedące 0 zapisują sie w macierzy L ponizej przekątnej i w macierzy U ponad przekątną. Jednak problem polega na tym ze nie ma wypelnienia calych macierzy, gdzyż kilka luk po wykonaniu skryptu zostaje pusta. I pytanie brzmi czemu ?
Czy ktoś mógłbym sprawdzić kod i pomóc mi w odnalezieniu błędu ?
Będę wdzięczny (IMG:style_emoticons/default/winksmiley.jpg)

Ponizej przesyłam link do kodu w php.


http://lew.wsinf.edu.pl/~blazej01/rozkladlu.txt

a o to wersja graficzna jak to wygląda :

http://lew.wsinf.edu.pl/~blazej01/rozkladlu2.php

a tutaj link do materiału na Wikipedia.pl :

http://pl.wikipedia.org/wiki/Metoda_LU

Ten post edytował kudlaty12 24.03.2010, 15:33:00
Go to the top of the page
+Quote Post
rzymek01
post
Post #2





Grupa: Zarejestrowani
Postów: 592
Pomógł: 62
Dołączył: 3.08.2006

Ostrzeżenie: (0%)
-----

Witaj,
dzięki za linka do LU, bo nie znałem wcześniej tej metody :]

odnośnie kodu - zawsze chciłem to powiedzieć - u mnie działa (IMG:style_emoticons/default/biggrin.gif)

dobra, było parę błędów w kodzie, ale jakoś poszło:
1. inicjalizuj zmienne, np. $suma = 0;
2. wg mnie sprawdzenie w liczeniu kolumn u[i][i]!=0 mija się z celem, bo wtedy cała metoda idzie sie j*** na łąkę i trzeba zastosować - z tego co tam jest napisane - metodę Crouta
3. zrób sobie z tego klase, zeby nie trzeba bylo przekazywac co chwilę macierzy, i,j,n (IMG:style_emoticons/default/smile.gif)
4. nie kamufluj błędów @-małpą, w ifie $i =1 (przypisanie!), akurat ten warunek był cały zbędny
5. natomiast główne błąd to logiczny w iteracji pętli, np. sumuje się włącznie do k-1, oraz w kolumnach bierzemy za $j o jeden więcej niż w wierszach

kod:
[PHP] pobierz, plaintext 
  1.     <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN"
  2. "http://www.w3.org/TR/html4/strict.dtd">
  3.  
  4. <html>
  5.  
  6. <head>
  7.  
  8. <title>Mnozenie Macierzy</title>
  9. <meta http-equiv="Content-type" content="text/html; charset=utf-8">
  10.  
  11. <style type="text/css">
  12. 	.pro{
  13. 		text-align: center;
  14. 	}
  15. 	.proo{
  16. 		margin-left: auto;
  17. 		margin-right: auto;
  18. 		width: 200px
  19.  
  20. 	}
  21. 	body { 
  22. 	background-color:#ff9900; 
  23.  
  24. 	}
  25.  
  26. </style>
  27. </head>
  28.  
  29.  
  30. 	 <?php       
  31. error_reporting(E_ALL);
  32.  
  33.       //Budujemy 3 tablice, A = L * U (przyklad ponizej to przyklad dla macierzy 3x3)
  34.  
  35.  
  36.  
  37. 	  //macierz A:	  
  38.  
  39.       $A = array( 'a11' => '5', 'a12' => '3', 'a13' => '2',
  40.  
  41.                   'a21' => '1', 'a22' => '2', 'a23' => '0',
  42.  
  43.                   'a31' => '3', 'a32' => '0', 'a33' => '4',   
  44.       );
  45.  
  46.       //kolumny macierzy:
  47.  
  48.       $L = array( 'l11' => '1', 'l12' => '0', 'l13' => '0',
  49.  
  50.                   'l21' => '' , 'l22' => '1', 'l23' => '0',
  51.  
  52.                   'l31' => '' , 'l32' => '' , 'l33' => '1',
  53.       );
  54.  
  55. 	  //wiersze macierzy:
  56.  
  57.       $U = array( 'u11' => '' , 'u12' => '' , 'u13' => '' ,
  58.  
  59.                   'u21' => '0', 'u22' => '' , 'u23' => '' ,
  60.  
  61.                   'u31' => '0', 'u32' => '0', 'u33' => '' ,
  62.       );
  63.  
  64.       //Aby wynik byl poprawny musimy obliczac na przemian wiersze macierzy U i kolumny macierzy L            
  65.  
  66.  
  67. 	  //funkcja liczaca $U czyli wiersze macierzy 
  68.  
  69.       function liczWiersz(&$A, &$L, &$U, $i, $j) {
  70.  
  71. 		$suma = 0;
  72.  
  73. 		for($k=1;$k<=$i-1;$k++)
  74. 			$suma += ($L['l'.$i.$k] * $U['u'.$k.$j]);
  75.  
  76. 		$U['u'.$i.$j] = round($A['a'.$i.$j] - $suma, 3);
  77.       }
  78.  
  79.       // funkcja liczaca $L czyli kolumny macierzy
  80.       function liczKolumne(&$A, &$L, &$U, $i, $j) {
  81. 		$suma = 0;
  82.  
  83.               if($U['u'.$i.$i] != 0 ) {
  84.  
  85.                       for($k=1;$k<=$i-1;$k++)
  86.  
  87.                               {
  88.  
  89.                                       $suma += ($L['l'.$j.$k] * $U['u'.$k.$i]);
  90.  
  91.                               }
  92.  
  93.               $L['l'.$j.$i] = round( (1 / $U['u'.$i.$i]) * ($A['a'.$j.$i] - $suma), 3);
  94.  
  95.               }
  96.  
  97.       }
  98.  
  99.  
  100.  
  101.       for($i=1;$i<=3;$i++) {
  102.  
  103.               for($j=$i;$j<=3;$j++) {
  104.  
  105. 			liczWiersz($A, $L, $U, $i, $j);
  106.  
  107.               }
  108.  
  109.               for($j=$i+1;$j<=3;$j++) {
  110.  
  111. 			liczKolumne($A, $L, $U, $i, $j);  
  112.               }
  113.  
  114.       }
  115.  
  116.     //funkcja wyswietlajaca macierz 
  117. 	// $M - dana macierz jaka chcemy wyswietlic
  118. 	// $p - przedrostek kolumny
  119.  
  120.       function wyswietlMacierz(&$M, $p) {
  121.  
  122.     // strtoupper - zwraca ciag znakow z wszystkimi literami alfabetu przekonwertowanymi na duże.
  123.  
  124.               echo "Macierz ".strtoupper($p);
  125.  
  126.               echo "<table style=\"border: 1px solid black; width: 150px; text-align: center;\">";
  127.  
  128.               for($i=1;$i<=3;$i++) {
  129.  
  130.                       echo "<tr>";
  131.  
  132.                       echo "<td>".$M[$p.$i.'1']."</td><td>".$M[$p.$i.'2']."</td><td>".$M[$p.$i.'3']."</td>";
  133.  
  134.                       echo "</tr>";
  135.  
  136.               }
  137.  
  138.               echo "</table>";
  139.  
  140.       }
  141.  
  142.  
  143.  
  144.       wyswietlMacierz($A, 'a');
  145.  
  146.       wyswietlMacierz($L, 'l');
  147.  
  148.       wyswietlMacierz($U, 'u');
  149.  
  150.  
  151.  
  152.  
  153.  
  154.       ?> 
  155. 	  <br>
  156. 	   <br>
  157. 	    <br>
  158. 		 <br>
  159. 		  <br>
  160. 		   <br>
  161.  
  162. 	  <footer>Stworzone na podstawie przykładu z Wikipedia.pl</footer>
[PHP] pobierz, plaintext
Go to the top of the page
+Quote Post
kudlaty12
post
Post #3





Grupa: Zarejestrowani
Postów: 2
Pomógł: 0
Dołączył: 24.03.2010

Ostrzeżenie: (0%)
-----

Wielkie dzieki. Mówiłem ze kilka błędów małych bo tak to reszta była dobrze.

Poprawiłem sobie jeszcze wedlug siebie ale i tak działa (IMG:style_emoticons/default/winksmiley.jpg)
Go to the top of the page
+Quote Post
« Następny starszy · Algorytmy, klasy, funkcje · Następny nowszy »
 

Reply to this topicStart new topic
2 Użytkowników czyta ten temat (2 Gości i 0 Anonimowych użytkowników)
0 Zarejestrowanych:

 

Tryb wyświetlania: Standardowy · Przełącz na: Linearny+ · Przełącz na: Drzewo

Śledź ten temat · Wyślij temat na e-mail · Wydrukuj ten temat · Subskrybuj to forum

RSS Aktualny czas: 24.12.2025 - 05:31
Powered By IP.Board © 2025  IPS, Inc.
All changes by PHP.pl Administrators
Hosting zapewnia NQ.pl hosting, trac, svn