| 9 | 1/1 | 返回列表 |
| 查看: 1327 | 回復(fù): 8 | ||||
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 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +65 |
|
Title "DIC-K"; Parameter kI(3),kII(3),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(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; Data; |

專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +65 |
|
Title "DIC-K"; Parameter kI(3),kII(3),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(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; Data; |

專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +65 |

專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +65 |

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +18 | rare12345 2026-03-18 | 18/900 |
|
|
[考研] 生物學(xué)調(diào)劑招人。。 +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|