| 查看: 1378 | 回复: 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 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有27人回复
-
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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