| 5 | 1/1 | 返回列表 |
| 查看: 1929 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
rhinoxixi金蟲(chóng) (正式寫(xiě)手)
|
[交流]
【求助】遞推公式的參數(shù)擬合 已有2人參與
|
||
|
自己建立的一個(gè)模型,是一個(gè)遞推公式: Yn=a+b*(Xn)^c+d*Yn-1 要擬合參數(shù),除了自己編程之外,有沒(méi)有什么軟件可以做到的? matlab行嗎?how? |
金蟲(chóng) (正式寫(xiě)手)
|
具體說(shuō)明如下 : 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í)我知道編程用最小二乘法可以求出來(lái),但是我不是很懂編程,所以問(wèn)有沒(méi)有軟件算。如果有人fortran很厲害,也可以幫我編個(gè)程序來(lái)擬合參數(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 ] |
鐵桿木蟲(chóng) (職業(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 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)354求調(diào)劑 +3 | 101次希望 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +6 | leaflight 2026-03-24 | 6/300 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑!。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|