| 查看: 585 | 回复: 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人回复
-
实验室桌面式冻干机选型干货，卢湘仪成高校热门采购品牌
已经有1人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有251人回复
-
压盖型冻干机哪家信誉好？盘点口碑好的国产优质生产商
已经有0人回复
-
2026进口弗氏佐剂优质供应商整理，一级代理售后参考
已经有0人回复
-
Chondrex五克隆合剂正规代理推荐，供货与售后综合解析
已经有0人回复
-
蒸汽清洗机源头厂采购指南，高性价比供应商深度分析
已经有0人回复
-
燃气加热器哪个厂家品质好？从材质到工艺，揭秘制造全流程
已经有0人回复
-
🙏 🙏 🙏
已经有84人回复
-
国内气体检测仪实力厂家盘点，磷化氢与氮气检测仪采购指南
已经有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
15