辛普生4等分程序.docx

1、辛普生4等分程序辛普生4等分程序 现在不少程序都要收费，有的虽好，但已知数据输入麻烦。本人现免费提供一个给各位测绘同行。一个十五个线形的匝道数据输入只要十分钟。本程序计算原理及程序已刊登在现代测绘2004年第三期上。欢迎试用后提出意见。本人还有更好的版本。另外本人还有可计算平曲线组的程序。（每个平曲线组内含直，一缓，园，二缓，直；一缓可以不等于二缓）.相当于电脑道路软件中的多交点组合曲线。方便测绘同行在野外工作。免得再回办公室去计算.CASIO-fx48匝道坐标放样程序XIANXINGZU（471步）Deg:Defm154:M“1SURU6YIZB7FY”=7=Goto2M=6=Goto5Z“

2、JIZUXIANXIN”:E=Z:T=1:Z5=T:L“1JISUAN2SHUSU”=2=Goto4Fixm:WGHD:Prog“D”:LbI1:Fixm:ABU:Prog“E:S=0:Prog“C”:Prog“A”:IszT:Prog“B”:DszZ:Goto1:Goto5:LbI2:Fixm:IszZ5:FS:F“ZHUANGHAO”S“D”:T=1:LbI3:FZ7T+5=IszT:TGoto3T=T-1:T1=T=1:TE+1=T=EProg“B”:L=F-W:Z520=Goto2:Prog“C”:“ENDOK!”:Z3=Exp-3:N=Z1:O=Z2:K=Z6:L=Z4:Prog“F

3、”:V=J:Prog“G”:J-V+ V :Prog“I”:Goto2:LbI4:Fixm:WGHDABU:Prog“D”:Prog“E”:Prog“A”:W=W+L:IszT:DszZ:Goto4:Goto5:LbI5:Fixm:NOKL:Prog“H”:Z1=N:Z2=O:Z6=K:Z4=L:Goto2(主程序结束)A（子程序）（132步）Z7T+1=G:Z7T+2=H:Z7T+3=U:Z7T+4=D:Z7T+5=W:Z7T+12=W+L:Z7T+6=A:Z7T+7=B:Z7T+11=V-9:Z7T+8=X:Z7T+9=YB（子程序）（87步）G=Z7T+1:H=Z7T+2:W=Z7T+5

4、:A=Z7T+6:B=Z7T+7:D=Z7T+4:C=Z7T+3C（子程序）（16步）P=AbsCC:C=AbsC:V=D+P(AL+L22C(B-A)r+9:Rad:X=L:Y=H+(Sin(D+P(AX+X22C(B-A，X，4:X=L:X=G+(Cos(D+P(AX+X22C(B-A，X，4:Deg:X=X+ScosVY=Y+SsinVD（子程序）（49步）W“QDZHUANGHAO”G“X1”H“Y1”D“A1”E（子程序）（54步）A“R1”：B“R2”：U“LS”:C=U:L=AbsC:A=A-1:B=B-1F（子程序）（3步）Pol(K-N，L-O+Z3G（子程序）（3步）Pol

5、(X-N，Y-O+Z3H（子程序）（53步）N“XYIQI”O“YYIQI”K“XHOUSI”L“YHOUSI”I(子程序）（65步）Ans=Ans+36V=Ans:6FracAns:V“JIAODU”=IntV+.1IntAns+6mFracAnsI“S=”说明：1.本程序根据微积分中牛顿-柯特斯（辛普生,分割指数n=4）法则编制。利用了4800p中的积分功能。可对任何线形（直、缓、园、复、卵、回头曲线等线形的任意组合）上的坐标及边桩坐标进行计算。并在输入置仪点和后视点坐标后及时给出放样所需的角度、边长。计算精度达到毫米级。本程序已知数据不需要人工判断放置在何处，输入简单。2.本程序一次最多

6、可输入二十组连续线形。若要增加或减少，请将主程序(1)中【Defm154:】的154改为（线形组最大数量+2）7即可。【CASIO-48扩展内存不能用字母输入是该机的一个设计缺陷。】3程序中数据自动存储在下列寄存器中：线形组数量Z,线形序号T,起算点坐标X1Z7T+1,Y1Z7T+2,曲线长度LSZ7T+3,起算点的切线方位角Z7T+4，起算点桩号Z7T+5,线形前半径R1Z7T+6,线形后半径R2Z7T+7,欲求点桩号F,边桩距离S，置仪点坐标XZ1,YZ2,后视点XZ6,YZ4，T为线形序号，程序自动计算。4.本程序直接输入各线形已知参数时需输入数据起算点桩号，起算点坐标X、Y;线形两端半

7、径R1、R2;起算点的切线方位角，线形长度。进行计算输入时第一次线形需输入上述-数据，以后的线形只需输入即可。本程序自动计算下一个线形的起算坐标，桩号，方位角。采用两种输入法主要是依据设计院给的参数。计算输入法比较方便。但有时设计院给定的起始数据不精确或各线形数据有矛盾（因为设计院往往是在CAD中作的，虽然误差很小，只有几毫米。但计算几公里后误差还是不能忽视的）。这时，用直接输入法比较准确。5.本程序直线半径输入时因为是无穷大，所以输入直线半径请按【EXP】【99】。子程序C(1)中V=D+P(AL+L22C(B-A)r的r用按键FUNCTION55输入。子程序C(2,3)中(D+的用按键FU

8、NCTION54输入。子程序R(2)中(6mFracAns)的m用按键FUNCTION11输入。6.本程序中边桩夹角为9度，如果要修改，请将子程序C(1)中V=D+P(AL+L22C(B-A)r+9的9改为你需要的数据。7.本程序中设置了直接进行计算功能，只要输入7可直接进入坐标求解、放样功能。还设置了修改置仪点、后视点坐标功能，只要输入6可在输入后进入放样功能。但要确保上述寄存器中数据未被破坏。8.本程序欲求桩号少于或大于起、终点桩号时，线形参数与首末端线形一致时坐标计算结果正确。否则，只能在起、终点桩号内计算。超出该范围桩号计算结果与实际不相符。9.程序输入中曲线右偏，曲线长度LS输入+号

9、。曲线左偏，曲线长度LS输入-号。（切记）1.程序中输入边桩距离S时，线路右边输入+号，线路左边输入-号。11.程序输完指定的几组参数后会自动提示输入置仪点和后视点坐标。按计算输入线形数据方法序号屏幕显示操作说明1XIANXINGZU【EXE】按【FILE】键，启动该文件。21SURU6YIZB7FY?.1【EXE】询问是输入数据。坐标或放样，本例输入1进入输入线形数据。3JIZUXIANXIN?3【EXE】询问是几个线形。本例按三个输入。41JISUAN2SHURU?1【EXE】询问是计算输入或直接输入各线形参数。本例按1计算输入。5QDZHUANGHAO?【EXE】询问起点桩号。输入K+。

10、第一线形按园曲线输入。6X1?29.656【EXE】询问起点坐标，输入坐标，下同。7Y1?751.155【EXE】8A1?7636【EXE】询问输入起点方位角，9R1?511.3【EXE】询问线形前半径，输入。园曲线两端半径相等。1R2?511.3【EXE】询问线形后半径，输入。11LS?-76.818【EXE】询问园曲线长。曲线左偏，输入曲线长时输入负号。切记！12X=286.31【EXE】显示园曲线终点坐标，即下一线形起点坐标。13Y=754.99【EXE】14R1?EXP99【EXE】询问输入第二线形前半径，第二线形按直线输入。15R2?EXP99【EXE】询问输入第二线形后半径，直线两

11、端半径相等且半径无穷大。16LS?-76.81858.489【EXE】询问直线长。曲线右偏，输入曲线长时输入正号。17X=2144.779【EXE】显示直线终点坐标，即下一线形起点坐标。18Y=753.375【EXE】19R1?EXP99【EXE】询问第三线形前半径，第二线形按直缓线输入。前半径无穷大。2R2?8【EXE】询问第三线形后半径，缓和曲线后半径与园曲线相接。19LS?58.4895【EXE】询问缓和曲线长。曲线右偏，输入曲线长时输入正号。2XYIQI?11【EXE】询问输入置仪点的坐标X。19YYIQI?12【EXE】询问输入置仪点的坐标Y。2XHOOSI?13【EXE】询问输入后

12、视点的坐标X。23YHOUSI?14【EXE】询问输入后视点的坐标Y。24ZHANGHAO?38.49【EXE】询问输入欲求点的桩号。25D?【EXE】询问输入欲求点距中线的距离。右边为+，左边为-。26X=247.912【EXE】显示欲求点的坐标X。27Y=754.473【EXE】显示欲求点的坐标Y。ENDOK!坐标计算结束。28JIAODU=289.4933842【EXE】显示置仪点自后视点到欲求点间顺时针方向夹角。已转换成六十进制角度2894933.8。29S=147.393【EXE】显示置仪点到欲求点间平距。3ZHANGHAO?16.5【EXE】继续重复24-29操作。求其他点数据。3

13、1D?-8【EXE】（线路前进方向左边，距中线8米）32X=2169.818【EXE】显示欲求点的坐标X。33Y=744.782【EXE】显示欲求点的坐标Y。34ENDOK!坐标计算结束。35JIAODU=291.5659266【EXE】显示置仪点自后视点右转到欲求点间顺时针方向夹角。已转换成六十进制角度。即2915659.336S=1162.641显示置仪点到欲求点间平距。37【AC】【AC】退出程序。关机休息后重新计算38XIANXINGZU【EXE】按【FILE】键，启动该文件。391SURU6YIZB7FY?167【EXE】询问是输入或直接放样，本例输入7进入放样程序。4ZHANGHA

14、O?76.818【EXE】询问输入欲求点的桩号。41D?【EXE】询问输入欲求点距中线的距离。右边为+，左边为-。42X=286.31【EXE】显示欲求点的坐标X。43Y=754.99【EXE】显示欲求点的坐标Y。ENDOK!坐标计算结束。44JIAODU=29.4242539【EXE】显示置仪点自后视点到欲求点间顺时针夹角。已转换成六十进制角度294253.9。45S=182.87【EXE】显示置仪点到欲求点间平距。继续重复24-29操作。求其他点数据【AC】退出程序。按直接输入方法序号屏幕显示操作说明1XIANXINGZU【EXE】按【FILE】键，启动该文件。21SURU6YIZB7FY

15、?1【EXE】询问是输入或直接计算坐标，本例输入1进入输入线形数据。3JIZUXIANXIN?3【EXE】询问是几个线形。本例按三个输入。41JISUAN2SHURU?2【EXE】询问是计算输入或直接输入各线形参数。本例按2直接输入。5QDZHUANGHAO?254.75【EXE】询问起点桩号。输入K+254.75。第一线形按缓和曲线输入。6X1?2263.472【EXE】询问起点坐标，输入坐标，下同。7Y1?755.364【EXE】8A1?32521【EXE】询问输入起点方位角。9R1?EXP99【EXE】询问线形前半径，输入。缓和曲线前半径无穷大。1R2?8.2【EXE】询问线形后半径，缓

16、和曲线后半径与下一园曲线半径相等。11LS?-5【EXE】询问缓和曲线长。曲线左偏，输入曲线长时输入负号。切记！12QDZHUANGHAO?34.75【EXE】显示第二起点桩号。程序已算好，不必输入。核对一下即可。13X1?【EXE】重复6-12步骤，直到指定的几组线形参数输入完毕。重复2-45步骤，直到按【AC】【AC】退出程序。修改置仪点和后视点坐标后放样方法1XIANXINGZU【EXE】按【FILE】键，启动该文件。21SURU6YIZB7FY?6【EXE】询问是输入起算数据。或输入仪器站、后视点坐标，或放样，本例输入6进入输入仪器站、后视点坐标功能。3XYIQI?1【EXE】询问输入

17、置仪点的坐标X。4YYIQI?2【EXE】询问输入置仪点的坐标Y。5XHOOSI?3【EXE】询问输入后视点的坐标X。6YHOUSI?4【EXE】询问输入后视点的坐标Y。7ZHANGHAO?38.49【EXE】询问输入欲求点的桩号。8D?【EXE】询问输入欲求点距中线的距离。右边为+，左边为-。9X=247.912【EXE】显示欲求点的坐标X。1Y=754.473【EXE】显示欲求点的坐标Y。ENDOK!坐标计算结束。11 L=1000K+M-20: V1=Prog “ZDGC” (如V取1,则将自动计算每隔20M)3Lbl 04 V=1 =L=L+20:Goto 1(要得到10M一个桩的高程

18、,把20改成10就可.或则在运行时取v=0即可，就是每次要输入桩号)5 K，M6 L=1000K+M7 Lbl 18 Prog “G”9Prog “SQX”10Prog “G1”11Prog “PQX-H”12 Goto 0（2为步号，在输入时，不要输，步号间用回车下同）主程序2SQX(程序名)1 T=Abs(R*(I-J)/2)2 R=(I-J)/Abs(I-J)*R3 N=A-T4 W=A+T5 U=(L-A)6 LN=Z=H+JU:Goto 1 7 LW=Z=H+JU+(N-L)2(平方)/(2R):Goto 1 8 LW=Z=H+IU 9 Lbl 1主程序3PQX-H(程序名)1 LS

19、=E=0.02:G=E:Goto 1 (其中为默认横披值)2 L(S+B)=E=0.02-(L-S)/B(0.02+C):G= -E:Abs (E) 0.02=G=0.02: Goto 1:3 L(S+Q-B)=E= -C:G= -E:Goto 1 4 L(S+Q)=E=0.02-(S+Q-L)/B(0.02+C):G= -E:Abs (E) 0.02=G=0.02: Goto 1:5 L(S+Q)=E=0.02:G=E6 Lbl 17 D（输入边桩距中桩距离）8 D=99=Prog “ZDGC” 9 F=1=O=Z-(D-1)E:P=Z-(D-1)G: F=-1=O=Z-(D-1)G:P=Z

20、-(D-1)E(D-1)中，把改成0则中央分隔带为M.默认中央分隔带为2M)10 Prog “XS1”子程序1（竖曲线引导程序）G(程序名)0 L31182.853=“INPUT LK31+182.853”:Prog“ZDGC”：1 L 31495=Prog“1”:Goto 1: （31495为第一竖曲线终点桩号）2 L31960=Prog“2”:Goto 1: 3 L32769.975=Prog“3”:Goto 1: 4 L33380.9=Prog“4”:Goto 1: 5 L33875=Prog“5”:Goto 1: 6 L34425=Prog“6”:Goto 1: 7 L35125=Prog