| 查看: 494 | 回复: 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个金币略表谢意。 |
回复此楼» 猜你喜欢
南方科技大学机械与能源工程系-邓辉研究员招收2026级博士生
已经有0人回复
-
南方科技大学机械与能源工程系-邓辉研究员招收2026级博士生
已经有0人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有238人回复
-
国内的超声波分散机源头厂家，制造商、供应商都有哪些？要质量好的奥
已经有0人回复
-
2026锡箔片推荐企业/供应商/哪家做得好/哪个品牌好？
已经有0人回复
-
【2026】氧化铜颗粒推荐供应商/哪家强/代理商,氧化铜颗粒供货商有哪些?哪家好？
已经有0人回复
-
2026年锡箔片推荐企业/供应商/哪家做得好/哪个品牌好
已经有1人回复
-
杭州精浩超声波搅拌机 —— 纳米材料搅拌 “神器”，20 分钟搞定之前 1 小时的活！
已经有1人回复
-
南方科技大学机械与能源工程系-邓辉研究员招收2026级博士生
已经有0人回复
-
2026锡制固体/锡制液体推荐企业/供应商/推荐品牌
已经有0人回复
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