版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

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

ahbiao

管理員

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

[交流] forcal 代碼運(yùn)行后出錯(cuò)

運(yùn)行了別人的Forcal代碼，出現(xiàn)了“Forcal表達(dá)式運(yùn)行錯(cuò)誤！
實(shí)數(shù)表達(dá)式名稱：main 所在模塊：-1
實(shí)數(shù)函數(shù)名：fcopt::SimOpt 錯(cuò)誤代碼：3”為什么？能不能幫忙調(diào)試下，指出問題所在�。≈x謝~~
代碼：!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的錯(cuò)誤輸出
  ERSET(0, 0, 0), //關(guān)閉IMSL錯(cuò)誤輸出
  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錯(cuò)誤輸出
};

[ Last edited by ahbiao on 2011-3-27 at 21:41 ]

回復(fù)此樓

» 猜你喜歡

天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有295人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有28人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有22人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)

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

求助用MS計(jì)算聲子譜出錯(cuò)的問題已經(jīng)有15人回復(fù)
fortran轉(zhuǎn)換CHGCAR出錯(cuò) 已經(jīng)有4人回復(fù)
VASP并行計(jì)算出錯(cuò) 已經(jīng)有9人回復(fù)
Gaussian 計(jì)算拉曼光譜最后時(shí)刻出錯(cuò) 已經(jīng)有14人回復(fù)
急！關(guān)于nsfc申請(qǐng)代碼的問題已經(jīng)有6人回復(fù)
得不到純產(chǎn)物，哪出錯(cuò)了，親們幫忙看看啊已經(jīng)有3人回復(fù)
高斯09中的L502出錯(cuò)的問題已經(jīng)有5人回復(fù)
關(guān)于CASTEP 中計(jì)算時(shí)出錯(cuò) 已經(jīng)有6人回復(fù)
出錯(cuò) warning#2 #180 已經(jīng)有7人回復(fù)
選錯(cuò)申請(qǐng)代碼了，影響大嗎已經(jīng)有29人回復(fù)
MS一打開就出錯(cuò)，任何任務(wù)都不能運(yùn)行已經(jīng)有12人回復(fù)
博士后基金申請(qǐng)?jiān)诖蛴∞D(zhuǎn)化成word時(shí)出錯(cuò) 已經(jīng)有13人回復(fù)

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

查看全部散金貼

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

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

ahbiao

專家顧問

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

引用回帖:

Originally posted by forcal at 2011-03-29 10:10:20:
由于fcopt庫更新，以上代碼已比較老，故出現(xiàn)錯(cuò)誤。
新的代碼為：
[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的回帖

查看全部 5 個(gè)回答

forcal

主管區(qū)長(zhǎng)

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的回帖

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的回帖

ahbiao

管理員

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

恩，謝謝！

回復(fù)此樓

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

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

查看全部 5 個(gè)回答

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 311求調(diào)劑 +8	李芷新1 2026-03-31	8/400	2026-04-01 01:39 by 1018329917
[考研] 一志愿武理材料工程302調(diào)劑環(huán)化或化工 +11	Doleres 2026-03-31	12/600	2026-04-01 01:29 by 1018329917
[考研] 352分-085602-一志愿985 +6	海納百川Ly 2026-03-29	6/300	2026-03-31 21:06 by yuq
[考研] 284求調(diào)劑 +9	小熊～～ 2026-03-31	9/450	2026-03-31 18:22 by 253863592
[考研] 070300化學(xué)354求調(diào)劑 +15	101次希望 2026-03-28	15/750	2026-03-31 17:58 by jp9609
[考研] 求調(diào)劑生物學(xué) 377分 +6	zzll03 2026-03-31	6/300	2026-03-31 17:33 by 唐沐兒
[考研] 315求調(diào)劑 +6	akie... 2026-03-28	7/350	2026-03-31 16:48 by asdfzly
[考研] 317求調(diào)劑 +12	蛋黃咸肉粽 2026-03-26	13/650	2026-03-31 16:44 by asdfzly
[考研] 土木304求調(diào)劑 +3	頂級(jí)擦擦 2026-03-31	3/150	2026-03-31 13:59 by 西京學(xué)院招辦
[考研] 22408 359分調(diào)劑 +4	Qshers 2026-03-27	8/400	2026-03-31 08:53 by Qshers
[有機(jī)交流] 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理 10+3	kaobao456 2026-03-29	4/200	2026-03-30 23:16 by nBu鋰
[考研] 調(diào)劑 +4	GK72 2026-03-30	4/200	2026-03-30 20:32 by dick_runner
[考研] 359求調(diào)劑 +5	王了個(gè)楠 2026-03-25	5/250	2026-03-30 19:36 by 源_2020
[考研] 283求調(diào)劑（080500） +14	A child 2026-03-27	14/700	2026-03-30 12:06 by 探123
[考研] 290求調(diào)劑 +3	dfffsar 2026-03-29	3/150	2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8	AZMK 2026-03-25	12/600	2026-03-28 18:15 by AZMK
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 265求調(diào)劑11408 +3	劉小鹿lu 2026-03-27	3/150	2026-03-27 20:53 by nihaoar
[考研] 08開頭275求調(diào)劑 +4	拉誰不重要 2026-03-26	4/200	2026-03-27 14:12 by Delta2012
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

ahbiao

[交流] forcal 代碼運(yùn)行后出錯(cuò)

» 猜你喜歡

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

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

ahbiao

forcal

forcal

ahbiao

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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

查看全部散金貼