| 6 | 1/1 | 返回列表 |
| 查看: 1920 | 回復(fù): 5 | |||
rhinoxixi金蟲 (正式寫手)
|
[交流]
【求助】遞推公式的參數(shù)擬合 已有2人參與
|
|
自己建立的一個(gè)模型,是一個(gè)遞推公式: Yn=a+b*(Xn)^c+d*Yn-1 要擬合參數(shù),除了自己編程之外,有沒有什么軟件可以做到的? matlab行嗎?how? |
金蟲 (正式寫手)
|
具體說(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),但是我不是很懂編程,所以問有沒有軟件算。如果有人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 ] |
鐵桿木蟲 (職業(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ā)表 | |
|---|---|---|---|---|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 070300化學(xué)319求調(diào)劑 +4 | 錦鯉0909 2026-03-17 | 4/200 |
|
|
[考研] 材料專碩306英一數(shù)二 +7 | z1z2z3879 2026-03-16 | 9/450 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +7 | 邱gl 2026-03-11 | 8/400 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +5 | ZY-05 2026-03-13 | 6/300 |
|
|
[考研] 285 求調(diào)劑 資源與環(huán)境 一志愿北京化工大學(xué) +3 | 未名考生 2026-03-10 | 3/150 |
|
|
[考研] 材料371求調(diào)劑 +9 | 鱷魚? 2026-03-11 | 11/550 |
|
|
[考研] 一志愿中科院,化學(xué)方向,295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
[考研] 求調(diào)劑 資源與環(huán)境 285 +3 | 未名考生 2026-03-10 | 3/150 |
|