フィルターのクリア

Solve first order system using bvp4c

5 ビュー (過去 30 日間)
Day Rosli
Day Rosli 2016 年 1 月 25 日
編集済み: Day Rosli 2016 年 1 月 26 日
Hi, I have problem to solve this in MATLAB. The equation already in first order system, yet I still confused to write the code and solve it. Or this equation cant be solved by using bvp4c perhaps?
Equation:
y(1)'=y(2);
((mu*y(2)')/R)=(((W-mu')*y(2))/R)+((i*R*(Abar*U+Bbar*V-Obar)+U)/R)*y(1)-(2*(1+V)/R)*y(5)+(Abarone*U'+Bbar*V')*y(3)+i*(kappa^2-((Abar*(i))/R))*y(4);
y(3)'=-i*y(1);
y(4)'=((i*(y(1)*W-(y(1)*mu)'))/R)-((i*R*(Abar*U+Bbar*V-Obar)+W')/R)*y(3);
y(5)'=y(6);
((mu*y(6)')/R)=(((W-mu')*y(6))/R)+((i*R*(Abar*U+Bbar*V-Obar)+U)/R)*y(5)+(2*(1+V)/R)*y(1)+(Abarone*V'-Bbar*U')*y(3)+i*(Bbar/R)*y(4);
with bc:
y(1)(0)=ya(2)(0)=ya(3)(0)=ya(4)(0)=ya(5)(0)=ya(6)(0)=0;
y(1)(infinity)=y(2)(infinity)=y(3)(infinity)=y(4)(infinity)=y(5)(infinity)=y(6)(infinity)=0;
where infinity = 20;
where;
i is imaginary number;
mu=((U'(x))^2+(V'(x))^2)^((n-1)/2);
Obar=(R^(n-1)/2)*omega;
Abarone=Abar-(i/R);
kappa=Abar^2+Bbar^2;
R=n=Abar=Bbar=1;
omega=0;
For U(x),V(x),W(x),U'(x),V'(x) and W'(x), the values I already got from previous computation. meaning each of them have values from x=0 until x=20.
The problem now is that,
1) How do I write the system in MATLAB since the mu' term makes it more confusing (for me).
2) Notice that there is U(x),V(x),W(x),U'(x),V'(x) and W'(x) which have their own values ranging from x=0 till x=20. How can I relate this
problem from my previous computation? Do I need to paste the values in the system for each U(x),V(x),W(x),U'(x),V'(x) and W'(x)?
Or is there other thing I can do?
3) Also, there is 12 boundary conditions where I only have 6 first order equations. Im so confused.
I hope someone could help me solving this. Im stuck with this almost 4 months.
Thank you in advance.

回答 (0 件)

カテゴリ

Help Center および File ExchangeOrdinary Differential Equations についてさらに検索

タグ

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by