 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2000 | 回復(fù): 5 | ||
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é)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有108人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有68人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
2樓2017-12-21 06:36:49
Thou_actor
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 358.8
- 散金: 18
- 帖子: 91
- 在線: 56.5小時
- 蟲號: 7621203
- 注冊: 2017-12-14
- 性別: GG
- 專業(yè): 機(jī)械測試?yán)碚撆c技術(shù)
3樓2019-04-03 16:03:43
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
xgp326265322
鐵蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 48.5
- 紅花: 3
- 帖子: 427
- 在線: 130.2小時
- 蟲號: 2409989
- 注冊: 2013-04-09
- 專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)
6樓2020-03-18 20:07:36
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 bustertang 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +8 | 簡單點(diǎn)0 2026-03-17 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 初試 317 +3 | 半拉月丙 2026-03-20 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-17 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|