应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » Kramers-Kronig Relations 推导介电常数
61/1返回列表
查看: 2045  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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楼 2016-02-16 03:38:15
已阅   回复此楼   关注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 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭