应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » matlab非线性拟合求助
51/1返回列表
查看: 2116  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

孤竹行者

金虫 (小有名气)

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

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

» 本帖附件资源列表

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

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

纳米

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
是你的终究是你的
1楼 2013-04-05 23:29:53
已阅   回复此楼   关注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;
有问题再多多交流！
赞一下(1人) 回复此楼
是你的终究是你的
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的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售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.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 03:55 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
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[教师之家] 跳槽后在研项目怎么办？ +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 淀粉搬运工
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
信息提示
关闭
请填处理意见
关闭