24小时热门版块排行榜    

查看: 2120  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

孤竹行者

金虫 (小有名气)

[求助] matlab非线性拟合求助

小弟是个matlab菜鸟,现在有一个方程想通过matlab进行拟合,得到其中的参数,文献中说是用Nelder-Mead simplex method,但是自己不知道怎么编写程序,望哪位大侠予以帮助,提供matlab程序,不胜感激!文献中的公式和实验数据如附件所示。
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 方程及文献数据.doc
  • 2013-04-05 23:28:45, 36.5 K

» 收录本帖的淘帖专辑推荐

纳米

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

是你的终究是你的
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujiguang

铁虫 (小有名气)

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

孤竹行者

金虫 (小有名气)

自己回复一下,千万别沉啊!
是你的终究是你的
2楼2013-04-07 15:03:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

littlepoint

木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孤竹行者

金虫 (小有名气)

引用回帖:
3楼: Originally posted by littlepoint at 2013-04-07 15:27:05
这种问题自己先试着解一下,看看基本原理。这种上来就求源代码的搞法不易于自己的提高哟。http://www.pudn.com/downloads62/sourcecode/math/detail213534.html

谢谢你的回复,不过我主要不是研究这方面的,而且我要急着用这个程序,一时间自己也没有时间和精力再从头学习爱,先希望高手们能帮我写出这个程序来,在下不胜感激!
是你的终究是你的
4楼2013-04-10 08:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[教师之家] 跳槽后在研项目怎么办? +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见