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

cheng1378653

應(yīng)助: 0 (幼兒園)
金幣: 85
帖子: 8
在線: 3.8小時(shí)
蟲(chóng)號(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ù)此樓

» 猜你喜歡

散金幣，求好運(yùn)，祝面上順利！已經(jīng)有80人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有211人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(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é)院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ù)

» 本主題相關(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ù)
如何用軟件畫(huà)出這個(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的回帖

cheng1378653

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 85
帖子: 8
在線: 3.8小時(shí)
蟲(chóng)號(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ù)，我忘改過(guò)來(lái)了，不好意思，不過(guò)和其它參數(shù)一樣也必須有一個(gè)合理的初始值才能得到好的結(jié)果。...

謝謝

回復(fù)此樓

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

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

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

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(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

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 35 (小學(xué)生)
金幣: 566.8
散金: 411
帖子: 371
在線: 176.9小時(shí)
蟲(chóng)號(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的回帖

lgycjpcqu

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 35 (小學(xué)生)
金幣: 566.8
散金: 411
帖子: 371
在線: 176.9小時(shí)
蟲(chóng)號(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的回帖

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

[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +3	星沉uu 2026-03-26	4/200	2026-03-26 16:24 by barness
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6	小小陳小小 2026-03-22	9/450	2026-03-26 15:58 by dick_runner
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +5	曼111 2026-03-24	6/300	2026-03-26 13:04 by 13756423260
[考研] 085600 材料與化工 329分求調(diào)劑 +9	Mr. Z 2026-03-25	9/450	2026-03-26 10:36 by baoball
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 求調(diào)劑323材料與化工 +4	1124361 2026-03-24	4/200	2026-03-25 11:19 by shulmg
[考研] 材料學(xué)碩，求調(diào)劑 6+4	糖葫蘆888ll 2026-03-22	9/450	2026-03-25 11:19 by greychen00
[考研] 考研化學(xué)308分求調(diào)劑 +10	你好明天你好 2026-03-23	11/550	2026-03-25 10:23 by userper
[考研] 0856材料專碩353求調(diào)劑 +6	NIFFFfff 2026-03-20	6/300	2026-03-24 21:02 by hello七七
[考研] 調(diào)劑 +4	13853210211 2026-03-24	4/200	2026-03-24 19:44 by ms629
[考研] 材料專碩找調(diào)劑 +5	哈哈哈吼吼吼哈 2026-03-23	5/250	2026-03-24 19:07 by 了了了了。。
[考研] 307求調(diào)劑 +3	余意卿 2026-03-21	6/300	2026-03-24 15:03 by 余意卿
[考研] 化學(xué)308分求調(diào)劑 +3	你好明天你好 2026-03-23	3/150	2026-03-23 20:11 by macy2011
[考研] 333求調(diào)劑 +3	ALULU4408 2026-03-23	3/150	2026-03-23 19:04 by macy2011
[考研] 298求調(diào)劑 +8	上岸6666@ 2026-03-20	8/400	2026-03-23 11:02 by laoshidan
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋?xiě)?/a>
[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky