版塊導(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í)熱門版塊排行榜

返回列表

wdxiao

銀蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 148.2
紅花: 1
帖子: 547
在線: 104.4小時(shí)
蟲號(hào): 882485
注冊(cè): 2009-10-24
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 最優(yōu)化計(jì)算求助！（1stopt或者M(jìn)atlab）

本人有一個(gè)最優(yōu)化計(jì)算的問題。用1stopt1.5注冊(cè)版編寫的代碼。但是運(yùn)行的時(shí)候總是沒有任何顯示。故意將代碼改錯(cuò)，編譯的時(shí)候的確能報(bào)告出錯(cuò)。

另外，貌似1stopt的pascal編程與標(biāo)準(zhǔn)pascal語言有很多不同，比如：
不能自定義結(jié)構(gòu)類型數(shù)據(jù)，不能自定義函數(shù)，自定義過程貌似也有問題。等等。

請(qǐng)大俠幫忙看看。另外，不知道MatLab能否處理這一問題？

先大致說說這個(gè)問題的物理本質(zhì)和程序的思路。

本質(zhì)就是計(jì)算三個(gè)分子間的范德華力相互作用能。
每個(gè)分子有24個(gè)原子。每個(gè)原子有x,y坐標(biāo)和范德華力相關(guān)參數(shù)兩個(gè)。
所有信息存放在三維數(shù)組里面。數(shù)組的第一指標(biāo)表示不同分子；第二指標(biāo)表示分子內(nèi)各個(gè)原子；第三指標(biāo)表示每個(gè)原子有x,y坐標(biāo)和范德華力相關(guān)參數(shù)兩個(gè)。

M0為分子標(biāo)準(zhǔn)的位置（0,0）和角度。
M1,M2,M3具有相同的分子角度，即將M0的標(biāo)準(zhǔn)角度繞(0,0)點(diǎn)旋轉(zhuǎn)phi角度,
M1,M2,M3分別為放置在(0,0), (b,0),和(a*cos(theta),a*sin(theta))位置。
然后分別計(jì)算不同分子各個(gè)原子對(duì)間的范德華力作用能。

程序是用Pascal語言編寫，如下：

Title "Type your title here"
Parameters a[1,4], b[1,4],theta[0,pi/2],phi[-pi/8,pi/8];
Minimum;
StartProgram;
const xx: array[1..24] of double
      =(0.703,1.698,1.698,0.703,-0.703,-1.698,-1.698,-0.703,1.246,3.008,3.008,1.246,-1.246,
      -3.008,-3.008,-1.246,0,-2.926,-4.138,-2.926,0,2.926,4.138,2.926);
   yy: array[1..24] of double
      =(1.698,0.703,-0.703,-1.698,-1.698,-0.703,0.703,1.698,3.008,1.246,-1.246,-3.008,-3.008,
      -1.246,1.246,3.008,4.138,2.926,0,-2.926,-4.138,-2.926,0,2.926);
var i,j: Integer;
Etot,r,angle,out1,out2,out3: double;
Ax,Ay,Ar,Ae,Bx,By,Br,Be,dist,rou: double;
M: array[0..3,1..24,1..4] of double;
Begin
for i:=1 to 24 do
begin
      M[0,i,1]:=xx;
      M[0,i,2]:=yy;
end;
for i:=1 to 8 do
begin
      M[0,i,3]:= 1.9;
      M[0,i,4]:= 0.044;
      M[0,i+8,3]:= 1.9;
      M[0,i+8,4]:= 0.044;
      M[0,i+16,3]:= 2.11;
      M[0,i+16,4]:= 0.202;
end;
M[1]:=M[0];
for i:=1 to 24 do
  begin
      r:=sqrt(M[1,i,1]*M[1,i,1]+M[1,i,2]*M[1,i,2]);
      if M[1,i,1]>=0 then
            angle:=arcsin(M[1,i,2]/r)
         else
            angle:= (M[1,i,2]/abs(M[1,i,2]))*(pi-abs(arcsin(M[1,i,2]/r)));
      angle:=angle+phi;
      M[1,i,1]:=r*cos(angle);
      M[1,i,2]:=r*sin(angle);
  end;
M[2]:=M[1];
M[3]:=M[1];
for i:=1 to 24 do
  begin
M[2,i,1]:=M[1,i,1]+b;
M[3,i,1]:=M[1,i,1]+a*cos(theta);
M[3,i,2]:=M[1,i,2]+a*sin(theta);
  end;
Etot:=0;
for i:=1 to 24 do
for j:=1 to 24 do
   begin
         Ax:=M[1,i,1];
         Ay:=M[1,i,2];
         Ar:=M[1,i,3];
         Ae:=M[1,i,4];
         Bx:=M[2,j,1];
         By:=M[2,j,2];
         Br:=M[2,j,3];
         Be:=M[2,j,4];
         dist:=sqrt((Ax-Bx)*(Ax-Bx)+(Ay-By)*(Ay-By));
         rou:=dist/(Ar+Br);
         if dist>3.311 then
            out1:= sqrt(Ae*Be)*(290000*exp(-12.5*rou)-2.25/exp(6*ln(rou)))
         else
            out1:= 336.176*sqrt(Ae*Be)/(rou*rou);
         Bx:=M[3,j,1];
         By:=M[3,j,2];
         Br:=M[3,j,3];
         Be:=M[3,j,4];
         dist:=sqrt((Ax-Bx)*(Ax-Bx)+(Ay-By)*(Ay-By));
         rou:=dist/(Ar+Br);
         if dist>3.311 then
            out2:= sqrt(Ae*Be)*(290000*exp(-12.5*rou)-2.25/exp(6*ln(rou)))
         else
            out2:= 336.176*sqrt(Ae*Be)/(rou*rou);
         Ax:=M[2,i,1];
         Ay:=M[2,i,2];
         Ar:=M[2,i,3];
         Ae:=M[2,i,4];
         dist:=sqrt((Ax-Bx)*(Ax-Bx)+(Ay-By)*(Ay-By));
         rou:=dist/(Ar+Br);
         if dist>3.311 then
            out3:= sqrt(Ae*Be)*(290000*exp(-12.5*rou)-2.25/exp(6*ln(rou)))
         else
            out3:= 336.176*sqrt(Ae*Be)/(rou*rou);
      Etot:=out1+out2+out3+Etot;
      end;
FunctionResult:= Etot;
End;
EndProgram;

[ Last edited by wdxiao on 2012-3-26 at 14:06 ]

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有85人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有29人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有53人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

1stOpt求解常微分方程邊值問題已經(jīng)有5人回復(fù)
哪位大哥能幫我用1stopt3.0以上版本跑一個(gè)擬合程序已經(jīng)有11人回復(fù)
1stOpt的函數(shù)作圖不靠譜嗎？已經(jīng)有4人回復(fù)
求助有關(guān)Matlab有約束非線性最優(yōu)化問題已經(jīng)有9人回復(fù)
用origin或1stopt擬合一條隱函數(shù)方程曲線已經(jīng)有4人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)
求助關(guān)于1stopt的問題已經(jīng)有6人回復(fù)
【求助】MUSIC的MD模塊計(jì)算甲烷在IRMOF-1上的擴(kuò)散系數(shù) 已經(jīng)有11人回復(fù)
用1stopt進(jìn)行曲線擬合每次給出結(jié)果不一樣已經(jīng)有6人回復(fù)
用1stopt進(jìn)行數(shù)據(jù)擬合時(shí)，無法運(yùn)行sharedmodel ,求助高手已經(jīng)有7人回復(fù)
求助關(guān)于1stopt擬合~··~~ 已經(jīng)有4人回復(fù)
【求助】精通matlab最優(yōu)化計(jì)算，一書中的funval函數(shù)是什么意思已經(jīng)有10人回復(fù)
【求助】如何計(jì)算有機(jī)小分子的三重態(tài)能級(jí)（T1-S0）？已經(jīng)有12人回復(fù)

1樓 2012-03-26 13:54:34

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

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

不知道你的M值是在哪賦值的？

贊一下

回復(fù)此樓

2樓2012-03-26 14:28:57

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

wdxiao

銀蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 148.2
紅花: 1
帖子: 547
在線: 104.4小時(shí)
蟲號(hào): 882485
注冊(cè): 2009-10-24
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

2樓: Originally posted by dingd at 2012-03-26 14:28:57:
不知道你的M值是在哪賦值的？

xx,yy數(shù)組分別賦值給M[0,1..24,1]和M[0,1..24,2]，表示標(biāo)準(zhǔn)位置的分子內(nèi)各個(gè)原子的x,y坐標(biāo)。對(duì)應(yīng)的表示范德華力的參數(shù)M[0,1..24,3]和M[0,1..24,4]由第一個(gè)for...do循環(huán)賦值。（分子的1~16號(hào)原子是同一種原子，17~24號(hào)原子是另一種原子）。

M0賦值給M1。然后M1旋轉(zhuǎn)phi角度后，賦值給M2和M3。然后M2和M3的各個(gè)原子坐標(biāo)經(jīng)過平移變換。

贊一下

回復(fù)此樓

3樓2012-03-26 14:47:22

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

wdxiao

銀蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 148.2
紅花: 1
帖子: 547
在線: 104.4小時(shí)
蟲號(hào): 882485
注冊(cè): 2009-10-24
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

2樓: Originally posted by dingd at 2012-03-26 14:28:57:
不知道你的M值是在哪賦值的？

贊一下

回復(fù)此樓

4樓2012-03-26 14:50:40

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

wdxiao

銀蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 148.2
紅花: 1
帖子: 547
在線: 104.4小時(shí)
蟲號(hào): 882485
注冊(cè): 2009-10-24
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

2樓: Originally posted by dingd at 2012-03-26 14:28:57:
不知道你的M值是在哪賦值的？

xx,yy數(shù)組分別賦值給M[0,1..24,1]和M[0,1..24,2]，表示標(biāo)準(zhǔn)位置的分子內(nèi)各個(gè)原子的x,y坐標(biāo)（第一個(gè)for....do循環(huán)）。對(duì)應(yīng)的表示范德華力的參數(shù)M[0,1..24,3]和M[0,1..24,4]由第二個(gè)for...do循環(huán)賦值。（分子的1~16號(hào)原子是同一種原子，17~24號(hào)原子是另一種原子）。

M0賦值給M1。然后M1旋轉(zhuǎn)phi角度后（第三個(gè)for...do循環(huán)），賦值給M2和M3。然后M2和M3的各個(gè)原子坐標(biāo)經(jīng)過平移變換（第四個(gè)for.....do 循環(huán)）。

另外，很奇怪的是，我的源程序上第一個(gè)for.....do 循環(huán)為：

for i:=1 to 24 do
begin
M[0,i,1]:=xx;
M[0,i,2]:=yy;
end;

發(fā)帖的時(shí)候自動(dòng)去掉了，不知道為啥。

贊一下

回復(fù)此樓

5樓2012-03-26 14:56:46

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

wdxiao

銀蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 148.2
紅花: 1
帖子: 547
在線: 104.4小時(shí)
蟲號(hào): 882485
注冊(cè): 2009-10-24
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

5樓: Originally posted by wdxiao at 2012-03-26 14:56:46:
xx,yy數(shù)組分別賦值給M和M，表示標(biāo)準(zhǔn)位置的分子內(nèi)各個(gè)原子的x,y坐標(biāo)（第一個(gè)for....do循環(huán)）。對(duì)應(yīng)的表示范德華力的參數(shù)M和M由第二個(gè)for...do循環(huán)賦值。（分子的1~16號(hào)原子是同一種原子，17~24號(hào)原子是另一種原子 ...

for i:=1 to 24 do
begin
M[1,i,1]:=xx ;
M[1,i,2]:=yy ;
end;

贊一下

回復(fù)此樓

6樓2012-03-26 15:10:32

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

wdxiao

銀蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 148.2
紅花: 1
帖子: 547
在線: 104.4小時(shí)
蟲號(hào): 882485
注冊(cè): 2009-10-24
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

奇怪，還是沒法寫成xx, yy

贊一下

回復(fù)此樓

7樓2012-03-26 15:11:09

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

wdxiao

銀蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 148.2
紅花: 1
帖子: 547
在線: 104.4小時(shí)
蟲號(hào): 882485
注冊(cè): 2009-10-24
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

"xx, yy"

回復(fù)此樓

8樓2012-03-26 15:11:36

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

wdxiao

銀蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 148.2
紅花: 1
帖子: 547
在線: 104.4小時(shí)
蟲號(hào): 882485
注冊(cè): 2009-10-24
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

xx【i】
yy【i】

回復(fù)此樓

9樓2012-03-26 15:12:17

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

wdxiao

銀蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 148.2
紅花: 1
帖子: 547
在線: 104.4小時(shí)
蟲號(hào): 882485
注冊(cè): 2009-10-24
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

真奇怪，發(fā)帖的時(shí)候沒法用英文的中括號(hào)寫xx i 。

贊一下

回復(fù)此樓

10樓2012-03-26 15:12:53

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 wdxiao 的主題更新

返回列表

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

[考研] 321求調(diào)劑 +3	何潤采123 2026-03-18	3/150	2026-03-18 21:27 by li123456789.
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +8	SUICHILD 2026-03-12	8/400	2026-03-18 20:55 by winsuccess
[考研] 材料專碩英一數(shù)二306 +4	z1z2z3879 2026-03-18	4/200	2026-03-18 20:00 by 楤哥
[考研] 085601材料工程專碩求調(diào)劑 +6	慕寒mio 2026-03-16	6/300	2026-03-18 14:26 by 007_lilei
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6	薛云鵬 2026-03-13	6/300	2026-03-18 14:14 by 脫穎而出
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[基金申請(qǐng)] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[論文投稿] 投稿問題 5+4	星光燦爛xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian

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

[求助] 最優(yōu)化計(jì)算求助！（1stopt或者M(jìn)atlab）

» 猜你喜歡

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

【答案】應(yīng)助回帖

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