 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 894 | 回復(fù): 6 | ||
hua99dushu鐵蟲 (小有名氣)
|
[求助]
求幫忙運行一下下面的1stopt的程序，微分方程-參數(shù)估計 已有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ù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有88人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團隊招收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)價值貼推薦，對您同樣有幫助:
- 求版主幫忙跑下面的1stopt的程序，參數(shù)估計。
已經(jīng)有6人回復(fù)
- 用1stopt軟件做微分方程組的參數(shù)估計
已經(jīng)有13人回復(fù)
- 怎么解這個復(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īng)有6人回復(fù)
- 用Matlab線性擬合求解兩個參數(shù)，十萬火急！
已經(jīng)有9人回復(fù)
- 1stopt中使用sharedmodel擬合曲線出現(xiàn)問題，求指示��！
已經(jīng)有6人回復(fù)
- 重金請教-如何用MAtlab對力學(xué)模型進行參數(shù)識別
已經(jīng)有12人回復(fù)
- matlab 擬合反應(yī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小時
- 蟲號: 2078718
- 注冊: 2012-10-22
- 性別: MM
- 專業(yè): 控制論中的數(shù)學(xué)方法
2樓2014-03-22 11:38:21
hua99dushu
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 74.1
- 帖子: 161
- 在線: 25.5小時
- 蟲號: 2078718
- 注冊: 2012-10-22
- 性別: MM
- 專業(yè): 控制論中的數(shù)學(xué)方法
3樓2014-03-22 11:39:41
hua99dushu
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 74.1
- 帖子: 161
- 在線: 25.5小時
- 蟲號: 2078718
- 注冊: 2012-10-22
- 性別: MM
- 專業(yè): 控制論中的數(shù)學(xué)方法
★
fegg7502: 金幣+1, 鼓勵交流 2014-03-24 08:33:10
fegg7502: 金幣+1, 鼓勵交流 2014-03-24 08:33:10
|
上面的錯誤挺多的，程序更正確定如下： 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]; 希望能夠得到估計的參數(shù)。。最好能有相關(guān)的圖 |
4樓2014-03-22 14:24:35
hua99dushu
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 74.1
- 帖子: 161
- 在線: 25.5小時
- 蟲號: 2078718
- 注冊: 2012-10-22
- 性別: MM
- 專業(yè): 控制論中的數(shù)學(xué)方法
|
請版主看此樓。。。。�？禳c幫忙運行一下吧。。 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è)作家)
- 計算強帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 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)計(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 計算 x 目標(biāo) y 計算 y 目標(biāo) z 計算 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小時
- 蟲號: 2078718
- 注冊: 2012-10-22
- 性別: MM
- 專業(yè): 控制論中的數(shù)學(xué)方法
7樓2014-03-24 20:49:49
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hua99dushu 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +9 | 阿斯蒂芬2004 2026-03-28 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-24 | 4/200 |
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +6 | 腳滑的守法公民 2026-03-27 | 6/300 |
|
|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
|
[考研]
 一志愿 西北大學(xué) 總分282 英語一62 求調(diào)劑
+7
|
18419759900 2026-03-25 | 8/400 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個好覺_16 2026-03-24 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5 | SYA！ 2026-03-23 | 5/250 |
|
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
|
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 8/400 |
|
|
|
[考研]
085602化學(xué)工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[有機交流]
 有機合成求助
 20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過 +6 | 幸運哩哩 2026-03-22 | 10/500 |
|