Ein System von ODEs in MATLAB lösen?

Question

Dies ist mein erstes Mal, also wenn du mehr Informationen brauchst oder etwas nicht richtig gemacht hast lass es mich wissen!

Ich möchte ein System von ODES plotten, das wäre kein Problem, wie ich das folgende tun würde;

F=@(t,x) [ 
.the RHS of my first order ODES. 
]; 


[t x]=ode45(F,[Range], [Initial conditions]); 

Wo ein Vektor x würde erstellt werden, die '(1), X' (2), ..., wobei X (1), x (2) würde die zeitabhängig variabel sein.

Allerdings ist mein System ein wenig komplizierter. Ich habe 6 Differentialgleichungen, die aus Gleichungen bestehen, die Differentialgleichungen und die zeitabhängigen Variablen enthalten. Zum Beispiel sind die RHS meiner ODES so etwas wie 2 * x (2) * x '(1) * f (a)

Wo f (a) könnte eine andere Funktion basierend auf Konstanten sein, a.

Ich habe meinen Code in der folgenden Reihenfolge eingerichtet;

Konstanten, dann Gleichungen der Form f (a), dann meine Differentialgleichungen gefolgt vom ode45 Solver und den Plotbefehlen. Allerdings bekomme ich mehrere Fehler "Undefinierte Funktion oder Variable", da frühe Gleichungen von Variablen/Gleichungen abhängen, die erst später definiert werden.

Vielen Dank schon mal so viel für Sie :)

Wie vorgeschlagen, hier ein Beispiel für die Art meines Codes helfen;

`%Constants 
a=34 
b=31 
c=20 
%Equations 
A=b*cos(2*pi) 
B=a*EQ1 
C=c*x(2) 
%DifferentialEquations 
EQ1=x(1)*A 
EQ2=(EQ3-EQ1)*(B-C) 
EQ3=x(2)*x(3) 
F=@(t,x) [EQ1;EQ2;EQ3;];[t x]=ode45(F,[0 10], [0 0 0 ]); 

Liefert den Fehler. Erzeugt die undefinierte Funktion oder Variable 'DEQ1A'.

Answer 1

Basierend auf dem Kommentar versuchen Sie eine Variable zu verwenden, bevor sie definiert wird. In MATLAB müssen Sie Ihre Variablen und Funktionen definieren, bevor Sie sie verwenden. Da die Gleichungen Funktionen voneinander sind, können Sie anonymous functions verwenden, um F zu definieren.

Beispiel:

a=34; 
b=31; 
c=20; 
A = b*cos(2*pi); 
EQ1 = @(x) x(1)*A; 
B = @(x) a*EQ1(x); 
C = @(x) c*x(2); 
EQ3 = @(x) x(2)*x(3); 
EQ2 = @(x) (EQ3(x) - EQ1(x))*(B(x)-C(x)); 
F = @(t,x) [EQ1(x);EQ2(x);EQ3(x)]; 
[t,x] = ode45(F,[0 10], [0 0 0]);