| 查看: 2753 | 回复: 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 ] |
回复此楼» 本帖@通知
» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有119人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高金求助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人回复
lgycjpcqu
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 566.8
- 散金: 411
- 帖子: 371
- 在线: 176.9小时
- 虫号: 1327328
- 注册: 2011-06-20
- 性别: GG
- 专业: 电化学
3楼2012-07-16 17:34:20