Numanal : Matlab
(http://numanal.inf.elte.hu/~krebsz/)
Mindent vektornak, mátrixnak lát
-> ami avval megoldható, azt avval oldja meg
(sorfolytonosan tárolja)
"diag(hova,mit)" -> megadható: milyen vektort a mátrix melyik
mellékátlójába tegye
A(1,:) = A(1,*)
"Epszilon" = "eps" ~ kb. 10 a mínusz6dikon (eltérések
számításáshoz)
"Végtelen" = "Inf"
"Nem-szám" = "NaN" = "Not-a-Number"
"Ciklus 1:0,1:5" = "For i:1 to 5 step 0,1 do"
Ha pontosvesszõt teszük az utasítás végére,
akkor nem írja ki az eredményt
"If feltétel utasítás else utasítás end"
"Save" -> mentés
"Clear" -> memória törlése
"%....%" -> megjegyzés
"diary on/off" -> naplózás ki/be
"disp('szöveg')" -> kiírja a szöveget
Rajzolás: kiteszi a pöttyöket (a vektor-értékeket),
és egyenesekkel köti össze õket
-> az x[i] értékekeknek monoton növekedõnek kell lennie!
-> a grafikon léptékét az elsõ kirajzolt függvény szerint határozza meg
--> "hold on" (ez az állapot), "hold off" (új grafikont kezd)
File-ok (adatfile, függvényfile -> *.m)
Az állománynévnek egyeznie kell a függvény nevével
Hibaszámítás
(a korábbi órai feladatok Matlabban)
Lineáris egyenletrendszer
(a korábbi órai feladatok Matlabban)
Példa Gauss eliminációra:
(ciklus nélkül,, "k" a ciklusváltozó)
"k=7; for j=(k+1):n, A(j,:)=A(j,:)-(A(j,k)/A(k,k))*A(k,:); end; A"
Részleges fõelemkiválasztás
(ciklus nélkül,, "k" a ciklusváltozó)
"[mert,mhely]=max(abs(A(sor(k:n),k)))
sor(k)=mhely+k-1; sor(mhely+k-1)=k;
for j=(k+1):n,
A(sor(j),:)=A(sor(j),:)-(A(sor(j),k)/A(sor(k),k))*A(sor(k),:);
end; A"
LU-felbontás
-> a Matlab általánosabb algoritmust használ:
szükség esetén oszlopot és sort is cserél
2002. november 8.
Szamosközi Péter