| 查看: 496 | 回复: 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个金币略表谢意。 |
回复此楼» 猜你喜欢
【26年入学】计算机视觉/电子信息/诚招博士生 光电教育部重点实验室
已经有2人回复
-
【26年入学】红外半导体/电子信息/嵌入式视觉/诚招博士生 光电教育部重点实验室
已经有0人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有56人回复
-
CSC难上岸，25年东京大学，早稻田，京都大学多Lab大量招募日本全奖博士生！
已经有6人回复
-
【26年入学】红外半导体/电子信息/嵌入式视觉/诚招博士生 光电教育部重点实验室
已经有0人回复
-
小论文投稿
已经有3人回复
-
谁让他附加功能艺术
已经有0人回复
-
很难说荣景园收入
已经有0人回复
-
好大太阳今年生日
已经有0人回复
-
读博入学申请
已经有0人回复
-
2026 博士自荐-机器人机构学方向
已经有12人回复
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