 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2015 | 回復: 5 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
bustertang新蟲 (初入文壇)
|
[求助]
Kramers-Kronig Relations 推導介電常數(shù) 已有2人參與
|
|
|
希望用epsilon2推導epsilon1，已經找到了別人的matlab程序，但是看不懂input中的alpha（the value of the moment considered）指的是什么意思，另外兩個參數(shù)omiga，imchi都有。不知道alpha該怎么設置。 %The program inputs are 1) omega, vector of the frequency %(or energy) components, 2) imchi, vector of the imaginary %part of the susceptibility under examination, and 3) alpha, %the value of the moment considered. The two vectors %1) and 2) must have the same length. %The output is the estimate of the real part as obtained %with K-K relations. %In order to use this program, save the whole text contained %in this section in a file and name it kkrebook.m if size(omega,1)>size(omega,2); omega=omega'; end; if size(imchi,1)>size(imchi,2); imchi=imchi'; end; %Here the program rearranges the two vectors so that, %whichever their initial shape, they become row vectors. g=size(omega,2); %Size of the vectors.% rechi=zeros(size(imchi)); %The output is initialized. a=zeros(size(imchi)); b=zeros(size(imchi)); %Two vectors for intermediate calculations are initialized deltaomega=omega(2)-omega(1); %Here we compute the frequency (or energy) interval j=1; beta1=0; for k=2:g; b(1)=beta1+imchi(k)*omega(k)^(2*alpha+1)/... (omega(k)^2-omega(1)^2); beta1=b(1); end; rechi(1)=2/pi*deltaomega*b(1)*omega(1)^(-2*alpha); %First element of the output: the principal part integration %is computed by excluding the first element of the input j=g; alpha1=0; for k=1:g-1; a(g)=alpha1+imchi(k)*omega(k)^(2*alpha+1)/... (omega(k)^2-omega(g)^2); alpha1=a(g); end; rechi(g)=2/pi*deltaomega*a(g)*omega(g)^(-2*alpha); %Last element of the output: the principal part integration %is computed by excluding the last element of the input for j=2:g-1; ; %Loop on the inner components of the output vector. alpha1=0; beta1=0; for k=1:j-1; a(j)=alpha1+imchi(k)*omega(k)^(2*alpha+1)/... (omega(k)^2-omega(j)^2); alpha1=a(j); end; for k=j+1:g; b(j)=beta1+imchi(k)*omega(k)^(2*alpha+1)/... (omega(k)^2-omega(j)^2); beta1=b(j); end; rechi(j)=2/pi*deltaomega*(a(j)+b(j))*omega(j)^(-2*alpha); end; %Last element of the output: the principal part integration %is computed by excluding the last element of the input |
回復此樓» 猜你喜歡
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8
已經有0人回復
-
0702一志愿吉大B區(qū)求調劑有論文
已經有4人回復
-
物理學I論文潤色/翻譯怎么收費?
已經有95人回復
-
歡迎調劑到贛南師范大學智能制造與未來能源學院
已經有28人回復
-
0702一志愿吉大B區(qū)求調劑 本科期間發(fā)表一篇Sci
已經有4人回復
-
歡迎調劑到贛南師范大學智能制造與未來能源學院
已經有26人回復
-
重慶交大26年碩士生招生擬調劑通知已出！歡迎加入光子學微結構與器件課題組。
已經有0人回復
-
UJN物理學專業(yè)調劑
已經有3人回復
-
0702一志愿吉大B區(qū)求調劑 本科期間發(fā)表一篇Sci
已經有0人回復
xgp326265322
鐵蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 48.5
- 紅花: 3
- 帖子: 427
- 在線: 130.2小時
- 蟲號: 2409989
- 注冊: 2013-04-09
- 專業(yè): 化工熱力學和基礎數(shù)據(jù)
4樓2020-03-18 17:24:22
2樓2017-12-21 06:36:49
Thou_actor
新蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 358.8
- 散金: 18
- 帖子: 91
- 在線: 56.5小時
- 蟲號: 7621203
- 注冊: 2017-12-14
- 性別: GG
- 專業(yè): 機械測試理論與技術
3樓2019-04-03 16:03:43
xgp326265322
鐵蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 48.5
- 紅花: 3
- 帖子: 427
- 在線: 130.2小時
- 蟲號: 2409989
- 注冊: 2013-04-09
- 專業(yè): 化工熱力學和基礎數(shù)據(jù)
5樓2020-03-18 19:45:27
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 385分 生物學（071000）求調劑 +3 | qf626 2026-04-01 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境工程297分求調劑一志愿杭高院 +10 | GENJIOW 2026-03-31 | 11/550 |
|
|
|
[考研] 070300化學279求調劑 +15 | 哈哈哈^_^ 2026-03-31 | 17/850 |
|
|
|
[考研] 261求調劑 +3 | 明仔· 2026-04-01 | 3/150 |
|
|
|
[考研] 一志愿鄭大材料工程290求調劑 +18 | Youth_ 2026-03-30 | 18/900 |
|
|
|
[考研] 生物學學碩，一志愿湖南大學，初試成績338 +8 | YYYYYNNNNN 2026-03-26 | 10/500 |
|
|
|
[考研] 一志愿中農0710生物學，微生物方向總分338求調劑 +3 | 柒xxxx. 2026-03-26 | 3/150 |
|
|
|
[考研] 一志愿北交材料工程總分358 +5 | cs0106 2026-04-01 | 7/350 |
|
|
|
[考研] 0855機械初試280求調劑 +3 | kazenotori 2026-03-31 | 3/150 |
|
|
|
[考研] 289求調劑 +7 | BrightLL 2026-03-29 | 7/350 |
|
|
|
[考研] 求調劑 +9 | akdhjs 2026-03-31 | 10/500 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
|
[考研] 083000學碩274求調劑 +12 | Li李魚 2026-03-26 | 12/600 |
|
|
|
[考研] 食品工程專碩一志愿中海洋309求調劑 +5 | 小張zxy張 2026-03-26 | 10/500 |
|
|
|
[考研] 327求調劑 +5 | 小卡不卡. 2026-03-29 | 5/250 |
|
|
|
[考研] 328求調劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
|
[考研] 0703 化學 求調劑，一志愿山東大學 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
|
[考研] 332求92調劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 一志愿南京航空航天大學材料學碩求調劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
|
[考研] 315調劑 +4 | 0860求調劑 2026-03-26 | 5/250 |
|