附件(点击下载) ： (在线预览)
程序源 2009.04.08校对件
2009.04.09附算例
已有朋友成功，请与成功的朋友多沟通 。2009.04.12
加桥台转轴计算，2009.04.19
加涵通转轴计算2009.05.05
附可替换子程序Q.4(5点法加速）2009.05.23
重新上转附件，附图，算例，及平分中矢法布桩坐标计算方法，2009.07.21 CASIO-----fx5800 (任意平曲线极坐标与直角坐标放样, 求里程及边距)
湖南津市 赵济汉
主程序PQX
Deg : Fix 3 : “XZ→0:YZ→1”？A : If A = 1: Then Goto 1 : IfEnd
If A = 0 : Then “ BS→0:XY→1:AND→2:DK→3:L(I)→4 ” ？O : IfEnd If O = 4: Then Goto 1 : IfEnd
If O = 3: Then Prog “ F.2 ” : If X= 0 : Then Goto 1 : IfEnd : IfEnd
If O ≠ 1: Then “ X1 ” ？X : “ Y1” ？Y : X→Z［11］: Y→Z［12］: “ X2 ” ？P : “ Y2” ？Q : Pol( P －X , Q －Y) : If J﹤0 : Then J + 360→J : IfEnd : Cls : “ S12= ” : Locate 6 ,1, I : “ B12= ” : J ►DMS◣
1→B : IfEnd
If O = 1: Then “XY(0) →0: ≠﹥1” ? B: IfEnd
Lbl 1 : If A = 1 And O = 3 :Then Prog “ F.2 ” : IfEnd : Prog “ Q.1 ”
子程序PPQX （给定综合曲线属性）
“ ZQX →0: *** →1: *** →2 ******* ” ？C
If C=0: Then 100→Z［98］: n0 →Z［99］: IfEnd
If C=1: Then m1→Z［98］: n1 →Z［99］: IfEnd
If C=2: Then m2→Z［98］: n2 →Z［99］: IfEnd
*******
Prog “ PQX ”
其中： mi为第i 条综合曲线矩阵变量起始单元，(第1条综合曲线矩阵变量起始单元应自100开始，以便其他程序临时使用100以前的矩阵变量，综合曲线中每一基本单元需8个矩阵存储单元) ； ni 为第i 条综合曲线中基本单元数，每增加一个条综合曲线则增加一行If C=*******→Z［99］: IfEnd 语句。
还应注特别意，原有程序中的矩阵变量定维语句DimZ 应取消，只要原有程序中的矩阵变量定维不大于100，原有程序的运行不受影响。
子程序Q . 5 (建立数据库)
Deg : Fix 5
0→I : “ N ” ? W : W →Z［99］: “Z［MIN ］”? V : V→Z［98］
Lbl 1 : “NO.I= ”: I ◣
Lbl C : V + 8 I→G
“→DK ”？A : A→Z［G ］
“→X ”？B : B→Z［G + 1］
“→Y ”？C : C→Z［G +2］
“→AT ”？D : D→Z［G +3］
“→(－1Y+1) ”？K : K→Z［G +4］
“→L ”？E : E→Z［G +5］
“→R1 ”？F : F→Z［G +6］
“→R2”？O : O →Z［G +7］
Lbl 2 : “JX→1: NO→0 : OK→2”? M : If M = 2 : Then Goto 3 : IfEnd : I
f M = 1: Then Goto C : IfEnd
I + 1→I : Goto 1
Lbl 3 : “END ”
子程序Q . 6
Lbl 1 : “S ” ? W : If W = 0 : Then Goto 2 : IfEnd
I －W→W : Cls : “⊿S= ” : Locate 5 , 1 , W : “ BP= ” : J ►DMS◣
Goto 1
Lbl 2 : Cls : “ OK ”
子程序Q.1
If X = 0 And O = 3 : Then Goto 2 : IfEnd
Lbl 1 : Z ［98］→ I : Z ［99］→ J : “ →DKP ” ? S
If S = －1: Then Prog “JH . ZJ ” : Goto 1 : IfEnd （加桥台转轴时） If S = －2: Then Prog “ZH . 2 ” : Goto 1 : IfEnd （加涵通转轴时） If S = 0: Then Goto 2 : IfEnd : If S﹤Z ［ I ］: Then “ DKP Z ［ I + 8 J ］: Then “ DKP>MAX ” ◣
Goto 1 : IfEnd
If O ≠ 4 And B ≠ 0 : Then “ →B ” ? V : V : Prog “ J ” : T→U : “→(－＋)D÷2 ” ? N : Else 0 →N : IfEnd
Prog “ Q.2 ” : Goto 1
Lbl 2 : “ END ”
子程序Q . 2
If O ≠3 : Then “ →” : IfEnd
0→I : Z ［98］+ 8 →M
If S ≥ Z ［M －8 ］And S ≤ Z ［M ］:Then Goto 1 : IfEnd
Lbl 2 : I +1 →I : M + 8 →M
If S ≤ Z ［M ］:Then Goto 1 : Else Goto 2 : IfEnd
Lbl 1 : If O = 3 : Then I →B :IfEnd : M→Z［27］: If O ≠ 4 : Then Prog “Q. 3” : Else “ L(I) , I = ” : I ◣
Cls : “ MIN = ” : Locate 6 , 1 , Z ［Z ［98］+8 I ］: “ MAX = ” : Locate 6 , 2 , Z ［Z ［98］+8 I +8］:I→B : IfEnd
子程序Q . 3
If O = 3 And U = 0 : Then B→I : IfEnd
Z ［98］+8 I →G
Z ［G ］→Z
Z ［G +1］→Z［1］: Z ［G +2］→Z［2］: Z ［G +3］→Z［3］
Z ［3］:Prog “ J ” : T→Z［3］
Z ［G +4］→Z［4］: Z ［G +5］→M : 1÷ Z ［G +6］→Z［5］: 1÷ Z ［G +7］→Z［6］
Z ［6］－Z ［5］→Z［6］: S －Z→D
If O = 3 : Then 0→N : IfEnd : Prog “Q . 4 ”
子程序Q .4 (正算)
Lbl 0 : 0→I : 0→E:0→F : 0→G : 0→H
Lbl 1 : I + 1→I : Z ［3］+ Z ［4］( ( I D÷12) ( Z ［5］+ I D Z ［6］÷24÷M ) ( 180÷л ) →W : If I ﹤12: Then If Frac ( I÷2 ) ﹥0 : Then E + cos ( W ) →E : F + sin ( W ) →F : Else G + cos ( W ) → G : H + sin ( W ) → H : If End : IfEnd
If I ﹤12 : Then Goto 1：IfEnd
Z ［1］+（D÷36）( cos ( Z ［3］) + cos ( W ) + 2G + 4E ) + N cos ( U + W
) →P
Z ［2］+（D÷36）( sin ( Z ［3］) + sin ( W ) + 2H + 4F ) + N sin ( U + W ) →Q
If O = 3 : Then Goto 2 : IfEnd
If O ≠ 0 : Then Cls : “ XP = ” : Locate 5 , 1 , P : “ YP = ” : Locate 5 , 2 , Q : “ AT = ” : W►DMS◣
IfEnd
Z ［27］→I : If Z ［I ］－ S = 0 And Z［I + 1］= 0 And N = 0 : Then P →Z［I + 1］: Q→Z［I + 2］: W : Prog “ H ” : T→Z［I +3］: IfEnd
If O ≠ 1: Then Pol( P －Z ［11］, Q －Z ［12］) : If J
Prog “ Q. 6 ” : IfEnd
Lbl 2
子程序F . 2 （反算）
“ ZH.D(LI) →0 : NO→1” ? U
Lbl 2 : If U = 0 : Then “ I ” ? B : IfEnd
Z ［98］+ 8 Z ［99］→M
If X = 0 : Then Z［13］→X : IfEnd : “ X ” ? X : If X = 0 : Then Goto 3 : IfEnd
“ Y ” ? Y : 1 →Z［8］
“=〉 ”
If U = 1 : Then －√ ( ( X －Z ［M + 1］) ² + (Y－Z ［M +2］) ² ) →Z［7］: Z ［7］+Z［M ］→S : Prog“Q . 2”: IfEnd
Z ［98］+ 8 B→M : Z ［M +3］: Prog“J ”: T －90→T
( Y －Z ［M +2］) cos (T) －( X －Z ［M + 1］) Sin (T) →Z［7］
If U =1 And Abs ( Z ［7］) ≤ 0.01: Then －0.01 →Z［7］: IfEnd
Z ［7］→D : Z ［M ］+ D →S :
If S > Z ［M +8］: Then Z［M +8］→S ; IfEnd
If U = 0 : Then Prog“Q . 3”: Else Prog“Q . 2”: IfEnd
“=〉 ”
Lbl 1 : ( Y －Q ) cos （W －90）－(X－P) sin (W－90) → Z ［8］
Z + Z ［7］+ Z ［8］→S
If Abs ( Z ［8］) ≥ 0.0001: Then Z ［7］+ Z ［8］→D : Prog“Q .4”: D→Z［7］: Goto 1: IfEnd
If U =1: Then If S + 0.01
If B
Cls : “ MIN= ” : Locate 6 , 1 , I : “ MAX= ” : Locate 6 , 2 , J : “ END, =〉ZH.D(LI)”◣
X →Z［13］: 0→X : 3→O : Goto 3 : IfEnd: IfEnd
If U =0: Then If S
Goto 2 : IfEnd : If S> Z ［Z ［98］+ 8(B+1) ］+ 0.01: Then “ DK>MAX , L→(I+1) ” ◣
Goto 2 : IfEnd : IfEnd
Cls : “ DKP=” : Locate 6 ,1, S : X －P→P : If P =0: Then ×10－9→P : IfEnd : Pol( P , Y －Q) : If sin（W －J ）>0: Then －I→I : IfEnd : “ D÷2= ”: Locate 6 , 2, I : Goto 2
Lbl 3
子程序 J
Ans ÷ . 36 －16 ( Int ( Ans ) ) ÷ 9－Int (100 (Ans－ Int (Ans ) ) ) ÷90→T
子程序 H
. 6 4 Int ( Ans )+. 36 Ans + . 004 Int ( ( 60 ( Ans － ( Int (Ans ) ) →T
为了实现桥台转轴坐标计算，现（2009.04.19）加一个子程序如下：
子程序 JH . ZJ (桥台转轴)
“=〉B ”? V : V : Prog“J : W + T→T
Lbl 1 : “=〉(－Y +) D ”? N : N = 0 =〉 Goto 2 : “=〉(－→ +) S ”? S Pol(N , S ) : Rec( I ,T －J )
Cls : “X = ”: Locate 6, 1 , I + P : “Y = ”: Locate 6 , 2 , J + Q ◣
Goto 1
Lbl 2 : “ OK ” ◣
在子程序Q.1 中语句 Lbl 1 : Z ［98］→ I : Z ［99］→ J : “ →DKP ” ? S 后加
If S = －1: Then Prog “JH . ZJ ” : Goto 1 : IfEnd （加桥台转轴时）就可以了，(2009.04.19以加，见绿色标注语句。
桥台（桩基）转轴直角坐标计算 运行提示符
直角坐标放样令DKP=-1时，转向桥台（桩基）转轴直角坐标计算
提 示 符 说 明
以上一测点中桩为中心作偏转
=〉 B ? 问转轴偏角（以小数点为度）
=〉(－Y + ) D? 问偏距（与直角坐标放样时意义一致）。赋值0时转向OK =〉(－→ + ) S? 问支距（小里程方向为－，大里程方向为+）
X= 给出点位纵坐标
Y= 给出点位横坐标
OK 本转轴各点计算结束，返回直角坐标放样。
子程序 ZH . 2 (涵通转轴)
“=〉B ”? V : V : Prog“J”: W + T→T
Lbl 1 : “=〉(－→ +) K ”? S : S = 0 =〉 G oto 2 : “=〉(－Y +) D ”? N N Cos ( T) + P → X : N Sin ( T) + Q→ Y
Cls : “X = ”: Locate 6, 1 ,S Cos( W ) + X : “Y = ”: Locate 6 , 2 , S Sin(W ) + Y ◣
Goto 1
Lbl 2 : “ OK ” ◣
在子程序Q.1 中语句If S = －1: Then Prog “JH . ZJ ” : Goto 1 : IfEnd 后加 If S = －2: Then Prog “ZH . 2 ” : Goto 1 : IfEnd （加涵通转轴时) 即可。
涵通转轴直角坐标计算 中桩直角坐标放样令DKP=－2时，
提 示 符 示 例 说 明
以上一测点中桩为中心作偏转
=〉 B ? 问涵通偏角（以小数点为度）
=〉(－→ + ) δ K? 问里程差（相对于涵通中桩，小里程为－、大里程为+），赋值0时转向OK
=〉(－Y + ) D? 问偏距（与直角坐标放样时意义一致）
X= 给出点位纵坐标
Y= 给出点位横坐标
OK 本转轴各点计算结束，返回直角坐标放样。
算例1
例某工程C 匝道:
根据算例数据确定综合曲线属性: 令名 C.ZD 令编号0 令（矩阵存起始储单元 ）Z ［MIN ］= 100 （综合曲线段数, 不含终点）N = 4
红色数据应在Q.5程序运行时按提示赋值
一 给矩阵存储单元定维:（COMP 状态下） 140→DimZ
二 编写程序PPQX （PROG 状态下给定综合曲线属性）
“ C. ZD →0 ” ？ C
If C=0: Then 起始单元100→Z［98］: 段数4 →Z［99］: IfEnd
Prog “ PQX ”
三 启动Q . 5 按下列数据建立矩阵存储表 （COMP 状态下按提示赋值，详见附件运行提示符）
自 然 段
名 称 L0 L1 L2 L3 L4
DK 0 190 355.927 472.168 561.791
X 20934.495 21066.119 21142.859 21076.993 20988.267
Y 89274.172 89411.182 89552.427 89640.286 89644.485
AT 47.08177 45.09322 93.33281 160.09323 185.50027
±1 -1 1 1 1 0
L 190 165.927 116.241 89.623 0
R1 5500 5500 100 100 0
R2 5500 100 100 ×10 20 0
备注： 起始矩阵存储单元 Z ［100］ 终了矩阵存储单元 Z ［139 ］。 (139 =100 +8 × 5-1)
四 启动程序PPQX （COMP 状态下指定需计算的综合曲线编号）
C. ZD →0 ？ 提示C 匝道赋值 0
回车
（自动转向PQX 进入各项计算）
为了提高计算速度，现提供Q.4 (5点法加速), 以该程序取代原有Q.4, 其它无需改变。取代后按指定区间反算时间大约4秒，有兴趣的朋友可以试试。
Q.4 (5点法加速)
If Z［4］= 0 : Then 1→Z［4］: IfEnd
Z ［4］Z ［G+6］→Z［9］: Z ［4］Z ［G+7］→Z［10］
( Z ［9］－Z ［10］) ÷ （ 2M Z［9］Z ［10］）→Z［13］: S －Z→ Z ［14］: Z ［13］Z ［14］→Z［13］: 1 ÷ Z ［9］→Z［22］
0.1739274226→ Z ［15］: 0.3260725774→ Z ［16］: 0.0694318442→ Z ［17］: 0. 3300094782→ Z ［18］
1－Z ［18］→Z［19］: 1－Z ［17］→Z［20］: 180 ÷ л → Z ［21］
Z ［3］+ Z ［17］Z ［14］(Z［22］+ Z ［17］Z ［13］) Z ［21］→Z［17］ Z ［3］+ Z ［18］Z ［14］(Z［22］+ Z ［18］Z ［13］) Z ［21］→Z［18］ Z ［3］+ Z ［19］Z ［14］(Z［22］+ Z ［19］Z ［13］) Z ［21］→Z［19］ Z ［3］+ Z ［20］Z ［14］(Z［22］+ Z ［20］Z ［13］) Z ［21］→Z［20］ Z ［3］+ Z ［14］ ( Z ［22］+ Z ［13］) Z ［21］→W
Z ［1］+ Z ［14］ ( Z ［15］cos( Z ［17］) + Z ［16］cos( Z ［18］) + Z ［16］c os( Z ［19］) + Z ［15］cos( Z ［20］) ) +N cos( U + W ) →P
Z ［2］+ Z ［14］ ( Z ［15］sin( Z ［17］) + Z ［16］sin ( Z ［18］) + Z ［16］s in ( Z ［19］) + Z ［15］sin ( Z ［20］) ) + N sin (U + W) →Q
If O = 3 : Then Goto 2 : IfEnd
If O ≠ 0 : Then Cls : “ XP= ” : Locate 5 , 1 , P : “ YP= ” : Locate 5 , 2 , Q : “ AT= ” : W►DMS◣
IfEnd
Z ［27］→I : If Z［I ］－ S = 0 And Z［I+1］= 0 And N = 0 : Then P→Z［I +1］: Q → Z ［I+2］: W : Prog “ H ” : T→ Z ［I+3］: IfEnd
If O ≠ 1: Then Pol( P －Z ［11］, Q －Z ［12］) : If J
Prog “ Q. 6 ” : IfEnd
Lbl 2