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

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Schreibweise ODE

 

Gast123456

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 06.05.2014, 16:27     Titel: Schreibweise ODE
  Antworten mit Zitat      
Hallo Leute,

ich würde gerne einen Solver mit Hilfe von Events steuern.

Dafür habe ich in meinem Programm bisher folgendes stehen, was sich auch ausführen lässt:
Code:
%main

tspan=[0 1];
u0=[0 0]
options=odeset('Events',@event_bew,'Refine',ref);

[t,y]=ode23(@f_bew,tspan,u0,options)


Nun würde ich aber gerne meiner Funktion f_bew noch ein paar Variablen mitgeben, z.B. a,b,c. Wie genau muss dazu der Code aussehen?

So funktioniert es bei mir nicht:

Code:
a=1;
b=2;
c=3;
tspan=[0 1];
u0=[0 0]
options=odeset('Events',@event_bew,'Refine',ref);

[t,y]=ode23(@f_bew,tspan,u0,options,[b]a,b,c[/b])


Ich danke euch schon mal im Voraus!!!!


Winkow
Moderator

Moderator



Beiträge: 3.842
Anmeldedatum: 04.11.11
Wohnort: Dresden
Version: R2014a 2015a
     Beitrag Verfasst am: 06.05.2014, 16:33     Titel:
  Antworten mit Zitat      
hi eine möglichkeit, wenn ich dich richtig verstehe, sind anonyme funktionen. so in der art
:
Code:
a=1
b=2
c=3
f=@(x) x*a*b*c

und das dann halt in den sover packen
_________________

richtig Fragen
Private Nachricht senden Benutzer-Profile anzeigen
 
Gast123456

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 06.05.2014, 16:40     Titel:
  Antworten mit Zitat      
Die Sache ist die, dass ich gerne in der Funktion noch ein paar Rechenoperationen ausführen möchte:

Code:
function [y]= f_bew(t, u ,a,b,c)

y=[u(2,:) ; a*b+c];

end
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.501
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 06.05.2014, 17:12     Titel:
  Antworten mit Zitat      
Hallo,

wenn etwas nicht funktioniert, bitte immer die Fehlermeldung angeben - insbesondere dann, wenn es nicht reproduzierbar ist, da man kein vollständiges Beispiel hat.

Aus der Doku von ode23 ist folgende Seite verlinkt, die zeigt, wie man mit anonymen Funktionen zusätzliche Daten übergeben kann
http://www.mathworks.com/help/matla.....-functions.html#bsgprpq-8

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Gast123456

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 10.05.2014, 20:33     Titel:
  Antworten mit Zitat      
Ich hätte gerne, dass der Parameter p, der sich bei jedem Schleifendurchgang um 1000 erhöht, an den Solver (bzw. f_test) weitergegeben wird.

Mein Main-Programm sieht folgendermaßen aus und lässt sich auch ausführen:

Code:
clear all; clc;


t0=0;
te=0.002;
t=t0

%leere Vektoren
tout=[];
yout=[];

%Anfangswerte
x0=[0; 0];

% Optionen für Integrator
options=odeset('Events',@event_test);

p=10^5;

while te > t

%f = @(t,x0) f_test(t,x0,p)

[t,x,tE,yE,iE]=ode23(@f_test,[t0 te],x0,options);

%[t,x,tE,yE,iE]=ode23(@f,[t0 te],x0,options);

tout=[tout;t]                      %hängt nach dem 1.Schritt den Vektor t an den leeren tout
yout=[yout;x]

x0=yout(end,:);
x0(2)=-0.3*x0(2)

%p=p+1000;

t0=tout(end)


end

figure;
plot(tout,yout(:,1))
grid on;

figure;
plot(tout,yout(:,2))
grid on;


Und hier die Funktion f_test:

Code:
function [y]= f_test(t,x)
%[y]= f_test(t,x,p)

A=1000;
m=3;      


F=1*10^6;

y=[x(2) ; F/m];
%y=[x(2) ; (F-p*A)/m];

end


Und hier die Event-Funktion:

Code:
function [q, isterminal, direction] = event_test(t, y)

q=[y(1)-0.046 y(2)];            


isterminal=[1 1];        
direction=[1 0];          

end





Die Überlegungen die ich angestellt habe, sind auskommentiert.

Wenn ich die Prozentzeichen entferne, erscheint diese Fehlermeldung:

Error: File: test.m Line: 24 Column: 23
"f" was previously used as a variable, conflicting
with its use here as the name of a function or
command.

Freue mich auf eure Rückmeldungen!
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.501
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 10.05.2014, 21:27     Titel:
  Antworten mit Zitat      
Hallo,

wenn du f als anonymous function handle definierst, musst du das direkt (d.h. ohne @ davor) übergeben - es ist ja schon ein Function Handle.

Code:
[t,x,tE,yE,iE]=ode23(f,[t0 te],x0,options);


Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



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 | goMatlab RSS Button 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.