版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

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ù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有59人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(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ù)擬合，通過全局優(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ù)擬合問題已經(jīng)有7人回復(fù)
采用matlab中nlinfit擬合，出現(xià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)過一個(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的回帖

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
不是說1sopt得到的是唯一解嗎？

贊一下

回復(fù)此樓

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

已閱回復(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樓一致，看來非線性最優(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的回帖

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的回帖

lgycjpcqu

金蟲 (正式寫手)

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

引用回帖:

5樓: Originally posted by dingd at 2013-03-29 11:15:58
如果你用的1stOpt老版本（如1.5），可是如下代碼：

Variable t,c;
Function c=if(t<T1,A*(1-exp(-K1*t)),A*(1-exp(-k1*T1))*exp(-k2*(t-T1)));
Data;
0 0
0.0833 0.0833
0.25 1.2491
0.5 1.7025
1 0. ...

這次對(duì)了，看來是我的版本太老了

贊一下

回復(fù)此樓

6樓2013-03-29 13:53:09

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

cheng1378653

新蟲 (初入文壇)

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

引用回帖:

4樓: Originally posted by lgycjpcqu at 2013-03-29 10:47:34
還是2樓的結(jié)果好，用剛才我用1sopt運(yùn)算的結(jié)果作為初始值，用MATLAB進(jìn)行擬合得到結(jié)果和1樓一致，看來非線性最優(yōu)化還是1sopt強(qiáng)大啊
A =
-0.0367
k1 =
-14.2282
k2 =
3.7015
T =
0.3179
你需要MATLAB 代碼的 ...

你好，你的這個(gè)代碼是默認(rèn)t=0.5時(shí)取得最大值，而我是想要根據(jù)擬合的結(jié)果，求出取得最大值得T值，謝謝你的幫助

贊一下

回復(fù)此樓

7樓2013-03-30 22:01:05

已閱回復(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)助回帖

★
fegg7502: 金幣+1, 應(yīng)助指數(shù)+1, 鼓勵(lì)交流 2013-04-02 09:25:26

引用回帖:

7樓: Originally posted by cheng1378653 at 2013-03-30 22:01:05
你好，你的這個(gè)代碼是默認(rèn)t=0.5時(shí)取得最大值，而我是想要根據(jù)擬合的結(jié)果，求出取得最大值得T值，謝謝你的幫助...

fun=@(beta,t)((t<0.5).*fun1(beta,t)+(t>=0.5).*fun2(beta,t)); 這一句中的0.5應(yīng)該換成beta(4)，這個(gè)是一個(gè)參數(shù)，我忘改過來了，不好意思，不過和其它參數(shù)一樣也必須有一個(gè)合理的初始值才能得到好的結(jié)果。

贊一下

回復(fù)此樓

8樓2013-03-31 10:34:12

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

cheng1378653

新蟲 (初入文壇)

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

引用回帖:

8樓: Originally posted by lgycjpcqu at 2013-03-31 10:34:12
fun=@(beta,t)((t<0.5).*fun1(beta,t)+(t>=0.5).*fun2(beta,t)); 這一句中的0.5應(yīng)該換成beta(4)，這個(gè)是一個(gè)參數(shù)，我忘改過來了，不好意思，不過和其它參數(shù)一樣也必須有一個(gè)合理的初始值才能得到好的結(jié)果。...

謝謝

回復(fù)此樓

9樓2013-03-31 12:25:54

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 cheng1378653 的主題更新

返回列表

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

[考研] 311求調(diào)劑 +8	26研0 2026-03-15	9/450	2026-03-22 00:09 by BruceLiu320
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 材料學(xué)碩333求調(diào)劑 +3	北道巷 2026-03-18	3/150	2026-03-21 18:17 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +3	沿岸有貝殼6 2026-03-21	3/150	2026-03-21 10:35 by 暮云清寒
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 化學(xué)求調(diào)劑 +4	臨澤境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 生物學(xué)調(diào)劑招人�。。� +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海

24小時(shí)熱門版塊排行榜

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

» 猜你喜歡

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

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

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

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

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

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

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