MATLABでODEのシステムを解く？

Question

これは私の初めての投稿ですので、もう情報が必要な場合、または何か正しいことをしていない場合はお知らせください！

ODESのシステムをプロットしたいと思います。以下のように問題はありません。ベクトルは、xを含む作成される

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


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

'（1）、X'（2）、...、X（1）、X（2）時間依存変数であろう場所。

私のシステムはもう少し複雑です。私は微分方程式と時間依存変数を含む方程式からなる6つの微分方程式を持っています。例えば、私のODESのRHSは、 のようなものです。2 * x（2）* x '（1）* f（a）

ここで、f（a）は定数に基づく別の関数です。

以下の順序でコードを設定しました。

定数、f（a）の形式の方程式、次に微分方程式の後にode45ソルバとplotコマンドが続きます。しかし、私は初期の方程式が後で定義されない変数/方程式に依存するため、いくつかのエラー "Undefined function or variable"を得ています。

は今まであなたが私のコードの種類の提案、ここでは一例として:)

を助けるために多くをいただきありがとうございます。

`%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 ]); 

エラーを出力します。 未定義の関数または変数 'DEQ1A'を生成します。

Answer 1

変数が定義される前に変数を使用しようとしているコメントに基づいています。 MATLABでは、使用する前に変数と関数を定義する必要があります。方程式はお互いの関数なのでanonymous functionsを使ってFを定義することができます。

例：

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]);