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

715211229

木蟲 (正式寫手)

應(yīng)助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時(shí)
蟲號(hào): 1357002
注冊(cè): 2011-07-27
性別: GG
專業(yè): 傳動(dòng)機(jī)械學(xué)

[求助] matlab擬合公式中未知系數(shù)

現(xiàn)在有2組數(shù)據(jù)，比較大，不方便貼出來，所以放附件了。
擬合給定的這2 個(gè)情況中f（x）中的系數(shù)ai。這個(gè)是我擬合的，效果不好，初值不知道該怎么調(diào)。

5、高速凸輪進(jìn)氣門上升段
f(x) = a6*(x^24)+a5*(x^22)+a4*(x^16)+a3*(x^12)+a2*(x^6)+a1*(x^2)+a0

   a0 =    5.675  (5.669, 5.68)
   a1 =    -9.946  (-9.978, -9.915)
   a2 =       6.59  (6.513, 6.667)
   a3 =    -3.32  (-3.421, -3.219)
   a4 =       1.42  (1.357, 1.482)
   a5 =    -0.2558  (-0.272, -0.2396)
   a6 =    0.08181  (0.07615, 0.08747)

6、高速凸輪進(jìn)氣門下降段  （效果不好）

   f(x) = a6*(x^26)+a5*(x^22)+a4*(x^20)+a3*(x^18)+a2*(x^5)+a1*(x^2)+a0

   a0 =    5.674  (5.625, 5.723)
   a1 =    -16.34  (-16.67, -16.01)
   a2 =    15.89  (15.25, 16.52)
   a3 =    -66.16  (-73.4, -58.92)
   a4 =       110  (97.04, 123)
   a5 =    -52.58  (-59.2, -45.96)
   a6 =    3.717  (3.193, 4.24)

[ Last edited by 715211229 on 2011-12-23 at 12:13 ]

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 數(shù)據(jù).rar

2011-12-23 12:12:11, 4.31 K

» 猜你喜歡

【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有300人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(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ù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)

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

在matlab中擬合曲面的表達(dá)式已經(jīng)有29人回復(fù)
matlab擬合擬均相動(dòng)力學(xué)參數(shù) 已經(jīng)有6人回復(fù)
需要用matlab擬合數(shù)據(jù)的外行人急需幫助已經(jīng)有10人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】matlab擬合二元線性方程已經(jīng)有3人回復(fù)
【求助】Matlab非線性最小二乘擬合活度系數(shù)模型（Willson、NRTL、UNIQUAC）已經(jīng)有17人回復(fù)
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復(fù)
【求助】matlab曲面擬合表達(dá)式已經(jīng)有10人回復(fù)
【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過一個(gè)點(diǎn) 已經(jīng)有8人回復(fù)
【求助】求MAtlab求解程序�。�！已經(jīng)有20人回復(fù)

我是蝸牛

1樓 2011-12-23 12:12:15

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

715211229

木蟲 (正式寫手)

應(yīng)助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時(shí)
蟲號(hào): 1357002
注冊(cè): 2011-07-27
性別: GG
專業(yè): 傳動(dòng)機(jī)械學(xué)

引用回帖:

2樓: Originally posted by dbb627 at 2011-12-23 15:09:47:
一定要選你的公式嗎
Linear model Poly9:
f(x) = p1*x^9 + p2*x^8 + p3*x^7 + p4*x^6 +
p5*x^5 + p6*x^4 + p7*x^3 + p8*x^2 + p9*x + p10
Coefficients (with 95% confidence boun ...

是的，求給定公式的系數(shù)

贊一下

回復(fù)此樓

我是蝸牛

4樓2011-12-23 20:03:48

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

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

一定要選你的公式嗎
Linear model Poly9:
   f(x) = p1*x^9 + p2*x^8 + p3*x^7 + p4*x^6 +
                  p5*x^5 + p6*x^4 + p7*x^3 + p8*x^2 + p9*x + p10
Coefficients (with 95% confidence bounds):
   p1 =    2.661  (-3.745, 9.066)
   p2 =    -40.35  (-80.31, -0.3958)
   p3 =    196.5  (91.87, 301.1)
   p4 =       -447  (-596.2, -297.8)
   p5 =       526  (400.1, 652)
   p6 =    -319.1  (-383, -255.1)
   p7 =    99.43  (80.45, 118.4)
   p8 =    -24.38  (-27.42, -21.34)
   p9 =    0.7499  (0.5288, 0.971)
   p10 =    5.671  (5.666, 5.677)

Goodness of fit:
  SSE: 0.003105
  R-square: 1
  Adjusted R-square: 1
  RMSE: 0.00413

Linear model Poly9:
   f(x) = p1*x^9 + p2*x^8 + p3*x^7 + p4*x^6 +
                  p5*x^5 + p6*x^4 + p7*x^3 + p8*x^2 + p9*x + p10
Coefficients (with 95% confidence bounds):
   p1 =    396.8  (308.6, 485)
   p2 =    -2180  (-2675, -1684)
   p3 =       4874  (3706, 6042)
   p4 =    -5600  (-7100, -4100)
   p5 =       3365  (2225, 4506)
   p6 =    -888.2  (-1410, -366.8)
   p7 =       18.5  (-120.8, 157.8)
   p8 =    9.332  (-10.71, 29.38)
   p9 =    -2.033  (-3.344, -0.721)
   p10 =    5.711  (5.683, 5.739)

Goodness of fit:
  SSE: 0.07685
  R-square: 0.9999
  Adjusted R-square: 0.9999
  RMSE: 0.02139

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2樓2011-12-23 15:09:47

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

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

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

感謝參與，應(yīng)助指數(shù) +1
715211229(金幣+1): ★有幫助嗯，是的，用多項(xiàng)式擬合可以擬合的很好，只是這個(gè)問題給定了方程，所以有點(diǎn)別扭 2011-12-23 20:05:44

就原數(shù)據(jù)和模型公式而言，除非改變公式，一樓的結(jié)果應(yīng)該就是最好的了。

贊一下

回復(fù)此樓

3樓2011-12-23 16:05:27

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

引用回帖:

4樓: Originally posted by 715211229 at 2011-12-23 20:03:48:
是的，求給定公式的系數(shù)

第一個(gè)不用在改了，第二個(gè)加了權(quán)重

CODE:

>>  x=x4;y=h4;

w=1./(y+0.00001);

opts=fitoptions('method','NonlinearLeastSquares','Weights',(w).^(1/8),'Lower',[-Inf -Inf -Inf -Inf -Inf  -Inf -Inf],'Upper',[Inf Inf Inf Inf Inf Inf Inf]);

opts.StartPoint= [ 3.7168  -52.5824  110.0003  -66.1606   15.8861  -16.3403    5.6742];

ft_ = fittype('a6*(x^26)+a5*(x^22)+a4*(x^20)+a3*(x^18)+a2*(x^5)+a1*(x^2)+a0','dependent',{'y'},...

    'independent',{'x'},'coefficients',{'a6', 'a5','a4','a3','a2','a1','a0'},'options',opts);

[cf_,good]= fit(x,y,ft_)

plot(x,cf_(x),'or-',x, y, '* ');

legend('擬合圖','原始數(shù)據(jù)')

cf_ =

   General model:
   cf_(x) = a6*(x^26)+a5*(x^22)+a4*(x^20)+a3*(x^18)+a2*(x^5)+a1*(x^2)+a0
   Coefficients (with 95% confidence bounds):
   a6 =    3.132  (2.697, 3.568)
   a5 =    -45.17  (-50.8, -39.54)
   a4 =    95.47  (84.3, 106.6)
   a3 =    -58.04  (-64.36, -51.72)
   a2 =       15.2  (14.57, 15.83)
   a1 =    -16.02  (-16.37, -15.67)
   a0 =    5.646  (5.59, 5.701)

good =

         sse: 4.3743
   rsquare: 0.9945
         dfe: 171
adjrsquare: 0.9943
      rmse: 0.1599

贊一下(2人)

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

5樓2011-12-24 11:02:12

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

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

[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級(jí)。 +4	XX小鄧 2026-03-29	4/200	2026-03-29 17:34 by 無際的草原
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5	hanamiko 2026-03-29	5/250	2026-03-29 16:46 by 學(xué)員8dgXkO
[考研] 291求調(diào)劑 +5	Y-cap 2026-03-29	6/300	2026-03-29 13:18 by mumin1990
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6	芝士雪baoo 2026-03-28	8/400	2026-03-29 08:16 by 松花缸1201
[考研] 346求調(diào)劑一志愿070303有機(jī)化學(xué) +3	蘿卜燉青菜 2026-03-28	3/150	2026-03-28 14:11 by 唐沐兒
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-27	6/300	2026-03-28 14:10 by 唐沐兒
[考研] 070300求調(diào)劑306分 +4	26要上岸 2026-03-27	4/200	2026-03-28 13:06 by 唐沐兒
[有機(jī)交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 305求調(diào)劑 +5	哇盧卡庫 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 286求調(diào)劑 +4	lim0922 2026-03-26	4/200	2026-03-27 10:28 by guoweigw
[考研] 求調(diào)劑 +6	林之夕 2026-03-24	6/300	2026-03-27 08:38 by hypershenger
[考研] 304材料求調(diào)劑 +4	鐘llll 2026-03-26	4/200	2026-03-27 03:42 by wxiongid
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚的貓
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 考研一志愿蘇州大學(xué)初始315（英一）求調(diào)劑 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu
[考研] 生物技術(shù)與工程 +3	1294608413 2026-03-25	4/200	2026-03-25 18:02 by 1294608413
[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 一志愿南航材料專317分求調(diào)劑 +5	炸呀炸呀炸薯?xiàng)l 2026-03-23	5/250	2026-03-24 16:52 by 星空星月
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky