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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

bloomi

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 76
帖子: 11
在線: 8.5小時
蟲號: 2974341
注冊: 2014-02-18
專業(yè): 遙感機(jī)理與方法

[求助] matlab非線性分段方程擬合急~謝謝~ 已有2人參與

初學(xué)者求幫忙！

如圖所示是一組非線性分段方程，其中待擬合參數(shù)有：T0，Ta，tm，δT，以及分段點ts。
我嘗試使用lsqcurvefit進(jìn)行擬合，但是代碼編寫有一些問題。。查閱相關(guān)資料也不太知道這種分段點就是未知參數(shù)的方程應(yīng)該如何編程，希望大家給予幫助

xdata = [ 5.5000 6.5000 7.5000 8.5000 9.5000 10.5000 11.5000 12.5000 13.5000 14.5000 15.5000 16.5000 17.5000 18.5000 19.5000 20.5000 21.5000 22.5000 23.5000 24.5000 1.5000 2.5000 3.5000 4.5000];
ydata = [299.0545  300.8932  304.5193  308.6715  311.6907  315.0223  315.7924  317.1786  315.0809  313.1319  312.8464  309.4712  305.4651  304.1893  302.5152  301.9055  303.4350  299.9959  299.2613  299.7677  299.4591  298.5113  298.7621  298.6650];
x0 = [299.05  40  14  0.5  17.4];

function F = fun(x,xdata)
if xdata > x(5)
w = 4/3*(x(3) - 5.4);
F = a + b*cos( 3.14/w.*( xdata -c ) );
elseif xdata <= x(5)
w = 4/3*(c - 5.4);
k = w/3.14 * ( atan( 3.14/w*(x(5)-x(3)) ) - d/b*asin( 3.14/w*( x(5)-x(3) ) ) );
F = a + d + ( b*cos( 3.14/w*(ee-c) ) - d ) * k/( k + xdata - ee );
end
end

[x,resnorm] = lsqcurvefit(@fun,x0,xdata,ydata);

這樣運(yùn)行之后會報錯：
Error using feval
Output argument "F" (and maybe others)
not assigned during call to
"F:\DTC\pros\fun1012.m>fun1012".

Error in lsqcurvefit (line 203)
         initVals.F =
         feval(funfcn_x_xdata{3},xCurrent,XDATA,varargin{:});

Caused by:
Failure in initial user-supplied
objective function evaluation.
LSQCURVEFIT cannot continue.

最后應(yīng)該就是沒有F的返回值了。。。我這個語句寫的應(yīng)該是肯定有問題但是我不太知道是什么問題。。希望能夠請高人指點一下謝謝�。。。�！

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有108人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有72人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

如何用Matlab擬合大量散點的邊界曲線已經(jīng)有5人回復(fù)
在用nlinfit對非線性函數(shù)擬合時遇到的一些問題請教已經(jīng)有5人回復(fù)
Matlab 二元線性函數(shù)擬合求助已經(jīng)有3人回復(fù)
不可線性化的非線性曲線方程擬合已經(jīng)有5人回復(fù)
求助：以下非線性擬合方程可以用1stopt求解么？用什么方法已經(jīng)有16人回復(fù)
求助用matlab 多元非線性擬合已經(jīng)有4人回復(fù)
求Origin或Matlab高手非線性擬合公式已經(jīng)有5人回復(fù)
matlab 擬合求助已經(jīng)有8人回復(fù)
關(guān)于matlab擬合已經(jīng)有4人回復(fù)
matlab非線性擬合求助已經(jīng)有7人回復(fù)
用matlab里的nlinfit如何擬合分段函數(shù) 已經(jīng)有8人回復(fù)
關(guān)于matlab微分、及非線性擬合問題已經(jīng)有5人回復(fù)
matlab數(shù)據(jù)擬合希望高手幫忙寫一下matlab的程序跪謝已經(jīng)有12人回復(fù)
matlab數(shù)據(jù)擬合問題求助已經(jīng)有7人回復(fù)
關(guān)于擬合：同一個方程，非線性擬合與轉(zhuǎn)換后線性擬合已經(jīng)有8人回復(fù)
等溫吸附線性非線性擬合結(jié)果不一致，希望大家多多指教已經(jīng)有15人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
用matlab擬合雙因變量實驗結(jié)果求助已經(jīng)有7人回復(fù)
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復(fù)
采用matlab中nlinfit擬合，出現(xiàn)問題了，不知道是哪錯了，該怎么解決，求高手指點一下已經(jīng)有9人回復(fù)

1樓 2014-10-12 22:47:18

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

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

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

感謝參與，應(yīng)助指數(shù) +1

推薦用1stopt軟件來做。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2014-10-13 09:56:31

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

bloomi

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 76
帖子: 11
在線: 8.5小時
蟲號: 2974341
注冊: 2014-02-18
專業(yè): 遙感機(jī)理與方法

引用回帖:

2樓: Originally posted by 月只藍(lán) at 2014-10-13 09:56:31
推薦用1stopt軟件來做。

因為要批量做，有2000多這樣的方程需要擬合，（xdata ydata值不一樣）所以希望寫到程序里去

贊一下

回復(fù)此樓

3樓2014-10-13 12:02:24

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

戴鋼盔的豬頭

木蟲 (知名作家)

應(yīng)助: 200 (高中生)
金幣: 5083.9
散金: 1310
紅花: 18
帖子: 6001
在線: 319.9小時
蟲號: 2854040
注冊: 2013-12-07
性別: GG
專業(yè): 波譜分析與成像分析

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
bloomi: 金幣+10, ★★★★★最佳答案 2014-10-15 10:09:56

分段函數(shù)最小二乘擬合目標(biāo)函數(shù)導(dǎo)數(shù)不連續(xù)，可采用兩段函數(shù)的加權(quán)加和作為擬合模型。加權(quán)因子假設(shè)分別為a和b，可設(shè)計連續(xù)光滑函數(shù)作為加權(quán)因子。這個函數(shù)需要表現(xiàn)出階越性，可采用某個三角函數(shù)。三角函數(shù)都是初等函數(shù)，連續(xù)可導(dǎo)，故可設(shè)計加權(quán)函數(shù)滿足：當(dāng)t小于ts時，a=1，b=0，否則a=0，b=1。具體哪個三角函數(shù)呢？樓主自己去研究下吧，授人以魚不如授人以漁。這個方法分兩段，分多段都可以用。

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

4樓2014-10-14 15:48:59

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

戴鋼盔的豬頭

木蟲 (知名作家)

應(yīng)助: 200 (高中生)
金幣: 5083.9
散金: 1310
紅花: 18
帖子: 6001
在線: 319.9小時
蟲號: 2854040
注冊: 2013-12-07
性別: GG
專業(yè): 波譜分析與成像分析

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

引用回帖:

4樓: Originally posted by 戴鋼盔的豬頭 at 2014-10-14 15:48:59
分段函數(shù)最小二乘擬合目標(biāo)函數(shù)導(dǎo)數(shù)不連續(xù)，可采用兩段函數(shù)的加權(quán)加和作為擬合模型。加權(quán)因子假設(shè)分別為a和b，可設(shè)計連續(xù)光滑函數(shù)作為加權(quán)因子。這個函數(shù)需要表現(xiàn)出階越性，可采用某個三角函數(shù)。三角函數(shù)都是初等函數(shù) ...

補(bǔ)充a和b不是同一個函數(shù)，但用到同一個三角函數(shù)。

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

5樓2014-10-14 15:50:21

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

bloomi

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 76
帖子: 11
在線: 8.5小時
蟲號: 2974341
注冊: 2014-02-18
專業(yè): 遙感機(jī)理與方法

引用回帖:

多謝您的幫助！我去研究下

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

6樓2014-10-15 10:09:39

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

bloomi

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 76
帖子: 11
在線: 8.5小時
蟲號: 2974341
注冊: 2014-02-18
專業(yè): 遙感機(jī)理與方法

引用回帖:

問個問題。。�！爱�(dāng)t小于ts時，a=1，b=0，否則a=0，b=1”
在這樣的條件下函數(shù)還能連續(xù)光滑嗎。。。

贊一下

回復(fù)此樓

7樓2014-10-28 15:53:28

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

戴鋼盔的豬頭

木蟲 (知名作家)

應(yīng)助: 200 (高中生)
金幣: 5083.9
散金: 1310
紅花: 18
帖子: 6001
在線: 319.9小時
蟲號: 2854040
注冊: 2013-12-07
性別: GG
專業(yè): 波譜分析與成像分析

引用回帖:

7樓: Originally posted by bloomi at 2014-10-28 15:53:28
問個問題。。�！爱�(dāng)t小于ts時，a=1，b=0，否則a=0，b=1”
在這樣的條件下函數(shù)還能連續(xù)光滑嗎。。。...

2/pi*arctan(w*x)，x是實驗參數(shù)（自變量），w是一個因子，越大階越性越強(qiáng)，可以自行給出，盡量縮短轉(zhuǎn)換帶。你自己去畫，再來說連不連續(xù)，可不可導(dǎo)。這個函數(shù)從-1階越到1（還是有轉(zhuǎn)換帶的，但可以人為使之縮短到遠(yuǎn)小于采樣間隔）。至于怎樣階越到0，自己去想，別再問我了。

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下(1人)

回復(fù)此樓

8樓2014-10-28 17:39:50

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

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

返回列表

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

[考研] 298求調(diào)劑 +6	上岸6666@ 2026-03-20	6/300	2026-03-22 20:21 by edmund7
[考研] 310求調(diào)劑 +4	baibai1314 2026-03-16	4/200	2026-03-22 20:19 by edmund7
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +3	jiajunX 2026-03-22	3/150	2026-03-22 19:32 by brblmd
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +3	niko- 2026-03-22	3/150	2026-03-22 15:15 by 楊楊楊紫
[考研] 070300化學(xué)求調(diào)劑 +3	苑豆豆 2026-03-20	3/150	2026-03-22 14:48 by ColorlessPI
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +3	酥酥魚.. 2026-03-21	3/150	2026-03-22 11:31 by 楊楊楊紫
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 求調(diào)劑 +3	白QF 2026-03-21	3/150	2026-03-21 13:12 by zhukairuo
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 二本跨考鄭大材料306英一數(shù)二 +3	z1z2z3879 2026-03-17	3/150	2026-03-21 02:29 by JourneyLucky
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 085600材料與化工 +8	安全上岸！ 2026-03-16	8/400	2026-03-20 22:13 by luoyongfeng
[考研] 08工學(xué)調(diào)劑 +5	用戶573181 2026-03-20	5/250	2026-03-20 15:47 by xia_2003
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong

24小時熱門版塊排行榜

[求助] matlab非線性分段方程擬合 急~謝謝~ 已有2人參與

» 猜你喜歡

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

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

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

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

[求助] matlab非線性分段方程擬合急~謝謝~ 已有2人參與