常用casio4500施工测量计算程序

１1、坐标反算程序程序:L1Fixm：{C，D}L2pol(C-A,D-B)◢L3ProgD-M-S原理:如下图,已知两点P1、P2的坐标分别为P1（A，B）、P2（C，D），则P1P2间的距离和方位角分别为：S=22B)-(DA)-(CT=arctg[(D-B)/(C-A)](C,D)(A,B)说明：起始点（A，B）应预置；调用子程序名称应与子程序名称严格一致。D-M-S：度·分·秒显示子程序L1W0W=W+360Fix4L2lntW+0.01lnt(60FracW)+0.006Frac(60FracW)◢２2、坐标正算程序程序：L1Fixm:{S,T}L2C=A+Rec(S,T)◢D=B+W◢原理：如下图，P1为已知点，测得未知点P2与P1间的距离S及方位角T，则P2点的坐标为：C=A+S×COS（T）D=B+S×SIN(T）(C,D)ST(A,B)说明：1、起始点（A，B）应预置。3、后方交会程序程序：L1Fixm:{OPQ}L2I=-O+P:J=Q-PL3K=(A-C)+(B-D)/tanIL4L=(D-B)+(A-C)/tanIL5M=(C-E)+(F-D)/tanJL6N=(F-D)+(E-C)/tanJL7U=(K+M)/(L+N)３L8X=C+(K-UL)/(1+U2)◢L9Y=D+U(K-UL)/(1+U2)◢原理：如下图，已知1、2、3点的坐标分别为（A，B）、（C，D）、（E，F），在未知点PP上设站，测得1、2、3点的方向值分别为O、P、Q，则PP点的坐标为：X=C+（K-U×L）/（1+U2）Y=D+U×（K-U×L）/（1+U2）其中U=（K+M）/（L+N）K=（A-C）+（B-D）/tgIL=（D-B）+（A-C）/tgIM=（C-E）+（F-D）/tgJN=（F-D）+（E-C）/tgJI=P-OJ=Q-P2(C,D)3(E,F)1(A,B)PQOPP(X,Y)说明：O、P、Q为顺时针观测三个合适已知点1、2、3的方向值。三个已知点的坐标值(A,B)、(C,D)、(E,F)应先输入。三个方向组成的两个夹角均不得为90°。４4、前方交会程序（1）程序：L1Fixm:{MN}L2X=(C/tanM+A/tanN+D-B)/(1/tanM+1/tanN)◢L3Y=(D/tanM+B/tanN+A-C)/(1/tanM+1/tanN)原理：如下图，已知P1、P2两点的坐标分别为（A,B）、(C,D)，在P1、P2上分别设站，测得未知点P与已知边所组成的夹角分别为M、N，则未知点P的坐标为：X=ctgNctgMctgNActgMCBDY=ctgNctgMctgNctgMDABCMN说明：两个已知点与交会点应按逆时针方向组成图形。两个已知点坐标（A,B）、(C,D)应先输入。M、N≠0°，90°及其整数倍。P1(A,B)P2(C,D)P(X,Y)５5、前方交会程序（2）程序：L1Fixm:{MN}L2pol(A-C,B-D)L3S=VSINM/SIN(M+N):T=W+NL4Rec(S,T):X=C+V◢Y=D+W◢原理：如下图，已知两点P1、P2坐标分别为（A，B）、（C，D），按前方交会法测得未知点P与已知边所组成的夹角分别为M、N，则P与P2间的距离S和方位角T分别为：S=P1P2×SIN（M）/SIN（M+N）T=N+TP2P1TP2P1为P2P1间的方位角则P点的坐标为：X=C+S×COS（T）Y=D+S×SIN（T）P(X,Y)P1(A,B)△MN△P2(C,D)说明：两个已知点与交会点应按逆时针方向组成图形。两个已知点坐标（A,B）、(C,D)应先输入。M、N可以为90°。若认为有时M、N为夹角不方便，可以用坐标方位角代替，这时只要将程序变为如下：L1Fixm:{MN}L2pol(A-C,B-D):M=W+180-M:M0M=M+360N0N=N+360L3N=N-W:L4S=VSINM/SIN(M+N):T=W+NL5Rec(S,T):X=C+V◢Y=D+W◢如下图：P(X,Y)MN(A,B)△△(C,D)６6、侧方交会程序程序：L1Fixm:{MNP}L2M=N:N=N:P=PL3pol(C-A,D-B)L4M≠0S=VSINM/SINP:T=W-M-PL5Rec(S,T):X=C+V◢Y=D+W◢L6N≠0S=VSINN/SINP:T=W+N+P-180L7Rec(S,T):X=A+V◢Y=B+W◢原理：如下图，三角形P1P2P3中P1、P2两已知点坐标分别为（A，B）、（C，D），测得角P及角M（或角N），可求得未知点P3的坐标X，Y。若测角P及角M，则P2P3间的距离S及方位角T分别为：S=P1P2×SIN（M）/SIN（P）T=180-M-P+Tp1p2Tp1p2为P1、P2间的方位角则有：X=C+S×COS（T），Y=D+S×SIN（T）若测角P及角N，则P1P2间的距离S及方位角T分别为：S=P1P2×SIN（N）/SIN（P）；T=M+P+Tp1p2-180Tp1p2为P1、P2间的方位角，则有：X=A+S×COS（T）；Y=B+S×SIN（T）P3(X,Y)PP1(A,B)△MN△P2(C,D)说明：两个已知点与交会点应按逆时针方向组成图形。未知点上测角为P，已知点P1或P2上没有测角时，M或N的值输0。两个已知点的坐标（A，B）、（C，D）应先输入。不采用条件语句时程序可缩改如下：L1Fixm:{MNP}L2pol(C-A,D-B)：S=V（SINM+SINN）/SINPL3T=W-M+N+(NM-MN)P-180•NML4X=Rec(S,T)+ANM+CMN◢L5Y=W+BNM+DMN◢7、边角后方交会程序程序：７L1Fixm:{CDKSP}L2pol(C-A,D-B)L3Q=90(1-K)+KSIN-1(SSINP/V)L4T=W+180-P-QL5Rec(S,T):X=A+V◢Y=B+W◢原理：P1(A,B)P3(X,Y)SPSPP3(X,Y)P2(C,D)说明：测边的已知点坐标（A,B）应先输入。а角、P角为以测边方向为起始方向，顺时针观测另外一个已知点的右角。测边对角为锐角时K=1，测边对角为钝角时K=-1。上述情况适于：只有两个已知点；两个已知点不通视或其中某一个很难设站的情况。8、计算面积程序程序：L1S=S:N=N:A=A:B=BL2Lbl2L3Fixm:{XY}L4S=S+(X+A)(Y-B)/2L5A=X:B=YL6N=1Goto8L7N=N-1:Goto2L8Lbl8:S◢８X(X2,Y2)(Xi,Yi)(X1,Y1)(Xi+1,Yi+1)Y(XN,YN)原理：如下图，多边形的顶点分别为(X1,Y1)、(X2,Y2)、(x3,y3)、...(xi,yi)、(xi+1,yi+1)、...(XN,YN)。则多边形的面积为：S=))((11121NiiiiiYYXX其中：XN+1=X1YN+1=Y1说明：此程序适于封闭多边形，利用多边形顶点坐标求面积；N为多边形的边数（即顶点数）；第一点（起始点）的坐标为（A，B）首先输入，按固定方向依次输入各顶点的坐标X、Y，最后一点的坐标仍输起始点坐标。9、方量表计算程序程序：L1V=V:I=IL2Lbl2L3Fixm:{DJ}L4D=0Goto7L5S=(I+J)/2◢G=SD◢９L6I=J:V=V+G:Goto2L7Lbl7:V◢原理：如下表，按断面法计算方量公式为：断面桩号断面面积平均面积断面间距体积IS◢DGJS◢DGJ说明：程序开始，对V清零或对上次方量保留，I为第一个断面的面积。D为下一个断面到上一个断面的距离，J为下一个断面面积，先D后J。最后以D=0结束，并显示总方量。10、圆曲线放样程序（1）程序：L1Fixm:{LS}L2T=P+KL×180/π/RL3Rec(S+R,T):C=X+V:D=Y+WL4Pol(C-A,D-B)◢L5ProgD-M-S)2(DJIV１０原理：如下图，已知圆曲线的圆心坐标为（X，Y），半径为R，起算方位角为P。在点（A，B）上设站，放样一点（C，D），该点距圆弧起算方向弧长为L，距弧线内侧（或外侧）S，则该点的坐标为：知道测站坐标可反算两点间的距离及方位角。PO（X，Y）LF（C，D）SE（A，B）说明：运行前预置：圆心坐标（X，Y），测站坐标（A，B），起算方位角（离心方向）P，弧长L的计算方向顺时针取K=1，逆时针时取K=-1，基准圆半径R。运行变量：L为从起算方向开始沿R半径的弧长，S为半径的伸（S0）,缩（S0）量。11、圆曲线放样程序（2）程序：L1Fixm:{ZS}L2T=P+KAbs(Z-G)×180/π/RL3Rec(S+R,T):C=X+V:D=Y+WL4Pol(C-A,D-B)◢L5ProgD-M-S原理：)180()()180()(RLKPSINSRYDRLKPCOSSRXC１１如下图，已知圆曲线的圆心坐标为（X，Y），半径为R，起算方位角为P，起算桩号为G。在点E（A，B）上设站，放样一点F（C，D），该点的桩号为Z，距弧线外内侧（或外侧）S，则该点的坐标为：用坐标反算可求出EF两点间的距离及方位角。O（X，Y）PRF（C，D）ZSE（A，B）G说明：运行前预置：圆心坐标（X，Y），测站坐标（A，B），起算方位角（离心方向）P，起算桩号G，基准圆半径R（确定桩号的圆半径），桩号增加方向为顺时针取K=1，逆时针时取K=-1；运行变量：Z为要放样的桩号线，S为对应的桩号线上相对于基准圆半径的伸（S0）,缩（S0）量；作为放样的测站点可以时任意的，但桩号只能对应于基准圆（半径为R）来归算。12、椭圆及其平行线放样程序程序：L1A：B：E：F：I：J：UL2lbl0:Fixm:{XKS}N=J+KSJ/AbsJL3X=0M=0：)180()()180()(RGZKSINSRYDRGZKCOSSRXC１２L4X≠0Y=J(1-X2/I2)：Q=I2Y/（J2X）L5T=tan-1Q+180-(K+X/AbsX)×90L6Rec(S,T):M=X+V:N=Y+WL7C=E+McosU-NsinUL8D=F+MsinU+NcosUL9Pol(C-A,D-B)◢L10ProgD-M-SL11Goto0原理：如图，P1（X，Y）为椭圆x2/a2+y2/b2=1上的一点，则距离椭圆S的椭圆平行线上的对应点P2的坐标（M，N）为:M=X+S×COS（T）N=Y+S×SIN（T）其中再将M，N转换为施工坐标系中的坐标即可。XP2（M，N）SP1（X，Y）Y说明：1、椭圆圆曲线标准方程为：x2/a2+y2/b2=1,程序运行前须预置的量为：（E,F）椭圆中心在施工坐标系中的坐标；（A，B）测站的施工坐标系坐标；U为椭圆坐标系相对于施工坐标系的旋转角；I=a;J=±b,第一、二象限曲线段J=b，第三、四象限曲线段J=-b；K=±1，平行曲线位于椭圆曲线外侧时取K=1，内侧取K=-1。)YXabarctg(90T22１３2、运行变量为：X∈[-a,a];S为平行曲线到基准曲线间的距离。13、非水平坡脚（顶）、扭曲坡面放样程序（1）程序：L1Fixm:{SH}L2Rec(S,N):C=E+V:D=F+WL3L=Abs(Ssin(N-W)):I=G+L/PL4J=Abs(H-I)×(O+QL):Rec(J,M)L5C=C+V:D=D+W:Pol(C-A,D-B)◢L6ProgD-M-S说明：(E，F，G)为起算点（起算主坡方向线与起算坡脚，坡顶线的交点）的三维坐标设计值。G时设计高程。M为起算主坡度线方向值，N为起算坡脚坡顶线的方向值，两者均由起算点离射。1：0为起算主坡度的坡比，无论上、下均取正值。1：P为起算坡脚（顶）线垂直于主坡方向的坡比，规定上升时P0,下降时P0。坡顶（脚）水平时取P≥S×103。Q