| 查看: 523 | 回复: 2 | ||
lanxiboy新虫 (初入文坛)
|
[求助]
关于七元非线性方程组求解MATLAB源程序的求助，求各位大神帮帮忙，不胜感激！
|
|
小弟现在做的是微合金元素复合析出计算，由于以前没有用过Matlab，而现在又比较着急用方程组的结果，来不及学习Matlab，因此在论坛里求助各位大神帮忙。参数方程组如下： R=8.314; L=-4260; t=1220; syms x a m Nb Ti C N f1=a*log(x*a*power(10,3.42-7900/(t+273.15)))/(Nb*Ti)+(1-a)*log(x*(1-a)*power(10,2.8-8500/(t+237.15)))/(Nb*N)+L*a*(1-a)/(R*(t+237.15)) f2=x*log(x*a*power(10,3.42-7900/(t+237.15)))/(Nb*C)+(1-x)*log(a*(1-x)*power(10,5.33-10475/(t+237.15)))/(Ti*C)+L*(1-a)^2/(R*(t+237.15)) f3=x*log(x*(1-a)*power(10,2.8-8500/(t+273.15)))/(Nb*N)+(1-x)*log((1-x)*(1-a)*power(10,0.32-8000/(t+237.15)))/(Ti*N)+L*a^2/(R*(t+237.15)) f4=m*(x/2)+(1-m)*Nb-0.0002118; f5=m*((1-x)/2)+(1-m)*Ti-0.0006646; f6=m*(a/2)+(1-m)*C-0.0037567; f7=m*((1-a)/2)+(1-m)*N-0.0001273; f=[f1;f2;f3;f4;f5;f6;f7]; x=[x;a;m;Nb;Ti;C;N]; x0=[0.010;0.010;0.010;0.0002118;0.0006646;0.0037567;0.0001273]; esp=[0.00000001;0.00000001;0.00000001;0.00000001;0.00000001;0.00000001;0.00000001]; N=10; X=newton(f,x,x0,esp,N); double(X); 要求出[x;a;m;Nb;Ti;C;N7个未知量，其他均为已知参量，自己在网上搜索了好多程序，牛顿下山法什么的，不是不收敛就是出现错误；求各位大神帮小弟写个原程序代码，由于小弟是新人，刚注册的号，金币不多，全部23个金币略表谢意。 |
回复此楼» 猜你喜欢
诚邀加盟！青岛大学郑晓钦教授团队诚招博士研究生、博士后及青年英才
已经有7人回复
-
北京信息科技大学联合清华大学某国家重点实验室，联合培养硕生、博生，并长期招博士后
已经有16人回复
-
电气科学与工程论文润色/翻译怎么收费?
已经有234人回复
-
北京信息科技大学智能机器人技术方向，联合清华培养硕生、博生，并长期招博士后
已经有21人回复
-
CSC-UTT : CNRS 法国国家科学研究中心/utt 特鲁瓦工程技术大学 CSC博士phd全奖招生
已经有8人回复
-
TPI层压板
已经有0人回复
-
MOF与高分子配位交换反应复现
已经有0人回复
-
低共熔溶剂
已经有0人回复
-
石河子大学（211、双一流）硕博研究生长期招生公告
已经有3人回复
-
本2一志愿京9求调剂能源动力方向
已经有1人回复
314977903 
版主 (职业作家)
小风子
- 应助: 161 (高中生)
- 贵宾: 0.152
- 金币: 14189.3
- 散金: 5787
- 红花: 39
- 沙发: 32
- 帖子: 4586
- 在线: 902.7小时
- 虫号: 1204391
- 注册: 2011-02-16
- 性别: GG
- 专业: 金属功能材料
- 管辖: 材料工程
2楼2016-06-04 08:31:27
314977903
版主 (职业作家)
小风子
- 应助: 161 (高中生)
- 贵宾: 0.152
- 金币: 14189.3
- 散金: 5787
- 红花: 39
- 沙发: 32
- 帖子: 4586
- 在线: 902.7小时
- 虫号: 1204391
- 注册: 2011-02-16
- 性别: GG
- 专业: 金属功能材料
- 管辖: 材料工程
3楼2016-06-07 10:39:15
