24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2008  |  回复: 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsh182

新虫 (初入文坛)

【答案】应助回帖

alpha设为0即可,该程序出自Kramers–Kronig Relations in Optical Materials Research一书,不知楼主后来正常使用没有?
2楼2017-12-21 06:36:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Thou_actor

新虫 (小有名气)

楼主现在知道alpha是什么意思,怎么设置了吗?我也有这个问题
3楼2019-04-03 16:03:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xgp326265322

铁虫 (正式写手)

您好,你的问题解决了吗?
4楼2020-03-18 17:24:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xgp326265322

铁虫 (正式写手)

您好,请问哪个代表的是反射光谱,如何使用反射光谱求解呢?
5楼2020-03-18 19:45:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xgp326265322

铁虫 (正式写手)

【答案】应助回帖

引用回帖:
2楼: Originally posted by lsh182 at 2017-12-21 06:36:49
alpha设为0即可,该程序出自Kramers–Kronig Relations in Optical Materials Research一书,不知楼主后来正常使用没有?

您好,它的这个代码哪个代表反射光谱,如何通过反射光谱计算介电常数
6楼2020-03-18 20:07:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bustertang 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂一志愿合肥工业大学 +3 秋二十二 2026-03-30 3/150 2026-03-30 23:59 by jp9609
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
[考研] 281求调剂 +5 亚克西good 2026-03-26 7/350 2026-03-30 20:42 by dophin1985
[考研] 一志愿:西北大学,英一数一408-284分求调剂 +5 12.27 2026-03-27 5/250 2026-03-30 19:31 by 源_2020
[考研] 0703 化学 求调剂,一志愿山东大学 342 分 +7 Shern—- 2026-03-28 7/350 2026-03-30 16:31 by nothing投稿中
[考研] 求调剂 +7 青春裁为三截 2026-03-29 7/350 2026-03-30 13:14 by laoshidan
[考研] 348求调剂 +6 小懒虫不懒了 2026-03-28 6/300 2026-03-30 10:29 by Evan_Liu
[考研] 一志愿211,335分,0856,求调剂院校和导师 +7 倾____萧 2026-03-27 8/400 2026-03-30 09:37 by longlotian
[考博] 26申博自荐 +6 whh869393 2026-03-24 6/300 2026-03-29 21:05 by nxgogo
[考研] 一志愿双一流机械285分求调剂 +4 幸运的三木 2026-03-29 5/250 2026-03-29 14:49 by Miko19
[考研] 一志愿北京理工大学本科211材料工程294求调剂 +8 mikasa的围巾 2026-03-28 8/400 2026-03-29 12:48 by 无际的草原
[考研] 298求调剂 +3 种圣赐 2026-03-29 3/150 2026-03-29 12:06 by longlotian
[考研] 330分求调剂 +5 qzenlc 2026-03-29 5/250 2026-03-29 07:37 by 无际的草原
[考研] 071000生物学求调剂,初试成绩343 +7 小小甜面团 2026-03-25 7/350 2026-03-28 20:25 by 唐沐儿
[考研] 085602 307分 求调剂 +7 不知道叫什么! 2026-03-26 7/350 2026-03-28 09:57 by 神马都不懂
[考研] 308求调剂 +7 墨墨漠 2026-03-25 7/350 2026-03-27 14:47 by 狂炫麦当当
[考研] 考研调剂 +10 呼呼?~+123456 2026-03-24 10/500 2026-03-27 11:46 by wangjy2002
[考研] 0703化学338求调剂! +6 Zuhui0306 2026-03-26 7/350 2026-03-27 10:35 by shangxh
[考研] 281求调剂 +6 Koxui 2026-03-24 7/350 2026-03-26 15:37 by 无际的草原
[考研] 296求调剂 +4 汪!?! 2026-03-25 7/350 2026-03-25 16:41 by 汪!?!
信息提示
请填处理意见