Matlab - Funktion < Lin. Gleich.-systeme < Numerik < Hochschule < Mathe < Vorhilfe
|
Status: |
(Frage) beantwortet | Datum: | 17:31 So 08.01.2006 | Autor: | alaffm |
Hallo!
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
Kann mir jemand helfen???
1)ich muss erstmal die Funktion beschreiben
function X = TriSol(L,R)
%X = TriSol(L,R)löst das GlSystem LX=R wobei L eine (n,n) nichtsingüläre untere Dreieckmatrix ist;
%die rechte Seite R ist eine (n,r) Matrix;
%X ist die gesuchte matrix.
[n,r] = size(R);
%Die Dimension Matrix R ist n Zeilen und r Spalten
X = zeros(n,r);
%Die Matrix X ist eine Nullmatrix mit n Zeilen und r Salten
for j=i:n-1
X(j,1:r) = R(j,1:r)/L(j,j);
R(j+1:n,j:r) = R(j+1:n,1:r) - L(j+1:n,j)*X(j,1:r);
end
%
X(n,1:r) = R(n,1:r)/L(n,n);
%X und R haben die gleiche Anzahl an Spalten
%r und L haben die gleiche Anzahl an Zeilen
2) Nutze die F. für die Berechnung der Inversen A = (aij)1<= i,j<=6 mit
aij := 0 , falls i < j
1 , falls i = j
-1 , falls i > j.
3) Berechne die Konditionszahl von A bei Nutzung der Zeilensummennorm
Mir fehlen noch ein paar Punkte,um zur Klausur zugelassen zu werden.
Ich wäre der Person,die mir dabei hilft sehr dankbar.
Ich muss die Aufgabe bis morgen 10 uhr schicken.Danke
|
|
|
|
Hallo alaffm,
> 1)ich muss erstmal die Funktion beschreiben
>
> function X = TriSol(L,R)
> %X = TriSol(L,R)löst das GlSystem LX=R wobei L eine (n,n)
> nichtsingüläre untere Dreieckmatrix ist;
> %die rechte Seite R ist eine (n,r) Matrix;
> %X ist die gesuchte matrix.
>
> [n,r] = size(R);
> %Die Dimension Matrix R ist n Zeilen und r Spalten
>
> X = zeros(n,r);
> %Die Matrix X ist eine Nullmatrix mit n Zeilen und r
> Salten
>
> for j=i:n-1
> X(j,1:r) = R(j,1:r)/L(j,j);
> R(j+1:n,j:r) = R(j+1:n,1:r) - L(j+1:n,j)*X(j,1:r);
> end
> %
>
> X(n,1:r) = R(n,1:r)/L(n,n);
Richtig bis auf die rot markierten Stellen da sollte wohl eine 1 stehen.
> %X und R haben die gleiche Anzahl an Spalten
> %r und L haben die gleiche Anzahl an Zeilen
>
> 2) Nutze die F. für die Berechnung der Inversen A =
> (aij)1<= i,j<=6 mit
>
> aij := 0 , falls i < j
> 1 , falls i = j
> -1 , falls i > j.
Die Funktion löst LX=R. Wenn man für R die Einheitsmatrix einsetzt ist X=....
> 3) Berechne die Konditionszahl von A bei Nutzung der
> Zeilensummennorm
[mm] cond(A)=||A||||A^{-1}||
[/mm]
> Mir fehlen noch ein paar Punkte,um zur Klausur zugelassen
> zu werden.
> Ich wäre der Person,die mir dabei hilft sehr dankbar.
> Ich muss die Aufgabe bis morgen 10 uhr schicken.Danke
Naja zu spät vllt. das nächste mal.
gruß
mathemaduenn
|
|
|
|