 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 741 | 回復: 4 | |||
[交流]
forcal 代碼運行后出錯
|
|
運行了別人的Forcal代碼，出現(xiàn)了“Forcal表達式運行錯誤！ 實數(shù)表達式名稱：main 所在模塊：-1 實數(shù)函數(shù)名：fcopt::SimOpt 錯誤代碼：3”為什么？能不能幫忙調(diào)試下，指出問題所在��！謝謝~~ 代碼：!using["sys","IMSL","fcopt"]; f(t,HCp,dHCp : R : k1,a,E1, Acid,T)= { R=8.3144, dHCp=-k1*Acid^a*exp[-E1/(R*T)]*HCp }; 優(yōu)化(_k1, _a, _E1 : i, j, s, t1, t2, HCp, IDO : k1,a,E1, Acid,T, hf, tArray, max, TempArray,AcidArray,maxHCp)= { k1=_k1, a=_a, E1=_E1, s=0, j=0, (j IDO=1, t1=tArray.A(0,0), HCp=tArray.A(0,j+1), i=1, (i IVPRK(&IDO, hf, &t1, t2, 1e-6, 0, &HCp), s=s+[HCp-tArray.A(i,j+1)]^2, i++ }, IVPRK(3, hf, &t1, t2+0.1, 1e-6, 0, &HCp), j++ }, s }; main(:i, k1, a, E1, min : hf, tArray, max, TempArray,AcidArray,maxHCp)= { ClearImslErr(), //清空IMSL的錯誤輸出 ERSET(0, 0, 0), //關閉IMSL錯誤輸出 hf=HFor("f"  ,max=10, maxHCp=8,//實驗數(shù)據(jù)組數(shù) tArray=new{rtoi(real_s),rtoi(max),rtoi(maxHCp+1)}.free().io::arrayns{ //存放實驗數(shù)據(jù) " 0 14.277 12.892 13.880 12.859 13.526 15.058 14.705 14.042 5 14.057 12.721 12.986 12.576 13.138 11.872 14.530 12.709 10 13.418 11.721 11.350 11.616 12.309 9.015 13.831 12.133 20 11.691 9.383 10.122 9.069 10.231 7.119 12.748 9.872 30 9.220 7.622 7.795 7.321 7.870 5.223 10.619 5.949 40 7.875 6.287 6.781 5.880 6.358 4.451 8.678 4.741 60 6.554 4.687 6.101 4.076 5.310 3.162 8.043 3.527 80 5.776 4.946 4.416 3.784 2.915 2.284 5.270 3.196 100 4.646 3.586 3.434 2.716 2.502 2.070 4.072 2.489 120 2.998 3.101 1.985 2.458 2.370 1.933 3.866 2.989 " }, TempArray=new{rtoi(real_s),rtoi(maxHCp)}.free().io::arrayns{"413.15 413.15 423.15 423.15 433.15 433.15 443.15 443.15"}, AcidArray=new{rtoi(real_s),rtoi(maxHCp)}.free().io::arrayns{"0.7 0.9 0.5 0.7 0.3 0.5 0.1 0.3"}, min=1.1, k1=0.3, a=1.2, E1=1e5, //提供一組合適的初值 i=SimOpt[HFor("優(yōu)化" , optstep,0, optexpand,1.5, optcontract,0.5, optmax,5000 : &k1, &a, &E1, &min], //求n維極值的單形調(diào)優(yōu)法printff{"\r\n實際迭代次數(shù)={1,i}, k1={2,r}, a={3,r}, E1={4,r}, 目標函數(shù)終值={5,r}\r\n",i, k1, a, E1, min}, ERSET(0, 2, 2), ERSET(0, 1, 0) //打開IMSL錯誤輸出 }; [ Last edited by ahbiao on 2011-3-27 at 21:41 ] |
回復此樓» 猜你喜歡
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有53人回復
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復
-
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復
-
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！
已經(jīng)有1人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復
-
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8
已經(jīng)有0人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求助用MS計算聲子譜出錯的問題
已經(jīng)有15人回復
- fortran轉(zhuǎn)換CHGCAR出錯
已經(jīng)有4人回復
- VASP并行計算出錯
已經(jīng)有9人回復
- Gaussian 計算拉曼光譜最后時刻出錯
已經(jīng)有14人回復
- 急！關于nsfc申請代碼的問題
已經(jīng)有6人回復
- 得不到純產(chǎn)物，哪出錯了，親們幫忙看看啊
已經(jīng)有3人回復
- 高斯09中的L502出錯的問題
已經(jīng)有5人回復
- 關于CASTEP 中計算時出錯
已經(jīng)有6人回復
- 出錯 warning#2 #180
已經(jīng)有7人回復
- 選錯申請代碼了，影響大嗎
已經(jīng)有29人回復
- MS一打開就出錯，任何任務都不能運行
已經(jīng)有12人回復
- 博士后基金申請在打印轉(zhuǎn)化成word時出錯
已經(jīng)有13人回復
» 搶金幣啦！回帖就可以得到:
-
青島大學功能紡織品與先進材料研究院2026碩士研究生招生調(diào)劑
+1/91
-
2026江西科技師范大學申亮組招碩士研究生
+1/88
-
南昌航空大學材料學院教授/博導團隊招收調(diào)劑研究生
+1/86
-
中山大學附屬醫(yī)院實驗和生信科研助理招聘
+1/77
-
新疆碳納米材料制備與應用重點實驗室2026年招聘啟事
+1/45
-
導師招生
+1/43
-
南京信息工程大學 化材院接受調(diào)劑研究生
+1/40
-
天津理工大學國家杰青王鐵團隊招聘
+1/29
-
上海大學 “生物有機電子材料及器件”團隊碩士、博士研究生招聘
+1/28
-
濟南大學水利與環(huán)境學院趙艷俠教授簡介
+1/21
-
西南科技大學李永升老師接收材料和化工類研究生調(diào)劑
+1/18
-
化工畢業(yè)三年轉(zhuǎn)電子信息可行嗎？？各位同行前輩
+1/11
-
寧夏大學化學化工學院魏逸彬研究員課題組招聘啟事
+1/7
-
22408一志愿211，329分，本科績優(yōu)，計算機科班
+1/6
-
2026年中國科學院城市環(huán)境所與福建農(nóng)林大學聯(lián)合培養(yǎng)碩士研究生項目介紹
+1/6
-
半導體光電子器件與集成
+1/3
-
環(huán)境微生物與污染控制研究所招收2026碩士研究生
+1/3
-
【博士招生】天津理工大學國家杰青王鐵課題組招收2026年博士研究生
+1/2
-
湖南農(nóng)業(yè)大學 環(huán)境科學與工程0830+資源與環(huán)境專碩0857 調(diào)劑名額
+1/2
-
【光學工程研究生調(diào)劑】科研資源豐厚，歡迎優(yōu)秀學子加入
+1/1
ahbiao(金幣+10): 謝謝~~ 2011-03-30 08:10:16
|
本帖內(nèi)容被屏蔽 |
2樓2011-03-29 10:10:20
3樓2011-03-30 08:12:07
ahbiao(金幣+40): 2011-03-31 08:34:11
|
本帖內(nèi)容被屏蔽 |
4樓2011-03-30 15:03:23
恩，謝謝！ 5樓2011-03-31 08:34:35
相關版塊跳轉(zhuǎn)
我要訂閱樓主 ahbiao 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 南京大學化學調(diào)劑 +9 | 景隨風 2026-03-29 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 283求調(diào)劑（080500） +7 | A child 2026-03-27 | 7/350 |
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+11
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
|
[考研] 296求調(diào)劑 +6 | 彼岸t 2026-03-29 | 6/300 |
|
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
|
[考研] 0856材料化工調(diào)劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +4 | wei-5 2026-03-26 | 4/200 |
|
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
|
[考研] 藥學105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 086502化學工程342求調(diào)劑 +6 | 阿姨復古不過 2026-03-27 | 6/300 |
|
|
|
[考研] 0703化學一志愿南京師范大學303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
|
[考研] 0703化學求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研] 生物學 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿山東大學藥學學碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|