Zuverlässigkeitsberechnung mit Monte Carlo Simulation

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

ponjio

Forum-Anfänger


	Beiträge: 11

	Anmeldedatum: 18.04.11

	Wohnort: ---

	Version: ---

Verfasst am: 22.11.2011, 15:12 Titel: Zuverlässigkeitsberechnung mit Monte Carlo Simulation

Hallo Member,

ich habe mal wieder ein Anliegen an euch, welches mich heute schon den ganzen Tag beschäftigt.
Und zwar möchte ich die Ausfallwahrscheinlichkeit eines elektrischen Systems mit Hilfe einer Monte-Carlo Simulation simulieren. Die Strukturfunktion für das System habe ich erstellt. Nun bin ich der Meinung, dass ich jeder Komponente über die Funktion rand eine Zufallszahl gebe und dann das ganze in einer for-Schleife 100.000x durchlaufen lasse, um einen annehmbaren Wert zu erhalten. So weit so gut. Problem ist jetzt, dass mein Programm die 100.000 Durchgänge immer mit ein und derselben Zufallszahl rechnet. Ich möchte aber 100.000 verschiedene Zufallszahlen zwischen 0 ... 1 und am Ende das ganze als Histogrammbalken ausgeben. Wie wird das realisiert? Ich blicke leider bei der Hilfe nicht mehr durch?! Als Anmerkung: Ich verwende die Version R2009a (7.8 ).

Hier mal meine Code-Schnipsel, die ich mir anhand des Internets und einem Buch versucht habe zusammenzustellen:

Code:

% MCSim = []; % Werte für die Zuverlässigkeit

n=100000;

stream=RandStream('mrg32k3a');
rand(stream,1,n)
% Versuch verschiedene Zufallszahlen zu generieren

x1=rand;
x2=rand;
x3=rand;
x4=rand;
x5=rand;
x6=rand;
x7=rand;

% MCSim = zeros(1, 100000);
for k=1:n % Zahl der Iterationen
X = [x1,x2,x3,x4,x5,x6,x7];
% MCSim = [MCSim; [X,0]];
end

P=((x1*x2*x3*x4*x5-x1*x2*x3*x4*x5*x6)*x7+x1*x2*x3*x4*x5*x6); %Strukturfunktion

% x = min(X):.2:max(X);
% hist(P,x)
% title('Zuverlässigkeitsverteilung')
% versuchte Ausgabe des ganzen

Funktion ohne Link?

Vielen Dank schon mal für die Hilfe.

soad

Forum-Century


	Beiträge: 150

	Anmeldedatum: 10.11.08

	Wohnort: ---

	Version: ---

Verfasst am: 22.11.2011, 15:39 Titel:

Hi,

das Problem ist, dass du die Zufallszahlen einmal generierst und anschließend die Iterationen durchführst. Daher musst du für jede Iteration neue Zufallszahlen generieren:

Code:

for k=1:n % Zahl der Iterationen
X = rand(1,7);
% MCSim = [MCSim; [X,0]];

P= [P; (prod(X(1:5))-prod(X(1:6)))*X(7)+prod(X(1:6))];
end

Funktion ohne Link?

ponjio

Themenstarter

Forum-Anfänger


	Beiträge: 11

	Anmeldedatum: 18.04.11

	Wohnort: ---

	Version: ---

Verfasst am: 22.11.2011, 17:45 Titel:

soad hat Folgendes geschrieben:

Hi,

das Problem ist, dass du die Zufallszahlen einmal generierst und anschließend die Iterationen durchführst. Daher musst du für jede Iteration neue Zufallszahlen generieren

So war es wohl. Vielen Dank für die Hilfe. Das Diagramm konnte ich jetzt auch generieren, nur das Ergebnis scheint noch nicht ganz zu stimmen, aber diesbezüglich werde ich nochmal die Ausgangsfunktion überprüfen.


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.