单面透镜光路和像差计算

function y=a2(W,M,J) %此函数文件作用：由物点W，方向向量M，透镜参数J得到 位于X=J(5)【J(5)代表像面距第一面透镜x轴方向的距离】处的像点坐标B,X=0处像点坐标K和方向向量E x=W(1); %计算实际光线通过透镜后的像点坐标及方向矢量 y=W(2); z=W(3); xi=M(1); yi=M(2); zi=M(3); d=J(1);%J(1)代表物点距第一面透镜的距离 n1=J(2);%J(2)代表物面所处介质折射率 n2=J(3);%J(3)代表像面所处介质折射率 r=J(4);%J(4)代表透镜的曲率半径 ci=1/r; e=d*xi-(x*xi+y*yi+z*zi); Mx=x+e*xi-d; M2=x*x+y*y+z*z-e*e+d*d-2*d*x; pp=M2*ci-2*Mx; cosI2=xi*xi-ci*pp; cosII2=1+(n1/n2)*(n1/n2)*(cosI2-1); g=sqrt(cosII2)-sqrt(cosI2)*(n1/n2); D=e+pp/(xi+sqrt(cosI2)); xx=x+D*xi-d;%像点坐标为（xx,yy,zz) yy=y+d*yi; zz=z+d*zi ; xxi=(n1/n2)*xi-g*(xx*ci-1);%方向矢量为（xxi,yyi,zzi) yyi=(n1/n2)*yi-g*yy*ci; zzi=(n1/n2)*zi-g*zz*ci; L=(J(5)-xx)/xxi; %计算位于x=J(5)处平面上的像点坐标 L=(x1-x')/xxi=(y1-y')/yyi=(z1-z,)/zzi y1=L*yyi+yy; z1=L*zzi+zz; global B; B(1)=J(5); B(2)=y1; B(3)=z1;%数组B用来存放x=60mm处平面上像点坐标（x1,y1,z1) global K; K(1)=xx; K(2)=yy; K(3)=zz;%数组K用来存放第一面透镜上像点坐标（x1,y1,z1) global E; E(1)=xxi; E(2)=yyi; E(3)=zzi;%数组K用来存放通过第一面透镜后光线方向向量