Gast123456
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am : 06.05.2014, 16:27
Titel : Schreibweise ODE
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:
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
Beiträge: 3.842
Anmeldedatum: 04.11.11
Wohnort: Dresden
Version: R2014a 2015a
Verfasst am : 06.05.2014, 16:33
Titel :
hi eine möglichkeit, wenn ich dich richtig verstehe, sind anonyme funktionen. so in der art
:
und das dann halt in den sover packen
_________________richtig Fragen
Gast123456
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am : 06.05.2014, 16:40
Titel :
Die Sache ist die, dass ich gerne in der Funktion noch ein paar Rechenoperationen ausführen möchte:
Harald
Forum-Meister
Beiträge: 24.501
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
Verfasst am : 06.05.2014, 17:12
Titel :
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
Gast123456
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am : 10.05.2014, 20:33
Titel :
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:
Und hier die Event-Funktion:
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
Beiträge: 24.501
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
Verfasst am : 10.05.2014, 21:27
Titel :
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.
Grüße,
Harald
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 - 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.