熔融石英折射率随波长变化用sellmeier公式Matlab
时间: 2024-04-27 15:23:29 浏览: 9
熔融石英的折射率与波长之间可以使用Sellmeier公式进行拟合,公式如下:
$n^2(\lambda) = 1 + \frac{B_1\lambda^2}{\lambda^2-C_1^2} + \frac{B_2\lambda^2}{\lambda^2-C_2^2} + \frac{B_3\lambda^2}{\lambda^2-C_3^2}$
其中,$n(\lambda)$表示熔融石英在波长为$\lambda$时的折射率,$B_1$、$B_2$、$B_3$、$C_1$、$C_2$和$C_3$是拟合系数。
在Matlab中,可以使用以下代码进行拟合和绘图:
```matlab
% 熔融石英的折射率数据
wavelength = [0.21, 0.23, 0.25, 0.27, 0.29, 0.31, 0.33, 0.35, 0.37, 0.39];
refractive_index = [1.4506, 1.4463, 1.4425, 1.4390, 1.4358, 1.4329, 1.4302, 1.4277, 1.4254, 1.4232];
% 使用Sellmeier公式进行拟合
lambda_sq = wavelength.^2;
B = [0.6961663, 0.4079426, 0.8974794];
C = [0.0684043^2, 0.1162414^2, 9.896161^2];
n_sq = 1 + B(1)*lambda_sq./(lambda_sq-C(1)) + B(2)*lambda_sq./(lambda_sq-C(2)) + B(3)*lambda_sq./(lambda_sq-C(3));
n = sqrt(n_sq);
% 绘制折射率随波长的变化曲线
plot(wavelength, refractive_index, 'o', wavelength, n, '-')
xlabel('波长(μm)')
ylabel('折射率')
legend('原始数据', '拟合曲线')
```
这段代码首先定义了熔融石英的折射率数据,然后使用Sellmeier公式进行拟合,得到折射率的取值$n$。最后,绘制出折射率随波长变化的曲线。