當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

ahbiao

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2037.5
帖子: 73
在線: 50.4小時
蟲號: 804318

[交流] 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   Acid=AcidArray.A(j), T=TempArray.A(j),
  IDO=1, t1=tArray.A(0,0), HCp=tArray.A(0,j+1),
  i=1, (i   t2=tArray.A(i,0),
  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ù)此樓

» 猜你喜歡

重慶交大材料學(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é)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有130人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有26人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有8人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

求助用MS計(jì)算聲子譜出錯的問題已經(jīng)有15人回復(fù)
fortran轉(zhuǎn)換CHGCAR出錯已經(jīng)有4人回復(fù)
VASP并行計(jì)算出錯已經(jīng)有9人回復(fù)
Gaussian 計(jì)算拉曼光譜最后時刻出錯已經(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ì)算時出錯已經(jīng)有6人回復(fù)
出錯 warning#2 #180 已經(jīng)有7人回復(fù)
選錯申請代碼了，影響大嗎已經(jīng)有29人回復(fù)
MS一打開就出錯，任何任務(wù)都不能運(yùn)行已經(jīng)有12人回復(fù)
博士后基金申請?jiān)诖蛴∞D(zhuǎn)化成word時出錯已經(jīng)有13人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-03-27 21:39:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

forcal

禁蟲 (初入文壇)

ahbiao(金幣+40): 2011-03-31 08:34:11

本帖內(nèi)容被屏蔽

4樓2011-03-30 15:03:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 5 個回答

forcal

禁蟲 (初入文壇)

ahbiao(金幣+10): 謝謝~~ 2011-03-30 08:10:16

本帖內(nèi)容被屏蔽

2樓2011-03-29 10:10:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ahbiao

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2037.5
帖子: 73
在線: 50.4小時
蟲號: 804318

引用回帖:

Originally posted by forcal at 2011-03-29 10:10:20:
由于fcopt庫更新，以上代碼已比較老，故出現(xiàn)錯誤。
新的代碼為：
[code]!using["math","IMSL","fcopt"];
f(t,HCp,dHCp : R : k1,a,E1, Acid,T)=
{ R=8.3144,
dHCp=-k1*Acid ...

使用Opt函數(shù)，如何用啊，麻煩您具體說下！

贊一下

回復(fù)此樓

3樓2011-03-30 08:12:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ahbiao

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2037.5
帖子: 73
在線: 50.4小時
蟲號: 804318

恩，謝謝！

回復(fù)此樓

5樓2011-03-31 08:34:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 274求調(diào)劑 +6	xiao愛同學(xué) 2026-03-30	6/300	2026-03-31 10:04 by cal0306
[考研] 本2一志愿C9-333分，材料科學(xué)與工程，求調(diào)劑 +5	升升不降 2026-03-31	5/250	2026-03-31 09:18 by 535743368
[考研] 一志愿食品科學(xué)與工程083200求調(diào)劑 +4	XQTJZ 2026-03-30	4/200	2026-03-31 04:10 by fmesaito
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +10	dongfang59 2026-03-27	10/500	2026-03-30 23:42 by 果果媽咪
[考研] 合肥區(qū)域性重點(diǎn)一本招收調(diào)劑 +4	6266jl 2026-03-30	4/200	2026-03-30 21:48 by zhuangyan123
[考研] 085602化工求調(diào)劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導(dǎo)導(dǎo)
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +12	辰沐5211314 2026-03-26	12/600	2026-03-30 19:28 by Wang200018
[考研] 材料化工340求調(diào)劑 +3	jhx777 2026-03-30	3/150	2026-03-30 17:54 by JourneyLucky
[碩博家園] 求調(diào)劑有機(jī)化學(xué)考研356分 +10	Nadiums 2026-03-25	11/550	2026-03-30 11:56 by yjolah
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6	王一一依依 2026-03-30	8/400	2026-03-30 11:52 by yjolah
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4	qianbai11 2026-03-29	4/200	2026-03-30 11:34 by 唐沐兒
[考研] 311求調(diào)劑 +6	冬十三 2026-03-24	6/300	2026-03-29 20:45 by 無際的草原
[考研] 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。 +11	nnnnnnn5 2026-03-25	11/550	2026-03-29 20:42 by 無際的草原
[考研] 總分293求調(diào)劑 +8	加一一九 2026-03-25	11/550	2026-03-29 19:53 by 無際的草原
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12	c297914 2026-03-28	12/600	2026-03-29 16:21 by Serene1974
[考研] 349求調(diào)劑 +6	李木子啊哈哈 2026-03-25	6/300	2026-03-29 12:47 by 無際的草原
[考研] 298求調(diào)劑 +3	種圣賜 2026-03-29	3/150	2026-03-29 12:06 by longlotian
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)
[考研] 309求調(diào)劑 +4	gajsj 2026-03-25	5/250	2026-03-26 00:27 by Dyhoer