| 5 | 1/1 | 返回列表 |
| 查看: 3637 | 回復(fù): 8 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
cheng1378653新蟲 (初入文壇)
|
[求助]
用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,(如果可以請劃一下圖形) |
金蟲 (正式寫手)
鐵桿木蟲 (職業(yè)作家)
金蟲 (正式寫手)
|
還是2樓的結(jié)果好,用剛才我用1sopt運算的結(jié)果作為初始值,用MATLAB進行擬合得到結(jié)果和1樓一致,看來非線性最優(yōu)化還是1sopt強大啊 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 |
鐵桿木蟲 (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求化學(xué)調(diào)劑 +5 | wulanna 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求收留 +5 | 1943443204 2026-03-28 | 5/250 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
[考研] 一志愿華北電力大學(xué)能動專碩,293,求調(diào)劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風 2026-03-23 | 3/150 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當牛 2026-03-24 | 4/200 |
|