當(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)行擬合，但是代碼編寫有一些問題。。查閱相關(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)該就是沒有F的返回值了。。。我這個(gè)語句寫的應(yīng)該是肯定有問題但是我不太知道是什么問題。。希望能夠請(qǐng)高人指點(diǎn)一下謝謝�。。。�！

回復(fù)此樓

» 猜你喜歡

【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有225人回復(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)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有10人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)

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

如何用Matlab擬合大量散點(diǎn)的邊界曲線已經(jīng)有5人回復(fù)
在用nlinfit對(duì)非線性函數(shù)擬合時(shí)遇到的一些問題請(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微分、及非線性擬合問題已經(jīng)有5人回復(fù)
matlab數(shù)據(jù)擬合希望高手幫忙寫一下matlab的程序跪謝已經(jīng)有12人回復(fù)
matlab數(shù)據(jù)擬合問題求助已經(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ù)擬合問題已經(jīng)有7人回復(fù)
采用matlab中nlinfit擬合，出現(xiàn)問題了，不知道是哪錯(cuò)了，該怎么解決，求高手指點(diǎn)一下已經(jīng)有9人回復(fù)

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

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

查看全部 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軟件來做。

贊一下

回復(fù)此樓

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

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

已閱回復(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軟件來做。

因?yàn)橐孔�，�?000多這樣的方程需要擬合，（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小時(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ā)表

[有機(jī)交流] 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理 10+3	kaobao456 2026-03-29	4/200	2026-03-30 23:16 by nBu鋰
[考研] 286分調(diào)劑 +10	Faune 2026-03-30	11/550	2026-03-30 20:35 by 啊李999
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級(jí)。 +9	XX小鄧 2026-03-29	9/450	2026-03-30 17:18 by limeifeng
[考研] 324求調(diào)劑 +9	hanamiko 2026-03-26	11/550	2026-03-30 14:27 by JourneyLucky
[考研] 070300化學(xué)354求調(diào)劑 +7	101次希望 2026-03-28	7/350	2026-03-30 13:22 by Iveryant
[考研] 求調(diào)劑 +7	青春裁為三截 2026-03-29	7/350	2026-03-30 13:14 by laoshidan
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4	qianbai11 2026-03-29	4/200	2026-03-30 11:34 by 唐沐兒
[考研] 一志愿南航 335分 \| 0856 \| GPA 4.07 \| 有科研經(jīng)歷 +8	cccchenso 2026-03-29	8/400	2026-03-29 23:53 by 我是小康
[考研] 299求調(diào)劑 +10	15188958825 2026-03-25	10/500	2026-03-29 17:51 by 王亮_大連醫(yī)科大
[考研] 340求調(diào)劑 +6	Amber00 2026-03-26	6/300	2026-03-29 12:06 by 無際的草原
[考研] 321求調(diào)劑 +7	璞玉~~ 2026-03-25	8/400	2026-03-29 06:41 by 544594351
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +9	天天好運(yùn)來上岸?/a> 2026-03-24	10/500	2026-03-28 22:17 by chemzp
[考研] 求調(diào)劑 +6	蘆lty 2026-03-25	7/350	2026-03-28 13:13 by 唐沐兒
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 085405 考的11408求各位老師帶走 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 09:19 by 樂呵呵的追夢(mèng)人
[考研] 274求調(diào)劑 +17	顧九笙要謙虛 2026-03-24	23/1150	2026-03-27 15:16 by caszguilin
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚的貓
[考研] 290分調(diào)劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky