| 查看: 2905 | 回复: 7 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
lgycjpcqu金虫 (正式写手)
|
[求助]
求助 MATLAB解方程组-fslove
|
|
|
求解x,y,z的数值，我用MATLAB 的fslove函数求解，解出来结果为 >> x=fsolve(@myfun1,[0.8 0.9 0.05],optimset('Display','off')) x = 0.800000000000194 0.900000000000932 0.050000000000305 这个结果似乎就是在初始值附件变动了一下，我改变初始值，结果也发生变化，是不是解法有问题，下面是我写的程序 function q=myfun1(p) K1=10^(-12.43); K2=10^(19.4); K3=5.6*10^(-20); Ccu=0.24; Coh=0.75; Csa=0.26; x(1)=p(1); x(2)=p(2); x(3)=p(3); q(1)=((1/Coh)*(x(1)/((1-x(1)-x(2))*(1-x(1)-3*x(2)-2*x(3)))))-K1; q(2)=((1/(Csa*Coh^3))*(x(2)/((1-x(2)-x(3))*(1-x(1)-x(2))*(1-x(1)-3*x(2)-2*x(3))^3)))-K2; q(3)=(Ccu*Coh^2*(1-x(2)-x(3))*(1-x(1)-3*x(2)-2*x(3))^2)-K3; %x=fsolve(@myfun1,[0.8 0.9 0.05],optimset('Display','off')); 要是该方程解出来就是这个这个结果，就说明我方程有问题，我就去重新考虑数学模型，现在求高手帮忙看看这些方程的求解是否有问题，多谢！或者哪位高手能够有其它方法准确的求解这个方程，还请指点一下。  方程组及参数 [ Last edited by lgycjpcqu on 2012-7-16 at 16:13 ] |
回复此楼» 本帖@通知
» 猜你喜欢
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有6人回复
-
物理学I论文润色/翻译怎么收费?
已经有149人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有30人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有62人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高金求助matlab解微分方程组
已经有12人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】如何解这样的方程组【有效期至2011年4月30日】
已经有9人回复
- 【求助】MATLAB 解方程组【已解决】
已经有7人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【求助】一个时变微分方程组的matlab实现
已经有8人回复
5楼2012-07-16 18:06:18
zt_chem
木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 2147.8
- 红花: 6
- 帖子: 872
- 在线: 981.6小时
- 虫号: 286920
- 注册: 2006-10-15
- 性别: GG
- 专业: 催化化学
2楼2012-07-16 17:00:58
lgycjpcqu
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 566.8
- 散金: 411
- 帖子: 371
- 在线: 176.9小时
- 虫号: 1327328
- 注册: 2011-06-20
- 性别: GG
- 专业: 电化学
3楼2012-07-16 17:34:20
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2012-07-16 18:00:32
