如图求主点里程及坐标。(线路起点里程0Km.docVIP

写如图求主点里程及坐标。（线路起点里程0Km） 程序如下： #includestdio.h #includemath.h #define PI 3.14159 main() { double x0,y0,K0,x2,y2,K2,x1,y1,R,l01,l02; double A,A12,A01; double b01,b02,p1,p2,m1,m2,b1; double T1,T2,L,r1,r2,E0; double D01,D02,KZH,KHY,KQZ,KYH,KHZ; double XZH,YZH,XHZ,YHZ,xHY,yHY,XHY,YHY,XYH,YYH; double KZ1=0.0,KH1=10.0,KY=0.0,KH2=0.0,KZ2=0.0; double X[60],Y[60],x[60],y[60],dr[60]; double l10,l20; int i=1,j,k; x0=0.0;y0=0.0;K0=0.0; x2=5156.0;y2=8316.0; x1=4000.0;y1=4000.0; R=800.0;l01=60.0;l02=60.0; printf(计算转角\n); A12=atan((y2-y1)/(x2-x1)); A01=atan((y1-y0)/(x1-x0)); A=(A12-A01)*180/PI;/*把A从弧度转为 °*/ if(A0) printf( 曲线右转\n); if(A0) printf( 曲线左转\n); printf(\n); /*计算曲线要素*/ printf(计算曲线要素\n); printf( 缓和曲线常数\n); b01=(l01/(2*R))*180/PI;/*把b01从弧度转为 °*/ b02=(l02/(2*R))*180/PI;/*把b02从弧度转为 °*/ p1=(l01*l01)/(24*R); p2=(l02*l02)/(24*R); m1=l01/2-(l01*l01*l01)/(240*R*R); m2=l02/2-(l02*l02*l02)/(240*R*R); printf( b01=%f b02=%f p1=%.3f p2=%.3f m1=%.3f m2=%.3f\n,b01,b02,p1,p2,m1,m2); printf( 曲线综合要素\n); T1=tan(A*PI/(2*180))*(R+p1)+m1-((p1-p2)/sin(A*PI/180));/*A用弧度*/ T2=tan(A*PI/(2*180))*(R+p2)+m1+((p1-p2)/sin(A*PI/180)); L=((A-b01-b02)*PI/180)*R+l01+l02; r1=(atan((R+p1)/(T1-m1)))*180/PI;/*把r1从弧度转为 °*/ r2=(atan((R+p2)/(T2-m2)))*180/PI;/*把r2从弧度转为 °*/ E0=((R+p1)/sin(r1*PI/180))-R;/*r1转为弧度*/ printf( T1=%.3f T2=%.3f L=%.3f r1=%f r2=%f E0=%.3f\n,T1,T2,L,r1,r2,E0); printf(\n); /*计算各主点里程*/ printf(计算各主点里程\n); printf( 点距离的计算\n); D01=sqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)); D02=sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1)); printf( D01=%.3f, D02=%.3f\n,D01,D02); printf( 曲线各主点里程\n); KZH=K0+(D01-T1); KHY=KZH+l01; KQZ=KHY+(PI/2-b01*PI/180-r1*PI/180); KYH=KZH+L-l01; KHZ=KZH+L; K2=(D01-T1)+L+(D02-T2); printf( KZH=%.3f KHY=%.3f KQZ=%.3f KYH=%.3f KHZ=%.3f K2=%.3f\n,KZH,KHY,KQZ,KYH,KHZ,K2); /*计算各主点坐标*/ XZH=x0+(KZH-K0)*cos(A01); YZH=y0+(