24小时热门版块排行榜    

查看: 1843  |  回复: 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的回帖

孤竹行者

金虫 (小有名气)

引用回帖:
6楼: Originally posted by sujiguang at 2013-07-08 11:30:39
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 ...

我这有一个自己编的程序,初步计算了一下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
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 成果系统访问量大,请15分钟后再尝试。由此给您造成的不便,敬请谅解。 +7 xhuama 2026-03-02 8/400 2026-03-02 22:33 by taidy
[考研] 调剂材料学硕 +3 词凝Y 2026-03-02 3/150 2026-03-02 21:51 by 杨杨杨紫
[考研] 求调剂 +3 Guo_yuxuan 2026-03-02 4/200 2026-03-02 21:18 by 无际的草原
[考研] 085600求调剂 +4 LRZZZZZZ 2026-03-02 5/250 2026-03-02 20:33 by 云旗yunqi
[考研] 材料284求调剂,一志愿郑州大学英一数二专硕 +15 想上岸的土拨鼠 2026-02-28 15/750 2026-03-02 20:13 by hypershenger
[考研] 中国科学技术大学材料与化工281求调剂,有科研和获奖经历 +5 wsxw 2026-03-02 5/250 2026-03-02 20:12 by hypershenger
[考研] 化学,材料,环境类求调剂 +7 考研版棒棒 2026-03-02 7/350 2026-03-02 19:56 by hypershenger
[考研] 306分材料调剂 +5 chuanzhu川烛 2026-03-01 6/300 2026-03-02 19:51 by 张晓芳0105
[考研] 275求调剂 +7 明远求学 2026-03-01 7/350 2026-03-02 19:22 by zhukairuo
[考博] 26超级电容器申博 +3 dhdjdjend 2026-02-25 3/150 2026-03-02 19:21 by 轻松不少随
[考研] 268求调剂 +3 好运连绵不绝 2026-03-02 3/150 2026-03-02 19:05 by zhukairuo
[考研] 一志愿华南理工大学材料与化工326分,求调剂 +3 wujinrui1 2026-02-28 3/150 2026-03-02 16:36 by chuocheng
[考研] 284求调剂 +10 天下熯 2026-02-28 11/550 2026-03-02 11:03 by 无际的草原
[考研] 材料学调剂 +10 提神豆沙包 2026-02-28 12/600 2026-03-02 09:26 by 李老师!
[考研] 279求调剂 +3 dua1 2026-03-01 4/200 2026-03-02 00:23 by 大脸蛋子
[基金申请] 成果系统访问量大,请一小时后再尝试。---NSFC啥时候好哦,已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 295复试调剂 +3 简木ChuFront 2026-03-01 3/150 2026-03-01 14:27 by zzxw520th
[考研] 寻找调剂 +4 LYidhsjabdj 2026-02-28 4/200 2026-03-01 10:56 by sunny81
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见