| 查看: 1841 | 回复: 7 | |||
孤竹行者金虫 (小有名气)
|
[求助]
matlab非线性拟合求助
|
小弟是个matlab菜鸟，现在有一个方程想通过matlab进行拟合，得到其中的参数，文献中说是用Nelder-Mead simplex method，但是自己不知道怎么编写程序，望哪位大侠予以帮助，提供matlab程序，不胜感激！文献中的公式和实验数据如附件所示。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 方程及文献数据.doc
2013-04-05 23:28:45, 36.5 K
» 收录本帖的淘帖专辑推荐
 纳米 |
» 猜你喜欢
理学，工学，农学调剂，少走弯路，这里欢迎您！
已经有8人回复
-
276求调剂
已经有6人回复
-
求调剂
已经有7人回复
-
085600求调剂
已经有5人回复
-
298求调剂
已经有13人回复
-
0856材料调剂
已经有5人回复
-
0856材料与化工，270求调剂
已经有10人回复
-
285求调剂
已经有9人回复
-
材料工程专硕283求调剂
已经有7人回复
-
环境调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于matlab微分、及非线性拟合问题
已经有5人回复
- matlab非线性拟合计算求助
已经有10人回复
- matlab数据拟合问题求助
已经有7人回复
- 诚恳求助如何做IC50拟合曲线
已经有7人回复
- matlab非线性参数拟合问题
已经有7人回复
- 急切求助用matlab曲线拟合色散方程
已经有10人回复
- 求助非线性拟合问题
已经有5人回复
- 【求助】Matlab非线性最小二乘拟合活度系数模型（Willson、NRTL、UNIQUAC）
已经有17人回复
- 【求助】急求关于matlab做非线性拟合的函数
已经有8人回复
- 【求助】如何使得matlab拟合的曲线强制经过一个点
已经有8人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
- 【求助】matlab 遗传算法工具箱绘图命令是什么？
已经有6人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【求助】matlab有公式的曲线拟合
已经有7人回复
- 【求助】如何把origin中曲线的拟合线延长？
已经有7人回复
孤竹行者
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1103.9
- 散金: 5
- 帖子: 137
- 在线: 52.7小时
- 虫号: 825902
- 注册: 2009-08-11
- 性别: GG
- 专业: 能源化工
2楼2013-04-07 15:03:47
littlepoint
木虫 (职业作家)
- 应助: 9 (幼儿园)
- 金币: 12030.1
- 散金: 8005
- 红花: 12
- 帖子: 3095
- 在线: 334.6小时
- 虫号: 1065217
- 注册: 2010-07-27
- 专业: 工程热力学
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢参与！ 2013-04-07 23:18:21
孤竹行者: 金币+5, ★★★很有帮助 2013-04-15 12:36:21
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢参与！ 2013-04-07 23:18:21
孤竹行者: 金币+5, ★★★很有帮助 2013-04-15 12:36:21
| 这种问题自己先试着解一下，看看基本原理。这种上来就求源代码的搞法不易于自己的提高哟。http://www.pudn.com/downloads62/sourcecode/math/detail213534.html |
3楼2013-04-07 15:27:05
孤竹行者
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1103.9
- 散金: 5
- 帖子: 137
- 在线: 52.7小时
- 虫号: 825902
- 注册: 2009-08-11
- 性别: GG
- 专业: 能源化工
4楼2013-04-10 08:25:17
孤竹行者
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1103.9
- 散金: 5
- 帖子: 137
- 在线: 52.7小时
- 虫号: 825902
- 注册: 2009-08-11
- 性别: GG
- 专业: 能源化工
5楼2013-04-15 12:35:30
sujiguang
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 292.3
- 散金: 78
- 帖子: 239
- 在线: 52.5小时
- 虫号: 1130581
- 注册: 2010-10-24
- 专业: 药物化学
|
NRTL模型的MATLAB程序： function f=nrtl(g) syms r1 r2 r3 r4 a t12 t21 R g21 g12 x2 x4 s1 s2 x1 x3 T; T=298.15; x1=O.999889456; x3=O.01617526; x2=l-x1; x4=l-x3; a=O.20; R=8.314; t21=g(1)/(R*T); t12=g(2)/(R*T); g21=exp(-a*t21); g12=exp(-a*t12); r1=exp(x2^2*(t12*(g12/(x1+x2*g21))^2+t12*g12/(x2+x1*g12)^2)); r2=exp(x1^2*(t12*(g12/(x2+x1*g12))^2+t21*g21/(x1+x2*g21)^2)); r3=exp(x4^2*(t12*(g12/(x3+x4*g21))^2+t12*g12/(x4+x3*g12)^2)); r4=exp(x3^2*(t12*(g12/(x4+x3*g12))^2+t21*g21/(x3+x2*g21)^2)); s1=x1*r1-x3*r3; s2=x2*r2-x4*r4; f=s1^2+s2^2 s1 s2 r1 r2 r3 r4 [x,fval]=fiminsearch(x1,x2) 不懂额，你有的话发我一份吧，我也是急用额，想要NRTL和UNIQUAC的计算活度系数的matlab程序 |
6楼2013-07-08 11:30:39
孤竹行者
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1103.9
- 散金: 5
- 帖子: 137
- 在线: 52.7小时
- 虫号: 825902
- 注册: 2009-08-11
- 性别: GG
- 专业: 能源化工
|
我这有一个自己编的程序，初步计算了一下wilson方程还算可以，你把里面的方程改成你的NRTL和UNIQUAC方程看看行不行！ function f = wilson1(Input) a12=Input(1); a13=Input(2); a21=Input(3); a23=Input(4); a31=Input(5); a32=Input(6); b12=Input(7); b13=Input(8); b21=Input(9); b23=Input(10); b31=Input(11); b32=Input(12); F=0; T=[ 298.45 302.7 308.05 313.35 317.75]; X1=[0.0031549 0.0037933 0.0048003 0.0061385 0.0074595]; X2=[0 0 0 0 0 0 0 0 0 ]; for i=1:length(X1) t=T(i); x1=X1(i); x2=X2(i); x3=1-x2; lnr1=1-log(x1+x2*(a12+b12/t)+x3*(a13+b13/t))-x1/(x1+x2*(a12+b12/t)+x3*(a13+b13/t))-x2*(a21+b21/t)/(x1*(a21+b21/t)+x2+x3*(a23+b23/t))-x3*(a31+b31/t)/(x1*(a31+b31/t)+x2*(a32+b32/t)+x3); A=21; B=9050; tcal=B/(A-lnr1-log(x1)); v=(tcal-t).^2; F=F+v; ff=(F/(length(X1)-1))^0.5; end f=ff; 有问题再多多交流！ |
7楼2013-07-10 14:10:04
sujiguang
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 292.3
- 散金: 78
- 帖子: 239
- 在线: 52.5小时
- 虫号: 1130581
- 注册: 2010-10-24
- 专业: 药物化学
，还是不会，哎，我再琢磨琢磨哈！为什么都是input未定义 8楼2013-07-10 19:21:17
5