| 查看: 1998 | 回复: 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 |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有162人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有26人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有34人回复
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
20