當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

Jounry-Jo

銀蟲 (小有名氣)

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

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

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

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(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

» 猜你喜歡

【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有137人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有10人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)

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

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

muzili2008

木蟲 (小有名氣)

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

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

引用回帖:

4樓: Originally posted by Jounry-Jo at 2015-06-28 20:18:52
希望R2最好能9.5以上，最好的是四個(gè)參數(shù)值都是幾千，初值可以設(shè)500-600。
實(shí)在不行的話，另一種情況是四個(gè)參數(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的回帖

查看全部 6 個(gè)回答

wuhanfeng

新蟲 (小有名氣)

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

有初始值嗎？或者大概有個(gè)范圍？

贊一下

回復(fù)此樓

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

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

dingd

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

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

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

★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
Jounry-Jo: 金幣+7, ★★★很有幫助, 謝謝 2015-06-28 20:21:16
fegg7502: 金幣+1, 鼓勵(lì)交流 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)計(jì)(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)計(jì)(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小時(shí)
蟲號: 2557609
注冊: 2013-07-21
性別: MM
專業(yè): 藥物化學(xué)

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

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

[考研] 083000環(huán)境科學(xué)與工程調(diào)劑，總分281 +4	橙子（勝意） 2026-03-30	4/200	2026-03-31 00:44 by Linzejun
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +5	小張zxy張 2026-03-26	10/500	2026-03-31 00:29 by jp9609
[考研] 材料與化工調(diào)劑一志愿大連海事085600，349 +5	吃的不少 2026-03-30	5/250	2026-03-30 23:31 by 果果媽咪
[考研] 291求調(diào)劑 +8	Y-cap 2026-03-29	11/550	2026-03-30 21:03 by Y-cap
[考研] 環(huán)境工程 085701，267求調(diào)劑 +7	minht 2026-03-29	7/350	2026-03-30 20:46 by dophin1985
[考研] 化工專碩求調(diào)劑 +5	question挽風(fēng) 2026-03-24	5/250	2026-03-30 18:34 by tianyyysss
[考研] 292求調(diào)劑 +13	是妍子也是研子 2026-03-30	13/650	2026-03-30 18:01 by 小徐0109
[考研] 332求調(diào)劑 +6	@MZB382400 2026-03-28	6/300	2026-03-30 16:57 by 無際的草原
[考研] 262求調(diào)劑 +4	ZZ..000 2026-03-30	4/200	2026-03-30 15:57 by wangjy2002
[考研] 化學(xué)308分調(diào)劑 +7	你好明天你好 2026-03-30	8/400	2026-03-30 13:02 by 你好明天你好
[考研] 340求調(diào)劑 +6	Amber00 2026-03-26	6/300	2026-03-29 12:06 by 無際的草原
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5	0857求調(diào)劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +3	AZMK 2026-03-27	5/250	2026-03-28 16:19 by xxxsssccc
[考研] 283求調(diào)劑 +7	A child 2026-03-28	7/350	2026-03-28 12:05 by zllcz
[考研] 274求調(diào)劑 +17	顧九笙要謙虛 2026-03-24	23/1150	2026-03-27 15:16 by caszguilin
[考研] 調(diào)劑推薦 +5	清酒714 2026-03-26	6/300	2026-03-27 11:12 by 不吃魚的貓
[考研] 336材料求調(diào)劑 +7	陳瀅瑩 2026-03-26	9/450	2026-03-27 00:20 by wxiongid
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5	軒小寧—— 2026-03-26	5/250	2026-03-26 12:05 by i_cooler
[考研] 334分一志愿武理-080500 材料求調(diào)劑 +4	李李不服輸 2026-03-25	4/200	2026-03-25 21:26 by 星空星月
[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325