导线计算平差

1、CASIO fx-5800P 程序集 导线近似平差计算程序DXPC程序名)Lb1 9:0 DimZ: 5N-3 DimZ: “ N ? N: “ A”？ A: “ B” ? B: “ C” ? C: “ D ?D: “E”？ E: “F”？ F: “G'? G: “H'? H:Pol(C-A,D-B):lfJ<0: ThenJ+360 J: IfEnd: JJ: J T:Pol(G-E,H-F): IfJ<0:T hen J+360 J: IfEnd： JJ:J M:0 L:0 U: 0 V:0 R: “？ T=” :DMST P: “ M= :“ DMS “ I”

2、？ I: Lb1 1: “Z”？ Z:R+1R: Z ZR: IfI=-1: ThenGoto C: IfEnd : Goto D:Lb1 C:P+Z-180 P:Goto E:Lb1 D:P-Z+180 P:Goto E:Lb1 E: IfP<0:Then P+36C P: Goto B:Else If P> 360: Then P-360 P:Goto B: IfEnd： PP:Goto B: Lb1 B: If R<N: Then Goto 1: ifEnd: Goto 4: Lb1 4:0 R: “P=” ：3600(P-M) P，IfI=-1: Then “Q= :

3、 -P - N QGoto N: IfEnd:“Q兰:p十 MGoto N:Lb1 N: T Z2N:Lb1 2: “S”？S:R+1 R:S ZN+R:L+S L: IfI=-1: T h e n Goto F: IfEnd: Goto G:Lb1 F:ZR+Z2N+R-1+Q - 3600 -180 Z2N+R:Goto H:Lb1 G:Z2N+R-1-ZR-Q -3600 +180 Z2N+R:Goto H:Lb1 H: IfZ2N+R<0: Th e n Z2N+R+360 Z2N+R:Goto A: Else If Z2N+R> 360: ThenZ2N+R-360 Z2

4、N+R:Goto A: IfEnd: Z2N+R Z2N+R:Goto A:Lb1 A: “ Z2N+R= ” :Z2N+R DMS“ Z3N+R-1=” :S cos (Z2N+R) Z3N+R-1 U+Z3N+R-1 U:“Z4N+R-2= ” :S sin (Z2N+R) Z4N+R-2 -V+Z4N+R-2 V:IfR+1<N: The n Goto 2: IfEnd: “L=” :L，“ U=” :1000(U+C-E) U，“ V=” : 1000(V+D-F) V，“K=” : Int(1000L -V(U2+V2) K，0 R:C X:D Y:Lb1 3:R+1 R:-U

5、ZN+R-L 0:-VZR+N“W= :W，X+Z3N+R-1+.001 0 X:Y+Z4N+R-2+. 001W Y: “X=” :X， “丫二” :Y，IfR+1<N: ThenGoto 3: IfEnd: Goto注:1.程序中5N-3在计算时应根据实际观测角数按5N-3式计算后输入实际数值 如观测角数N为13,则“ 5N-3”改为“ 62”第128頁符合导线计算表工程名称：韶关市曲江至南雄段高速公路合同段：C02标占号八、观测角(右)O!"角度改正/ ff 值()坐标方位角o!边长(m)坐标增量值(m)改正数(mm)坐标(m)备注 X Y X XXXGPS0748 58

6、 43.3342802.30063406.609已知GPS08208 25 24-0.49443134.81763788.839已知20 33 19.82470.964440.9789165.3623-8.69.8SG035131 25 02-0.49443575.78763954.21169 08 18.32292.529104.1729273.3518-5.46.1SG036207 16 18-0.49443679.95564227.56941 52 0.81497.719370.6500332.1788-9.110.4SG037135 01 02-0.49444050.59664559.

7、75886 50 59.31377.22620.7299376.6560-6.97.9SG037-1263 02 24-0.49444071.31964936.4223 48 3 5.8260.675260.098917.3211-4.85.4SG03890 19 18-0.49444331.41364953.74993 29 18.29403.306-24.5398402.5587-7.48.4SG039242 45 10-0.49444306.86665356.31630 44 8.79475.05408.3213242.7884-8.79.9SG040152 18 34-0.494447

8、15.17865599.11458 25 35.28474.029248.1980403.8580-8.79.9SG041168 19 58-0.49444963.36866002.98270 05 37.77407.383138.7061383.0425-7.58.5SG042209 58 35-0.49445102.06666386.03340 07 3.27481.730368.3904310.4067-8.810.0SG043201 06 17-0.49445470.44866696.45019 00 46.76375.488355.0032122.3274-6.97.8SG044-1

9、156 49 30-0.49445825.44466818.78542 11 17.26570.944423.0374383.4272-10.511.9GPS09199 30 30-0.49446248.47167202.224已知22 40 47.75GPS1046803.15067434.023已知辅助计算右=2366° 18' 02"刀 D=5087.043X=3113.7472Y=3413.2789f B = -6.42 " fx=-0.0932 fy=-0.1061f=0.1412f B 容 = + 36 ”fB <f B 容 符合精度要求

10、K=1/（刀D/f） 1/36022<1/1500符合精度要求示意 图第129頁说明:第一步：N?13观测角数.第一步：A?:42802.300起始设站点后视点已知X坐标.第二步:B?:63406.609起始设站点后视点已知Y坐标第三步：C?:43134.817起始设站点已知X坐标第四步:D?63788.839）起始设站点已知Y坐标第六步：E.?46248.471终点设站点已知X坐标第七步：F?67202.224终点设站点已知Y坐标第八步：G.?46803.150终点设站点前视点已知X坐标.第九步：H?67434.023终点设站点前视点已知 Y坐标第130頁第十步:1:1观测角为左角时“

11、 |=-1 ” ,观测角为右角时“ 1=1第 一步：P: 48 ° 58' 43.33 "第一次显示起始已知方位角.（以度分秒显示）第十二步:M:22 ° 40' 47.75 "显第十三步：Z?208°25'24"观测角值.第十四步：Z?131° 25'02"观测角值第十五步：Z? 207° 16'18"观测角值第十六步：Z?135°01 '02 "观测角值第十七步：Z?263° 02'24"观测角值

12、第十八步：Z?90 °19'18 "观测角值第十九步：Z?242° 45'10"观测角值第二十步：Z? 152° 18'34"观测角值示终边已知方位角第二十一步:Z? 168 ° 19' 58O58第二十二步：Z? 209"观测角值35 "观测角值第二十三步:Z? 201 ° 06' 17"观测角值第二十四步：Z? 1564930 "观测角值第二十五步:Z? 199° 30' 30 "观测角值第二十六步：P:

13、 -6.42 第二次显示导线角度闭合差.（以秒显示）第二十七步:Q=-0.494显示角度改正值.（以秒计）第二十八步:S?470.964平距观测值.第二十九步：Z2N+R= 20 ° 3319.82 "显示改正后方位角.（以度分秒显示）第三十步：Z3N+R-1=440.9789 显示 X 增量.（以 M计）Z4N+R-2=165.3623显示 Y增量.（以M计）第三一步：S?292.529平距观测值第三十二步：Z2N+R= 69 ° 0818.32 "显示改正后方位角.（以度分秒显示）第三十三步：Z3N+R-1=104.1729 显示 X增量.（以M计）

14、Z4N+R-2=273.3518显示 Y增量.（以M第三十四步：S?497.719平距观测值第三十五步：Z2N+R=41 ° 52' 0.81 "显示改正后方位角.（以度分秒显示第三十六步：Z3N+R-1=370.6500 显示 X增量.（以M计）Z4N+R-2=332.1788显示 Y增量.（以M计）第三十七步：S?377.226平距观测值第三十八步：Z2N+R= 86 ° 50' 59.31 "显示改正后方位角.（以度分秒显示）第三十九步：Z3N+R-1=20.7299 显示 X 增量.（以 M计）Z4N+R-2=376.6560显示

15、 Y 增量.（以 M计）第四十步：S?260.675平距观测值第四一步：Z2N+R= 3 ° 48' 35.8 "显示改正后方位角.（以度分秒显示）第四十二步：Z3N+R-1=260.0989 显示 X 增量.（以 M 计）Z4N+R-2=17.3211显示 Y 增量.（以 M计）第四十三步：S?403.306平距观测值第四十四步：Z2N+R= 93 ° 29' 18.29 "显示改正后方位角.（以度分秒显示）第四十五步：Z3N+R-1=-24.5398 显示 X 增量.（以 M计）Z4N+R-2=402.5587显示 Y 增量.（以 M

16、计）第四十六步：S?475.050平距观测值第四十七步：Z2N+R= 30° 44' 8.79 "显示改正后方位角.（以度分秒显示）第四十八步：Z3N+R-1=408.3213 显示 X 增量.（以 M计） Z4N+R-2=242.7884显示 Y增量.（以M计）第四十九步：S?474.029平距观测值第五十步：Z2N+R= 58° 25' 35.28 "显示改正后方位角.（以度分秒显示）第五十一步：Z3N+R-1=248.198 显示 X 增量.（以 M计）Z4N+R-2=403.858显示 Y 增量.（以 M计）第五十二步：S?407

17、.383平距观测值第五十三步：Z2N+R= 70 ° 05' 37.77 "显示改正后方位角.（以度分秒显示）第五十四步：Z3N+R-1=138.7061 显示 X 增量.（以 M计）Z4N+R-2=383.0425显示 Y增量.（以M计）第五十五步：S?481.730平距观测值第五十六步：Z2N+R= 40 ° 07' 3.27 "显示改正后方位角.（以度分秒显示第五十七步：Z3N+R-1=368.3904 显示 X增量.（以M计） Z4N+R-2=310.4067显示 Y增量.（以M计）第五十八步：S?375.488平距观测值第十九步

18、：Z2N+R= 19° 00' 46.76 "显示改正后方位角.（以度分秒显示）第六十步：Z3N+R-1=355.0032 显示 X增量.（以M计）Z4N+R-2=122.3274显示 Y 增量.（以 M第六一步：S?570.944平距观测值第六十二步:Z2N+R= 42 ° 11 ' 17.26 "显示改正后方位角.（以度分秒显示）第六十三步:Z3N+R-1=423.0374 显示 X 增量.（以 M计）Z4N+R-2=383.4272 显示 Y增量.（以第六十四步:L=5087.043显示距离总长度.（以M计）第六十五步:U=93.2

19、纵坐标X闭合差.（以毫米计）第六十六步：V=-106.1横坐标Y闭合差.（以毫米计）第六十七步:K=1/36022边长相对精度.第六十八步:0=-8.6纵坐标X改正值.（以毫米计）W=9.8横坐标Y改正值.（以毫米计）第六十九步:X=43375.787平差后坐标值.Y=63954.211平差后坐标值第七十步:O=-5.4 纵坐标X改正值.（以毫米计）W=6.1横坐标Y改正值.（以毫米计）第七十一步:X=43679.955平差后坐标值.Y=64559.758平差后坐标值第七十二步:O=-9.1纵坐标X改正值.（以毫米计）W=7.9横坐标Y改正值.（以毫米计）第七十三步:X=44050.596平差

20、后坐标值.Y=64559.758平差后坐标值第七十四步:O=-6.9纵坐标X改正值.（以毫米计）W=7.9横坐标Y改正值.（以毫米计）O第七十五步:X=44071.319平差后坐标值.Y=64936.422平差后坐标值第七十六步O=-4.8纵坐标X改正值.（以毫米计）W=5.4横坐标Y改正值.（以毫米计）第七十七步:X=44331.413平差后坐标值.Y=64953.749平差后坐标值第七十八步:O=-7.4纵坐标X改正值.（以毫米计）W=8.4横坐标Y改正值.（以毫米计）第七十九步:X=44306.866平差后坐标值.Y=65356.316平差后坐标值第吧十步：0=8.7纵坐标X改正值.（以

21、毫米计）W=9.9横坐标Y改正值.（以毫米计）O第吧一步:X=44713.178平差后坐标值.Y=65599.114平差后坐标值第八十二步:O=-8.7纵坐标X改正值.（以毫米计）W=9.9横坐标Y改正值.（以毫米计）第八十三步:X=44963.368平差后坐标值.Y=66002.982平差后坐标值第八十四步:O=-7.5纵坐标X改正值.（以毫米计）W=8.5横坐标Y改正值.（以毫米计）第131頁平差后坐标值横坐标丫改正值.（以毫米计）平差后坐标值横坐标丫改正值.（以毫米计）0平差后坐标值横坐标丫改正值.（以毫米计）平差后坐标值第八十五步：X=45102.066平差后坐标值.Y=66386.0

22、33第八十六步：0=-8.8纵坐标X改正值.（以毫米计）W=10.0第八十七步：X=45470.448平差后坐标值.丫=66696.450第八十八步：0=-6.9纵坐标X改正值.（以毫米计）W=7.8第八十九步：X=45825.444平差后坐标值.丫=66818.785第九十步：0=-10.5纵坐标X改正值.（以毫米计）W=11.9第九一步：X=46248.471平差后坐标值.丫=67202.224CASIO -fx-5800P 程序集锥坡坐标计算程序ZPJS（锥坡计算程序）LbI 9:“A”?A:“B”？B:“C”？C:“M'?M: “N” ？N:“I” ?I:N-Abs (N) J

23、: IfM<AbsN : The n Goto E:IfE nd: Goto F: LbI E: MAbs sin (I)E： N-JMcos(I) F: IfE> AbsF: Then E W Abs (F) V: Goto A: IfEnd: i V : Abs( F) W Goto A: LbI F: Abs(Nsin(l) )t E : M-N Jcos(I) F : IfE>F : ThenE W： F V： Goto A : IfEnd: E t V : F t W: Goto A :LbI A : W-V2 Z:V( n 2 V2+4Z) D: 2W- (V+W

24、) T: (W-2V) -Zt S: IfV>.5W : Then D+V(T+V 2S- W) n V+4Z) t L:Goto B： Else IfV<.5W : Then vumvuv-wvA>A>rwrD+V (T+w"s- V2) n 2v"+4Z ) t L: Goto B : IfEnd: D+V(T n 2v"+4Z ) t L: Goto B : LbIB: L - 4t L: “ L=” :L 丄 If M< Abs (N ): ThenA+Mcos(l)cos (C+I) t G: B+Mcos(I)sin (C

25、+I) t H: Goto C : IfEnd: A+NJcos(I) cos(C)t G: B+Jn cos(I) si n(C)t H: Goto C:LbI C: -1 t R: int L t K: LbI 0: R+1tr：VF- K t q： w/(1-R -K ) tP：|fM<Abs(N)Then G+Pcos(C+I)+JQsin(C+I) tx：H+Psin(C+l)-JQcos(C+l) ty： “X=” ：X丄“Y=” :YGoto2: IfEnd: G+Pcos(C)-QJsin (C)tx： H+Psin(C)+QJcos(C)ty： “X=” :X丄 “Y仝

26、：丫 丄 Goto 2: LbI 2: IfR<K: ThenGoto 0: IfEnd:Goto 9切线方位角说明：1 .A? 0锥顶点X坐标B? 0锥顶点Y坐标2 .C? 0锥顶点纵轴（锥顶点至轴端点方向）线方位角。（以桥台背墙线中心点处切线方位角为准）3. M? 12锥坡纵向 M轴长N ? 5锥坡横向N轴长（左“-”右“ +”）4. I ? 60锥坡纵向轴M与横向轴N之间的夹角（左“-”右“ +”）5. L=11.2926596锥坡曲线计算长度8 . X=12 .Y=0显示锥坡长轴等分点坐标，等分数取L的整数部分9. X=11.9607.Y=0.3936显示锥坡长轴等分点坐标10.

27、 X=11.8417.Y=0.7873 显示锥坡长轴等分点坐标11. X=11.6399.Y=1.1809 显示锥坡长轴等分点坐标12. X=11.3496.Y=1.5746 显示锥坡长轴等分点坐标13. X=10.9619 .Y=1.9682显示锥坡长轴等分点坐14. X=10.4623 .Y=2.3619 显示锥坡长轴等分点坐标15. X=9.8282 . Y=2.7555显示锥坡长轴等分点坐标16. X=9.0203.Y=3.1492 显示锥坡长轴等分点坐标17. X=7.9621.Y=3.5428 显示锥坡长轴等分点坐标18. X=6.4577.Y=3.9365 显示锥坡长轴等分点坐标

28、19. X=2.5.Y=4.3301显示锥坡长轴等分点坐标CASIO fx-5800P 程序集直线、圆曲线、缓和曲线坐标计算程序(单独通用)DHQXJS主程序名)0 Dimz: 2Dimz: “ 1.HS说Y”：“ 2.XYHS : “ N ?N:“ A” ?A:“B” ?B“C ?C:“F” ?F:“E” ?E:“R ?R:“Z” ?Z:“T” ?T: NZ1: Z 2 -24- R-ZA (4) *2688 -RA (3)+ZA (6) - 506880- RA (5)-ZA(8) - 154828800 - RA (7) Y:T2- 24 - R-T A (4) - 2688 - RA

29、(3)+T A (6) - 506880 - RA (5) -TA (8)-154828800 - RA (7) I:.5Z-ZA(3) - 240 -氏+ZA (5) - 34560 - RA (4)- ZA (7) - 8386560 -RA (6)D:.5T-T A (3) - 240 - R+T A (5) - 34560 - R A (4)- T A (7) - 8386560 - R A (6) G:(R+Y) tan(.5 Abs(E)+D-(Y-I)- sin Abs( (E) Q:(R+I) tan(.5 Abs (E)+G+(Y-I)- sin (Abs (E) V: Ab

30、s( E) n R-180-.5(Z+T) Z2: Int(E) - Abs(Int (E) D: If N=1: Th-n Goto A:lfEnd:Goto B:Lb1 A:“H”？ H:“S”？ S:Prog “ DHQXJS1 : Z1 N: “U=” :U DMS “X=” :X “ Y=” :Y “ D”？ D:IfD 工0:ThnGoto A : IfEnd:“P”？ P: “ K”？ K:XOY I: U+K U:L bl 8: “ J”？ J: “ W ? W: “ M' ? M: C+.5MJ cos( U).r r .rtr r -sin( P)+W cos (U

31、+P)-sin( P) L:I+.5MJ sin( U) - sin (P)+W sin (U+P)-sin( P) N: “L=” 丄“N=” : N / IfW 工 0:Then Goto 8:Else If J 工 0:Then Goto 8:IfEnd: Z1 N: Goto A : Lb1 B:“ J”？ J:“P”？ P: Prog “ DHQXJS” ： Z1 N: F-Q+W H: “ H=' :H “S=” :S Z1 N:Goto BDHQXJS1子 程序名 1)If H < F-Q: Then Goto 5: Else IfH> F-Q+Z+Z2+T:

32、 Then Goto 6: Else If H < F-Q+Z: Then Goto 1 : Else.-iinrvi"-rwirvwwjwrr-j-v-v-x-a-a-urvB-s-jHnnirM-Hs-sr2If H > F-Q+Z+Z2 : Then Goto 3 : IfEnd: Goto 2:Lb1 1: H-F+Q I: 90- n *R*Z Y: C+DYI U:K cos (YX 2),0,I) G: j( sin (YX 2), 0,I) C A-Q cos(C)+G cos(C)-D Csin (C) N: B-Q sin (C)+G sin(C)+D

33、Ccos (C) I:Goto 4:Lb1 2: H-F+Q I: 180(I-.5Z)* n * F Y: C+DY U: R sin (Y)+.5Z-ZA (3) - 240- R2 G: R(1-cos (Y)+Z 2- 24- R C A-Q cos (C)+G cos( C)-DCsin( C) N: B-Q sin( C)+G sin (C)+DCcos( C) I:Goto 4:Lb1 3: F-Q+Z+Z2+T-H I: 90 *n - R- T Y: C+E-DYI2 U: Kcos (YX2),20,I) G: /(sin (YX),0,I) C A+V cos(C+E)-

34、G cos(C+E)-D O sin(C+E) N: B+V sin(C+E)-Gsin(C+E)+D O cos(C+E) I:Goto 4:Lb1 5: H-F+Q I: A+(I-Q)cos(C) N: B+(I-Q)sin(C) I: C U:Goto 4:L b1 6: H+V-F+Q-Z-Z2-T I: A+I cos(C+E) N: B+I sin (C+E ) I: C+E U:Goto 4:Lb14: IfU v 0: Then U+360 U:Goto D: Else If U >360: Then U-360 U:Goto D: IfEnd:U U:Goto D:L

35、b1D: N-Ssin( U)X: I+ S cos( U )DHQXJS2子程序名 2）(J-A+Qcos(C)cos(C)+(P-B+Qsin(C)sin(C)W: F-Q+W H:0 S:Lb1 0:Prog“ DHQXJS” : (P-Y)sinW: F-Q+W H: 0 S:Goto(U)+(J-X)cos( U) S: IfAbs( S)<10 A (-4) : ThenGoto 1: IfEnd: W+S0:Lb1 1: 0 S:Prog “DHQXJS” ： (P-Y)-cos (U) S程序运行示例及说明运行主程序“ DHQXJS、任意点坐标计算（道路、桥、涵）第一步：

36、N? 1（进入道路、桥、涵任意点坐标计算）.第二步：A? 39720.624 (输入交点的“ X”坐标).第三步：B? 61232.290 (输入交点的“ Y”坐标).第四步：C? 5 ° 55' 35.76 ” （输入交点的第一方位角）.第五步：F=11180.245 （输入交点的里程桩号）.第六步：E? 39 ° 20' 39.3 ” （输入交点的转角左“-”、右“ +”）.第七步：R? 2200 （输入曲线半径）.第八步：Z? 300 （输入第一缓和曲线长度）.（独立圆曲线时“ Z=0”）第九步：T? 300 （输入第二缓和曲线长度）.（独立圆曲线时“

37、 T=0”）第90页第十步：H? 11000 （输入待求点桩号）.第一步：S? 0 （输入待求点至中桩的距离（中心桩为“ 0”、左为“-” 右为“ +”）. 第十二步：U=21° 43' 51.8 "（待求点所计算的切线方位角）.第十三步：X=39524.931 （显示待求点所计算的“ X”坐标）.Y=61297.303 （显示待求点所计算的“ Y”坐标）.第十四步：D? 1（输入“1”时返回第十步，再重复步骤第十一步十四步计算各点的道路坐标；当输入“0”时进入桥、涵坐标计算，再继续下面步骤）.第十五步：P? 135（输入桥、涵中心点处纵向轴线（路线前进方向）与横向

38、轴线（从左到右方向）之间的夹角.如图）.第十六步：K? 0 （输入桥、涵中心点处纵向轴线（路线前进方向）偏离切线方向的偏角(左“-”右“ +”）如图）.第十七步：J? 8（输入桥、涵的垂直跨距（桥、涵各点至中心处横向轴线的垂直距离的2 倍).第十八步：W? 20（输入桥、涵上各点至纵向轴线（路线前进方向）的垂直距离（左“-”右“ +”）如图）.第十九步：M? 1（计算以桥涵中心处横轴线为界在前面各点的坐标;（M=1为前面、M=-1为后面、M=0为横轴线）如图）.第二十步：L=39504.202（显示待求点所计算的“X”坐标）.N=61310.571（显示待求点所计算的“ Y”坐标）.第二十一步

39、：返回第十七步，再根据需要计算各点坐标.第二十二步：J、W伺时输入“ 0”则可以退岀桥、涵坐标计算，并返回到程序第十步状态 二、坐标反算桩号、距离FZ7纵向轴线(L、N)交点(A BN=2时(J、PN=1时(X Y横向轴线M=OM=-1Z8ZH第一步：N? 2（进入坐标反算桩号、距离）.第二步：J? 39524.931（输入已知点X坐标）.第三步：P? 61297.303（输入已知点 Y坐标）.第四步：H? 15000 （输入在本程序数据库计算范围之内的任意桩号）.第五步：H=11000 （显示待求点所计算的里程桩号）.第六步：S=0 （显示待求点所计算的距离左为“-”、右为“ +”）.CAS

40、IO fx-5800P 程序集独立(直线、圆曲线、缓和曲线、卵形曲线、回头曲线坐标计算程序)LQXJS主程序名)“1.HS 说Y” ： “ 2.XYTS” ： “Z”？ Z: “ A”？ A: “ B”？ B: “ C”？ C: “D'? D: “ E”？ E:“F”？ F: “ G'? G:Fix 4: (E-D)- Abs(G-F)宀 I: 90I +n R: 180D - n V: If Z=1:Then Goto 1: IfEnd: Goto 2:Lb1 1: “ H'? H: “S” ？ S: Abs(H-F ) W：Prog“ LQX1 ：“U='：

41、U DMS “X= ”：X/“丫二” :Y/ “Q？ Q:IfQ 工 0: Then Goto 1: IfEnd:“P”？P: “K”？ K: U+K U:Lb1 8:“J”？ J: “W ？ W: “M'? M: O+.5MJ cos( U)-sin( P)+W cos (U+P)- sin( P) L:T+.5MJ sin( U) - sin (P)+W sin (U+P)-sin( P) N: “L=” :L / “N=” : NdfW工 0:Then Goto 8:Else If J 工0:Then Goto 8:IfEnd: Goto 1 : Lb1 2:“J”？ J: “

42、P” ？ P: Prog “LQX2 : F+W H: “H=:出 “ S= ”：S / Goto 2LQX1子程序名1)C+RWVW U: If Uv0:Ten U+36( U:Goto 7: Else If U>360: Ten U-360 U:Goto 7: IfEnd: I U: Goto 7: Lb1 7: A+ X cos(C+RX+VX,0,W)-S sin (U) Q C X:2B+J(si n(C+ RX +VX,0,W)+S cos (U) T:T YLQX2子程序名2)(J-A)cos( C)+(P-B)sin( C) W: 0 S: Lb1 0:Prog“LQX1

43、' : (P-T)sin(U)+(J- C)cos (U) S: If Abs (S)<10 a(-4): ThenGoto 1 : IfEnd: W+SW: 0 S:Goto 0:Lb1 1: 0 S:Prog “LQX1: (P-T) - cos( U) S程序运行示例及说明运行主程序“ DHQXJS三、任意点坐标计算(道路、桥、涵)第一步：Z? 1(进入道路、桥、涵任意点坐标计算).第二步：A? 742589.425(输入起点的“ X”坐标).第三步：B? 463404.387(输入起点的“ 丫坐标).第四步：C? 15° 23' 31.7 (输入起点的方位角).第五步：D?1十340 (输入起点曲率).第六步