 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 740 | 回復(fù): 4 | |||
[交流]
forcal 代碼運(yùn)行后出錯
|
|
運(yùn)行了別人的Forcal代碼，出現(xiàn)了“Forcal表達(dá)式運(yùn)行錯誤！ 實(shí)數(shù)表達(dá)式名稱：main 所在模塊：-1 實(shí)數(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), //關(guān)閉IMSL錯誤輸出 hf=HFor("f"  ,max=10, maxHCp=8,//實(shí)驗(yàn)數(shù)據(jù)組數(shù) tArray=new{rtoi(real_s),rtoi(max),rtoi(maxHCp+1)}.free().io::arrayns{ //存放實(shí)驗(yàn)數(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í)際迭代次數(shù)={1,i}, k1={2,r}, a={3,r}, E1={4,r}, 目標(biāo)函數(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 ] |
回復(fù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有128人回復(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ù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊招收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)價值貼推薦，對您同樣有幫助:
- 求助用MS計算聲子譜出錯的問題
已經(jīng)有15人回復(fù)
- fortran轉(zhuǎn)換CHGCAR出錯
已經(jīng)有4人回復(fù)
- VASP并行計算出錯
已經(jīng)有9人回復(fù)
- Gaussian 計算拉曼光譜最后時刻出錯
已經(jīng)有14人回復(fù)
- 急！關(guān)于nsfc申請代碼的問題
已經(jīng)有6人回復(fù)
- 得不到純產(chǎn)物，哪出錯了，親們幫忙看看啊
已經(jīng)有3人回復(fù)
- 高斯09中的L502出錯的問題
已經(jīng)有5人回復(fù)
- 關(guān)于CASTEP 中計算時出錯
已經(jīng)有6人回復(fù)
- 出錯 warning#2 #180
已經(jīng)有7人回復(fù)
- 選錯申請代碼了，影響大嗎
已經(jīng)有29人回復(fù)
- MS一打開就出錯，任何任務(wù)都不能運(yùn)行
已經(jīng)有12人回復(fù)
- 博士后基金申請在打印轉(zhuǎn)化成word時出錯
已經(jīng)有13人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
下一站，欽州
+1/80
-
講師身份申請面上項目中標(biāo)率有影響嗎
+1/74
-
青島農(nóng)業(yè)大學(xué)海洋學(xué)院招調(diào)劑碩士研究生
+1/73
-
北京-89175-事業(yè)單位-誠征女友
+1/63
-
歡迎報考贛南師范大學(xué)智能制造與未來能源學(xué)院
+2/62
-
歡迎化學(xué)、生物與醫(yī)藥、藥學(xué)，生物、食品等相關(guān)專業(yè)的同學(xué)
+1/39
-
歡迎化學(xué)、化工、材料、環(huán)境、藥學(xué)、納米、能源等專業(yè)優(yōu)秀同學(xué)加入！名額充足
+1/39
-
西華大學(xué)材料學(xué)院表面科學(xué)與工程技術(shù)科研團(tuán)隊2026年招收研究生
+1/33
-
哈爾濱工業(yè)大學(xué)（深圳）-何自開教授團(tuán)隊誠招化學(xué)工程與技術(shù)專業(yè)博士研究生
+1/30
-
濟(jì)南大學(xué)水利與環(huán)境學(xué)院趙艷俠教授簡介
+1/29
-
上海大學(xué) “生物有機(jī)電子材料及器件”團(tuán)隊碩士、博士研究生招聘
+1/28
-
【青島大學(xué)】2026年生物與醫(yī)藥申請考核制博士生招生（含少數(shù)民族骨干人才）
+1/21
-
西華大學(xué)材料學(xué)院表面科學(xué)與工程技術(shù)科研團(tuán)隊2026年招收研究生
+1/14
-
食品學(xué)碩考研調(diào)劑
+1/12
-
港科大（廣州）招收光學(xué)傳感器，集成電路設(shè)計，醫(yī)學(xué)影像設(shè)備以及可穿戴設(shè)備方向的博士
+1/9
-
寧夏大學(xué)化學(xué)化工學(xué)院魏逸彬研究員課題組招聘啟事
+1/7
-
5科大招收材料專業(yè)碩士
+1/3
-
【博士招生】天津理工大學(xué)國家杰青王鐵課題組招收2026年博士研究生
+1/2
-
上海師范大學(xué)有機(jī)化學(xué)專業(yè)研究生招生
+1/2
-
青島大學(xué)物理科學(xué)學(xué)院國家級人才團(tuán)隊招收材料、物理、能源等相關(guān)專業(yè)調(diào)劑碩士研究生
+1/2
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
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ahbiao 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +8 | leaflight 2026-03-24 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 348求調(diào)劑 +3 | 小懶蟲不懶了 2026-03-28 | 3/150 |
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +4 | 大饞小子 2026-03-28 | 4/200 |
|
|
|
[考研] 070300化學(xué)354求調(diào)劑 +3 | 101次希望 2026-03-28 | 3/150 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|