1程序(TYQXjs)
1SZ > XY:2XY > SZ:N:UX0:VY0:OS0:GF0:HLS:PR0:RRN:Q:C1÷P:D(PR)÷(2HPR):E180÷π:N1>Goto 1:≠>Go to 2Δ←┘
Lbl 1:{SZ}:SZ:WAbs(SO):Prog SUB1:XXSX◢
YYSY◢
Goto 1←┘
Lbl 2:{XY}:XY:IX:JY:Prog SUB2:SSO+W◢
ZZZ◢
Goto 2
2 正算子程序(SUB1)
A01739274226:B03260725774:K00694318442:L03300094782:F1L:M1K:XU+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD))):YV+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))):
FG+QEW(C+WD)+90:
XX+ZcosF:YY+ZsinF
3 反算子程序(SUB2)
TG90:WAbs((YV)cosT(XU)sinT):Z0:Lbl 0:Prog SUB1:LT+QEW(C+WD):Z(JY)cosL(IX)sinL:AbsZ<1E6>Goto1:≠>WW+Z:Goto 0Δ←┘
Lbl 1:Z0:Prog SUB1:Z(JY)÷sinF
二说明
1规定
(1) 道路中线前进方(里程增方)区分左右
线元左偏时
Q1线元右偏时Q1线元直线时Q0
(2) 求点位中线时Z0位中线左铡时Z取负值
位中线中线右侧时Z取正值
(3) 线元直线时起点止点曲率半径穷
1045次代
(4) 线元圆曲线时起点止点什线元相接曲率
半径均等圆弧半径
(5) 线元完整缓曲线时起点直线相接时曲率半
径穷1045次代圆曲线相接时曲率半
径等圆曲线半径止点直线相接时曲率半径穷
1045次代圆曲线相接时曲率半径等圆曲线半径
(6) 线元非完整缓曲线时起点直线相接时曲率半径等设计规定值
圆曲线相 接时曲率半径等
圆曲线半径止点直线相接时曲率半径等设计规定值圆曲
线相接时曲率半径等圆曲线半径
2输入显示说明
输入部分:
1 SZ > XY
2 XY > SZ
N 选择计算方式输入1表示进行里程边距计算坐标 输入2
表示坐标反算里程边距
X0 ?线元起点X坐标
Y0 ?线元起点Y坐标
S0 ?线元起点里程
F0 ?线元起点切线方位角
LS ?线元长度
R0 ?线元起点曲率半径
RN ?线元止点曲率半径
Q ? 线 元左右偏标志(左偏Q1右偏Q1直线段Q0)
S ? 正算时求点里程
Z ?正算时求点距中线边距(左侧取负值右侧取正值
中线取零)
X ?反算时求点X坐标
Y ?反算时求点Y坐标
显示部分:
XS××× 正算时计算出求点X坐标
YS××× 正算时计算出求点Y坐标
S××× 反算时计算出求点里程
Z××× 反算时计算出求点边距
四算例
某匝道五段线元(直线+完整缓曲线+圆曲线+非完整缓曲线+直线)
组成段线元素(起点里程S0起点坐标X0 Y0起点切线方位角
F0线元长度LS起点曲率半径R0止点曲率半径RN线 元左右偏标志Q)
:
S0 X0 Y0 F0 LS R0 RN Q
500000 19942837 28343561 125 16 3100 269256 1E45 1E45 0
769256 19787340 28563378 125 16 3100 37492 1E45 22175 1
806748 19766566 28594574 120 25 5407 112779 22175 22175 1
919527 19736072 28701893 91 17 3063 80285 22175 9579228 1
999812 19744038 28781659 80 40 5000 100000 1E45 1E45 0
1正算
(注意:略计算方式线元素输入请行根求点线元输入线元素)
S700 Z5 计算 XS1983141785 YS2850972590
S700 Z0 计算 XS1982733592 YS2850683837
S700 Z 5 计算 XS1982325398 YS2850395084
S780 Z5 计算 XS1978525749 YS2857502270
S780 Z0 计算 XS1978115561 YS2857216358
S780 Z 5 计算 XS1977705373 YS2856930446
S870 Z5 计算 XS1974753609 YS2865413091
S870 Z0 计算 XS1974268648 YS2865291379
S870 Z 5 计算 XS1973783688 YS2865169668
S940 Z5123 计算 XS19741 59118 YS2872205802
S940 Z0 计算 XS1973647687 YS2872235642
S940 Z 3009 计算 XS1973347298 YS2872253168
2 反算
X19831418 Y28509726 计算 S6999999974 Z 5 00018164
X19827336 Y28506838 计算 S6999996493 Z 0000145136
X1982325398 Y2850395084 计算 S6999999985 Z 5000003137
X1978525749 Y2857502270计算 S7800000035 Z 5 000001663
X1978115561 Y2857216358 计算 S7800000025 Z 0000002979
X1977705373 Y2856930446 计算 S7800000016 Z 499999578
X19747536 Y28654131计算 S8700001137 Z 499941049
X19742686 Y28652914 计算 S8700003175 Z 000041814
X9737837 Y28651697 计算 S8700002748 Z 4999808656
X197415912 Y287220580 计算 S9399999786 Z 5123024937
X197364769 Y287223564 计算 S9399999862 Z 0000027710
X197334730 Y287225317 计算 S940000 0238 Z 300898694
3 修改意见:
程序正算反算似正确满足测量精度求明白三点:
<1>意曲线正反算达精度求?
<2>字母代表意思?推导程?会会出错?
<3>应该补充什功?提高精度求?
修改程序没条件限制运行正常完全满足客户精度求实愿家分享
4 例题
圆曲线半径R700米圆曲线长350米曲线右转起点里程DK1+200坐标(200456) 切线方位角120°15′266〃计算DK1+520右100米点坐标反算桩号左右距
读者算源程序计算结果:
正算 X7303846 Y6033210
反算 W 1512069875 Z954193 (明显超出范围)
W 1512069875 Z954193计算坐标坐标肯定错
文香网httpwwwxiangdangnet
《香当网》用户分享的内容,不代表《香当网》观点或立场,请自行判断内容的真实性和可靠性!
该内容是文档的文本内容,更好的格式请下载文档