| 查看: 4742 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
acnes实习版主
|
[求助]
matlab数值求解边界条件微分方程组
|
||
期望高手指点一下我的疑问，指出一个大方向也好，我现在就是处于无从下手的阶段，求点拨 解一个方程组，文献上给出了如下的形式： du/dx+dv/dx=0 u=dn/dx v=-dm/dx du/dx=-exp(-x)+n*m/(m+n) 边界条件和初值条件 u(p)=0,u(0)=v(p),v(0)=0 n(0)=q 其中函数自变量为x，边界条件和初值条件中p为常量，q为变量 按照matlab求解边界条件微分方程bvp4c方法，我进行了如下处理： y(1)=n,y(2)=m,y(3)=u,y(4)=v y'(1)=y(3),y'(2)=-y(4),y'(3)=-exp(-x)+y(1)*y(2)/(y(1)+y(2)),y'(4)=exp(-x)-y(1)*y(2)/(y(1)+y(2)) 方程M文件 function dydx=ivode(x,y) dydx=[y(3);-y(4);-exp(-x)+2*y(1)*y(2)/(2*y(2)+y(1));exp(-x)-2*y(1)*y(2)/(2*y(2)+y(1))]; 边界条件M文件 function res=ivbc(ya,yb) res=[ya(1)-q;ya(4);yb(3);ya(3)-yb(4)]; command: solinit=bvpinit(linspace(0,1,10),[1 1 0 0]); sol=bvp4c(@ivode,@ivbc,solinit); ??? Error using ==> bvp4c Unable to solve the collocation equations -- a singular Jacobian encountered 请问这个问题怎么解决呢？还有边界条件中u(0)=v(p)是否必须需要等于一个常量，方程组才能给出数值解。手工处理du/dx=-exp(-x)+n*m/(m+n)方程，分离变量也没弄明白。。。曾经学的高数都还给老师了。。。 |
回复此楼» 猜你喜欢
309分085801求调剂
已经有7人回复
-
生物学308分求调剂（一志愿华东师大）
已经有3人回复
-
求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过
已经有3人回复
-
环境工程 085701，267求调剂
已经有14人回复
-
308求调剂
已经有12人回复
-
求调剂：085600材料与化工，考材科基，总分319
已经有21人回复
-
311求调剂
已经有8人回复
-
294分080500材料科学与工程求调剂
已经有13人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有19人回复
-
288资源与环境专硕求调剂，不限专业，有学上就行
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 如何用Matlab求时变微分方程的解？
已经有5人回复
- matlab求解高阶方程组
已经有3人回复
- MATLAB求解非线性方程组
已经有5人回复
- 关于matlab求解非线性指数方程组出现问题
已经有3人回复
- 有偿求助：matlab怎么解联立方程组？
已经有11人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 帮忙用Mathematica求解一个微分方程组
已经有6人回复
- MATLAB求解代数方程组问题，请帮忙
已经有9人回复
- Matlab求解偏微分方程组
已经有22人回复
- 请问如何求解二元一阶微分方程组
已经有24人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 一个偏微分方程组（包括边界条件）如何求解得到自振频率和振型？
已经有4人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】请教非齐次常微分方程组的解析解法
已经有4人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】求教这种微分方程的详细求解过程
已经有14人回复
- 【分享】KS方程求解器-matlab版
已经有10人回复
dbb627
新虫
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
6楼2011-12-21 11:17:32
wangyuewu08
兑换贵宾
- 应助: 13 (小学生)
- 金币: 2794
- 红花: 4
- 帖子: 1141
- 在线: 279.1小时
- 虫号: 1290073
- 注册: 2011-05-09
- 性别: GG
- 专业: 制造系统与自动化
2楼2011-12-20 18:38:48
acnes
专家顾问
- 应助: 3 (幼儿园)
- 金币: 897.3
- 散金: 328
- 红花: 6
- 帖子: 370
- 在线: 158.9小时
- 虫号: 440096
- 注册: 2007-10-26
- 性别: GG
- 专业: 理论和计算化学
3楼2011-12-20 19:08:24
wangyuewu08
版主
- 应助: 13 (小学生)
- 金币: 2794
- 红花: 4
- 帖子: 1141
- 在线: 279.1小时
- 虫号: 1290073
- 注册: 2011-05-09
- 性别: GG
- 专业: 制造系统与自动化
4楼2011-12-20 20:10:30
