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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 急~~~~~~~求利用MATLAB或1STOPT非线性拟合两组数据
61/1返回列表
查看: 1341  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Jounry-Jo

银虫 (小有名气)

[求助] 急~~~~~~~求利用MATLAB或1STOPT非线性拟合两组数据 已有1人参与

1stopt5.0是试用版不能拟合三个参数以上的非线性拟合，matlab试了好久都没能出正确数据，不知是否编程有问题。着急得数据分析，希望有软件的大侠帮帮忙。不胜感激
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 拟合数据和公式.xlsx
    • 2015-06-27 22:14:21, 14.02 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-06-27 22:14:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuhanfeng

新虫 (小有名气)
有初始值吗？或者大概有个范围？
赞一下 回复此楼
2楼2015-06-28 12:16:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

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

Jounry-Jo

银虫 (小有名气)
引用回帖:
2楼: Originally posted by wuhanfeng at 2015-06-28 12:16:22
有初始值吗？或者大概有个范围？

希望R2最好能9.5以上，最好的是四个参数值都是几千，初值可以设500-600。
实在不行的话，另一种情况是四个参数值都很小。谢谢
赞一下 回复此楼
4楼2015-06-28 20:18:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

muzili2008

木虫 (小有名气)
★
fegg7502: 金币+1, 鼓励交流 2015-07-07 09:13:30
引用回帖:
4楼: Originally posted by Jounry-Jo at 2015-06-28 20:18:52
希望R2最好能9.5以上，最好的是四个参数值都是几千，初值可以设500-600。
实在不行的话，另一种情况是四个参数值都很小。谢谢...

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
赞一下(1人) 回复此楼
5楼2015-07-06 22:43:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jounry-Jo

银虫 (小有名气)
引用回帖:
5楼: Originally posted by muzili2008 at 2015-07-06 22:43:22
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 =  ...

谢谢，不过不是应助回帖，无法给你金币了
赞一下 回复此楼
6楼2015-07-11 21:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Jounry-Jo 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博26年 +4 八6八68 2026-03-19 4/200 2026-03-24 15:49 by 小Ben呵呵
[考研] 305分求调剂（食品工程） +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 284求调剂 +10 Zhao anqi 2026-03-22 10/500 2026-03-24 00:08 by Equinoxhua
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 一志愿上海交大生物与医药专硕324分，求调剂 +5 jiajunX 2026-03-22 5/250 2026-03-23 18:07 by YMU施老师
[考研] 263求调剂 +6 yqdszhdap－ 2026-03-22 9/450 2026-03-23 12:57 by yqdszhdap－
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
信息提示
关闭
请填处理意见
关闭