| 9 | 1/1 | 返回列表 |
| 查看: 1331 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +4 | Mqqqqqq 2026-03-19 | 4/200 |
|
|
[考研] 274求調(diào)劑 +8 | S.H1 2026-03-18 | 8/400 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 8/400 |
|
|
[考研] 317求調(diào)劑 +4 | 申子申申 2026-03-19 | 8/400 |
|
|
[考研] 329求調(diào)劑 +6 | 想上學(xué)吖吖 2026-03-19 | 6/300 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|