 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2793 | 回復: 15 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[交流]
【求助】使用MATLAB怎么實現(xiàn)擬合力場參數(shù)的程序？【已解決】 已有6人參與
|
|||
|
可以將根據(jù)VDW公式計算出的E_VDW_total視為各個原子間距離r（r是向量）以及原子VDW參數(shù){sigma}、{epsilon}變量集的函數(shù)。r (i)代表第i個掃描點的原子間距離向量，定義誤差函數(shù)ErrF=∑(r(i)下量化計算的能量-r(i)下的E_VDW_total)^2，然后對ErrF對{sigma}和{epsilon}變量集中的每個變量求導得0，如果寫出來是線性形式可以用矩陣方程來解，如果不是的話可以用非線性優(yōu)化方法解，解出{sigma}、{epsilon}，用這樣參數(shù)計算的VDW作用可以對所有掃描的點的能量都能較好描述�？赡苷f得比較抽象，這與擬合ESP電荷的方法比較像，可看相關(guān)文獻。 [ Last edited by nono2009 on 2010-12-1 at 08:20 ] |
回復此樓» 猜你喜歡
08工科275分求調(diào)劑
已經(jīng)有3人回復
-
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有14人回復
-
安徽大學專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑
已經(jīng)有11人回復
-
070300化學求調(diào)劑
已經(jīng)有12人回復
-
化學308分調(diào)劑
已經(jīng)有12人回復
-
359求調(diào)劑
已經(jīng)有3人回復
-
一志愿武理材料工程302調(diào)劑環(huán)化或化工
已經(jīng)有8人回復
-
材料工程085601數(shù)二英一335求調(diào)劑
已經(jīng)有5人回復
-
一志愿北京科技大學材料學碩328分求調(diào)劑
已經(jīng)有3人回復
-
0856求調(diào)劑
已經(jīng)有9人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab擬合擬均相動力學參數(shù)
已經(jīng)有6人回復
- 需要用matlab擬合數(shù)據(jù)的外行人急需幫助
已經(jīng)有10人回復
- 求助��！求一程序，用matlab程序做，用蒙特卡洛方法模擬
已經(jīng)有10人回復
- matlab非線性參數(shù)擬合問題
已經(jīng)有7人回復
- 如何用MATLAB 實現(xiàn)化學反應方程式（寫程序代碼）?
已經(jīng)有20人回復
- Matlab人工神經(jīng)網(wǎng)絡(luò)工具箱的位置及使用
已經(jīng)有5人回復
- 看不懂matlab程序該怎么辦？
已經(jīng)有7人回復
- 【求助】如何根據(jù)勢能曲線擬合力場參數(shù)
已經(jīng)有19人回復
9樓2010-06-02 08:19:23
2樓2010-05-31 19:05:46
lijinfeng042
木蟲 (小有名氣)
Matlab
- 仿真EPI: 2
- 應助: 1 (幼兒園)
- 金幣: 2156.1
- 散金: 115
- 帖子: 291
- 在線: 31.5小時
- 蟲號: 1019062
- 注冊: 2010-05-15
- 性別: GG
- 專業(yè): 通信理論與系統(tǒng)
3樓2010-05-31 20:01:07
|
我是在問問題，這個是我寫的程序，其實跟最小二乘法有關(guān)： function [a1,b1,a2,b2]=forcefieldfitting(r,e) %使用DFT計算的數(shù)據(jù)擬合力場參數(shù),e是使用量化計算出現(xiàn)的相互作用能，r是兩個分子片段之間的距離 %LJ-12-06公式中的EPSILON(C):a1 %LJ-12-06公式中的SIGMA(C):b1 %LJ-12-06公式中的EPSILON(O):a2 %LJ-12-06公式中的SIGMA(O):b2 %單個的CO2的能量可以認為是兩個O的能量和一個C的能量的加和 if(length(e)==length(r)) n=length(e); else disp('e和r的維數(shù)不相等！'); return; end %維數(shù)檢查 A=zeros(4,3) B=zeros(4,2); for i=1:n A(1,1)=A(1,1)+(384*a1^2*b1^23)/r(i)^24; A(1,2)=A(1,2)+(192*a1^2*b1^11)/r(i)^12; A(1,3)=A(1,3)+(48*a1*e*b1^5)/r(i)^6; A(2,1)=A(2,1)+(32*a1*b1^24)/r(i)^24; A(2,2)=A(2,2)+(32*a1*b1^12)/r(i)^12 A(2,3)=A(2,3)+(8*e(i)*b1^6)/r(i)^6; A(3,1)=A(3,1)+(1536*a2^2*b2^23)/r(i)^24; A(3,2)=A(3,2)+(768*a2^2*b2^23)/r(i)^12; A(3,3)=A(3,3)+(96*a2*e(i)*b2^5)/r(i)^6; A(4,1)=A(4,1)+(128*a2*b2^24)/r(i)^24; A(4,2)=A(4,2)+(128*a2*b2^12)/r(i)^12; A(4,3)=A(4,3)+(16*e(i)*b2^6)/r(i)^6; B(1,1)=B(1,1)+(576*a1^2*b1^17)/r(i)^18; B(1,2)=B(1,2)+(96*a1*e(i)*b1^11)/r(i)^12; B(2,1)=B(2,1)+(64*a1*b1^18)/r(i)^18; B(2,2)=B(2,2)+(8*e(i)*a1^12)/r(i)^12; B(3,1)=B(3,1)+(2304*a2^2*b2^17)/r(i)^18; B(3,2)=B(3,2)+(192*a2*e(i)*b2^11)/r(i)^12; B(4,1)=B(4,1)+(256*a2*b2^18)/r(i)^18; B(4,3)=B(4,3)+(16*e(i)*b2^12)/r(i)^12; end s=A\B; a1=s(1); b1=s(2); a2=s(3); b2=s(4); |
4樓2010-05-31 22:50:07
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑：085600材料與化工，考材科基，總分319 +16 | 678lucky 2026-03-31 | 20/1000 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境工程 085701，267求調(diào)劑 +13 | minht 2026-03-29 | 13/650 |
|
|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +4 | 余航航 2026-03-31 | 4/200 |
|
|
|
[考研] 生醫(yī)工0831調(diào)劑求推薦 +4 | 小熊睿睿_s 2026-03-27 | 6/300 |
|
|
|
[考研] 材料科學與工程調(diào)劑 +13 | 深V宿舍吧 2026-03-30 | 14/700 |
|
|
|
[考研] 271求調(diào)劑 +5 | 勒布朗@ 2026-03-31 | 7/350 |
|
|
|
[考研] 學碩274求調(diào)劑 +17 | Li李魚 2026-03-26 | 17/850 |
|
|
|
[考研] 土木304求調(diào)劑 +4 | 兔突突突， 2026-03-31 | 5/250 |
|
|
|
[考研]
 一志愿大連理工大學材料求調(diào)劑
+6
|
Gymno 2026-03-30 | 6/300 |
|
|
|
[考研] 11408總分309，一志愿東南大學求調(diào)劑，不挑專業(yè) +5 | 天賦帶到THU 2026-03-29 | 6/300 |
|
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
|
[考研] 0703化學/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +13 | 丹青奶蓋 2026-03-26 | 15/750 |
|
|
|
[考研] 332求調(diào)劑 +6 | Lyy930824@ 2026-03-29 | 6/300 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 331環(huán)境科學與工程求調(diào)劑 +3 | 熠然好運氣 2026-03-27 | 3/150 |
|
|
|
[考博] 26申博 +3 | 加油沖啊！ 2026-03-26 | 3/150 |
|
|
|
[考研] 中國科學院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|