第二近轴光的实际参数.cpp

第二近轴光的实际参数.cpp 近轴光线计算

#include <math.h>
#include <stdio.h>
#define P 3.1415926
double r[]={4.699,-855,-9.117,3.77,104.5,3.77,-6.406},
n[]={1,1.78830877,1,1.68893168,1,1.61659181,1.7844303,1},
d[]={1.33373,0.588139,0.6059186,0.3354896,0.6085918,1.503606};
double PA[20],qiucha1,qiucha2,qiucha3,qiucha4,wsecha1,wsecha2,wsecha3,wsecha4;
double Ip[20],IIp[20],UUp[20],hmax,umax,LLF,LLC,YY02,YY03;
double L,LL,U,x[12],LLp,Up,Lp,nF[20],nC[20],nD[20];
double UU,I,II,h,t2,YY0,x2,D,Xt,Xs,Xts,ff=1.0,Kt[10],YYF,YYC,YFC;
double lz,l,lp,ll,llp,l2,i,ii,u,up,uu,uup,wsecha,U1[30],L1[30],UU1[30],LL1[30];
double W,yy0,LLL,UUU,qiucha,PA1,PA2,x1,tt,k,ss,s,syy,llll,YY;
int j,t,jj,c,cm;
char ch;

void LU(double L,double U)
{
	int i,h;
double a,b,L2,I,II,UU,LL;
h=1;
U=0;
for(i=0;i<7;i++)
{
if(U==0) a=h/r[i];
else a=(L-r[i])*(sin(U*P/180))/r[i];
I=asin(a)/P*180;
b=a*n[i]/n[i+1];
II=asin(b)/P*180;
UU=U+I-II;
LL=b*r[i]/(sin(UU*P/180))+r[i];
L2=LL-d[i];
L=L2;U=UU;

}
printf("\n第一近轴光的实际参数如下：U'=%lf,L'=%lf\n",UU,LL);
}
void lup(double lp,double up)
{
	int i;
double a,b,l2,uup,llp,W;
lp=0.3222;
W=20.556;
up=sin(W*3.1415926/180);
for(i=0;i<7;i++)
{
a=(lp-r[i])*up/r[i];
b=a*n[i]/n[i+1];
uup=up+a-b;
llp=b*r[i]/uup+r[i];
l2=llp-d[i];
lp=l2;up=uup;
}
printf("lp=0.3222,W=20.556");
printf("\n第二近轴光的理想参数如下：up'=%lf,lp'=%lf\n",uup,llp);
}

void LUp(double Lp,double Up)
{
double a,b,L2;
Lp=0.3222;
W=20.556;
up=sin(W*3.1415926/180);
for(jj=0;jj<7;jj++)
{
a=(Lp-r[jj])*(sin(up*P/180))/r[jj];
I=asin(a)/P*180;
b=a*n[jj]/n[jj+1];
II=asin(b)/P*180;
UU=up+I-II;
LLp=b*r[jj]/(sin(UU*P/180))+r[jj];
L2=LLp-d[jj];
PA[jj]=Lp*sin(up*P/180)/cos((I-up)/2*P/180);
Ip[jj]=I;IIp[jj]=II;UUp[jj]=UU;
Lp=L2;up=UU;
}

printf("\n第二近轴光的实际参数如下：\nIp'=%lf,Up'=%lf,Lp'=%lf\n",II,UU,LLp);
}



void main()
{
	int i,h1=1;
  double U,L,F;
  double j[10];
  double jj[10];
  double l[10]={exp(15)};
  double u[10]={0};
  for(i=0;i<7;i++){
      if(i==0)
		   j[i]=h1/r[i];
	   else
	   j[i]=(l[i]-r[i])*u[i]/r[i];
	   jj[i]=(n[i]/n[i+1])*j[i];
	   u[i+1]=u[i]+j[i]-jj[i];
	   l[i+1]=r[i]+jj[i]*r[i]/u[i+1]-d[i];
	   U=u[i+1];
	   L=l[i+1];
	   F=h1/u[i+1];
  }
  printf("系统焦距F=%f\n",F);
  printf("第一近轴光的理想参数如下:u'=%f,像距l'=%f\n",U,L);
  double LL,UU;
LU(LL,UU);
   double lz,W,lp,up;
  lup(lp,up);
LUp(lz,W);


}