单面透镜光路和像差计算

计算光线经过单面透镜折射后的实际光线的光路和像差，并绘制像差曲线图

function y=a1（Wn） %此文件作用：由一个物点W和光阑放大倍数n计算实际方向向量F
%物点坐标为Wn为光阑半径放大倍数
%数组X，YZ分别存放以1为半径的光阑上33个点的坐标
X=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];   
Y=[0 0 0 0 sqrt（2）/8 2*sqrt（2）/8 3*sqrt（2）/8 4*sqrt（2）/8 2/8 4/8 6/8 8/8 sqrt（2）/8 2*sqrt（2）/8 3*sqrt（2）/8 4*sqrt（2）/8 0 0 0 0 -1*sqrt（2）/8 -2*sqrt（2）/8 -3*sqrt（2）/8 -4*sqrt（2）/8 -2/8 -4/8 -6/8 -8/8 -1*sqrt（2）/8 -2*sqrt（2）/8 -3*sqrt（2）/8 -4*sqrt（2）/8 0];
Z=[2/8 4/8 6/8 8/8 1*sqrt（2）/8 2*sqrt（2）/8 3*sqrt（2）/8 4*sqrt（2）/8 0 0 0 0 -1*sqrt（2）/8 -2*sqrt（2）/8 -3*sqrt（2）/8 -4*sqrt（2）/8 -2/8 -4/8 -6/8 -8/8 -1*sqrt（2）/8 -2*sqrt（2）/8 -3*sqrt（2）/8 -4*sqrt（2）/8 0 0 0 0 1*sqrt（2）/8 2*sqrt（2）/8 3*sqrt（2）/8 4*sqrt（2）/8 0];
XX=n*X;%数组XXYYZZ分别放置实际光阑上点的xyz坐标
YY=n*Y;
ZZ=n*Z;
for i=1:33
    D（i）=sqrt（（XX（i）-W（1）+30）^2+（YY（i）-W（2））^2+（ZZ（i）-W（3））^2）; %光阑上的点到物点的距离为D，并考虑物面对于光阑处坐标系而言，处于X=-30mm处
    global FX;%数组FX，FY，FZ分别存放方向向量F的xyz坐标
    global FY;
    global FZ;
    FX（i）=（XX（i）-W（1）+30）/D（i）;
    FY（i）=（YY（i）-W（2））/D（i）; 
    FZ（i）=（ZZ（i）-W（3））/D（i）;     
    i=i+1;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1179  2012-05-23 09:59  单面透镜\a1.m

     文件        621  2012-05-23 10:21  单面透镜\a11.m

     文件        579  2012-05-23 10:24  单面透镜\a111.m

     文件        728  2012-05-23 10:25  单面透镜\a113.m

     文件        726  2012-05-23 10:24  单面透镜\a13.m

     文件       1214  2012-05-23 09:50  单面透镜\a2.m

     文件        736  2012-05-23 10:14  单面透镜\a3.m

     文件        393  2012-05-16 15:08  单面透镜\H0.m

     文件        112  2012-05-16 19:43  单面透镜\H00.m

     文件        296  2012-05-23 10:05  单面透镜\H1.m

     文件        146  2012-05-25 17:00  单面透镜\H11.m

     文件         76  2012-05-23 10:43  单面透镜\T1.m

     文件        701  2012-05-23 10:44  单面透镜\X1.m

     文件        758  2012-05-23 09:58  单面透镜\XX1.m

     目录          0  2013-07-30 10:27  单面透镜

----------- ---------  ---------- -----  ----

                 8265                    15