global epsilon0 epsilon11 epsilon33 k11 k22 k33 theta phi v h d p0 q0 dim lamda no ne
clear

%设定光学常数
lamda = [650e-9, 550e-9, 450e-9];
beta = pi/4;

%设定材料常数
k11=1.37e-12;k22=1.0e-12;k33=1.68e-12;
epsilon0=8.85e-12;
epsilon11=10.3epsilon0;
epsilon33=3.6epsilon0;
vth=pi*sqrt(k33/(epsilon33-epsilon0)*epsilon0);

%设定结构常数
d=6.0e-6; %厚度
dim=201; %网络节点数目
h=d/dim; %网络节点间距

theta(1)=(pi/180)*89; %预倾角1
theta(dim)=(pi/180)*89; %预倾角2
…（完整代码见文件，下载链接见下方）

指向矢仿真matlab代码下载链接

透过率仿真matlab代码下载链接