Systemmatrix Neumann Randbedingung

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

fossifoh

Forum-Newbie


	Beiträge: 1

	Anmeldedatum: 02.01.11

	Wohnort: Wolfenbüttel

	Version: ---

Verfasst am: 02.01.2011, 18:53 Titel: Systemmatrix Neumann Randbedingung

Hallo Leute,
ich habe folgendes Problem. Ich muss eine elliptische DGL in einem beliebigen rechteckigem Gebiet mit Neumann Randbedingungen lösen.
Nur mit meiner programmierten Systemmatrix funktioniert das nicht und ich bekomm unendlich große Werte. Habe das Problem auch schon zum Teil herausbekommen (glaube ich).

Hier erstmal der Programmabschnitt für die Matrix:

Code:

h=1/3; % Festlegung der Schrittweite h
k=1/4; % Festlegung der Schrittweite k
hk=k/h; % k faktoriell in h abgebildet
%---Anfang Systemmatrix Neumann---%
%---Anfang Blockdiagramm---%
d1=diag((2+2/hk^2)*ones(n,1));
r1=-diag(1*ones(n-1,1),1);
l1=-diag(1*ones(n-1,1),-1);
d1=d1+r1+l1;

d1(1,2)=-2;
d1(n,n-1)=-2;
%---Ende Blockdiagramm---%

%---Anfang Seitendiagonale rechts---%
r2=-diag((1/hk^2)*ones(nm-n,1),n);

for kk=1:n
r2(0+kk,n+kk)=-1/hk^2*2;
end
%---Ende Seitendiagonale rechts---%

%---Anfang Seitendiagonale links---%
l2=-diag((1/hk^2)*ones(nm-n,1),-n);

for kk=1:n
l2(nm+1-kk,nm-n+1-kk)=-1/hk^2*2;
end
%---Ende Seitendiagonale links---%

mat=d1;

for kk=1:m-1
mat=blkdiag(mat,d1);
end

mat1=mat+r2+l2;

clear mat;
%---Ende Systemmatrix Neumann---%
disp(mat1)
mat1=sparse(mat1);
%;

Funktion ohne Link?

Die Matrix, die erstellt wird, sieht genauso aus, wie die, die wir in den Vorlesungen auf dem Papier erstellt haben.
Wenn ich die zwei "-2" im Block einen anderen Wert gebe, bekomme ich vernünftige Zahlen.

5.5556 -2.0000 0 0
-1.0000 5.5556 -1.0000 0
0 -1.0000 5.5556 -1.0000
0 0 -2.0000 5.5556

Ich hoffe, es ist nur ein kleiner Fehler, den mir ein gübter Matlab'er/erin beantworten kann.
Erstmal vielen Dank, ich werde morgen wieder vorbeischauen.


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.