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

bloomi

新蟲 (初入文壇)

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

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

初學(xué)者求幫忙！

如圖所示是一組非線性分段方程，其中待擬合參數(shù)有：T0，Ta，tm，δT，以及分段點(diǎn)ts。
我嘗試使用lsqcurvefit進(jìn)行擬合，但是代碼編寫有一些問(wèn)題。。查閱相關(guān)資料也不太知道這種分段點(diǎ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)行之后會(huì)報(bào)錯(cuò)：
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)該就是沒(méi)有F的返回值了。。。我這個(gè)語(yǔ)句寫的應(yīng)該是肯定有問(wèn)題但是我不太知道是什么問(wèn)題。。希望能夠請(qǐng)高人指點(diǎn)一下謝謝�。。。。�

回復(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ù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有74人回復(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ù)
天津城建大學(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ù)

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

如何用Matlab擬合大量散點(diǎn)的邊界曲線已經(jīng)有5人回復(fù)
在用nlinfit對(duì)非線性函數(shù)擬合時(shí)遇到的一些問(wèn)題請(qǐng)教已經(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微分、及非線性擬合問(wèn)題已經(jīng)有5人回復(fù)
matlab數(shù)據(jù)擬合希望高手幫忙寫一下matlab的程序跪謝已經(jīng)有12人回復(fù)
matlab數(shù)據(jù)擬合問(wèn)題求助已經(jīng)有7人回復(fù)
關(guān)于擬合：同一個(gè)方程，非線性擬合與轉(zhuǎn)換后線性擬合已經(jīng)有8人回復(fù)
等溫吸附線性非線性擬合結(jié)果不一致，希望大家多多指教已經(jīng)有15人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(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ù)

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

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

bloomi

新蟲 (初入文壇)

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

引用回帖:

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

因?yàn)橐孔�，�?000多這樣的方程需要擬合，（xdata ydata值不一樣）所以希望寫到程序里去

贊一下

回復(fù)此樓

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

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

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

月只藍(lán)

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

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

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

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

推薦用1stopt軟件來(lái)做。

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

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

戴鋼盔的豬頭

木蟲 (知名作家)

應(yīng)助: 200 (高中生)
金幣: 5083.9
散金: 1310
紅花: 18
帖子: 6001
在線: 319.9小時(shí)
蟲號(hào): 2854040
注冊(cè): 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è)計(jì)連續(xù)光滑函數(shù)作為加權(quán)因子。這個(gè)函數(shù)需要表現(xiàn)出階越性，可采用某個(gè)三角函數(shù)。三角函數(shù)都是初等函數(shù)，連續(xù)可導(dǎo)，故可設(shè)計(jì)加權(quán)函數(shù)滿足：當(dāng)t小于ts時(shí)，a=1，b=0，否則a=0，b=1。具體哪個(gè)三角函數(shù)呢？樓主自己去研究下吧，授人以魚不如授人以漁。這個(gè)方法分兩段，分多段都可以用。

[ 發(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小時(shí)
蟲號(hào): 2854040
注冊(cè): 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è)計(jì)連續(xù)光滑函數(shù)作為加權(quán)因子。這個(gè)函數(shù)需要表現(xiàn)出階越性，可采用某個(gè)三角函數(shù)。三角函數(shù)都是初等函數(shù) ...

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

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

贊一下

回復(fù)此樓

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

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

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

[考研] 322求調(diào)劑 +8	宋明欣 2026-03-27	8/400	2026-03-30 13:38 by chemdavid
[考研] 0703化學(xué) +14	妮妮ninicgb 2026-03-27	14/700	2026-03-30 12:46 by fangnagu
[考研] 材料與化工272求調(diào)劑 +21	阿斯蒂芬2004 2026-03-28	21/1050	2026-03-30 10:52 by 晴空210210
[考研] 290求調(diào)劑 +3	dfffsar 2026-03-29	3/150	2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +9	Koxui 2026-03-28	9/450	2026-03-29 20:43 by 無(wú)際的草原
[考研] 085600 材料與化工 329分求調(diào)劑 +14	Mr. Z 2026-03-25	14/700	2026-03-29 17:27 by 學(xué)員8dgXkO
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3	llss0711 2026-03-28	6/300	2026-03-29 10:26 by llss0711
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 085602 化工專碩 338分求調(diào)劑 +12	路癡小琪 2026-03-27	12/600	2026-03-28 15:41 by L135790
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[考研] 材料求調(diào)劑一志愿哈工大324 +7	閆旭東 2026-03-28	9/450	2026-03-28 08:51 by Xu de nuo
[考研] 286求調(diào)劑 +4	丟掉懶惰 2026-03-27	7/350	2026-03-28 08:07 by baoball
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[有機(jī)交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 08開(kāi)頭275求調(diào)劑 +4	拉誰(shuí)不重要 2026-03-26	4/200	2026-03-27 14:12 by Delta2012
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 309求調(diào)劑 +4	gajsj 2026-03-25	5/250	2026-03-26 00:27 by Dyhoer
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5	超贊的 2026-03-24	5/250	2026-03-24 21:07 by 星空星月
[考研] 269求調(diào)劑 +4	我想讀研11 2026-03-23	4/200	2026-03-23 21:25 by pswait