| 查看: 514 | 回复: 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个金币略表谢意。 |
回复此楼» 猜你喜欢
国产vs进口半导体制冷模块：实测对比4大本土品牌，差距缩小了吗？
已经有0人回复
-
催化剂焙烧炉供应商怎么选？优质厂商深度解析与推荐
已经有0人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有254人回复
-
南方科技大学机械与能源工程系-邓辉研究员招收2026级博士生
已经有0人回复
-
西北工业大学民航学院招博士与硕士复合材料方向
已经有30人回复
-
PEG-HMDI预聚体接枝环氧树脂以及LNBR的讨论
已经有11人回复
-
广东工业大学精密电子制造技术与装备国家重点实验室刘江文教授课题组招收博士生、博后
已经有0人回复
-
一站式解决工业环境气象监测难题：优选综合服务商型防爆气象站厂家推荐
已经有0人回复
-
南方科技大学机械与能源工程系-邓辉研究员招收2026级博士生
已经有0人回复
-
IEEE会议Ei收录高校主办-IEEE第三届能源与电气工程国际学术会议-EEE 2026
已经有0人回复
-
Ei会议-第六届流体与化学工程国际学术会议-ICFCE2026
已经有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
