写一下sellmeier公式并用Matlab表达出来

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}+...$$

其中$n(\lambda)$表示波长为$\lambda$的光在该介质中的折射率，$B_i$和$C_i$是与该介质的性质有关的常数。

在MATLAB中，可以通过编写一个函数来计算Sellmeier公式的值。以下是一个简单的MATLAB代码示例：

function n = sellmeier(B, C, lambda)
% 计算Sellmeier公式的值
% B和C是常数向量，lambda是波长向量
n2 = 1 + B(1).*lambda.^2./(lambda.^2-C(1).^2) + ...
     B(2).*lambda.^2./(lambda.^2-C(2).^2) + ...
     B(3).*lambda.^2./(lambda.^2-C(3).^2);
n = sqrt(n2);
end

使用该函数，可以计算出在给定波长下介质的折射率。例如，计算玻璃的折射率在波长为550nm时的值，可以这样做：

B = [1.03961212, 0.231792344, 1.01046945];
C = [0.00600069867, 0.0200179144, 103.560653];
lambda = 550e-9;
n = sellmeier(B, C, lambda);
disp(n)  % 输出结果为1.5157