|
|
| monchhichi |
Forum-Newbie
|
 |
Beiträge: 3
|
|
|
|
Anmeldedatum: 06.12.08
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 07.12.2008, 15:02
Titel: fminsearch
|
 |
| |
 |
 |
Meine Aufgabe ist es, Parameter mit fminsearch zu optimieren und jeden Integrationsschritt auszugeben.
Mein Programm läuf jetzt...allerdings stürzt mein matlab irgenwann ab- und ich hab keine Ahnung woran das legt:(
Velleicht kann mir ja irgendwer helfen
meine Werte für die Startwerte sind mir bekannt:
"Startwerte=[my_max_1 Y_XS K_S];
%Parameter enthält gesuchten Parameter
options = optimset;
Parameter=fminsearch(@odefun_aufg_fmin,Startwerte,options)
function ausgabe= odefun_aufg_fmin (Startwerte)
%Messwerte
S=[499.5, 442.8, 440.9, 362.1, 305.4, 276.3, 223.2, 148.9, 93.5, 41.8,...
4.6, 0.1, 0.0]; %[mg/L]
X=[38.8, 54.1, 108.6, 142.8, 209.7, 234.7, 287.7, 314.5, 355.7, 394.1,...
374.6, 390.4, 392.6];%[mg/L]
%Integrationsintervall
zeit=[0, 1, 2, 3, 4, 4.5, 5, 5.5, 6, 6.5, 7, 8, 9];%[h]
%Startwert sind X(t=0), S(t=0)
Start=[38.8; 499.5];
p.my_max=Startwerte(1);
p.Y_XS=Startwerte(2);
p.K_S=Startwerte(3);
%ode45
options = odeset;
[t_ode,x]=ode45(@odefun_aufg2,zeit,Start,options,p);
%Zerlegung der Ergebnismatrix:
x_1 = x(:,1);
s_1 = x(:,2);
figure(9);
plot(t_ode,x_1,'r', t_ode,s_1,'g', zeit,X,'r.', zeit,S,'g.')
title('Aufgabe 2b')
xlabel('Zeit in h')
ylabel('Substart bzw. Biomasse [mg/L]')
legend('Biomasse nach Optimierung','Substrat nach Optimierung','Biomasse- Messwerte','Substrat- Messwerte)')
ausgabe = sum(sum([X;S]-x').^2); % Berechnen der Fehlerquadratsumme
return
function[u]= odefun_aufg2 (t_ode,x,p)
%Zustandsvektor wird zerlegt
Biomasse=x(1);
Substrat=x(2);
%Hilfsgroeße berechnen
my=(p.my_max*Substrat)/((p.K_S)+Substrat);
%zeitliche Ableitung bilden
dBiomasse_dt=my*Biomasse;
dSubstrat_dt=-(my*Biomasse)/(p.Y_XS);
%Ausgabe zusammensetzen
u =[dBiomasse_dt dSubstrat_dt]';
"
wie lass ich mir denn alle Integrationsschritte ausgeben...und warum schafft es mein Computer einfach nicht, das Programm richtig auszuführen???
Danke für eure Hilfe!!
|
|
|
|
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
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 - 2024
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.
|
|
