| 查看: 1280 | 回复: 4 | ||
驿边桥CQU新虫 (小有名气)
|
[求助]
边界条件微分方程用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 |
回复此楼» 本帖@通知
» 猜你喜欢
310求调剂
已经有14人回复
-
290调剂生物0860
已经有28人回复
-
山东双非院校考核超级无底线，领导幸灾乐祸，教师遭殃恐
已经有4人回复
-
279学硕食品专业求调剂院校
已经有14人回复
-
0831生医工第一轮调剂失败求助
已经有16人回复
-
电子信息270求调剂
已经有18人回复
-
食品与营养（0955）271求调剂
已经有5人回复
-
293求调剂
已经有11人回复
-
322求调剂
已经有13人回复
-
272分材料子求调剂
已经有50人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程初值问题和边值问题的区别？
已经有5人回复
- 关于二阶偏微分方程求解
已经有21人回复
- 用fortran解一个二阶微分方程
已经有14人回复
- 有限差分法 第二类边界条件
已经有12人回复
- 请问这个偏微分方程组有一般的解法么？
已经有12人回复
- 求解一个微分方程，另请教微分方程如何确定其有无解？
已经有7人回复
- matlab 非线性微分方程求解
已经有3人回复
- 一个椭圆型偏微分方程的数值解
已经有19人回复
- MATLAB 解偏微分方程时的边界条件传递问题
已经有4人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 如何用Matlab求时变微分方程的解？
已经有5人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 抛物型偏微分方程的数值解求解问题，谢谢！
已经有14人回复
- 数学物理方程求解！（带边界条件的偏微分方程）
已经有13人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- Matlab: 用bvp解二阶微分方程（边界值）遇到困难，请高手指点
已经有5人回复
- 1stOpt求解常微分方程边值问题
已经有5人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 【求助】用pdepe求解偏微分方程
已经有10人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
 本帖仅楼主可见 |
2楼2015-04-10 10:26:55
 |
3楼2015-04-10 10:30:10
驿边桥CQU
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 80
- 散金: 5
- 帖子: 60
- 在线: 36.1小时
- 虫号: 2838211
- 注册: 2013-11-30
- 专业: 可再生与替代能源利用中的
帮忙指点下吧 4楼2015-04-10 11:59:54
驿边桥CQU
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 80
- 散金: 5
- 帖子: 60
- 在线: 36.1小时
- 虫号: 2838211
- 注册: 2013-11-30
- 专业: 可再生与替代能源利用中的
5楼2015-04-17 22:38:38
10