 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 893 | 回復(fù): 6 | ||
hua99dushu鐵蟲 (小有名氣)
|
[求助]
求幫忙運(yùn)行一下下面的1stopt的程序，微分方程-參數(shù)估計(jì) 已有1人參與
|
|
Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1],c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,1],C=[0,1],E=[0,1],N=[0,1]; Sharedmodel; Variable t, x,y,z; Function x'=a1*x*(y/M-1)-a2*y+a3*z; y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E); z'=c1*x*(x/N-1)-c2*y-c3*z; RowData： x=[1.0000,1.098,1.226,1.426,1.709,2.102,2.442,2.906,3.438,3.796,4.532,5.189]; y=[1.0000,1.017,1.081,1.219,1.470,1.838,2.064,2.247,2.397,2.449,2.742,3.209]; z=[1.0000,1.106,1.240,1.455,1.754,2.174,2.542,3.042,3.622,4.028,4.843,5.741]; 最好能給出誤差圖和相圖呢。。 |
回復(fù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有206人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 求版主幫忙跑下面的1stopt的程序，參數(shù)估計(jì)。
已經(jīng)有6人回復(fù)
- 用1stopt軟件做微分方程組的參數(shù)估計(jì)
已經(jīng)有13人回復(fù)
- 怎么解這個(gè)復(fù)數(shù)方程
已經(jīng)有3人回復(fù)
- 求用高版本的1stOpt擬合方程
已經(jīng)有8人回復(fù)
- 求大神用matlab給我處理一組數(shù)據(jù)，有人沒有�。�
已經(jīng)有15人回復(fù)
- matlab如何求解含有bessel方程中的參數(shù)
已經(jīng)有5人回復(fù)
- 1stopt擬合問題求助
已經(jīng)有7人回復(fù)
- 1stopt 擬合計(jì)算
已經(jīng)有6人回復(fù)
- 用Matlab線性擬合求解兩個(gè)參數(shù)，十萬火急！
已經(jīng)有9人回復(fù)
- 1stopt中使用sharedmodel擬合曲線出現(xiàn)問題，求指示�。�
已經(jīng)有6人回復(fù)
- 重金請教-如何用MAtlab對力學(xué)模型進(jìn)行參數(shù)識別
已經(jīng)有12人回復(fù)
- matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- 跪求1stOpt2.5以上版本，或者請高人幫忙擬合包含復(fù)數(shù)的函數(shù)的參數(shù)
已經(jīng)有8人回復(fù)
- 需要用matlab擬合數(shù)據(jù)的外行人急需幫助
已經(jīng)有10人回復(fù)
hua99dushu
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 74.1
- 帖子: 161
- 在線: 25.5小時(shí)
- 蟲號: 2078718
- 注冊: 2012-10-22
- 性別: MM
- 專業(yè): 控制論中的數(shù)學(xué)方法
2樓2014-03-22 11:38:21
hua99dushu
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 74.1
- 帖子: 161
- 在線: 25.5小時(shí)
- 蟲號: 2078718
- 注冊: 2012-10-22
- 性別: MM
- 專業(yè): 控制論中的數(shù)學(xué)方法
3樓2014-03-22 11:39:41
hua99dushu
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 74.1
- 帖子: 161
- 在線: 25.5小時(shí)
- 蟲號: 2078718
- 注冊: 2012-10-22
- 性別: MM
- 專業(yè): 控制論中的數(shù)學(xué)方法
★
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-03-24 08:33:10
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-03-24 08:33:10
|
上面的錯(cuò)誤挺多的，程序更正確定如下： Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1], c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,1],C=[0,1],E=[0,1],N=[0,1]; Sharedmodel; Variable t,x,y,z; Function x'=a1*x*(y/M-1)-a2*y+a3*z; y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E); z'=c1*x*(x/N-1)-c2*y-c3*z; RowData： t=[0,1,2,3,4,5,6,7,8,9,10,11]; x=[1.0000,1.098,1.226,1.426,1.709,2.102,2.442,2.906,3.438,3.796,4.532,5.189]; y=[1.0000,1.017,1.081,1.219,1.470,1.838,2.064,2.247,2.397,2.449,2.742,3.209]; z=[1.0000,1.106,1.240,1.455,1.754,2.174,2.542,3.042,3.622,4.028,4.843,5.741]; 希望能夠得到估計(jì)的參數(shù)。。最好能有相關(guān)的圖 |
4樓2014-03-22 14:24:35
hua99dushu
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 74.1
- 帖子: 161
- 在線: 25.5小時(shí)
- 蟲號: 2078718
- 注冊: 2012-10-22
- 性別: MM
- 專業(yè): 控制論中的數(shù)學(xué)方法
|
請版主看此樓。。。。�？禳c(diǎn)幫忙運(yùn)行一下吧。。 Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1], c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,5],C=[0,10],E=[0,10],N=[0,5]; Sharedmodel; Variable t,x,y,z; Function x'=a1*x*(y/M-1)-a2*y+a3*z; y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E); z'=c1*x*(x/N-1)-c2*y-c3*z; RowData： t=[0,1,2,3,4,5,6,7,8,9,10,11]; x=[1.0000,1.0777,1.1443,1.3497,1.5881,1.9814,2.1905,2.5799,2.9763,3.2445,3.5765,3.9315]; y=[1.0000,1.0597,1.0982,1.2035,1.3566,1.4610,1.4771,1.4862,1.6052,1.6094,1.8892,2.0334]; z=[1.0000,1.0920,1.2033,1.4031,1.6920,1.9382,2.2159,2.6186,2.9489,3.1536,3.5979,4.0233]; |
5樓2014-03-23 20:18:44
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應(yīng)助回帖
★ ★
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+2, 應(yīng)助指數(shù)+1, 3ks 2014-03-24 08:33:25
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+2, 應(yīng)助指數(shù)+1, 3ks 2014-03-24 08:33:25
|
只大概算了下，供參考： 均方差(RMSE):0.227821606236476 殘差平方和(SSE):1.71278858084954 相關(guān)系數(shù)(R): 0.968347715467601 相關(guān)系數(shù)之平方(R^2): 0.937697298051321 確定系數(shù)(DC): 0.810472079088207 F統(tǒng)計(jì)(F-Statistic): -6.59531364009876 參數(shù) 最佳估算 -------------------- ------------- a1 3.45198672383651E-11 a2 5.83733964128734E-8 a3 0.138472304076658 b1 1.44029340345734E-6 b2 0.0518147236689287 b3 0.0406931861755383 c1 3.70541984860404E-6 c2 1.37133289578244E-7 c3 8.59791920071737E-10 m 3.16205888404511 c 2.08529659127593 e 3.70210463181138 n 6.20898311572885E-5 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) x 計(jì)算 x 目標(biāo) y 計(jì)算 y 目標(biāo) z 計(jì)算 z 1 1 1.0777 1.14300503145253 1.0597 1.05729662499709 1.092 1.06851434515071 2 2 1.1443 1.29685305954089 1.0982 1.11588828762999 1.2033 1.1572932896264 3 3 1.3497 1.46465905303486 1.2035 1.17578665477138 1.4031 1.2709880836245 4 4 1.5881 1.65027065718425 1.3566 1.23690130994621 1.692 1.41560520068449 5 5 1.9814 1.85849139764935 1.461 1.29894832240681 1.9382 1.59905260863551 6 6 2.1905 2.09539536339523 1.4771 1.36129294495917 2.2159 1.8319321926032 7 7 2.5799 2.36877609710511 1.4862 1.42267433005819 2.6186 2.12871429124705 8 8 2.9763 2.68879681981474 1.6052 1.48071370354723 2.9489 2.50951876009061 9 9 3.2445 3.06894995110712 1.6094 1.53101337771447 3.1536 3.00288297139041 10 10 3.5765 3.52750422345532 1.8892 1.56545537275185 3.5979 3.65018117220115 11 11 3.9315 4.08974278597234 2.0334 1.56886778417041 4.0233 4.51289544111373 |
6樓2014-03-23 21:52:53
hua99dushu
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 74.1
- 帖子: 161
- 在線: 25.5小時(shí)
- 蟲號: 2078718
- 注冊: 2012-10-22
- 性別: MM
- 專業(yè): 控制論中的數(shù)學(xué)方法
7樓2014-03-24 20:49:49
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hua99dushu 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +5 | 芝士雪baoo 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)求調(diào)劑 +9 | 奶油草莓. 2026-03-22 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢 2026-03-23 | 16/800 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 總分293求調(diào)劑 +6 | 加一一九 2026-03-25 | 8/400 |
|
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +9 | Mr. Z 2026-03-25 | 9/450 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +5 | plum李子 2026-03-21 | 8/400 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+4
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|