715211229

木蟲 (正式寫手)

應助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時
蟲號: 1357002
注冊: 2011-07-27
性別: GG
專業(yè): 傳動機械學

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

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

5、高速凸輪進氣門上升段
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、高速凸輪進氣門下降段  （效果不好）

   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 ]

回復此樓

» 本帖附件資源列表

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

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

» 猜你喜歡

調劑已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有95人回復
中德博士后2026項目5月30日截止- 國內博士后入站、德國兩年全職已經(jīng)有0人回復
面向全球申請人的中德博士后2026項目5月30日截止- 國內博后入站德國兩年全職已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。已經(jīng)有0人回復
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！已經(jīng)有1人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調劑，10人已經(jīng)有0人回復

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

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

我是蝸牛

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

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

dingd

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

計算強帖: 4
應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應助回帖

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

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

贊一下

回復此樓

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

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

查看全部 6 個回答

dbb627

榮譽版主 (著名寫手)

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

一定要選你的公式嗎
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

贊一下

回復此樓

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

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

715211229

木蟲 (正式寫手)

應助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時
蟲號: 1357002
注冊: 2011-07-27
性別: GG
專業(yè): 傳動機械學

引用回帖:

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ù)

贊一下

回復此樓

我是蝸牛

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

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

dbb627

榮譽版主 (著名寫手)

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

引用回帖:

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

第一個不用在改了，第二個加了權重

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人)

回復此樓

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

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

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

[考研] 壓國家一區(qū)線，求導師收留，有恩必謝！ +3	迷人的哈哈 2026-03-28	3/150	2026-03-28 07:34 by ACS Nano——
[考研] 083000學碩274求調劑 +6	Li李魚 2026-03-26	6/300	2026-03-28 07:25 by barnett0632
[考研] 085701環(huán)境工程，267求調劑 +13	minht 2026-03-26	13/650	2026-03-28 07:14 by mmm just
[考研] 070300化學求調劑 +4	起個名咋這么難 2026-03-27	4/200	2026-03-27 21:39 by 83503孫老師
[考研] 考研化學308分求調劑 +10	你好明天你好 2026-03-23	12/600	2026-03-27 14:43 by shangxh
[考研] 085600，材料與化工321分，求調劑 +9	大饞小子 2026-03-27	9/450	2026-03-27 14:30 by mmm just
[考研] 求調劑 +6	林之夕 2026-03-24	6/300	2026-03-27 08:38 by hypershenger
[考研] 081200-11408-276學碩求調劑 +4	崔wj 2026-03-26	4/200	2026-03-27 08:04 by chemisry
[考研] 321求調劑 +5	材料cailiao 2026-03-21	5/250	2026-03-26 20:41 by fmesaito
[考研] 324求調劑 +4	wysyjs25 2026-03-21	4/200	2026-03-26 20:38 by fmesaito
[考研] 286求調劑 +13	Faune 2026-03-21	13/650	2026-03-26 19:52 by peike
[考研] 生物學 296 求調劑 +4	朵朵- 2026-03-26	6/300	2026-03-26 19:01 by 不吃魚的貓
[考研] 340求調劑 +3	Amber00 2026-03-26	3/150	2026-03-26 18:57 by 不吃魚的貓
[考研] 一志愿南京郵電大學 288分材料考研求調劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 一志愿中南大學化學學碩0703總分337求調劑 +7	niko- 2026-03-22	7/350	2026-03-25 20:14 by qingfeng258
[考研] 生物技術與工程 +3	1294608413 2026-03-25	4/200	2026-03-25 18:02 by 1294608413
[考研] 296求調劑 +4	汪�。�！ 2026-03-25	7/350	2026-03-25 16:41 by 汪�。�！
[考研] 340求調劑 +5	話梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin
[考研] 一志愿南航材料專317分求調劑 +5	炸呀炸呀炸薯條 2026-03-23	5/250	2026-03-24 16:52 by 星空星月
[考研] 308求調劑 +3	墨墨漠 2026-03-21	3/150	2026-03-22 16:54 by i_cooler