| 5 | 1/1 | 返回列表 |
| 查看: 1930 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
rhinoxixi金蟲 (正式寫手)
|
[交流]
【求助】遞推公式的參數(shù)擬合 已有2人參與
|
||
|
自己建立的一個(gè)模型,是一個(gè)遞推公式: Yn=a+b*(Xn)^c+d*Yn-1 要擬合參數(shù),除了自己編程之外,有沒(méi)有什么軟件可以做到的? matlab行嗎?how? |
|
用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 ] |
金蟲 (正式寫手)
|
具體說(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ù),謝謝啦! |
鐵桿木蟲 (職業(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ā)表 | |
|---|---|---|---|---|
|
[考研] 283求調(diào)劑(080500) +7 | A child 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +7 | 李李不服輸 2026-03-26 | 7/350 |
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +8 | 谷雨上岸 2026-03-23 | 9/450 |
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 0856,材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒(méi)煩惱 2026-03-23 | 4/200 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|