使用Mathematica评估Kramers-Kronig关系

问题描述

我正在尝试使用Kramers-Kronig关系从Rb的透射数据计算折射率，但是积分无法收敛，我尝试了其他多个问题的代码建议，但都无济于事，谁能帮忙我吗？

我的代码如下：

ClearAll["Global`*"];
c = 3*10^8;
Subscript[\[Lambda],0] = 780.24*10^-9;
Subscript[\[Nu],0] = c/Subscript[\[Lambda],0]; 
RbT24L75 = Import["RbT24L75.txt",{"Data",All,{1,2}}];
ListLinePlot[{RbT24L75},AxesOrigin -> {-5*10^9,0},Frame -> True,FrameLabel -> {"Frequency detuning (Hz)","Transmission"}]
Freq = 2 \[Pi] (RbT24L75[[All,1]] + Subscript[\[Nu],0]);
FreqDt = 2 \[Pi] RbT24L75[[All,1]];
RbT24L75[[All,1]] = Freq;
Trans = RbT24L75[[All,2]];
Im\[Chi]1 = -((Log[Trans] c )/(0.75 Freq));
RbT24L75[[All,2]] = Im\[Chi]1;
func = Interpolation[RbT24L75];
Plot[func[x],{x,Min[Freq],Max[Freq]},PlotRange -> All]
Re\[Chi]1[\[Omega]_] := (2 \[Omega]/\[Pi])*NIntegrate[func[x]/(x^2 - \[Omega]^2),Method -> "PrincipalValue",Exclusions -> x^2 - \[Omega]^2 == 0] //Quiet
Plot[Re\[Chi]1[\[Omega]],{\[Omega],Max[Freq]}]

我导入的数据给出了失谐频率与传输之间的关系： Transmission of Rb

非常感谢！

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）