我想要一个代码解决一个耦合的三阶和二阶颂歌使用拍摄方法和RK-4数值技术,请如果任何人都可以帮助
15的观点(30天)
显示旧的评论
(1 + 2米* eta) f ' + 2 m * f + * f - f的^ 2 - k1 * f +λ*θ= 0 - - - - - - - - - - - (1)
(1 + 2米* eta)θ“+ 2 m *θ+公关(f *θ“- f”*θ= 0 - - - - - - - (2)
“f”和“θ”是“埃塔”的功能,“埃塔”是一个独立的变量
3给出初始条件:η= 0,f (0) = 0, f(0) = 1,θ(0)= 1
说我减少这些方程(1)和(2)五颂歌(射击方法)
f = z;f (0) = 0 - - - - - - (3)
z ' = p;z (0) = 1 - - - - - - (4)
p ' =(2米*“- f * f + f ' ^ 2 + k1 * f 'λ*θ)/ (1 + 2 * M *η);p(0) =(猜值)——(5)
θ' =问;θ(0)= 1 - - - - - - (6)
q =(2米*θ的公关(f *θ“- f”*θ))/ (1 + 2 * M *η);问(0)=(猜值)- - - - - - (7)
需要满足的边界条件是:f(η= 10)= 0和θ(η= 10)= 0η= 10
考虑到:
M = 1
k1 = 0.1
λ= 0.1
公关= 0.7
以步长:h = 0.01
接受的答案
Torsten
2017年11月22日
试一试
函数主要
全球公关k1 Mλ
公关= 0.7;k1 = 0.1;M = 1;λ= 0.1;
rlow = 0;
rhigh = 10;
N = 1000;
选择= bvpset (“统计数据”,“上”,“RelTol”1 e-5);
solinit = bvpinit (linspace (rlow rhigh N), [0 1 0 1 0]);%
索尔= bvp4c (@projode @mybcs、solinit选项);
函数dy = projode (n, y)
全球公关k1 Mλ
dy = [y (2);y (3);(2 * M * y (3) - y (1) * (3) + y (2) ^ 2 + k1 * y(2)λ* y (4)) / (1 + 2 * M * n);y (5);(2 * M * y(5)公关* y(1) *(5) +公关* y (2) * y (4)) / (1 + 2 * M * n)];
函数res = mybcs (ya,)
res =[丫(1);丫(2)-1.0;丫(4)-1.0;yb (2);yb (4)];
最好的祝愿
Torsten。