| 查看: 1647 | 回复: 5 | ||
[求助]
matlab下bvp4c求解非线方程问题
|
|
请问各位大侠，我自己编了个用bvp4c求解器求解非线性边值条件方程组的程序。当方程中系数n=1时，即系数简化的情况下，程序能较好的运行并算出较好的数值结果。当n不等于1，程序提示出现奇异性，不能求解方程。请问这是给定初值不合适的原因吗还是有别的问题，能帮忙解决一下吗？谢谢！ global b t E EI eta H l r Es EIz Q EIA m n ha; b=100/10^9;l=1000/10^9;t=50/10^9;r=0.89; E=76*10^9;Es=1.22; EI=E*b*t^3/12; H=2*(t+b)*r; EIz=EI+0.5*Es*b*t^2+Es*t^3/6; EIA=E*b*t+2*Es*(b+t); m=H*l^2/EIz; %n=EIA*l^2/EIz; n=1; ha=H/EIA; eta=t/l; Q=-1; %以上是给定常数参数 myinit=[0 1 0 1 0.001 -0.1];%[0 1 0 0 7 -5],[0 1 0 1 0.1 -0.1] solint=bvpinit([0 0.25 0.5 0.75 1],myinit); sol=bvp4c(@ss01,@ss02,solint); function odefun1=ss01(x,y); %给出控制方程 global eta m n; odefun1=[y(2);-eta^2*y(4)*y(5);y(4);y(5);y(6);-m*y(5)-n*(y(2)+0.5*eta^2*y(4)^2)*y(5)]; return function bcfun1=ss02(ya,yb);%给出边界条件 global Q n eta ha; bcfun1=[ya(1);ya(4);ya(3);yb(5);yb(2)+0.5*eta^2*yb(4)^2+ha;eta*(yb(2)+0.5*eta^2*yb(4)^2)*yb(4)+eta*yb(6)/n+ha*eta*yb(4)-Q]; return |
回复此楼» 猜你喜欢
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
今天务委会开完了，明天出结果吗
已经有21人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab用多元线性最小二乘法求解多元线性方程的问题
已经有11人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 一个微分方程组，求解
已经有18人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- MATLAB非线性方程组求解~~~反应釜模拟
已经有9人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2012-01-17 23:14:07
3楼2012-01-18 11:00:31
4楼2012-01-18 11:01:15
5楼2012-01-18 11:01:42
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
6楼2012-01-18 13:52:28
28