| 查看: 741 | 回复: 3 | ||
中将木虫 (小有名气)
|
[求助]
求高手帮我分析下我这个matlab求解非线性微分方程组解的情况
|
|
我这个是一组非线性方程组，可是每次给不同的初值都达不到我要求结果，不知道怎么办，求大神指教，不知道是方法不对还是怎么的，我希望结果输出的都是正数在0和1之间，而且最后的total要求最好小于1。。。。。。。。这样的非线性方程组怎么求好呢，有什么好的方法呢，求教求教 这个是我函数方程 function F=fun(x) ya0=0.05;ye0=0.2;yf0=0.75; Kp1=6.6090e9;Kp2=568.35;Kp3=2.4474e5;Kp5=6.6946e16;Kp6=3.3131e-3;Kp8=1.4006e13; F=[Kp1*((1+2/3*x(2)-2*x(5)-2*x(4))*ya0-1/3*(3*x(1)+2*x(2)+x(3)+x(5)+x(4)+x(6)))*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))-x(2)*((1+2/3*x(2)-2*x(5)-2*x(4))*yf0+x(1)-x(5)+x(6)); Kp2*((1+2/3*x(2)-2*x(5)-2*x(4))*ya0-1/3*(3*x(1)+2*x(2)+x(3)+x(5)+x(4)+x(6)))*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))- x(2)*x(3); Kp3*((1+2/3*x(2)-2*x(5)-2*x(4))*ya0-1/3*(3*x(1)+2*x(2)+x(3)+x(5)+x(4)+x(6)))*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))^2-x(3)^3; Kp5*((1+2/3*x(2)-2*x(5)-2*x(4))*ya0-1/3*(3*x(1)+2*x(2)+x(3)+x(5)+x(4)+x(6)))*((1+2/3*x(2)-2*x(5)-2*x(4))*yf0+x(1)-x(5)+x(6))^3-10e6*x(5)^3*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))^7; Kp6*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))*x(5)-x(4)*((1+2/3*x(2)-2*x(5)-2*x(4))*yf0+x(1)-x(5)+x(6)); Kp8*x(3)*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))-x(6)*((1+2/3*x(2)-2*x(5)-2*x(4))*yf0+x(1)-x(5)+x(6)); ]; 下面是调用 x0=[0.1;0.1;0.1;0.1;0.1;0.1]; options=optimset('Display','iter'); y=fsolve(@fun,x0,optimset('fsolve')) yA=((1+2/3*y(2)-2*y(5)-2*y(4))*0.1-1/3*(3*y(1)+2*y(2)+y(3)+y(5)+y(4)+y(6))) yE=((1+2/3*y(2)-2*y(5)-2*y(4))*0.5-1/3*(3*y(1)+y(2)+2*y(3)-7*y(5)+5*y(6)-4*y(4))) yF=((1+2/3*y(2)-2*y(5)-2*y(4))*0.4+y(1)-y(5)+y(6)) > total=y(1)+y(2)+y(3)+y(4)+y(5)+y(6)+yA+yE+yF |
回复此楼» 猜你喜欢
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有32人回复
-
物理学I论文润色/翻译怎么收费?
已经有72人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有78人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- matlab解偏微分方程求助
已经有12人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 二阶非线性微分方程求解
已经有14人回复
- 请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。
已经有17人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 请教用matlab求解一个非线性偏微分方程组的数值解
已经有8人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab求解非线性方程组
已经有16人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab解微分方程组
已经有15人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
中将
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1958.7
- 散金: 86
- 红花: 1
- 帖子: 221
- 在线: 76小时
- 虫号: 886327
- 注册: 2009-10-28
- 性别: GG
- 专业: 分离过程
2楼2013-06-22 16:52:34
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★
感谢参与，应助指数 +1
csgt0: 金币+1, 谢谢 2013-06-24 17:16:47
感谢参与，应助指数 +1
csgt0: 金币+1, 谢谢 2013-06-24 17:16:47
|
与初值选取不当有个吧，试试1stOpt，不用初值，不过本题似乎有多组解： 1： x1: 0.452159092308927 x2: -3.5576820815459E-17 x3: -1.6802621629643E-17 x4: 0.70583278260955 x5: 0.289090165134331 x6: 0.579315494431816 2： x1: 2.3247555351878 x2: -3.51606195970296 x3: 2.54582708579748E-8 x4: -1.05916256056018 x5: 1.08911300900033 x6: -0.182685684997559 |
3楼2013-06-24 09:06:15
中将
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1958.7
- 散金: 86
- 红花: 1
- 帖子: 221
- 在线: 76小时
- 虫号: 886327
- 注册: 2009-10-28
- 性别: GG
- 专业: 分离过程
4楼2013-06-27 12:24:25
40