Verfasst am: 14.04.2015, 17:33
Titel: Mehrere Graphen in einem plot
Moin!
Bei dem Code erstelle ich ein Chloridprofil, wobei ich e immer ändern kann.
Wie bekomme ich das hin, dass ich mir e=0,2 ausgeben lasse und dann meinetwegen e=1 berechne und mir das in einer Graphik anzeigen lasse. Also die zwei graphen in einm plot.
Und wie lege ich eine Legende an mit den zwei Graphen?
Vielen Dank für die Antwort!
Code:
clearall
%% Implementieren der Paramter
D1=10^-8; % Diffusionskoeffizient [cm^2]
D2=10^-9; % Diffusionskoeffizient [cm^2]
dt=86400; %[s]
dx=0.05; %[cm]
n=400; %Anzahl der Spalten(Betrachtungspunkte)
laengeges=dx*n; %Gesamtlänge [cm]
x=0:dx:laengeges-15
e=0.2; %Größe der Randzone [cm]
k=e/dx %Punkte bis Schnittstelle
Ne=(D1*dt)/dx^2; %Neumann-Zahl für Randzone [-]
Ke=(D2*dt)/dx^2; %Ke-Zahl für Betonkern
abst=10*365*24*60*60;
steps=round(abst/dt);
%% Iteration bis zum Punkt 4
C=zeros(steps,n+1); % Array
C(:,1)=0.5% 1. Spalte = 0.5 -> Randbedingung
%C(:,n+1)=1%
for L=1:steps-1 for i=2:k
C(L+1,i)=Ne*C(L,i-1)+C(L,i)-2*Ne*C(L,i)+Ne*C(L,i+1);
end
for i=k+2:n
C(L+1,i)=Ke*C(L,i-1)+C(L,i)-2*Ke*C(L,i)+Ke*C(L,i+1);
end
% C(L+1,n+1)= 0.5*C(L+1,n); end
%% Plotten von C plot(x,C(steps,1:101),'+','color','blue');
hold on
plot(x,C(steps,1:101),'color','blue');
holf off
grid on
einfach ne forschleife über e laufen lassen.
das was du mit 2 plot befehlen machst kann man auch mit einem machen.
das clear all würde ich weglassen. warum wurde hier schon ausgiebig im forum besprochen
Code:
for e=[0.2,1]
%% Implementieren der Paramter
D1=10^-8; % Diffusionskoeffizient [cm^2]
D2=10^-9; % Diffusionskoeffizient [cm^2]
dt=86400; %[s]
dx=0.05; %[cm]
n=400; %Anzahl der Spalten(Betrachtungspunkte)
laengeges=dx*n; %Gesamtlänge [cm]
x=0:dx:laengeges-15;
%Größe der Randzone [cm]
k=e/dx; %Punkte bis Schnittstelle
Ne=(D1*dt)/dx^2; %Neumann-Zahl für Randzone [-]
Ke=(D2*dt)/dx^2; %Ke-Zahl für Betonkern
abst=10*365*24*60*60;
steps=round(abst/dt);
%% Iteration bis zum Punkt 4
C=zeros(steps,n+1); % Array
C(:,1)=0.5; % 1. Spalte = 0.5 -> Randbedingung
%C(:,n+1)=1%
for L=1:steps-1 for i=2:k
C(L+1,i)=Ne*C(L,i-1)+C(L,i)-2*Ne*C(L,i)+Ne*C(L,i+1);
end
die color angabe kann auch in einem RGB vektor erfolgen. die kann man natürlich vorher festlegen und dann durchlaufen lassen wenn man das nicht automatisch machen will.
_________________
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
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.