来自 Kramers Kronig 关系的实际折射率

信息处理 Python 频率响应 转换功能 希尔伯特变换 scipy
2022-02-23 03:34:08

我测量了折射率的复杂部分k(其中折射率为m=n+ik) 在跨越几个数量级的波长或频率的非线性网格上测量。

给定k,n应该从 Kramers-Kronig 关系推导出来。

但是,使用测试数据集n众所周知,我没有恢复正确的解决方案。知道我做错了什么(nu=频率阵列)吗?

import scipy.fftpack as ft

plt.loglog(nu, n, label='n')
plt.loglog(nu, k, label='k')
plt.loglog(nu, 1.0 + ft.hilbert(k) / np.pi, label='n from k (orig. grid)')
plt.xlabel('frequency [Hz]')
plt.legend()

在此处输入图像描述

0个回答
没有发现任何回复~