Jounry-Jo

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 976
散金: 1
帖子: 53
在線: 22.1小時
蟲號: 2557609
注冊: 2013-07-21
性別: MM
專業(yè): 藥物化學(xué)

[求助] 急~~~~~~~求利用MATLAB或1STOPT非線性擬合兩組數(shù)據(jù) 已有1人參與

1stopt5.0是試用版不能擬合三個參數(shù)以上的非線性擬合，matlab試了好久都沒能出正確數(shù)據(jù)，不知是否編程有問題。著急得數(shù)據(jù)分析，希望有軟件的大俠幫幫忙。不勝感激

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 擬合數(shù)據(jù)和公式.xlsx

2015-06-27 22:14:21, 14.02 K

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有55人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

求matlab大神指點，關(guān)于logistic擬合已經(jīng)有13人回復(fù)
matlab 非線性擬合已經(jīng)有16人回復(fù)
求大神跑下1stOpt，一個sharemodel，謝謝！已經(jīng)有7人回復(fù)
求高手給做一道題，用MATLAB或1stOpt，請把結(jié)果直接給出來，我的軟件就不行，謝謝已經(jīng)有5人回復(fù)
求助1stopt擬合動力學(xué)參數(shù) 已經(jīng)有4人回復(fù)
求大神幫忙擬合一個非線性方程，求出模型參數(shù) 已經(jīng)有15人回復(fù)
用1stOpt進行多模型多變量非線性擬合得不到結(jié)果已經(jīng)有4人回復(fù)
在進行多元非線性擬合的時候，用matlab擬合的值殘差都大于0 已經(jīng)有14人回復(fù)
求助用matlab 多元非線性擬合已經(jīng)有4人回復(fù)
matlab 擬合求助已經(jīng)有8人回復(fù)
求助Matlab最小二乘法非線性擬合已經(jīng)有4人回復(fù)
緊急求助，利用Matlab對實驗數(shù)據(jù)進行擬合求解參數(shù)。已經(jīng)有27人回復(fù)
大家在使用1stOpt進行數(shù)據(jù)擬合的時候有沒有遇到過這種情況？已經(jīng)有8人回復(fù)
matlab數(shù)據(jù)擬合問題求助已經(jīng)有7人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)

1樓 2015-06-27 22:14:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wuhanfeng

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 113.6
紅花: 2
帖子: 64
在線: 94.9小時
蟲號: 1612956
注冊: 2012-02-12
專業(yè): 傳熱傳質(zhì)學(xué)

有初始值嗎？或者大概有個范圍？

贊一下

回復(fù)此樓

2樓2015-06-28 12:16:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +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
相關(guān)系數(shù)(R): 0.939291118069207
相關(guān)系數(shù)之平方(R^2): 0.8822678044837
確定系數(shù)(DC): 0.882147234885844
卡方系數(shù)(Chi-Square): 1.7441087621366
F統(tǒng)計(F-Statistic): 61.0931783656144

參數(shù) 最佳估算
---------- -------------
a21 11.9497208389958
a31 17.4598494080349
a12 -6.19190121843902
a32 0.00659902374601305
a13 4.12126008280832
a23 -0.296670034331323

第二組：
均方差(RMSE):0.973759323354068
殘差平方和(SSE):34.135459913483
相關(guān)系數(shù)(R): 0.97164462582794
相關(guān)系數(shù)之平方(R^2): 0.944093278900318
確定系數(shù)(DC): 0.944088417392726
卡方系數(shù)(Chi-Square): 1.25185088936219
F統(tǒng)計(F-Statistic): 107.312290588236

參數(shù) 最佳估算
---------- -------------
a21 7.4558436152277
a31 13.1441376631054
a12 6.26895600341922
a32 3.27183196434597
a13 -3.59060566306702
a23 3.09901050411686

贊一下

回復(fù)此樓

3樓2015-06-28 16:11:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Jounry-Jo

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 976
散金: 1
帖子: 53
在線: 22.1小時
蟲號: 2557609
注冊: 2013-07-21
性別: MM
專業(yè): 藥物化學(xué)

引用回帖:

2樓: Originally posted by wuhanfeng at 2015-06-28 12:16:22
有初始值嗎？或者大概有個范圍？

希望R2最好能9.5以上，最好的是四個參數(shù)值都是幾千，初值可以設(shè)500-600。
實在不行的話，另一種情況是四個參數(shù)值都很小。謝謝

贊一下

回復(fù)此樓

4樓2015-06-28 20:18:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

muzili2008

木蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1855.2
紅花: 1
帖子: 189
在線: 138.5小時
蟲號: 1089618
注冊: 2010-09-04
專業(yè): 化工系統(tǒng)工程

★
fegg7502: 金幣+1, 鼓勵交流 2015-07-07 09:13:30

引用回帖:

4樓: Originally posted by Jounry-Jo at 2015-06-28 20:18:52
希望R2最好能9.5以上，最好的是四個參數(shù)值都是幾千，初值可以設(shè)500-600。
實在不行的話，另一種情況是四個參數(shù)值都很小。謝謝...

clc
clear
% http://www.gaoyang168.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人)

回復(fù)此樓

5樓2015-07-06 22:43:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Jounry-Jo

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 976
散金: 1
帖子: 53
在線: 22.1小時
蟲號: 2557609
注冊: 2013-07-21
性別: MM
專業(yè): 藥物化學(xué)

引用回帖:

5樓: Originally posted by muzili2008 at 2015-07-06 22:43:22
clc
clear
% http://www.gaoyang168.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 = ...

謝謝，不過不是應(yīng)助回帖，無法給你金幣了

贊一下

回復(fù)此樓

6樓2015-07-11 21:19:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 生物學(xué)一志愿985，分數(shù)349求調(diào)劑 +4	zxts12 2026-03-21	4/200	2026-03-22 07:20 by ilovexiaobin
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +5	生物工程調(diào)劑 2026-03-17	9/450	2026-03-21 23:32 by zhujy1982
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7	墨墨漠 2026-03-20	7/350	2026-03-21 16:36 by barlinike
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 311求調(diào)劑 +5	冬十三 2026-03-18	5/250	2026-03-21 00:16 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 350求調(diào)劑 +5	weudhdk 2026-03-19	5/250	2026-03-20 22:04 by luoyongfeng
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11	adaie 2026-03-19	11/550	2026-03-20 21:30 by laoshidan
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 環(huán)境工程調(diào)劑 +9	大可digkids 2026-03-16	9/450	2026-03-20 17:38 by 醉在風(fēng)里
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678