| 查看: 2043 | 回复: 0 | |||
[求助]
求助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)是否必须需要等于一个常量，方程组才能给出数值解。 |
回复此楼» 猜你喜欢
人气不行了
已经有10人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有6人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
matlab解微分方程
已经有10人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】急请微分方程高手看下这个方程组！
已经有14人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
20