| 24小時(shí)熱門版塊排行榜 |
| 7 | 1/1 | 返回列表 |
| 查看: 1858 | 回復(fù): 6 | ||
wljabgf銀蟲 (小有名氣)
|
[求助]
我用1stopt對(duì)一個(gè)函數(shù)中的2個(gè)參數(shù)進(jìn)行雙參數(shù)擬合,為何每次計(jì)算結(jié)果不同呢?求解答? 已有3人參與
|
|
各位大神們,我用1stopt 1.5版 對(duì)一個(gè)函數(shù)中的2個(gè)參數(shù)進(jìn)行雙參數(shù)擬合,為何每次計(jì)算結(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)的。 第一次計(jì)算結(jié)果: "Type your title here" 迭代數(shù): 17 計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:00:01:782 優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法 計(jì)算結(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)計(jì)(F-Statistic): 356.755375369784 參數(shù) 最佳估算 ---------- ------------- a 56.9724281749416 h 0.0308130218572859 ====== 結(jié)果輸出 ===== No 實(shí)測(cè)值f 計(jì)算值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 第二次計(jì)算結(jié)果: "Type your title here" 迭代數(shù): 23 計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:00:02:434 優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法 計(jì)算結(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)計(jì)(F-Statistic): 356.755375369253 參數(shù) 最佳估算 ---------- ------------- a 67.5875462999094 h 0.0259161681817663 ====== 結(jié)果輸出 ===== No 實(shí)測(cè)值f 計(jì)算值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 第三次計(jì)算結(jié)果: "Type your title here" 迭代數(shù): 16 計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:00:01:669 優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法 計(jì)算結(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)計(jì)(F-Statistic): 356.755375369213 參數(shù) 最佳估算 ---------- ------------- a 34.9984214789949 h 0.0506078021383304 ====== 結(jié)果輸出 ===== No 實(shí)測(cè)值f 計(jì)算值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 后面也算過(guò)很多次,但每一次結(jié)果都不同,請(qǐng)問(wèn)大神們,這個(gè)能不能通過(guò)什么算法設(shè)置或者約束條件,讓這兩個(gè)參數(shù)的結(jié)果能夠確定下來(lái)得到唯一固定的值啊,比如:可否設(shè)置另相關(guān)系數(shù)R^2最大來(lái)回歸得到這兩個(gè)參數(shù)的最有結(jié)果,跪求各位大神的指導(dǎo)啊,謝謝! |
鐵桿木蟲 (職業(yè)作家)
銀蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
|
樓主的公式就是過(guò)擬合: 公式中:(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)) 可以用一個(gè)參數(shù)來(lái)代替,比如h1,原公式等價(jià)變?yōu)椋?br /> f=(1-v)^1.5+(a-1)*v*h1 上面公式中的:(a-1)*v*h1就會(huì)導(dǎo)致過(guò)擬合,多解 |
銀蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
銀蟲 (小有名氣)
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +3 | 匹克i 2026-03-23 | 3/150 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
[考研] 招08考數(shù)學(xué) +6 | laoshidan 2026-03-20 | 14/700 |
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 生物學(xué)調(diào)劑招人!! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|