| 6 | 1/1 | 返回列表 |
| 查看: 1921 | 回復(fù): 5 | |||
rhinoxixi金蟲 (正式寫手)
|
[交流]
【求助】遞推公式的參數(shù)擬合 已有2人參與
|
|
自己建立的一個(gè)模型,是一個(gè)遞推公式: Yn=a+b*(Xn)^c+d*Yn-1 要擬合參數(shù),除了自己編程之外,有沒有什么軟件可以做到的? matlab行嗎?how? |
金蟲 (正式寫手)
|
具體說明如下 : n Xn Yn 1 0.82900 0.12500 2 0.81600 0.14600 3 0.64300 0.24300 4 0.57000 0.24700 5 0.47200 0.26800 6 0.36500 0.28400 7 0.24500 0.30700 8 0.11800 0.27100 9 0.37200 0.30000 10 0.55800 0.32400 Y0=0 Y1=a+b*X1^c+d*Y0 Y2=a+b*X2^c+d*Y1 ....... 我要擬合出a、b、c、d 其實(shí)我知道編程用最小二乘法可以求出來,但是我不是很懂編程,所以問有沒有軟件算。如果有人fortran很厲害,也可以幫我編個(gè)程序來擬合參數(shù),謝謝啦! |
|
用MATLAB做的,beta為結(jié)果,不知是否符合要求,可能換別的初值會(huì)更好,樓主可以多試試。 function [beta r]=datafit(beta) close all hold on load data.txt X=data(:,2); Y=data(:,3); beta=nlinfit(X,Y,'fun',beta); Y1=fun(beta,X); r=corrcoef(Y,Y1); beta=beta'; r=r(1,2); function Y=fun(beta,X) a=beta(1); b=beta(2); c=beta(3); d=beta(4); n=length(X); Y(1)=a+b*X(1)^c; for i=2:n Y(i)=a+b*X(i)^c+d*Y(i-1); end Y=Y'; >> [beta r]=datafit([0 1 0 1]) Warning: Iteration limit exceeded. Returning results from final iteration. > In nlinfit at 220 In datafit at 7 beta = -4.4118 4.5129 0.0035 0.7410 r = 0.9689 >> [ Last edited by fspdlh on 2009-4-11 at 23:29 ] |
鐵桿木蟲 (職業(yè)作家)
|
用1stOpt處理,結(jié)果更好更穩(wěn)定,也不需要猜初值: 均方差(RMSE): 0.0154322799500792 殘差平方和(SSE): 0.00238155264457617 相關(guān)系數(shù)(R): 0.971104117667667 相關(guān)系數(shù)之平方(R^2): 0.943043207351097 參數(shù) 最佳估算 -------------------- ------------- a 0.100663466550055 b -2.12112269614125E-24 c -23.8930576866545 d 0.69639916154168 |
木蟲 (正式寫手)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +5 | yy要上岸呀 2026-03-17 | 5/250 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 312求調(diào)劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
[考研] 274求調(diào)劑0856材料化工 +13 | z2839474511 2026-03-11 | 14/700 |
|
|
[考研] 求調(diào)劑,總分315,考的生物醫(yī)藥,一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進(jìn)步 2026-03-11 | 5/250 |
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
[考研] 求調(diào)劑 資源與環(huán)境 285 +3 | 未名考生 2026-03-10 | 3/150 |
|