 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2014 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
bustertang新蟲 (初入文壇)
|
[求助]
Kramers-Kronig Relations 推導(dǎo)介電常數(shù) 已有2人參與
|
|
|
希望用epsilon2推導(dǎo)epsilon1，已經(jīng)找到了別人的matlab程序，但是看不懂input中的alpha（the value of the moment considered）指的是什么意思，另外兩個參數(shù)omiga，imchi都有。不知道alpha該怎么設(shè)置。 %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 |
回復(fù)此樓» 猜你喜歡
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有90人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
計算與實驗科學(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)有28人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有20人回復(fù)
Thou_actor
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 358.8
- 散金: 18
- 帖子: 91
- 在線: 56.5小時
- 蟲號: 7621203
- 注冊: 2017-12-14
- 性別: GG
- 專業(yè): 機械測試?yán)碚撆c技術(shù)
3樓2019-04-03 16:03:43
2樓2017-12-21 06:36:49
xgp326265322
鐵蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 48.5
- 紅花: 3
- 帖子: 427
- 在線: 130.2小時
- 蟲號: 2409989
- 注冊: 2013-04-09
- 專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)
4樓2020-03-18 17:24:22
xgp326265322
鐵蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 48.5
- 紅花: 3
- 帖子: 427
- 在線: 130.2小時
- 蟲號: 2409989
- 注冊: 2013-04-09
- 專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)
5樓2020-03-18 19:45:27
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 309分085801求調(diào)劑 +3 | 學(xué)員Gtwj7W 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +10 | 美味可樂雞翅 2026-03-26 | 11/550 |
|
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
|
[考研] 315求調(diào)劑 +6 | akie... 2026-03-28 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +9 | 藍月亮亮 2026-03-30 | 9/450 |
|
|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 8/400 |
|
|
|
[考研] 調(diào)劑 +4 | GK72 2026-03-30 | 4/200 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 +7 | YYYYYNNNNN 2026-03-26 | 9/450 |
|
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +7 | JKSOIID 2026-03-26 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +5 | 小卡不卡. 2026-03-29 | 5/250 |
|
|
|
[考研] 0703化學(xué)321分求調(diào)劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
|
[考研] 總分293求調(diào)劑 +8 | 加一一九 2026-03-25 | 11/550 |
|
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 一志愿上海理工能源動力（085800）310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|