| 查看: 1363 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
Jounry-Jo银虫 (小有名气)
|
[求助]
急~~~~~~~求利用MATLAB或1STOPT非线性拟合两组数据 已有1人参与
|
|
| 1stopt5.0是试用版不能拟合三个参数以上的非线性拟合，matlab试了好久都没能出正确数据，不知是否编程有问题。着急得数据分析，希望有软件的大侠帮帮忙。不胜感激 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 拟合数据和公式.xlsx
2015-06-27 22:14:21, 14.02 K
» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有218人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求matlab大神指点，关于logistic拟合
已经有13人回复
- matlab 非线性拟合
已经有16人回复
- 求大神跑下1stOpt，一个sharemodel，谢谢！
已经有7人回复
- 求高手给做一道题，用MATLAB或1stOpt，请把结果直接给出来，我的软件就不行，谢谢
已经有5人回复
- 求助1stopt拟合动力学参数
已经有4人回复
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- 用1stOpt进行多模型多变量非线性拟合得不到结果
已经有4人回复
- 在进行多元非线性拟合的时候，用matlab拟合的值残差都大于0
已经有14人回复
- 求助用matlab 多元非线性拟合
已经有4人回复
- matlab 拟合求助
已经有8人回复
- 求助Matlab最小二乘法非线性拟合
已经有4人回复
- 紧急求助，利用Matlab对实验数据进行拟合求解参数。
已经有27人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- matlab数据拟合问题求助
已经有7人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
Jounry-Jo
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 976
- 散金: 1
- 帖子: 53
- 在线: 22.1小时
- 虫号: 2557609
- 注册: 2013-07-21
- 性别: MM
- 专业: 药物化学
4楼2015-06-28 20:18:52
2楼2015-06-28 12:16:22
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Jounry-Jo: 金币+7, ★★★很有帮助, 谢谢 2015-06-28 20:21:16
fegg7502: 金币+1, 鼓励交流 2015-06-29 08:36:56
感谢参与，应助指数 +1
Jounry-Jo: 金币+7, ★★★很有帮助, 谢谢 2015-06-28 20:21:16
fegg7502: 金币+1, 鼓励交流 2015-06-29 08:36:56
|
用1.5老版本的都行吧： 第一组： 均方差(RMSE):1.07205511441145 残差平方和(SSE):48.2706910701012 相关系数(R): 0.939291118069207 相关系数之平方(R^2): 0.8822678044837 确定系数(DC): 0.882147234885844 卡方系数(Chi-Square): 1.7441087621366 F统计(F-Statistic): 61.0931783656144 参数 最佳估算 ---------- ------------- a21 11.9497208389958 a31 17.4598494080349 a12 -6.19190121843902 a32 0.00659902374601305 a13 4.12126008280832 a23 -0.296670034331323 第二组： 均方差(RMSE):0.973759323354068 残差平方和(SSE):34.135459913483 相关系数(R): 0.97164462582794 相关系数之平方(R^2): 0.944093278900318 确定系数(DC): 0.944088417392726 卡方系数(Chi-Square): 1.25185088936219 F统计(F-Statistic): 107.312290588236 参数 最佳估算 ---------- ------------- a21 7.4558436152277 a31 13.1441376631054 a12 6.26895600341922 a32 3.27183196434597 a13 -3.59060566306702 a23 3.09901050411686 |
3楼2015-06-28 16:11:34
muzili2008
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1855.2
- 红花: 1
- 帖子: 189
- 在线: 138.5小时
- 虫号: 1089618
- 注册: 2010-09-04
- 专业: 化工系统工程
★
fegg7502: 金币+1, 鼓励交流 2015-07-07 09:13:30
fegg7502: 金币+1, 鼓励交流 2015-07-07 09:13:30
|
clc clear % http://muchong.com/bbs/viewthread.php?tid=9102854&fpage=4&target=blank data = xlsread('data.xlsx'); x1 = data(:,1); x2 = data(:,2); x3 = data(:,3); y = data(:,4); beta0 = ones(6,1)*1e-5; options = optimset('TolFun' , 1e-15, ... 'TolX' , 1e-15, ... 'MaxFunEvals', 1e5, ... 'MaxIter' , 1e4 ); [X,RESNORM,RESIDUAL,EXITFLAG,OUTPUT,LAMBDA,JACOBIAN]= ... lsqnonlin(@OptObjFunc,beta0,[],[],options,x1,x2,x3,y); y_err = OptObjFunc(X,x1,x2,x3,y); y_cal = y-y_err; plot(y_cal,y,'r*') axis([0 20 0 20]) % =================================== function ans = OptObjFunc(param,x1,x2,x3,y) % [a12,a21,a13,a31,a23,a32] = beta; m = size(x1,1); y_cal = zeros(m,1); for i = 1:m y_cal(i) = nrtl3comp(x1(i),x2(i),x3(i),param); end ans = y-y_cal; % ============================================ function y = nrtl3comp(x1,x2,x3,param) a12 = param(1); a21 = param(2); a13 = param(3); a31 = param(4); a23 = param(5); a32 = param(6); y=(a21*exp(-0.3*a21)*x2+a31*exp(-0.3*a31)*x3)/(x1+exp(-0.3*a21)*x2+exp(-0.3*a31)*x3) ... -(a21*exp(-0.3*a21)*x2+a31*exp(-0.3*a31)*x3)/(x1+exp(-0.3*a21)*x2+exp(-0.3*a31)*x3) ... *x1/(x1+exp(-0.3*a21)*x2+exp(-0.3*a31)*x3)+ ... a12-(a12*exp(-0.3*a12)*x1+a32*exp(-0.3*a32)*x3)/ ... (x2+exp(-0.3*a12)*x1+exp(-0.3*a32)*x3)*x2* ... exp(-0.3*a12)/(x2+exp(-0.3*a12)*x1+exp(-0.3*a32)*x3)+a13 ... -(a13*exp(-0.3*a13)*x1+a23*exp(-0.3*a23)*x2)/ ... (x3+exp(-0.3*a13)*x1+exp(-0.3*a23)*x2)*exp(-0.3*a13)*x3/ ... (x3+exp(-0.3*a13)*x1+exp(-0.3*a23)*x2) ; 5.4478；6.184；-3.5576；13.645；-14.644；2.5835 |
5楼2015-07-06 22:43:22
