倍频晶体KTP的允许角随波长的变化matlab代码

倍频晶体KTP的允许角随波长的变化可以通过以下代码实现：

% 倍频晶体KTP的允许角随波长的变化

% 声明常量
lambda = linspace(0.5, 1.5, 1000); % 波长范围
L = length(lambda);
d = 3.5; % KTP 晶体厚度
n1 = sqrt(3.0066 + 0.03901./(lambda.^2 - 0.04251)); % 基波折射率
n2 = sqrt(3.0333 + 0.04154./(lambda.^2 - 0.04547)); % 二倍频折射率
theta = zeros(1, L); % 允许角

% 计算允许角
for i = 1:L
    k1 = 2*pi/lambda(i)*n1(i);
    k2 = 2*pi/lambda(i)*n2(i);
    theta(i) = asin(sqrt((k1^2 - k2^2)/(2*k1^2)));
end

% 绘制允许角随波长的变化曲线
plot(lambda, theta*180/pi, 'linewidth', 2);
xlabel('波长 (μm)');
ylabel('允许角 (°)');
title('倍频晶体KTP的允许角随波长的变化');

代码中，使用 `linspace` 函数定义了波长范围，使用 `sqrt` 函数计算了基波和二倍频的折射率，并使用 `asin` 函数计算了允许角。最后，使用 `plot` 函数绘制了允许角随波长的变化曲线。