| 5 | 1/1 | 返回列表 |
| 查看: 1330 | 回復(fù): 8 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
wangdeqiang金蟲 (小有名氣)
|
[求助]
1stopt擬合中出現(xiàn)的問題 已有1人參與
|
|
|
Title "DIC-K"; Parameters kI,kII,b; constant n=3,g=8.1*10^4,u=0.26,k=2.17; Variable x,y,z; Function z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2))) -sum((kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b; Data; 376.78243 -112.14027 9.73371 374.92666 -111.4325 9.61182 373.12866 -110.7178 9.49729 371.38928 -109.99629 9.39441 369.70935 -109.26811 9.3051 368.08966 -108.5334 9.22794 366.53104 -107.79231 9.15739 365.03425 -107.04503 9.08779 363.60006 -106.29172 9.0149 本人想通過1stopt求解kI,kII,但是運行時系統(tǒng)提示參數(shù)kI數(shù)值范圍定義有誤! 請問該如何寫程序? |

金蟲 (小有名氣)

鐵桿木蟲 (職業(yè)作家)
|
公式: z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2))) -sum((kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b; 改為: z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2))) -sum(i=1:n)(kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b; 結(jié)果: 均方差(RMSE):0.0334465479380339 殘差平方和(SSE):0.0100680441207408 相關(guān)系數(shù)(R): 0.989298144852019 相關(guān)系數(shù)之平方(R^2): 0.978710819407646 確定系數(shù)(DC): 0.978710819407646 卡方系數(shù)(Chi-Square): 0.000543327950684942 F統(tǒng)計(F-Statistic): 140.916649515266 參數(shù) 最佳估算 ---------- ------------- ki 0.245068724509857 kii -5.08642699323206 b 9.4238457540924 |
金蟲 (小有名氣)

鐵桿木蟲 (職業(yè)作家)
|
公式寫成如下: Function z=sum(i=1:n, kI)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2))) -sum(i=1:n, KII)(kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b; 均方差(RMSE):0.000319000884196416 殘差平方和(SSE):9.15854077062857E-7 相關(guān)系數(shù)(R): 0.999999031699104 相關(guān)系數(shù)之平方(R^2): 0.999998063399146 確定系數(shù)(DC): 0.999998063399146 卡方系數(shù)(Chi-Square): 4.91879921856916E-8 F統(tǒng)計(F-Statistic): 0 參數(shù) 最佳估算 ---------- ------------- b 9.37438411982016 ki1 240.991105375969 ki2 45.2897205090565 ki3 -0.00759010072910402 kii1 1285.30644313843 kii2 -40.4189608273874 kii3 -4.56724263105769 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +4 | adaie 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 3/150 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +3 | sbdksD 2026-03-19 | 3/150 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊 2026-03-18 | 5/250 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +8 | zlingli 2026-03-13 | 8/400 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|