[求助] 边界条件微分方程用bvp4c求解，猜测值不同有不同曲线解的原因以及边界条件明显不符

如题所示，具体如下： 程序： function dydx=tketi(x,y) Ja=0.1784;Ra=126.5;Pr=0.0052; R=0;P=3.81; A=Ja*R/Ra/Pr;B=Ja/Ra/Pr; dydx=[y(2);B/y(1)^2-(y(2)^2+A)/y(1)-y(2)/x;0]; function bc=tketiBCjiashe(ya,yb) Ja=0.1784;Ra=126.5;Pr=0.0052;R=0;P=38.1; bc=[ya(2)-P/Ra/ya(3);ya(2)-1;ya(3)-ya(1)]; clear all solinit=bvpinit(linspace(20,170,4000),[40 0 10]); sol=bvp4c(@tketi,@tketiBCjiashe,solinit); x=linspace(20,170,4000); y=deval(sol,x);figure(9); plot(x,y(1,,\'-*r\');title(\'40 0 10\') 从所得图形中可以看到，边界条件y(170)=1明显不符，还有 当把solinit=bvpinit(linspace(20,170,4000),[40 0 10]);中的假设换为[10 0 10]时会得到不同的曲线解 对上面的问题还请大家指点！！！ 微分方程及边界条件.png 变形后方程.png 所得图形.jpg

回复此楼