| 5 | 1/1 | 返回列表 |
| 查看: 1333 | 回復(fù): 8 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
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,但是運(yùn)行時(shí)系統(tǒng)提示參數(shù)kI數(shù)值范圍定義有誤! 請(qǐng)問該如何寫程序? |

金蟲 (小有名氣)

鐵桿木蟲 (職業(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)計(jì)(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)計(jì)(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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料080500調(diào)劑求收留 +7 | 一顆meteor 2026-03-13 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +4 | Charlieyq 2026-03-19 | 4/200 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[考研] 274求調(diào)劑 +8 | S.H1 2026-03-18 | 8/400 |
|
|
[考研] 329求調(diào)劑 +6 | 想上學(xué)吖吖 2026-03-19 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|