| 7 | 1/1 | 返回列表 |
| 查看: 1857 | 回復(fù): 6 | ||
wljabgf銀蟲 (小有名氣)
|
[求助]
我用1stopt對一個函數(shù)中的2個參數(shù)進(jìn)行雙參數(shù)擬合,為何每次計算結(jié)果不同呢?求解答? 已有3人參與
|
|
各位大神們,我用1stopt 1.5版 對一個函數(shù)中的2個參數(shù)進(jìn)行雙參數(shù)擬合,為何每次計算結(jié)果不同呢?求解答?代碼如下: Title "Type your title here"; //Parameters ; //Variable ; //Function ; Parameter a=0 [0,100], h=0 [0,0.1]; Variable v,f; Function f=(1-v)^1.5+(a-1)*v*(6.0881e-7*((34.5+h)^3-34.5^3)+0.0000170824*((25.25+h)^3-25.25^3)+0.0001727927*((14.25+h)^3-14.25^3)+0.0004147025*((7.125+h)^3-7.125^3)+0.0011128873*((3.555+h)^3-3.555^3)); Data; 0 1 0.1 0.934033162 0.2 0.808121285 0.3 0.720435573 0.4 0.649447618 0.5 0.588298163 0.6 0.467070538 選擇的算法是:1. 麥夸特法:Levenberg-Marquardt(LM),其中參數(shù)都是默認(rèn)的。 第一次計算結(jié)果: "Type your title here" 迭代數(shù): 17 計算用時(時:分:秒:微秒): 00:00:01:782 優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法 計算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn) 均方差(RMSE): 0.0236564369215793 殘差平方和(SSE): 0.00391738905477264 相關(guān)系數(shù)(R): 0.993065206977526 相關(guān)系數(shù)之平方(R^2): 0.986178505309316 決定系數(shù)(DC): 0.9818600058285 卡方系數(shù)(Chi-Square): 0.00316216817997876 F統(tǒng)計(F-Statistic): 356.755375369784 參數(shù) 最佳估算 ---------- ------------- a 56.9724281749416 h 0.0308130218572859 ====== 結(jié)果輸出 ===== No 實(shí)測值f 計算值f 1 1 1.0000000 2 0.934033162 0.8963063 3 0.808121285 0.8005244 4 0.720435573 0.7131360 5 0.649447618 0.6347233 6 0.588298163 0.5660100 7 0.467070538 0.5079301 第二次計算結(jié)果: "Type your title here" 迭代數(shù): 23 計算用時(時:分:秒:微秒): 00:00:02:434 優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法 計算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn) 均方差(RMSE): 0.0236564369215793 殘差平方和(SSE): 0.00391738905477264 相關(guān)系數(shù)(R): 0.993065206977515 相關(guān)系數(shù)之平方(R^2): 0.986178505309295 決定系數(shù)(DC): 0.9818600058285 卡方系數(shù)(Chi-Square): 0.00316216817999158 F統(tǒng)計(F-Statistic): 356.755375369253 參數(shù) 最佳估算 ---------- ------------- a 67.5875462999094 h 0.0259161681817663 ====== 結(jié)果輸出 ===== No 實(shí)測值f 計算值f 1 1 1.0000000 2 0.934033162 0.8963063 3 0.808121285 0.8005244 4 0.720435573 0.7131360 5 0.649447618 0.6347233 6 0.588298163 0.5660100 7 0.467070538 0.5079301 第三次計算結(jié)果: "Type your title here" 迭代數(shù): 16 計算用時(時:分:秒:微秒): 00:00:01:669 優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法 計算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn) 均方差(RMSE): 0.0236564369215794 殘差平方和(SSE): 0.00391738905477264 相關(guān)系數(shù)(R): 0.993065206977515 相關(guān)系數(shù)之平方(R^2): 0.986178505309294 決定系數(shù)(DC): 0.9818600058285 卡方系數(shù)(Chi-Square): 0.00316216817999256 F統(tǒng)計(F-Statistic): 356.755375369213 參數(shù) 最佳估算 ---------- ------------- a 34.9984214789949 h 0.0506078021383304 ====== 結(jié)果輸出 ===== No 實(shí)測值f 計算值f 1 1 1.0000000 2 0.934033162 0.8963063 3 0.808121285 0.8005244 4 0.720435573 0.7131360 5 0.649447618 0.6347233 6 0.588298163 0.5660100 7 0.467070538 0.5079301 后面也算過很多次,但每一次結(jié)果都不同,請問大神們,這個能不能通過什么算法設(shè)置或者約束條件,讓這兩個參數(shù)的結(jié)果能夠確定下來得到唯一固定的值啊,比如:可否設(shè)置另相關(guān)系數(shù)R^2最大來回歸得到這兩個參數(shù)的最有結(jié)果,跪求各位大神的指導(dǎo)啊,謝謝! |
鐵桿木蟲 (職業(yè)作家)
銀蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
|
樓主的公式就是過擬合: 公式中:(6.0881e-7*((34.5+h)^3-34.5^3)+0.0000170824*((25.25+h)^3-25.25^3)+0.0001727927*((14.25+h)^3-14.25^3)+0.0004147025*((7.125+h)^3-7.125^3)+0.0011128873*((3.555+h)^3-3.555^3)) 可以用一個參數(shù)來代替,比如h1,原公式等價變?yōu)椋?br /> f=(1-v)^1.5+(a-1)*v*h1 上面公式中的:(a-1)*v*h1就會導(dǎo)致過擬合,多解 |
銀蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
銀蟲 (小有名氣)
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +6 | vv迷 2026-03-22 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 0854電子信息求調(diào)劑 324 +3 | Promise-jyl 2026-03-23 | 3/150 |
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 初始318分求調(diào)劑(有工作經(jīng)驗(yàn)) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|