| 查看: 2143 | 回复: 5 | ||
bustertang新虫 (初入文坛)
|
[求助]
Kramers-Kronig Relations 推导介电常数 已有2人参与
|
|
希望用epsilon2推导epsilon1，已经找到了别人的matlab程序，但是看不懂input中的alpha（the value of the moment considered）指的是什么意思，另外两个参数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 |
回复此楼» 猜你喜欢
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
物理学I论文润色/翻译怎么收费?
已经有259人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有109人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有64人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
2楼2017-12-21 06:36:49
Thou_actor
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 358.8
- 散金: 18
- 帖子: 91
- 在线: 56.5小时
- 虫号: 7621203
- 注册: 2017-12-14
- 性别: GG
- 专业: 机械测试理论与技术
3楼2019-04-03 16:03:43
xgp326265322
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 48.5
- 红花: 3
- 帖子: 427
- 在线: 130.2小时
- 虫号: 2409989
- 注册: 2013-04-09
- 专业: 化工热力学和基础数据
4楼2020-03-18 17:24:22
xgp326265322
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 48.5
- 红花: 3
- 帖子: 427
- 在线: 130.2小时
- 虫号: 2409989
- 注册: 2013-04-09
- 专业: 化工热力学和基础数据
5楼2020-03-18 19:45:27
xgp326265322
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 48.5
- 红花: 3
- 帖子: 427
- 在线: 130.2小时
- 虫号: 2409989
- 注册: 2013-04-09
- 专业: 化工热力学和基础数据
6楼2020-03-18 20:07:36
10