clc
clear all
close all
a=0.633*10.^(-6);
b=4000;
c=2.0*10.^(-6);%线性膨胀系数
c0=10*10.^(-3);%长度
d=2.2*10.^(-6);%电场方向系数
d0=5*10.^(-3);%电场方向的厚度
T0=298;%室温
T=233:0.1:333;%温度范围
n=sqrt(4.913-2.78*10.^(-2)*a.^2+(0.1173+1.65*10.^(-8)*T.^2)./(a.^2-(0.212+2.7*10.^(-8)*T.^2).^2));%寻常光折射率与温度的关系
G=(a/(2*b))   *   (     (  d0+d*(T-T0) )./ (   c0+c*(T-T0))     ).*(1./(n.^3)   );%电光系数与温度的关系
figure('color',[1,1,1])
plot(T,G)
ylim([0.11,0.112]*1e-10);
title('T-G(T)温度与电光系数关系图')
xlabel('T/温度')
ylabel('G(T)/电光系数')

grid on
figure('color',[1,1,1])
plot(T,n)
ylim([1.522,1.529]);
title('T-n(T)温度与折射率关系图')
xlabel('T/温度')
ylabel('n(T)/电光系数')

grid on