Azzera filtri
Azzera filtri

Solve first order system using bvp4c

5 visualizzazioni (ultimi 30 giorni)
Day Rosli
Day Rosli il 25 Gen 2016
Modificato: Day Rosli il 26 Gen 2016
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.

Risposte (0)

Community Treasure Hunt

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

Start Hunting!

Translated by