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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » matlab非线性拟合求助
81/1返回列表
查看: 1898  |  回复: 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的回帖

孤竹行者

金虫 (小有名气)
自己回复一下，千万别沉啊！
赞一下 回复此楼
是你的终究是你的
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的回帖

孤竹行者

金虫 (小有名气)
各位大侠，帮帮忙吧！小弟拜谢了！
赞一下 回复此楼
是你的终究是你的
5楼2013-04-15 12:35:30
已阅   回复此楼   关注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程序
赞一下(1人) 回复此楼
6楼2013-07-08 11:30:39
已阅   回复此楼   关注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的回帖

sujiguang

铁虫 (小有名气)
引用回帖:
7楼: Originally posted by 孤竹行者 at 2013-07-10 14:10:04
我这有一个自己编的程序，初步计算了一下wilson方程还算可以，你把里面的方程改成你的NRTL和UNIQUAC方程看看行不行！
function f = wilson1(Input)
a12=Input(1);
a13=Input(2);
a21=Input(3);
a23=Input(4); ...

，还是不会，哎，我再琢磨琢磨哈！为什么都是input未定义
赞一下 回复此楼
8楼2013-07-10 19:21:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 孤竹行者 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂 +10 月@163.com 2026-04-11 10/500 2026-04-12 09:14 by zhouyuwinner
[考研] 305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-12 07:35 by zhouxiaoyu
[考研] 化学工程调剂289 +44 yang婷 2026-04-07 50/2500 2026-04-12 02:36 by 秋豆菜芽
[考研] 305求调剂 +6 77Qi 2026-04-07 6/300 2026-04-12 02:30 by 秋豆菜芽
[考研] 11408。358求调剂 +3 TMYzds 2026-04-07 3/150 2026-04-11 17:10 by 氮气气气
[考研] 机械专硕270求调剂，接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 080500求调剂 +17 黄宇博 2026-04-06 17/850 2026-04-11 08:36 by zhq0425
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 282，电气工程专业，求调剂，不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[考研] 296求调剂 +6 汪！？！ 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
[考研] 生物学求调剂 一志愿沪9，326分 +7 刘墨墨 2026-04-06 7/350 2026-04-10 08:11 by kangsm
[考研] 求调剂希望还是希望在山河四省附近 +3 快乐的小白鸽 2026-04-05 3/150 2026-04-09 17:36 by wp06
[考研] 368化学求调剂 +13 wwwwabcde 2026-04-07 14/700 2026-04-09 14:47 by heaven_jay
[考研] 一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药 +7 超人不会飞@ 2026-04-08 7/350 2026-04-09 09:35 by gong120082
[考研] 316求调剂 +4 15318418673 2026-04-07 4/200 2026-04-07 22:12 by hemengdong
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 调剂 +4 mcbbc 2026-04-06 5/250 2026-04-07 12:33 by upczlm1989
[考研] 材料调剂 +5 小刘同学吖吖 2026-04-06 5/250 2026-04-06 18:34 by sherry_1901
信息提示
关闭
请填处理意见
关闭