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

cheng1378653

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 85
帖子: 8
在線: 3.8小時(shí)
蟲號(hào): 2246356
注冊(cè): 2013-01-17
專業(yè): 計(jì)算機(jī)軟件

[求助] 用matlab里的nlinfit如何擬合分段函數(shù)

Sample Text
t=[0 0.0833 0.25 0.5 1 2];
c=[0 0.0833 1.2491 1.7025 0.2623 0.0757];
函數(shù)形式為c=A*(1-exp(-K1*t)) 0<=t<=T
c=A*(1-exp(-k1*T))*exp(-k2*(t-T)) t>=T
如何用腳本式或是函數(shù)式M文件擬合，求出A，k1,k2,(如果可以請(qǐng)劃一下圖形)

回復(fù)此樓

» 猜你喜歡

面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有166人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(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é)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有2人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有7人回復(fù)

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

對(duì)兩個(gè)自變量，一個(gè)因變量進(jìn)行非線性擬合回歸，求出目標(biāo)函數(shù)的兩個(gè)參數(shù)。已經(jīng)有5人回復(fù)
原創(chuàng)：用MATLAB做的GUI做數(shù)據(jù)擬合，通過(guò)全局優(yōu)化算法（含原文件，code，教程）已經(jīng)有55人回復(fù)
matlab非線性擬合計(jì)算求助已經(jīng)有10人回復(fù)
matlab中關(guān)于非線性擬合，用nlinfit.程序總出錯(cuò)。謝謝各位大俠已經(jīng)有3人回復(fù)
用MATLAB或者ORIGIN怎么擬合復(fù)數(shù)數(shù)據(jù) 已經(jīng)有10人回復(fù)
matlab擬合求助已經(jīng)有11人回復(fù)
matlab多元非線性擬合已經(jīng)有9人回復(fù)
如何用軟件畫出這個(gè)方程的圖形已經(jīng)有13人回復(fù)
MATLAB用于三元回歸參數(shù) 已經(jīng)有3人回復(fù)
誠(chéng)懇求助如何做IC50擬合曲線已經(jīng)有7人回復(fù)
用matlab擬合雙因變量實(shí)驗(yàn)結(jié)果求助已經(jīng)有7人回復(fù)
matlab非線性參數(shù)擬合問(wèn)題已經(jīng)有7人回復(fù)
采用matlab中nlinfit擬合，出現(xiàn)問(wèn)題了，不知道是哪錯(cuò)了，該怎么解決，求高手指點(diǎn)一下已經(jīng)有9人回復(fù)
matlab的簡(jiǎn)單函數(shù)擬合已經(jīng)有6人回復(fù)
求助matlab隱函數(shù)參數(shù)擬合已經(jīng)有16人回復(fù)
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復(fù)
【求助】求助matlab數(shù)值擬合已經(jīng)有8人回復(fù)
如何用MATLAB擬合指數(shù)函數(shù) 已經(jīng)有3人回復(fù)
【求助】matlab函數(shù)擬合，擬合結(jié)果并保存，方便調(diào)用已經(jīng)有4人回復(fù)
【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過(guò)一個(gè)點(diǎn) 已經(jīng)有8人回復(fù)
【求助】使用MATLAB怎么實(shí)現(xiàn)擬合力場(chǎng)參數(shù)的程序？【已解決】已經(jīng)有15人回復(fù)
Matlab 的自定義函數(shù)擬合已經(jīng)有3人回復(fù)
【其他】自定義的二元函數(shù)擬合已經(jīng)有6人回復(fù)

1樓 2013-03-28 22:01:09

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

lgycjpcqu

金蟲 (正式寫手)

應(yīng)助: 35 (小學(xué)生)
金幣: 566.8
散金: 411
帖子: 371
在線: 176.9小時(shí)
蟲號(hào): 1327328
注冊(cè): 2011-06-20
性別: GG
專業(yè): 電化學(xué)

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

★ ★ ★ ★
cheng1378653: 回帖置頂 2013-03-29 18:09:29
cheng1378653: 金幣+2, ★有幫助 2013-03-30 22:02:11
fegg7502: 金幣+2, 應(yīng)助指數(shù)+1, 鼓勵(lì)交流 2013-04-02 09:24:40

還是2樓的結(jié)果好，用剛才我用1sopt運(yùn)算的結(jié)果作為初始值，用MATLAB進(jìn)行擬合得到結(jié)果和1樓一致，看來(lái)非線性最優(yōu)化還是1sopt強(qiáng)大啊
A =
-0.0367
k1 =
-14.2282
k2 =
3.7015
T =
0.3179
你需要MATLAB 代碼的話如下
clear
tic
fun1=inline('beta(1).*(1-exp(-beta(2).*t))','beta','t');
fun2=inline('beta(1).*(1-exp(-beta(2).*beta(4)))*exp(-beta(3)*(t-beta(4)))','beta','t');
fun=@(beta,t)((t<0.5).*fun1(beta,t)+(t>=0.5).*fun2(beta,t));
t=[0 0.0833 0.25 0.5 1 2];
c=[0 0.0833 1.2491 1.7025 0.2623 0.0757];
beta=nlinfit(t,c,fun,[-0.05 -6.8 3.7 0.5]);
A=beta(1);
k1=beta(2);
k2=beta(3);
T=beta(4);
warning off all
c1=beta(1).*(1-exp(-beta(2).*t(1:4)));
c2=beta(1).*(1-exp(-beta(2).*beta(4))).*exp(-beta(3).*(t(4:6)-beta(4)));
cfit(1:4)=c1;
cfit(4:6)=c2;
plot(t,c,'*')
hold on;
plot(t,cfit,'r-')
toc

贊一下(2人)

回復(fù)此樓

4樓2013-03-29 10:47:34

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

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

dingd

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

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

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

★
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+1, 鼓勵(lì)交流 2013-04-02 09:24:49

1stOpt更簡(jiǎn)單：

CODE:

Function c=if(t
Data;

t=[0 0.0833 0.25 0.5 1 2];

c=[0 0.0833 1.2491 1.7025 0.2623 0.0757];

均方差(RMSE): 0.0282872879644411
殘差平方和(SSE): 0.00480102396229929
相關(guān)系數(shù)(R): 0.999311354215736
相關(guān)系數(shù)之平方(R^2): 0.998623182664489
決定系數(shù)(DC): 0.998184285936899
卡方系數(shù)(Chi-Square): 0.0580621164226486
F統(tǒng)計(jì)(F-Statistic): 367.165006987895

參數(shù) 最佳估算
---------- -------------
t1 0.317859946002189
a -0.0366738585259589
k1 -14.2282073650289
k2 3.70145325521581

jv1.jpg

贊一下

回復(fù)此樓

2樓2013-03-29 08:50:59

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

lgycjpcqu

金蟲 (正式寫手)

應(yīng)助: 35 (小學(xué)生)
金幣: 566.8
散金: 411
帖子: 371
在線: 176.9小時(shí)
蟲號(hào): 1327328
注冊(cè): 2011-06-20
性別: GG
專業(yè): 電化學(xué)

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

★
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+1, 鼓勵(lì)交流 2013-04-02 09:24:59

引用回帖:

2樓: Originally posted by dingd at 2013-03-29 08:50:59
1stOpt更簡(jiǎn)單：

Function c=if(t<T1,A*(1-exp(-K1*t)),A*(1-exp(-k1*T1))*exp(-k2*(t-T1)));
Data;
t=;
c=;

均方差(RMSE): 0.0282872879644411
殘差平方和(SSE): 0.00480102396229929
相關(guān)系數(shù)(R): ...

為啥我的1sopt運(yùn)算結(jié)果不一致呢？
參數(shù) 最佳估算
---------- -------------
t1          0.500000049987681
a       -0.0568673833778535
k1          -6.86274581292278
k2          3.70079073587775
不是說(shuō)1sopt得到的是唯一解嗎？

贊一下

回復(fù)此樓

3樓2013-03-29 10:38:11

已閱回復(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í)
蟲號(hào): 291104
注冊(cè): 2006-10-28

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

★
fegg7502: 金幣+1, 應(yīng)助指數(shù)+1, 3ks 2013-04-02 09:25:13

引用回帖:

3樓: Originally posted by lgycjpcqu at 2013-03-29 10:38:11
為啥我的1sopt運(yùn)算結(jié)果不一致呢？
參數(shù) 最佳估算
---------- -------------
t1          0.500000049987681
a       -0.0568673833778535
k1          -6.86274581292278
k2          3.70079073587775
...

如果你用的1stOpt老版本（如1.5），可是如下代碼：

CODE:

Variable t,c;

Function c=if(t
Data;

0        0

0.0833        0.0833

0.25        1.2491

0.5        1.7025

1        0.2623

2        0.0757

贊一下

回復(fù)此樓

5樓2013-03-29 11:15:58

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

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

[考研] 085600材料與化工調(diào)劑 +8	kikiki7 2026-03-30	8/400	2026-03-30 20:24 by tianyyysss
[考研] 材料科學(xué)與工程 317求調(diào)劑 +7	JKSOIID 2026-03-26	7/350	2026-03-30 20:10 by 滴滴上岸呀
[考研] 266分，求材料相關(guān)專業(yè)調(diào)劑 +5	哇呼哼呼哼 2026-03-30	6/300	2026-03-30 19:40 by 2026材料調(diào)劑
[考研] 311求調(diào)劑 +7	藍(lán)月亮亮 2026-03-30	7/350	2026-03-30 19:24 by 無(wú)際的草原
[考研] 0703化學(xué)321分求調(diào)劑 +10	三dd. 2026-03-30	11/550	2026-03-30 19:24 by markhwc
[碩博家園] 求調(diào)劑有機(jī)化學(xué)考研356分 +10	Nadiums 2026-03-25	11/550	2026-03-30 11:56 by yjolah
[考研] 285求調(diào)劑 +5	AZMK 2026-03-29	8/400	2026-03-30 10:41 by AZMK
[考研] 材料科學(xué)與工程求調(diào)劑 +6	深V宿舍吧 2026-03-29	6/300	2026-03-29 20:51 by 唐沐兒
[考研] 求收留 +5	1943443204 2026-03-28	5/250	2026-03-29 12:11 by 無(wú)際的草原
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4	小張zxy張 2026-03-26	8/400	2026-03-28 19:25 by lbsjt
[考研] 347求調(diào)劑 +3	山頂見α 2026-03-25	3/150	2026-03-28 14:13 by 唐沐兒
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當(dāng)當(dāng)
[考研] 315調(diào)劑 +4	0860求調(diào)劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002
[考研] 292求調(diào)劑 +4	求求了收下我吧?/a> 2026-03-26	4/200	2026-03-27 10:37 by zhshch
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002
[考研] 286求調(diào)劑 +4	lim0922 2026-03-26	4/200	2026-03-27 10:28 by guoweigw
[考研] 材料調(diào)劑 5+4	想要一壺桃花水 2026-03-25	10/500	2026-03-26 19:56 by 不吃魚的貓
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機(jī) 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4	w320357296 2026-03-25	6/300	2026-03-25 17:59 by 255671
[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325