辛普生4等分程序.docx

辛普生4等分程序.docx

《辛普生4等分程序.docx》由会员分享，可在线阅读，更多相关《辛普生4等分程序.docx（28页珍藏版）》请在冰豆网上搜索。

辛普生4等分程序

辛普生4等分程序

现在不少程序都要收费，有的虽好，但已知数据输入麻烦。

本人现免费提供一个给各位测绘同行。

一个十五个线形的匝道数据输入只要十分钟。

本程序计算原理及程序已刊登在《现代测绘》2004年第三期上。

欢迎试用后提出意见。

本人还有更好的版本。

另外本人还有可计算平曲线组的程序。

（每个平曲线组内含直，一缓，园，二缓，直；一缓可以不等于二缓）.相当于电脑道路软件中的多交点组合曲线。

方便测绘同行在野外工作。

免得再回办公室去计算.

 CASIO-fx48ØØ 匝道坐标放样程序                                    

 XIAN XING ZU  （471步）

 Deg:

 Defm 154:

 M“1SURU 6YIZB 7FY”=7 => Goto 2⊿ M=6 => Goto 5⊿ Z“JIZU XIANXIN”:

 E=Z:

 T=1:

 Z［5］= T :

 L“1JISUAN 2SHUSU”=2 => Goto 4⊿ Fixm :

｛WGHD｝:

 Prog“D”:

 LbI 1:

 Fixm :

｛ABU｝:

Prog“E:

 S=0:

  Prog“C”:

Prog“A”:

Isz T:

 Prog“B”:

Dsz Z:

 Goto 1:

  Goto 5 :

LbI 2:

 Fixm :

 Isz Z［5］:

｛FS｝:

 F“ZHUANG HAO”S“D”:

 T=1:

 LbI 3:

 F≥Z［7T+5］=> Isz T:

 TGoto 3 ⊿⊿ T=T-1:

 T≤1 => T=1:

 T＞E+1=> T=E ⊿⊿ Prog“B”:

 L=F-W:

 Z［5］>20=>Goto 2:

 Prog“C”:

 “END OK!

 ”:

Z［3］= Exp -3Ø:

 N=Z［1］:

 O=Z［2］:

 K=Z［6］:

 L=Z［4］:

Prog“F”:

 V=J:

 Prog“G”:

J-V+V:

 Prog“I”:

 Goto 2:

LbI 4:

Fixm :

｛WGHDABU｝:

 Prog“D”:

 Prog“E”:

 Prog“A”:

 W=W+L:

 Isz T:

 Dsz Z:

 Goto4:

Goto 5:

LbI 5:

Fixm :

｛NOKL｝:

 Prog“H”:

  Z［1］=N:

 Z［2］=O:

 Z［6］=K:

 Z［4］=L:

 Goto 2

 （主程序结束）

A     （子程序）（132步）

 Z［7T+1］=G:

 Z［7T+2］=H:

 Z［7T+3］=U:

 Z［7T+4］=D:

 Z［7T+5］=W:

 Z［7T+12］=W+L:

 Z［7T+6］=A:

 Z［7T+7］=B:

 Z［7T+11］=V-9Ø:

 Z［7T+8］=X:

 Z［7T+9］=Y 

B     （子程序）（87步）

 G =Z［7T+1］:

 H=Z［7T+2］:

 W=Z［7T+5］:

 A=Z［7T+6］:

 B=Z［7T+7］:

 D=Z［7T+4］:

 C=Z［7T+3］

C      （子程序）（16Ø步）

 P= Abs C÷C :

 C= Abs C:

 V =D + P（ AL + L2 ÷2C×（B-A））r +9Ø :

 Rad :

 X=L:

 Y = H + ∫（ （Sin （ Dº+ P（ AX + X2÷2C×（B-A ，Ø ，X ，4 :

 X=L:

 X = G + ∫（ （Cos （ Dº+ P（ AX + X2÷2C×（B-A ，Ø ，X ，4 :

 Deg :

 X = X + S cos V ◢ Y = Y + S sin V ◢ 

D     （子程序）（49步）

 W“QD ZHUANG HAO”G“X1” H“Y1” D“A1 >”

E     （子程序）（54步）

 A “R1”：

 B “R2”：

 U “LS”:

 C=U :

 L= Abs C:

 A = A-1 :

 B = B-1 

F     （子程序）（3Ø步）

 Pol（ K-N ，L-O + Z[3] 

G     （子程序）（3Ø步）

 Pol（ X-N，Y-O + Z[3]

H     （子程序）（53步）

 N“X YIQI”O“Y YIQI”K“X HOUSI”L“Y HOUSI”

I     （子程序）（65步）

 Ans＜Ø => Ans + 36Ø ⊿

 V= Ans:

 6Ø Frac Ans:

V “＜JIAO DU”=Int V +.Ø1 Int Ans + 6m Frac Ans ◢ I“S=” ◢

 说明：

1. 本程序根据微积分中牛顿-柯特斯（辛普生,分割指数n=4）法则编制。

利用了4800p中的积分功能。

可对任何线形（直、缓、园、复、卵、回头曲线等线形的任意组合）上的坐标及边桩坐标进行计算。

并在输入置仪点和后视点坐标后及时给出放样所需的角度、边长。

计算精度达到毫米级。

本程序已知数据不需要人工判断放置在何处，输入简单。

2. 本程序一次最多可输入二十组连续线形。

若要增加或减少，请将主程序

（1）中【Defm 154:

】的154改为（线形组最大数量+2）×7即可。

【CASIO-48ØØ扩展内存不能用字母输入是该机的一个设计缺陷。

】

3． 程序中数据自动存储在下列寄存器中：

线形组数量→Z, 线形序号→T, 起算点坐标X1→Z［7T+1］,Y1→Z［7T+2］, 曲线长度LS→Z［7T+3］, 起算点的切线方位角→Z［7T+4］，起算点桩号→Z［7T+5］,线形前半径R1→Z［7T+6］, 线形后半径R2→Z［7T+7］, 欲求点桩号→F , 边桩距离→S ，置仪点坐标X→Z［1］, Y→Z［2］, 后视点X→Z［6］, Y→Z［4］，T为线形序号，程序自动计算。

4. 本程序直接输入各线形已知参数时需输入数据 ①起算点桩号，②起算点坐标X、Y; ③线形两端半径R1、R2; ④起算点的切线方位角，⑤线形长度。

       进行计算输入时第一次线形需输入上述①-⑤数据，以后的线形只需输入③⑤即可。

本程序自动计算下一个线形的起算坐标，桩号，方位角。

       采用两种输入法主要是依据设计院给的参数。

计算输入法比较方便。

但有时设计院给定的起始数据不精确或各线形数据有矛盾（因为设计院往往是在CAD中作的，虽然误差很小，只有几毫米。

但计算几公里后误差还是不能忽视的）。

这时，用直接输入法比较准确。

5. 本程序直线半径输入时因为是无穷大，所以输入直线半径请按【EXP】【99】。

 子程序C

（1）中V =D + P（ AL + L2 ÷2C×（B-A））r 的 r   用按键[FUNCTION] [5] [5]输入。

 子程序C （2,3）中（ Dº+  的 º  用按键[FUNCTION] [5] [4]输入。

 子程序R 

（2）中（ 6m Frac Ans ◢ ） 的 m  用按键[FUNCTION] [1][﹀][﹀][﹀] [1]输入。

6. 本程序中边桩夹角为9Ø度，如果要修改，请将子程序C

（1）中V =D + P（ AL + L2 ÷2C×（B-A））r +9Ø 的9Ø改为你需要的数据。

7. 本程序中设置了直接进行计算功能，只要输入7可直接进入坐标求解、放样功能。

还设置了修改置仪点、后视点坐标功能，只要输入6可在输入后进入放样功能。

但要确保上述寄存器中数据未被破坏。

8. 本程序欲求桩号少于或大于起、终点桩号时，线形参数与首末端线形一致时坐标计算结果正确。

否则，只能在起、终点桩号内计算。

超出该范围桩号计算结果与实际不相符。

9. 程序输入中曲线右偏，曲线长度LS输入+号。

曲线左偏，曲线长度LS输入-号。

（切记）

1Ø. 程序中输入边桩距离S时，线路右边输入+号，线路左边输入-号。

11. 程序输完指定的几组参数后会自动提示输入置仪点和后视点坐标。

             按计算输入线形数据方法

序号       屏 幕 显 示          操    作                说     明

1 XIAN XING ZU                   【EXE】 按【FILE】键，启动该文件。

2 1SURU 6YIZB 7FY ?

  Ø.ØØØ      1  【EXE】 询问是输入数据。

坐标或放样，本例输入1进入输入线形数据。

3 JIZU XIANXIN?

               3  【EXE】 询问是几个线形。

本例按三个输入。

4 1JISUAN 2SHURU?

    Ø       1  【EXE】 询问是计算输入或直接输入各线形参数。

本例按1计算输入。

5 QD ZHUANG HAO?

     Ø       Ø  【EXE】 询问起点桩号。

输入KØ+ØØØ。

第一线形按园曲线输入。

6 X1?

          Ø         2ØØ9.656 【EXE】 询问起点坐标，输入坐标，下同。

7 Y1?

          Ø          751.155  【EXE】 

8 A1＞?

                  7ºØ6´36 【EXE】 询问输入起点方位角，

9 R1?

                  511.Ø3  【EXE】 询问线形前半径，输入。

园曲线两端半径相等。

1Ø R2?

                  511.Ø3  【EXE】 询问线形后半径，输入。

11 LS?

                  -76.818 【EXE】 询问园曲线长。

曲线左偏，输入曲线长时输入负号。

切记！

12 X=                       2Ø86.31Ø【EXE】 显示园曲线终点坐标，即下一线形起点坐标。

13 Y=                      754.9Ø9  【EXE】 

14 R1?

     Ø        EXP  99   【EXE】 询问输入第二线形前半径，第二线形按直线输入。

15 R2?

     Ø        EXP  99   【EXE】 询问输入第二线形后半径，直线两端半径相等且半径无穷大。

16 LS?

   -76.818         58.489   【EXE】 询问直线长。

曲线右偏，输入曲线长时输入正号。

17 X=                    2144.779   【EXE】 显示直线终点坐标，即下一线形起点坐标。

18 Y=                      753.375  【EXE】 

19 R1?

    Ø        EXP 99    【EXE】 询问第三线形前半径，第二线形按直缓线输入。

前半径无穷大。

2Ø R2?

    Ø            8ØØ   【EXE】 询问第三线形后半径，缓和曲线后半径与园曲线相接。

19 LS?

   58.489              5Ø  【EXE】 询问缓和曲线长。

曲线右偏，输入曲线长时输入正号。

2Ø X YIQI?

                   11ØØ   【EXE】 询问输入置仪点的坐标X。

19 Y YIQI?

                   12ØØ   【EXE】 询问输入置仪点的坐标Y。

2Ø X HOOSI?

           13ØØ   【EXE】 询问输入后视点的坐标X。

23 Y HOUSI?

           14ØØ   【EXE】 询问输入后视点的坐标Y。

24 ZHANG HAO?

         38.4Ø9   【EXE】 询问输入欲求点的桩号。

25 D?

                      Ø   【EXE】 询问输入欲求点距中线的距离。

右边为 +，左边为 -。

26 X=     2Ø47.912                  【EXE】 显示欲求点的坐标X。

27 Y=     754.473                  【EXE】 显示欲求点的坐标Y。

         END OK!

                         坐标计算结束。

28 

已转换成六十进制角度289º49´33.8"。

29 S=           1Ø47.393          【EXE】 显示置仪点到欲求点间平距。

3Ø ZHANG HAO?

 16Ø.5            【EXE】 继续重复24---29操作。

求其他点数据。

31 D?

                    -8    【EXE】 （线路前进方向左边，距中线8米）

32 X=           2169.818          【EXE】 显示欲求点的坐标X。

33 Y=            744.782          【EXE】 显示欲求点的坐标Y。

34 END OK!

  坐标计算结束。

35 

已转换成六十进制角度。

即291º56´59.3"

36 S=           1162.641                   显示置仪点到欲求点间平距。

37                【AC】   【AC】                 退出程序。

         关机休息后重新计算  

38 XIAN XING ZU                   【EXE】 按【FILE】键，启动该文件。

39 1SURU 6YIZB 7FY ?

  16      7   【EXE】 询问是输入或直接放样，本例输入7进入放样程序。

4Ø ZHANG HAO?

         76.818   【EXE】 询问输入欲求点的桩号。

41 D?

                    Ø     【EXE】 询问输入欲求点距中线的距离。

右边为 +，左边为 -。

42 X=     2Ø86.31Ø                  【EXE】 显示欲求点的坐标X。

43 Y=     754.9Ø9                  【EXE】 显示欲求点的坐标Y。

 END OK!

  坐标计算结束。

44 

已转换成六十进制角度29غ42´53.9"。

45 S=           1Ø82.Ø87          【EXE】 显示置仪点到欲求点间平距。

   继续重复24---29操作。

求其他点数据

                【AC】               退出程序。

 按直接输入方法  

序号      屏 幕 显 示            操    作                      说     明

1 XIAN XING ZU                   【EXE】 按【FILE】键，启动该文件。

2 1SURU 6YIZB 7FY ?

    Ø       1  【EXE】 询问是输入或直接计算坐标，本例输入1进入输入线形数据。

3 JIZU XIANXIN?

               3  【EXE】 询问是几个线形。

本例按三个输入。

4 1JISUAN 2SHURU?

    Ø       2  【EXE】 询问是计算输入或直接输入各线形参数。

本例按2直接输入。

5 QD ZHUANG HAO?

     Ø 254.Ø75  【EXE】 询问起点桩号。

输入KØ+254.Ø75。

第一线形按缓和曲线输入。

6 X1?

                Ø 2263.472 【EXE】 询问起点坐标，输入坐标，下同。

7 Y1?

                Ø  755.364  【EXE】 

8 A1＞?

              3º25´21"    【EXE】 询问输入起点方位角。

9 R1?

                EXP  99   【EXE】 询问线形前半径，输入。

缓和曲线前半径无穷大。

1Ø R2?

                8ØØ.Ø2    【EXE】 询问线形后半径，缓和曲线后半径与下一园曲线半径相等。

11 LS?

                     -5Ø  【EXE】 询问缓和曲线长。

曲线左偏，输入曲线长时输入负号。

切记！

12 QD ZHUANG HAO?

         3Ø4.Ø75   【EXE】 显示第二起点桩号。

程序已算好，不必输入。

核对一下即可。

13 X1?

                          Ø   【EXE】 重复6-12步骤，直到指定的几组线形参数输入完毕。

   重复2Ø-45步骤，直到按【AC】【AC】退出程序。

 修改置仪点和后视点坐标后放样方法 

1 XIAN XING ZU                   【EXE】 按【FILE】键，启动该文件。

2 1SURU 6YIZB 7FY ?

    Ø      6   【EXE】 询问是输入起算数据。

或输入仪器站、后视点坐标，或放样， 本例输入6进入输入仪器站、后视点坐标功能。

3 X YIQI?

                    1ØØ   【EXE】 询问输入置仪点的坐标X。

4 Y YIQI?

                    2ØØ   【EXE】 询问输入置仪点的坐标Y。

5 X HOOSI?

            3ØØ   【EXE】 询问输入后视点的坐标X。

6 Y HOUSI?

            4ØØ   【EXE】 询问输入后视点的坐标Y。

7 ZHANG HAO?

         38.4Ø9   【EXE】 询问输入欲求点的桩号。

8 D?

                      Ø   【EXE】 询问输入欲求点距中线的距离。

右边为 +，左边为 -。

9 X=     2Ø47.912                  【EXE】 显示欲求点的坐标X。

1Ø Y=      754.473                  【EXE】 显示欲求点的坐标Y。

 END OK!

  坐标计算结束。

11 

已转换成六十进制角度。

为33غ53´2Ø.4"

12 S=           2Ø25.291          【EXE】 显示置仪点到欲求点间平距。

13   重复7-11步骤，直到按【AC】【AC】退出程序。

       现在4850p只要380元。

程序步骤28500，增加了十进制转换六十进制功能，不必再为此编一个程序了。

公路高程计算程序计算任一点中边桩高程（CASIO4800）

            文章来源：

中华工程网收集点击数：

4431更新时间：

2007-7-20【字体：

小大】

            ZDGC（程序名）

            １V“Change0-OnlyOne:

1-AutoNextOne”（回车）

            ２V=1=>L=1000K+M-20:

≠>V﹥1=>Prog“ZDGC”⊿⊿（如V取1,则将自动计算每隔20M）

            3　Lbl0

            4V=1=>L=L+20:

Goto

1⊿（要得到10M一个桩的高程,把20改成10就可.或则在运行时取v=0即可，就是每次要输入桩号）

            5{K，M}

            6L=1000K+M

            7Lbl1

            8Prog“G”

            9  Prog“SQX”

            10  Prog“G1”

            11  Prog“PQX-H”

            12Goto0

            （１－１2为步号，在输入４８００时，不要输，步号间用回车．下同）

            主程序2

            SQX（程序名）

            1T=Abs（R*（I-J）/2）

            2R=（（I-J）/Abs（I-J））*R

            3N=A-T

            4W=A+T

            5U=（L-A）

            6L≤N=>Z=H+JU:

Goto1⊿

            7L≤W=>Z=H+JU+（N-L）2（平方）/（2R）:

Goto1⊿

            8L>W=>Z=H+IU⊿

            9Lbl1

            主程序3

            PQX-H（程序名）

            1L≤S=>E=0.02:

G=E:

Goto1⊿（其中Ｅ＝．０２为默认横披值）

            2L≤（S+B）=>E=0.02-（L-S）/B×（0.02+C）:

G=-E:

Abs（E）≤0.02=>G=0.02:

            ⊿Goto1:

⊿

            3L≤（S+Q-B）=>E=-C:

G=-E:

Goto1⊿

            4L≤（S+Q）=>E=0.02-（S+Q-L）/B×（0.02+C）:

G=-E:

Abs（E）≤0.02=>G=0.02:

            ⊿Goto1:

⊿

            5L＞（S+Q）=>E=0.02:

G=E⊿

            6Lbl1

            7{D}（输入边桩距中桩距离）

            8D=99=>Prog“ZDGC”⊿

            9F=1=>O=Z-（D-1）E:

P=Z-（D-1）G:

            ≠>F=-1=>O=Z-（D-1）G:

P=Z-（D-1）E⊿⊿（（D-1）中，把１改成0．５则中央分隔带为１M.默认中央分隔带为2M）

            10Prog“XS1”

            子程序1（竖曲线引导程序）

            G（程序名）

            0L﹤31182.853=>“INPUTL≥K31+182.853”:

Prog“ZDGC”：

⊿

            1L≤31495=>Prog“1”:

Goto1:

⊿（31495为第一竖曲线终点桩号）

            2L≤31960=>Prog“2”:

Goto1:

⊿

            3L≤32769.975=>Prog“3”:

Goto1:

⊿

            4L≤33380.9=>Prog“4”:

Goto1:

⊿

            5L≤33875=>Prog“5”:

Goto1:

⊿

            6L≤34425=>Prog“6”:

Goto1:

⊿

            7L≤35125=>Prog