曲线任意里程中边桩坐标正反算（4850P计算器）程序（第四次修改

曲线任意里程中边桩坐标正反算(4850P计算器)程序(第四次修改版）

曲线任意里程中边桩坐标正反算(CASIO fx-4850P计算器)程序 （ 第四次修改版 ）

一、程序功能及原理

1.功能说明：本程序由一个主程序(TYQXjs)和四个子程——正算子程序(SUB1)、反算子程序(SUB2)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。本程序也可以在CASIO fx-4800P计算器运行。

2．计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。 二、源程序 1.主程序(TYQXjs)

Lbl 0：Defm 50:\=> XY,2.XY => SZ\：N： N=1=>Goto 1：≠>Goto 2Δ←

Lbl 1：{SZ}：SZ

S≤Z[18]=>J=1：Prog “DAT1”：Goto 3Δ S≤Z[26]=>J=2：Prog “DAT1”：Goto 3Δ S≤Z[34]=>J=3：Prog “DAT1”：Goto 3Δ S≤Z[42]=>J=4：Prog “DAT1”：Goto 3Δ ……………

S≤Z[8（N+1）+2]=>J=n：Prog “DAT1”：Goto 3Δ Lbl 3：W=Abs(S-O)：Prog \：\◢ \◢

”FWI=”: F=F-M :F→DMS◢

Goto 4←┘

Lbl 2：{XY}：XY：Z[4]=X：Z[5] =Y：N=0

Lbl A：Isz N：A=Z[8N+3]-M：B=Z[8（N+1）+3]-M：Prog \：Z[6]×Z[7]<0=>J=N：Prog DAT1\：Goto BΔA=Z[8N+3]+M：B=Z[8（N+1）+3]+M：Prog \：Z[6]×Z[7]<0=>J=N：Prog DAT1\：Goto B：≠> Goto AΔ Lbl B：Prog \：\O+W◢ \◢ Goto 2

Lbl 4:J=0:I=Pol(X-Z[1],Y-Z[2]):F=J:F<0=>F=F+360Δ ”DIST=”: I◢”FW=”:F→DMS◢ Goto 1

2. 正算子程序(SUB1)

A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:Z[3]=0.5:X=U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Ncos(G+QEZ[3]W(C+Z[3]WD))+Bcos(G+QE(1-L)W(C+(1-L)WD))+Acos(G+QE(1-K)W(C+(1-K)WD))): Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Nsin(G+QEZ[3]W(C+Z[3]WD))+Bsin (G+QE(1-L)W(C+(1-L)WD))+Asin (G+QE(1-K)W(C+(1-K)WD)))：M”ANG=”: F=G+QEW(C+WD)+M：X=X+ZcosF：Y=Y+ZsinF

3. 反算子程序(SUB2)

M”ANG=”:T=G-M：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog \：L=T+QEW(C+

WD)：Z=(Z[5]-Y)cosL-(Z[4]-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘

Lbl 1：Z=0：Prog \：Z=(Z[5]-Y)÷sinF 4. 垂距计算子程序(ZX1)

Z[6]=(Z[5]-Z[8N+1])COS A-(Z[4]-Z[8N])SIN A Z[7]=(Z[5]-Z[8(N+1)+1])COS B-(Z[4]-Z[8(N+1)])SIN B 5.曲线元要素数据库：DAT1

Lbl1:J=1=>U=Z[8]:V=Z[9]:O=Z[10]:G=Z[11]:H=Z[12]:P=Z[13]:R=Z[14]:Q=Z[15]Δ

J=2=>U=Z[16]:V=Z[17]:O=Z[18]:G=Z[19]:H=Z[20]:P=Z[21]:R=Z[22]:Q=Z[23]Δ J=3=> U=Z[24]:V=Z[25]:O=Z[26]:G=Z[27]:H=Z[28]:P=Z[29]:R=Z[30]:Q=Z[31]Δ J=4=> U=Z[32]:V=Z[33]:O=Z[34]:G=Z[35]:H=Z[36]:P=Z[37]:R=Z[38]:Q=Z[39]Δ ……………………………

J=N=> U=Z[8N]:V=Z[8N+1]:O=Z[8N+2]:G=Z[8N+3]:H=Z[8N+4]:P=Z[8N+5]:R=Z[8N+6]:Q=Z[8N+7]Δ

(注：如有多个曲线元要素继续添加入数据库DAT1中)

Lbl 2 ：U\：V\：O\：G\：H\：P\：R\：Q：C=1÷P：D=(P-R)÷(2HPR)：E=180÷π 三、使用说明 1、规定

(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，

Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右

侧时，Z取正值。

(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。

(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆 弧的半径。

(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45

次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半

径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的

值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等

于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

(7)曲线元要素数据库（DAT1）可根据线型不同分为各个线元段输入到DAT1中，即分为直线段、缓和曲线、圆曲线等。

(8)正算时可仅输入里程Lp和边距Dp及右交角ANG全线计算，反算时通过输入测点的X、Y坐标和右交角ANG后计算器自动判断该点所属曲线元并利用该线元的曲线要素执行反算中桩里程Lp及支距Dp。 2、输入与显示说明 （一）、 输入部分： 1. SZ => XY 2. XY = > SZ

1、 N ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算 里程和边距。

2、X0 ？线元起点的X坐标 3、Y0 ？线元起点的Y坐标 4、S0 ？线元起点里程 5、F0 ？线元起点切线方位角 6、LS ？线元长度 7、R0 ？线元起点曲率半径 8、RN ？线元止点曲率半径

9、Q ？ 线 元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0) 10、S ？ 正算时所求点的里程

11、Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)

12、ANG？正算边桩时左右边桩连线与线路中线的右交角 13、J？曲线元数据库曲线段判断系数（J=1、2…..n） 14、X ？反算时所求点的X坐标 15、Y ？反算时所求点的Y坐标 16、M?斜交右角

17、Z[1]-----测站点X坐标 18、Z[2]----- 测站点Y坐标

19、A、B、N是Gauss-Legendre求积公式中的插值系数 20、K 、L、Z[3] 是Gauss-Legendre求积公式中的求积节点

21、N=1时：Z[8]、Z[9]、Z[10]、Z[11]、Z[12]、Z[13]、Z[14]、Z[15]分别是各曲线元的X、Y坐标、起始里程S0、初始方位角F0、线元长度Ls、起点半径R0、终点半径Rn、曲线偏向系数Q

……………………………………………………

N=n时：Z[8N]、Z[8N+1]、Z[8N+2]、Z[8N+3]、Z[8N+4]、Z[8N+5]、Z[8N+6]、Z[8N+7]分别是各22、曲线元的X、Y坐标、起始里程S0、初始方位角F0、线元长度Ls、起点半径R0、终点半径Rn、曲线偏向系数Q

23、正算中的Z[18]、Z[26]、Z[34]、Z[42]……Z[8（N+1）+2]分别为各线元中的终点里程

24、主程序中的“Defm 50”应视具体线元的数量，适当增加扩展变量。 25、本次修改版不考虑CASIO fx-4500PA，因其容量太小，扩展变量实在太少而不能胜任。 （二）、显示部分：

XS=××× 正算时，计算得出的所求点的X坐标 YS=××× 正算时，计算得出的所求点的Y坐标 S=××× 反算时，计算得出的所求点的里程 Z=××× 反算时，计算得出的所求点的边距 四、算例

某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组

成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度

LS、起点曲率半径R0、止点曲率半径RN、线 元左右偏标志Q）如下： S0 X0 Y0 F0 LS R0 RN Q 500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0

769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.7

5 -1

806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.75 -1

919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228 -1

999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0 1、正算

（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元 要素）

S=700 Z=-5 计算得 XS=19831.41785 YS=28509.72590 S=700 Z=0 计算得 XS=19827.33592 YS=28506.83837 S=700 Z= 5 计算得 XS=19823.25398 YS=28503.95084

S=780 Z=-5 计算得 XS=19785.25749 YS=28575.02270 S=780 Z=0 计算得 XS=19781.15561 YS=28572.16358 S=780 Z= 5 计算得 XS=19777.05373 YS=28569.30446

S=870 Z=-5 计算得 XS=19747.53609 YS=28654.13091 S=870 Z=0 计算得 XS=19742.68648 YS=28652.91379 S=870 Z= 5 计算得 XS=19737.83688 YS=28651.69668

S=940 Z=-5.123 计算得 XS=19741. 59118 YS=28722.05802 S=940 Z=0 计算得 XS=19736.47687 YS=28722.35642 S=940 Z= 3.009 计算得 XS=19733.47298 YS=28722.53168 2、 反算

X=19831.418 Y=28509.726 计算得 S=699.9999974 Z= -5 .00018164 X=19827.336 Y=28506.838 计算得 S=699.9996493 Z= 0.000145136 X=19823.25398 Y=28503.95084 计算得 S=699.9999985 Z= 5.000003

137

X=19785.25749 Y=28575.02270 计算得 S=780.0000035 Z= -5 .000001663

X=19781.15561 Y=28572.16358 计算得 S=780.0000025 Z=- 0.000002979

X=19777.05373 Y=28569.30446 计算得 S=780.0000016 Z= 4.99999578

X=19747.536 Y=28654.131 计算得 S=870.0001137 Z= -4.99941049 X=19742.686 Y=28652.914 计算得 S=870.0003175 Z=- 0.00041814 X=19737.837 Y=28651.697 计算得 S=870.0002748 Z= 4.999808656

X=19741.5912 Y=28722.0580 计算得 S=939.9999786 Z= -5.123024937 X=19736.4769 Y=28722.3564 计算得 S=939.9999862 Z=- 0.000027710 X=19733.4730 Y=28722.5317 计算得 S=940.0000238 Z= 3.00898694

注：1、承蒙yshf的源程序：《曲线任意里程中边桩坐标正反算程序》的共享，这里特意感谢！

2、希望广大网友同行支持程序共享并积极扩展，多多网上交流！

查看楼主其它文章 <<返回

→ 回复内容 只查看楼主与本人回复 1 楼：测量员 回 hanniu_1973 :老师您才是真正的大师.中国有您这样如此康概的测量工作员,中国的测量事业有前途了． 2 楼：jinyong0813 各位大虾要是能加上高程的程序,共用此程序输入数据就太棒了!一次计算出三维坐标. 3 楼：hanniu_1973楼主 回 评：+4 分 回 评：+5 分 关于公路三维动态控制测量系统，我正在不断完善中，现在并不是很成熟，故未发布，如有各位同行有更好的程序或编程思路可发布网上进行深度交流！ 4 楼：hanniu_1973楼主 回 我已完成了曲线平行线坐标及其与直线等交点坐标计算的严密解算程序，如有这方面想法、意见可与我交流。 hanniu_1973@163.com 5 楼：280050876 hanniu_1973 : 你的主程序中\F=F-M F→DMS◢ \无法运行,我改为\-M →DMS◢ \后可行。 另外有个问题想请教你：就是关于扩展变量“Defm”的取值是否取决于“J=n”呢？拿上面算例来说，“Defm”的取值该取多少为宜？？我发现在上面算例中正算时“Defm”的取值问题不是很大，但反算时“Defm”的取值不正确程序运行时会出现“Mem ERROR”，导致反算无法运行。请指教！！！！！ 6 楼：hanniu_1973楼主 280050876: 关于上面第一个问题，不好意思！因为网络上传或文字处理软件的问题，故我已在上面的文章里修改了，请在上面查看。而第二个问题:Defm 8N+8 7 楼：280050876 hanniu_1973 : 在主程序的Lbl A中B=Z[8（N+1）+3]-M及B=Z[8（N+1）+3]+M和垂距计算子程序(ZX1)中Z[8(N+1)+1]等扩展变量在算例中最后一线元中找不到这些变量,故程序运行时总显示“Mem ERROR”,这是我的个人观点,不知对否?请指教!!!! 8 楼：hanniu_1973楼主 我在上面已说过了，即Defm 8N+8 9 楼：280050876 hanniu_1973 : 程序中Defm 8N+8取值不错,但程序运行时如上算例中Z[8（N+1）+3]是不是就是Z[51]了,但变量中无Z[51]呀?????? 10 楼：hanniu_1973楼主 准确地说是：Defm 8（N+1）+8 11 楼：111 hanniu_1973 ： 这个程序我看了一下发现你这次改编的思路非常好，但我有个问题想问问你；分离式路基这个里程怎么处理啊，会重复了啊！ 12 楼：ws 回 回 评：+5 分 回 回 回 回 回 回 程序已经很完美，不过，你们是高手，我想供广大网友提供思路，能改进就好了。 修改的程序对于已知数据区还不如按照第三种修改方案，反正要输入，免得扩展。但是单个数据源的求解不如按照原程序，所以不知道高手能否按照以下思路再改改。 第一，建立主数据程序，即选择不同的匝道数据库进入。如加一个条件关系就可以了。把不同的匝道数据作子程序。 第二，选择S,即任意里程，进入不同的线元。 第三，进入运算主程序。已知附值还是编进去，虽然已知数据会根据数据库程序重复出现一次，但是也没有多大关系，为了熊掌和鱼兼得，可以从该主程序进入，就对已知数据库以外的线元能单独计算。 一点小的建议，冒昧！ 13 楼：HANNIU_1973楼主 111: 对于分离式路基，可对数据库和条件判断进行修改即可，因其比较简单和使用不多，请各位自行扩展！ 14 楼：HANNIU_1973楼主 ws： 首先说明一下，你提的思路不错，希望各位网友或同行积极加入共享行列！其实编程最重要的是数学模型和解算思路，其余的只是程序员的份内事了。当然要问到具体的实现步骤，这有点像学书法一样，对一位尚未入门或刚入门的学者来说，的确不易，能把学写工整就已不错了，更何况章法和韵律呢？言归正传，把各条线的线元数据作为数据库编入子程序，然后在主程序的正算条件判断中加入路由并在反算中设置N的初值即可，因其比较简单，请各位自行扩展，我在这里就不再赘述了。另外，在这里我向大家提二点思路：1、可加入平行线解算功能及其与曲线求交功能；2、可加入方向坐标法或称虚拟测站法放样功能。 15 楼：qw-197312 请教几个问题： １＼Defm 后面的变量如何设置，好象用８（ｎ＋１）＋８坐标反算时老时出问题，设置其它也不可以． ２＼坐标正算时曲线在ＨＹ＼ＨＺ两点坐标计算均有问题． ３＼程序＂ＤＡＴ１＂中的＂ＬＢ１＂和 ＂ＬＢ２＂好象无ＧＯＴＯ语句对应． 我用的是ＦＸ４８００计算器， 特此请教！ 16 楼：280050876 hanniu_1973 : 我不知道这个程序你是否已调试通过,反正我是对照你的程序一字一句的输入的.但我的Defm取值不管是Defm 8N+8或是Defm 8（N+1）+8,反算时均未能通过,总是出现“Mem ERROR”.我个人认为是不是你主程序的Lbl A中B=Z[8（N+1）+3]-M及B=Z[8（N+1）+3]+M和垂距计算子程序(ZX1)中Z[8(N+1)+1]等扩展变量在子程序DAT1中的变量不存在的回 回 回 回 缘故?????拿上面算例来说,当J=5时,Z[8（N+1）+3]=Z[51],Z[8(N+1)+1]=Z[49],而Z[51]、Z[49]这两个变量在主程序和子程序中根本找不到相对应的数值。 17 楼：ws hanniu_1973:你好！ 我对你复合SIMPSON公式计算的曲线元坐标程序感兴趣，这两天怎么进不去。再说，我调用你的这个程序，和“谢了”网友说的一样。能不能你再检查检查！谢了 18 楼：NetBee FW=”:F→DMS◢ Goto 1 中的F→DMS怎么实现？找不到呀？ 19 楼：111 hanniu_1973: 你好我是个新手，你能告诉我再程序中如：Z［10］如何给它赋值啊！！！！！！！！！！ 20 楼：hanniu_1973楼主 直接输入Z[10]=1793.258 21 楼：hanniu_1973楼主 NetBee: 你的计算器是4850P吗？你可直接输入即可，但如是4800P那需编程解决！ 22 楼：yshf hanniu_1973老兄： 你的4800P是什么样的厂生产，只要定义了扩展变量（Defm N , N>=10且在内存容许的情况下），就可以直接在4800P上输入“Z[10]=1793.258”，你恐怕是没试用过此吧！ 23 楼：hanniu_1973楼主 回 回 回 回 回 回 回 我用的是4850P，不过这一层面，4800P亦应该如此，即为扩展变量赋值时就用Z[10]=1793.285。如果不是这样，我真有点纳闷了，难道还有别的输入途径，请大家贴出来推广，我愿意赐教！ 24 楼：111 hanniu_1973： 是不是在子程序“DAT1”中直接这样输入啊，O=Z[10]=1793.258: 25 楼：hanniu_1973楼主 111: 那倒不是，是在运行程序前，直接在计算器里输入Z[10]=1793.258回车即可。 回 回

26 楼：wyqzm 呵呵,此次修改比之以前的要强了一些了,提点意思行吗?(高手看了不要见笑了) 回 1,建议取消所有扩展变量,因为有部分用户用的还是4800P,或4500,并且可能已经存储了很多自己的小程序,如果使用了Defm 50用来扩展50个定长字节串,那就很有可能产生存储器超界的.注意说的是\字节串\不是字符串,由于这种变量无论你是否使用,它都要占用固定长度的字节空间,并且每使用一个变量,他就会比默认变量多好几个字节,一旦程序变长,特别是数据文件的变长,这种多余字节的影响对4800的用户还是很大的,对于4500的用户那就不要想了,它的存储容量太小了, 2,子程序SUB1使用的常量应是可以代换的,个人认为设置这样的常量并不可取,虽然说在一定程序上节省了点空间,但相对你的整个SUB1来说,这点空间实在是微不足道的,并且这个常数大多数人并不好理解,我认为还是改成计算表达式更为好一点. 3,你的程序极尽计算器编程之能事啊,用一段子程序SUB1就实现了,说实话,我看起来真是有点费劲,不知道其他同仁有没有同感啊?对于此式我也有一句要说,希望你不要介意才是: 诚然,在程序中使用一定的技巧可以达到很好的效果,但现代程序设计的感点认为,条理清晰,简单易懂的程序才是标准哦,你的计算表达写得过于繁杂,一般人是难以看懂的,而且算法体现不出来啊,要不你另外提供一个这个程序的计算表达式?毕竟坐标并不闹着玩的,它很重要. 4,既然是线元曲率半径法,那当然就是直线和圆曲线元两端曲率半径相等,回旋线元不相等了,你在程序中规定直线元曲率半径输入10的45次方,我觉得这一点设计不可取啊,虽然10的45次方已经够大了,但有个更简单的方法怎么不用呢?无穷大的倒数就是无穷小,此外需入0就行了,没有必要输入那么个大数吧? 5,回旋线的精度可以通过循环迭代获得提高,但这个设置值在哪里呢?我看了半天也没有看明白,请指教 [img]http://www.cz88.net/ip/pic_5.aspx[/IMG] 27 楼：njy135njy 回 在主程序的Lbl A中B=Z[8（N+1）+3]-M及B=Z[8（N+1）+3]+M程序运行时总显示“Mem ERROR”, 28 楼：njy135njy 回 我输了16个线元的数据，反算时在主程序的Lbl A中B=Z[8（N+1）+3]-M及B=Z[8（N+1）+3]+M程序运行时总显示“Mem ERROR”,我的defm设定为200，请指点一下啊！谢了 29 楼：yshf njy135njy： 定义一个扩展变量要耗有10字节，(200+1)×8×10=16080，如果是4800P，已远远大于4500字节了；如果是4850P，还回 勉强可以（机内所有程序耗用的字节数不得超过12420情况下）。个人认为要计算16个线元，使用PC—1500、PC—E500，或者掌上电脑，可以不当心出“Mem ERROR”的问题。 30 楼：yshf wyqz: (1). 第一个观点略有同感； (2). SUB1中的常量乃积分得出的真数，怎样用计算表达式代替？那样做岂不是再将简单的问题复杂化，无异于站起来翻身； (3). 载重量是几百千克的面包车，能期望像一列火车或是万吨巨轮那样的载重量吗？不问具体的运用环境，乱用标准，对微分几何知之甚至少，给在别人的感觉就成了：“井底之蛙、赶面杖吹火”；要是人人都看懂了，hanniu_1973岂不吃饱了没事干，费尽心想写这样的程序逗大伙玩吗？ (4).10的45次方小了，你再试试10的100次方、1000次方、......，那肯定会味道不错的； 照你的想法，在4800P或4850P上运行一下：“1÷0”，说不定你会有新的收获的！ (5).回旋线什么样的精度，怎样通过循环迭代来获得提高呢？还是那句话：“你对微分几何知道的不多”，当然看不明白那设置值啦，明确地告诉你：正算中就没有循环迭代。建议你再找有关资料看看，信口雌雄，盲加评论，这可不是好事，这对一个人是很重要的。 31 楼：njy135njy 嗯，现在我的程序可以用行了，我是这样处理的： 把DAT1程序中的 Lbl1:J=1=>U=Z[8]:V=Z[9]:O=Z[10]:G=Z[11]:H=Z[12]:P=Z[13]:R=Z[14]:Q=Z[15]Δ J=2=>U=Z[16]:V=Z[17]:O=Z[18]:G=Z[19]:H=Z[20]:P=Z[21]:R=Z[22]:Q=Z[23]Δ J=3=> U=Z[24]:V=Z[25]:O=Z[26]:G=Z[27]:H=Z[28]:P=Z[29]:R=Z[30]:Q=Z[31]Δ J=4=> U=Z[32]:V=Z[33]:O=Z[34]:G=Z[35]:H=Z[36]:P=Z[37]:R=Z[38]:Q=Z[39]Δ …………………………… J=N=> U=Z[8N]:V=Z[8N+1]:O=Z[8N+2]:G=Z[8N+3]:H=Z[8N+4]:P=Z[8N+5]:R=Z[8N+6]:Q=Z[8N+7]Δ 换为Lbl 1: U=Z[8J]:V=Z[8J+1]:O=Z[8J+2]:G=Z[8J+3]:H=Z[8J+4]:P=Z[8J+5]:R=Z[8J+6]:Q=Z[8J+7]Δ就不出现哪个“Mem ERROR”了。 32 楼：njy135njy yshf 定义一个扩展变量要耗10字节，defm 200不是200X10=2000字节吗？不是200个线元啊？16线元应该是(16+1)×8×10=1360。没有你说的哪么多啊！？ 回 回 回 感谢你的程序，真的是very good！！！ 33 楼：wyqzm yshf (看了你的几句话,我有些不吐不快)我真是不明白你怎么会发那么大的火,我之前已经声明了,高手看了不要见笑,我只不过是个低低手,甚至更低,言词有所冒犯还请你原谅! 1,估计没有看过windows的极限编程,像这样的数据存储,只怕是4800P也存储不了几个,留给有4850P的朋友使用吧?由于定义的数据过多,此类程序用手工输入方法倒还不如常规程序方法; 2,你用的什么方法我不太清楚,但表达式写得如此之长,凭这一点就可以肯定并不是什么好的表达方法,也许我的思路与你的并不一样,编写的程序也就不相同,但说回来,难道一个程序写到人家都看不懂或是不易看懂,那能算是好程序吗?我只是建议,提点意思罢了,用得着发那么大的火吗?朋友,脾气不用那么爆燥吧?呵 3,条理清晰,简单易懂的程序难道就差么?非得写那么复杂?4800P不也可以计算积分吗?就算是写不了,另外给坛中的朋友提供一个相关资料也行吧? 4,我并不是说10的45次方太小了,你大概是没有看清楚吧?我是说太大了,用于计算直线还用得着这个,那是高射炮打蚊子----大材小用了,大概是你自己把简单问题复杂化了吧,请你看看<<全站仪与高等级公路测量>>这本老书吧,它后面的几页也讲述了这个论题,你的经典理论很精湛,我这个低低手是及不上你的了,那本书的计算式也是积分得到的,好像也没有你的这么复杂吧? 5,只能说算法不够灵活,有时取简化公式的前两项就能满足要求,但有些时候仅仅两项是不够的,我好像看到哪里说这个程序是计算前6项,以我认为,一些时候计算6项都是多余,还浪费不少时间,消耗电能;从另外一个方面来说,也许你睦砺鄹哂谒撕芏?但这个算法描述也绝不是一个好的算法描述,之所以就算你有好的理论,却无法用算法表达出来,那也只是枉然,说实话,这个程序我也就扫了一眼,知道有更好的算法来实现,从你的上一回复也能看得出来,计算器你并不了解,程序你更不懂......所以,我也并不是信口雌雄，盲加评论,再说了,我也不想跟你再争论什么,在我看来,既然你懂得这么高深的理论,那么你应该也是一个脾气温和的人,肯请你注意!这里是论坛,是\论\的地方,我并不想为这么一个小东西争论,相信你也是吧! 温馨提示你:在现代道桥路线理论中,除经典的交点导线法之处,还有线元积木法,模式法,关于此类理论,<<测绘通报>>很早以前就有很多相关描述,但使用最多的还是经典的交点导线法和线元积木法,这在国内有名的路线软件如中交一公院的<<纬地道路辅助设计系统>>,南京李方的EICAD,国外著名路线软件德国IB&T公司的CARD/1中已经被广泛使用,如你的理论(也是线元法吧?)得以重用,也许可以开辟路线设计的新篇章了 34 楼：yshf wyqzm： 正如你所说的，“windows的极限编程”的确没看过，它的理论和数据结构能否引用到计算器中本人真不得而知；对计算器回 回 嘛我不了解，程序嘛更谈不上懂，至于什么导线法、线元法、模具式法、拟合法、积木法、综合法、弦切法、圆弧移动法、端点受限法、......，还有什么设计系统等啦更是闻所未闻了，跟你比起来差得不只是一点点，惭愧啊。不过，“windows”操作系统大家都在用，它的程序假如不编译，看得懂的不是大多数人，按你的的观点那是好程序吗？再如，500万像素的数码像机，真懂它原理的人有多少，按你的观点它岂能称为照像机或者好像机吗？我不爆燥也没有发火，既然有人提意见，就得有回复，不然怎称之为“论谈”吗，何况像你这样没看清没看明白的，我更不会那样了！本人要表达只是：没有最好，只有更好，能运行、快速准确得出所要求的结果即可，提意见或建议请透过现象看本质、请看清楚、请看明白了后再进行。是的，仅凭“扫了一眼”就看出半径为10的45次的曲率可以用0代替这简单的问题被复杂化，真神了，真服了你，上次回复不是让你试一下：“1÷0”，试了吗、你知道其中的奥秘吗？这怎能不让人联想起猪八戒吃人生果嘛！是的，计算直线用了那么长的表达式，真是杀鸡用了宰牛刀，为什么要那样做你知道吗？你期望的是用一个以上的子程序来解决坐标计算，这样会思路清晰、条理清楚、简单而易于理解，我想再笨的人也不会如此的，源程序的说明岂不成了画蛇添足，看明白了吗？你所期望的灵活是不是可能以样理解：将面包车的车箱造成可拆卸的，拉货时换上一个货箱，拉人时再换上车箱？程序已从使用范围、精度、速度、储存空间上作了统筹兼顾提炼而成的，难道这一点也没看出来？若论积分计算，4500P前的就能进行，何况4800P乎？用真数进行积分计算意义和速度你知道吗？乌呼，不知道这程序设计的概念是不是在远古的恐龙时代就已产生了，以至于在文明的现代，这样的程序就像4800P计算器本身能进行各种计算的原理一样让人难以看懂，看来开辟新篇章时代已到来！ 35 楼：lshy 能编出这样的程序的确不简单，使用方便，计算准确。到目前为止，还未发现比这更好的呢！ 36 楼：njy135njy 同意lshy，这个程序也是我用过中的最好的！ 37 楼：njy135njy 忘记评分了 38 楼：wyqzm 回 回 评：+5 分 回 回 评：+5 分 你这个人真是不可理喻,我已说了,不想为这个东西争论什么,你是高手高手高高手,行了吧????好像这个世界上就你行一样! 简直无聊之极! 39 楼：yshf 回 是的，我这个人真是不可“理喻”，你不是要“看了半天也没有看明白,请指教”吗？，怎么才指教了几句，就受不啦，不理不论，怎么能辩清可取不可取？这程序对你来说太复杂了，对你不适用，因为“现代程序设计的感点认为,条理清晰,简单易懂的程序才是标准”。你不是“我也就扫了一眼,知道有更好的算法来实现”了，拿出来嘛，见见太阳，别捂着发霉了！哦，忘记了，听你的口气，《测绘通报》你不是经常看吗，这程序的核心算法来源于那里，难道说这也没看出来？天下不是我最行，但天下就不行！我还是要再次提醒你：请看清楚看明白了，再提意见或建议，否则别之乎者也。本人非在取笑你，只是认为太草率了，明白么？ 40 楼：yshf 回 更正一下：”天下不是我最行，但天下就不行！”应为：“天下不是我最行，但天下就你不行！”。还有一点：此程序不是本人写的，但它是采取本人程序的核心算法扩展而成的。 41 楼：wyqzm 回 你不是说人家扯大旗么?那你这又是什么?难道不是更大的大旗?看楼上的,还有人尊称你为\老师\那就是敬重你,既是这样,难道不应该做做表率?而你都做了什么?这样的话只会让人看你没有斗量,为了如此小事就大为光火,不要怪别人说你没有素质,如此这样,我真是为你感到惭愧,我都不在乎你如此说了,你竟然还...... 另外,我也看了其他一些贴子,难道你不觉得自己就是在扯大旗?也许那些发贴的朋友水平是不如你,但你犯不着态度如此之差吧?你说你这不是扯大旗是什么?还对其他发贴的朋友如此强硬,小心风大了闪了舌头,就算你水平如此之高,但做人如此不厚道,态度不谦逊,那只也能说明你这个人俗不可耐,还说什么\但天下就你不行\之类,.......不要那么狂妄自大了,借你的话:\井底之蛙、撼面杖吹火\不要自以为了不起!迟早你会让你的无知自负付出代价的,既然我已经在第一贴就说了\高手看了不要见笑了)\难道你还不明白 42 楼：wyqzm 回 真是搞不明白,自己才来这个地方,小心发点小小的意思(还敬请高手不要见笑),居然还惹来如此苦笑不得的事,真不知道得罪了哪路神仙了,如此这般的话,还是敬请管理员删除我的贴子和ID吧,以免影响其他发贴看贴的朋友! 在此,我向其他的发贴看贴的朋友道歉了. 43 楼：wyqzm 更正一下：”小小的意思”应为：“小小的意见”。 44 楼：wuheqi 回 回 我刚刚在研究你的程序，我对你的程序的数学模型没接触过！ 计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。 我去查了Gauss-Legendre 没查到，我不知道这个是什么含义，具体是什么样的一个原理， 还有Lbl 0：Defm 50:\=> XY\：\=> SZ\：N： N=1=>Goto 1：≠>Goto 2Δ← 这句中，我对1.SZ => XY\：\=> SZ\不理解 45 楼：wuheqi 希望各位朋友能相互学习，我的QQ49619996 回

46 楼：yshf 回 我真是不明白：你这是眼高高、手低低、腔油油、调滑滑、自谓具有X光眼的家伙，知错了怎么还不思悔改，强词夺理、狡辩干啥，见个相似的标准你就套，别迷信是放之四海而皆准，殊不知那是凝似病例，你身体中每一个细胞就像程序中的一行，它们是如何分离、怎样获取营养、又是如何新陈代谢的？你看得明白、你看得懂吗？如果看不懂，依据你的观点，你这人还能称得上是人或者好人吗？狂妄自大的是你，扯大旗的也是你，态度如此之差的也是你，无聊之极的还是你，如此不厚道、态度不谦逊的仍然是你，为自己无知且自负而付出代价的同样会是你，蚊子咬称杆—称嘴硬的不是你还会是谁？在家里、在学校，难道就没有人给你说过：知之为知之、不知为不知；没有金刚钻不揽磁器活，是只乌鸦别扮凤凰，是头猪干嘛鼻子上插根葱装象、是条狗嘴里为何吐象牙、要了贞洁立牌坊干什么呢？对你这样的，真不相信他们没有教导过！乌呼，吾闻池鱼，不知海水深。最后，又不得不再一次提醒你：请看清楚看明白了，再提意见或建议，别再玩什么“小小的意思”了，“坐标并不闹着玩的,它很重要”！ 47 楼：wuheqi 请不要吵了，大家都是高手，请能解释下，我刚才提的意见吗？ 48 楼：yshf 回 回 \=> XY\：\=> SZ\作用是显示两个提示以便于进行操作，即提示：输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算。至于“Gauss-Legendre ”请将电子箱号码贴出来，我给你传一份过来。“Dca×Dcb<0=>该测点在其线元内”，我也没弄明白，还是请教 hanniu_1973吧。 49 楼：wuheqi 谢谢，我的油箱是wu_he_qi@163.com 还有如果你有QQ可以加入我的QQ49619996 50 楼：lshy 请将“Gauss-Legendre”发一份到 lsy53041@sohu.com 好吗。谢谢了。 51 楼：yshf 上述二位，邮件已发，请注意查收。 52 楼：wuheqi 收到了！！谢谢 53 楼：wuheqi 回 回 回 回 回 正算程序中 A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:Z[3]=0.5 是什么意思？ 54 楼：qw-197312 请发一份“Gauss-Legendre”到 qwh-7@163.com 好吗。非常谢谢。 55 楼：njy135njy 回 回 yshy 请给我也发一份“Gauss-Legendre”到njy135njy@163.com中，非常感谢！ 56 楼：qw-197312 请问yshf： 主程序“Lbl A：Isz N：A=Z[8N+3]-M：B=Z[8（N+1）+3]-M：”中的M值是代表什么？好像在此处M没有赋于它含义。 另外，坐标正算程序(包括匝道)已有较成熟的程序且比此程序输入方便，我个人建议只保留坐标反算程序较合适。 第三次修改版我用过，可以计算，但四次修改版我怎么也没有正确使用，不知各位测友是否能正确使用？ 57 楼：njy135njy qw-197312 我的可以正确使用了,我把dat1z改了一下 M值是反算时输入斜交角度,不过好像只能输90,不然会算错,我的就是这样的,不知是什么原因? 58 楼：njy135njy dat1 59 楼：yshf 回 回 回 回 “M”跟本人贴出的、由佚名扩展的“曲线任意里程中边桩坐标正反算(CASIO fx-4850P计算器)程序（修改版）”中的“T”是一回事，在说明中已经指明了：“（5）由于加入了T值，反算时，会发生迭代不收敛现象，也就不可能 得到正确计算结果。”，“M”值的含意是“12、ANG？正算边桩时左右边桩连线与线路中线的右交角”。在第四修改版中，启动程序选择反算后至执行“SUB1”或“SUB2”之前是没有让用户对其进行赋值的，其值为进行反算前保留在“M”中的值。如此处理：将程序“Lbl 2：{XY}：XY：Z[4]=X：Z[5] =Y：N=0”改为“Lbl 2：{XY}：XY：Z[4]=X：Z[5] =Y：N=0：M=90”，反算过程中，当要求输入 “NAG ？”时直接回车即可。“SUB1”的那几个常数是Gauss-Legendre求积公式中的插值系数和节点。大家应根据自己的习惯，选用较优的输入部分来组合程序中的核心算法部分，以获得满意的效果。 60 楼：yshf qwh-7及njy135njy邮件已发，请查收。 61 楼：njy135njy yshf 邮件是收到了，可是没有内容啊？麻烦你重发一下，谢谢！！！ 62 楼：njy135njy 收到了，谢谢 63 楼：qw-197312 回 回 回 回 谢谢njy135njy ，我再按你修改的程序再试试! 谢谢yshf，邮件已收到! 64 楼：njy135njy yshf 发现一个新问题,反算时如是直线段的桩号就会算错,算的结果显示的是后一段曲线的桩号.为什么呢? 65 楼：njy135njy yshf 反算时如直线后是缓和曲线就算错,如直线后是圆曲线就计算无误,why? 66 楼：yshf 回 回 回 \算的结果显示的是后一段曲线的桩号\出错，这可能是主程序中“Lbl A：”至“Lbl B：”之间的问题。至于“反算时如直线后是缓和曲线就算错,如直线后是圆曲线就计算无误”的问题，因没有将该程序输入到计算器中进行测试，故没有发觉，但可能还是与主程序中“Lbl A：”至“Lbl B：”之间的有关 正算及反算的核心程序“SUB1“、“SUB2”除这个M”ANG=”之外，是没问题的。 。 67 楼：yshf \算的结果显示的是后一段曲线的桩号\出错，这可能是主程序中“Lbl A：”至 “Lbl B：”之间的问题。至于“反算时如直线后是缓和曲线就算错,如直线后是 圆曲线就计算无误”的问题，因没有将该程序输入到计算器中进行测试，故没有发 觉，但可能还是与主程序中“Lbl A：”至“Lbl B：”之间的有关 正算及反算的 核心程序“SUB1“、“SUB2”除这个M”ANG=”之外，是没问题的。 68 楼：njy135njy 期望高手们给检查一下 69 楼：hanniu_1973楼主 回 回 回 大家好！国庆放假七天，故本人未能上网不能及时答复，请见凉！我看了wyqzm和yshf好友的争论，我感觉很好！我并不是看别人的热闹，其实网上论坛，就是要发挥其“论”的特点，无争则无论嘛！wyqzm网友提的第一条意见是可取的，但其它的就有点定论过早，我感觉应先看清或先作一番分析后再定论则可取。Yshf网友的分析虽语言有点尖锐但很到位、专业素质高并且底蕴很好，我感觉网上论坛就应如此，否则网上生源令人堪忧！本人有如下几点需说明的：1、本人修改的程序中的数据库使用扩展变量的确占用闪存资源，如果只为了正算全线贯通可取消扩展变量改成真数表示，可为了在反算时亦能全线贯通之不得已的情况下才使用扩展变量的，因为在反算时为了让计算器能自动判断测点所在的线元区域才加入扩展变量的，如大伙有更好的解决办法，请上传交流，本人愿意赐教！2、“Dca×Dcb<0=>该测点在其线元内”的意思是判断该点至该线元始终点的法线的垂距符号是否相异，如是则该点在其中，否则不在该线元中，接着继续判断下一个线元。。。。。3、关于“M”，在只进行反算时可取消直接改成90，因为本人在编程序时另外考虑了平行线的功能故未删除。其实一个程序只要其核心没有问题就算通过，其它一些输入错误或小错误对一个内行的人是不难发现并更正的。故此我希望各位热爱CASIO系列计算器编程的网友，不要局限于拿来主义----现成的程序使用，理应不断探讨更优化的数学模型、测量理论和更多的编程思路，只有这样才能挖掘编程的乐趣！让我们一起共勉吧！ 70 楼：njy135njy hanniu_1973 这个程序我实在是看不明白,所以请你解决一下上边我的问题:帮助一下!! 反算直线时如直线后是缓和曲线就算错,算出的里程和支距为下一段能缓和曲线的,如直线后是圆曲线就计算无误,why? 71 楼：280050876 njy135njy老兄: 你输入的程序能反算?该程序我正算还可以,反算时总显示“Mem ERROR”.真搞不懂什么原因,各位大虾能否指点一二??? 72 楼：hanniu_1973楼主 280050876: 程序本身没问题，也许是你设置扩展变量的总量不对所致，请你好好检查再说！ 73 楼：280050876 hanniu_1973: 依上面你的算例来说,你设置扩展变量的总量是多少?我是按你的程序输的,反算也没通过呀???? 74 楼：280050876 回 回 回 回 回 我不会编程, 建议各位大虾能否集思广义,将上面程序完善一下,使其成为仅输入任一点的三维坐标,就能反算出对应的中桩及边距,同时根据计算出来的边距与理论边桩相比较,求出其是否位于边线上,这就解决了施工现场高挖(填)边坡顶(脚)线放线的问题. 75 楼：hanniu_1973楼主 280050876: 关于你的第一个问题，我在上面的答复中已述，请你自行斟酌，而第二个问题，其实就是线路三维坐标动态控制系统 ，原理地球人都知道，可要编程成集成度高的三维坐标控制系统并不易。如果只是确定平面坐标点相对于线路中线的位置其实很简单，但加入Z坐标（高程）后就不易了，因为每一线元中纵横向设计（竖曲线、横坡、超高、边坡形状等回 等）的数学模型较多并较杂（当然具体线路具体分析亦不是很难，但做为通用型程序就不易了）。当然在专业的通用软件（如VB、C++、Dephi等）或PDA上编程较容易一点，但如果在CASIO系列的计算器上编程（受其容量、性能和速度的限制）就不是那么容易了。同时我亦相信，包括我在内的广大同行亦在做这方面的研究，也许有些人编出的程序不尽完善，正在不断调试中，亦有些人已编成功，但要有偿使用。不管别人编的程序再多完善，终究不是自己的，只有加入自己的东西才算其拥有该程序的内涵，亦可在其中品味编程的乐趣！在此，我建议各位初入行的朋友们，不要只“拿来主义”，应主动出击，多搜集一些测绘先进的理论，建立优化的数学模型和参考别人优秀的例程、学习先进的算法并多动手编程。让我们共勉共进吧！（注：关于测绘先进的论文可登陆WWW.CNKI.NET网站下载） 76 楼：haoji 回 程序还算不错，就觉得正算子程序是不是太笼统和复杂了点，还有几处的显示应跟个人的需要来显示，有些个用常数代替，如M，：“DIST\显示的是什么？ 数据库就个人觉得用扩展变量来保持有点不是很好，因为我们不单只用这条程序，如运行过其他的程序，那里面也有扩展，拉在用这条是不是出错了， 另外程序的运行速度很慢，如果线元多的话，我想在里面循环，判断，那样会更慢！ 77 楼：hanniu_1973楼主 haoji: 至于你提的意见很好，但第一条我不敢苟同，亦许你对上述正算子程序的编制原理理解并不是很到位，其实正算子程序是根据Gauss-Legendre五点公式计算近似插值的，只是公式罗列罢了，使用的变量较多，但亦只能这样呀！你只要弄明白了五点公式就明白了，请你仔细查阅相关数值计算方法的书籍。而对其它的问题，其实任何一种程序都具有优劣性的，我在这里所要表达的是计算器编程的思路，至于具体程序中的编程细节，应依个人喜好而定，如显示、输入和输出等环节。扩展变量型与常量型数据库，各有优劣，应依个人取舍！即扩展变量的条理清晰、使用方便，但其占用容量大、维护性差；而常量型的占用容量相对较小、维护性稍为好一点，但其条理性不如前者。故一个程序不管如何好，总是有局限性的，需要我们进一步跟踪测绘先进理论、优化数学模型并不断扩展， 78 楼：测量员 hanniu_1973: 网上都说这个程序方便，一直没来得急看，冒昧的问一句：那个数据库全部是扩展（Lbl1:J=1=>U=Z[8]:V=Z[9]:O=Z[10]:G=Z[11]:H=Z[12]:P=Z[13]:R=Z[14]:Q=Z[15]Δ）那真正的数据难道还要另一数库与所有的扩展变量对应吗？？？ 79 楼：hanniu_1973楼主 测量员: 如此看来你对CASIO系列的计算器的扩展变量并不是很清楚，其实具体的数据只要程序运行前在COMP状态下输入即可，一次输入终身享用，除非你用 Defm 0清除或重新输入Defm N，如：Z[10]=123.456回车即可。 回 回 回

80 楼：测量员 hanniu_1973: 回 谢谢真还不太清楚,那整个计算器里如果有其它程序用同样的扩展变量,启不 乱套,comp状态下用扩展也得注意了啊.以前真还没么用过....还有那扩展变量的所有值都得一个个在comp状态下赋值吗???? 81 楼：KENT 回 我有一个PC-E500S计算机程序，没有说明书，不知怎么用，我传上来，请各位大师帮个忙，告诉我怎么用， 82 楼：KENT PC-E500S计算机全中文综合曲线计算坐标程序 1000 CLS :WAIT &:TT=4321:C$=\1E64C080FFFE00C03C1E64C080000010D0FFFED0B022A222E2E2A2233220\(&,&7):GPRINT C$;\0\1005 D$=\848486C58000000601FF7F20190E4743FF7F01073C18\(&1,&F):GPRINT D$;\1010 LOCATE &C,&:INPUT \&C,&1:INPUT \1015 GCURSOR (&2,&17):GPRINT \4C080FFFE00C03C1E64C080000010D0FFFED0B022A222E2E2A2233220\(&14,&17):GPRINT A$ 1020 GCURSOR (&1,&1F):GPRINT \000818484848584FFFF8584848486C58000000601FF7F20190E4743FF7F01073C18\(&10,&1F):GPRINT B$ 1025 LOCATE &C,&2:INPUT \&C,&3:INPUT \1030 XCHA=X2-X1:YCHA=Y2-Y1+1E-55 1031 S=SQR (XCHA*XCHA+YCHA*YCHA):T=&B4-&5A*SGN YCHA-ATN (XCHA/YCHA) 1035 GCURSOR (&90,&7):GPRINT C$;\0010D0FFFF901000FEFC0407F6E40604 1040 GCURSOR (&91,&F):GPRINT D$;\7F018340BF8F8181FFFF8181\&1A,&2:PRINT \&1A,&3:INPUT V 1045 CLS :WAIT &:A$=\(&10,&7):GPRINT \2F8F0101010F0F81\1050 GCURSOR (&32,&7):GPRINT \ \1055 GCURSOR (&10,&F):GPRINT \1060 GCURSOR (&32,&F):GPRINT \回 1065 LOCATE &E,&:INPUT \&E,&1:INPUT \A:LOCATE &1A,&:INPUT \1066 X0=LS-LS^3/40/R/R:Y0=LS^2/6/R-LS^4/336/R^3 1067 P=LS*LS/24/R:Q=LS/2-LS^3/240/R/R:B0=90*LS/R/PI 1070 TH=(R+P)*TAN(A/2)+Q:LH=R*(A-2*B0)*PI/180+2*LS:LY=R*(A-2*B0)*PI/180:EH=(R+P)/COS(A/2)-R:DH=2*TH-LH 1071 ZH=V:JD=ZH+TH:HY=ZH+LS:YH=HY+LY:HZ=YH+LS:QZ=HZ-LH/2 1072 LOCATE &1A,&1:PRINT \&,&2:WAIT :PRINT \H;A$ 1076 CLS :WAIT &:GCURSOR (&10,&7):GPRINT 83 楼：KENT 接上 \(&32,&7):GPRINT \010FFFE10FFFE1010F0F810006058CCE772209090FFFE909197D280\\1077 GCURSOR (&10,&F):GPRINT \OR (&32,&F):GPRINT \1079 WAIT &:LOCATE &18,&3:PRINT \&18,&:PRINT \&18,&1:PRINT \HY:LOCATE &5,&2:PRINT \1080 LOCATE &18,&2:PRINT \&5,&3:WAIT :PRINT \& 1085 USING:CLS 1090 WAIT &:GCURSOR (&1,&F):GPRINT \40001010101010F1F69E949090D09018100000F8F010D89C979290D09010F0F810 1095 GCURSOR (&4A,&7):GPRINT \(&49,&17):GPRINT \1F1210101010901810 1100 GCURSOR (&49,&F):GPRINT \(&49,&1F):GPRINT \7F2121212121FF7F01\&10,&1:PRINT \&10,&3:PRINT \1105 GCURSOR (&1,&17):GPRINT \0703C0F030040C07F3F0000000000FF7F001F0F08081F0F4040FF7F\&15,&2:INPUT \O 1110 GCURSOR (&A1,&F):GPRINT \1115 GCURSOR (&A1,&17):GPRINT \ATE &22,&2:INPUT W:CLS 1120 LG=INT (V*&3E8)/&3E8 回 1125 N$=STR$ (INT (LG/&3E8)):N$=N$+\(LG-INT (LG/&3E8)*&3E8):LOCATE &,&1:PRINT N$; 1130 USING \1140 LOCATE &10,&1:PRINT \&10,&2:PRINT \&1:WAIT:PRINT \ 1150 G=INT ((LG+W)/W)*W:CLS:WAIT & 1155 N$=STR$ (INT (G/&3E8)):N$=N$+\(G-INT (G/&3E8)*&3E8):LOCATE &,&1:PRINT N$; 1160 L=G-LG:X=L-L^5/40/R^2/LS^2:Y=O*L^3/6/R/LS 1165 X3=X1-Y*SIN T+X*COS T:Y3=Y1+Y*COS T+X*SIN T 1170 USING \1180 LOCATE &10,&1:PRINT \&10,&2:PRINT \&1:WAIT:PRINT \& 1185 CLS:G=G+W:IF G>HY THEN GOTO 1500 1190 N$=STR$ (INT (G/&3E8)):N$=N$+\(G-INT (G/&3E8)*&3E8):LOCATE &,&1:PRINT N$; 1195 L=G-LG:X=L-L^5/40/R^2/LS^2:Y=O*L^3/6/R/LS 1200 X3=X1-Y*SIN T+X*COS T:Y3=Y1+Y*COS T+X*SIN T 1205 USING \1210 LOCATE &10,&1:PRINT \&10,&2:PRINT \&1:WAIT:PRINT \& 1215 GOTO 1185 1500 LG=INT(HY*&3E8)/&3E8 1501 N$=STR$ (INT (LG/&3E8)):N$=N$+\(LG-INT (LG/&3E8)*&3E8):LOCATE &,&1:PRINT N$; 84 楼：KENT 接上 1503 X=LS-LS^3/40/R^2:Y=O*LS^2/6/R 1505 X3=X1-Y*SIN T+X*COS T:Y3=Y1+Y*COS T+X*SIN T 1510 LOCATE &10,&1:PRINT \&10,&2:PRINT \&1:WAIT:PRINT \& 1515 IF (HY-ZH)>=LH THEN END 1520 G=INT ((LG+W)/W)*W:CLS:WAIT &:IF G>YH GOSUB 1521:GOTO 2000 ELSE GOTO 1530 1521 G=INT(QZ*&3E8)/&3E8:CLS 1522 N$=STR$ (INT (G/&3E8)):N$=N$+\(G-INT (G/&3E8)*&3E8):LOCATE &,&1:PRINT N$; 1523 L=LY/2:C=L*180/R/PI+B0:X=R*SIN C+Q:Y=O*R*(1-COS C)+P:TT=1234 1524 X3=X1-Y*SIN T+X*COS T:Y3=Y1+Y*COS T+X*SIN T 1525 LOCATE &10,&1:PRINT \&10,&2:PRINT \&1:WAIT:PRINT \& 1526 RETURN 回 1530 N$=STR$ (INT (G/&3E8)):N$=N$+\(G-INT (G/&3E8)*&3E8):LOCATE &,&1:PRINT N$; 1531 L=G-LG:C=L*180/R/PI+B0:X=R*SIN C+Q:Y=O*R*(1-COS C)+P 1535 X3=X1-Y*SIN T+X*COS T:Y3=Y1+Y*COS T+X*SIN T 1545 LOCATE &10,&1:PRINT \&10,&2:PRINT \&1:WAIT:PRINT \& 1585 CLS:G=G+W:TRG=G:IF G>QZ AND TT=4321 THEN GOSUB 1521:G=TRG 1589 IF G>YH THEN GOTO 2000 1590 CLS:N$=STR$ (INT (G/&3E8)):N$=N$+\(G-INT (G/&3E8)*&3E8):LOCATE &,&1:PRINT N$; 1595 L=G-LG:C=L*180/R/PI+B0:X=R*SIN C+Q:Y=O*R*(1-COS C)+P 1600 X3=X1-Y*SIN T+X*COS T:Y3=Y1+Y*COS T+X*SIN T 1610 LOCATE &10,&1:PRINT \&10,&2:PRINT \&1:WAIT:PRINT \& 1615 GOTO 1585 2000 CLS:LG=INT(YH*&3E8)/&3E8 2005 N$=STR$ (INT (LG/&3E8)):N$=N$+\(LG-INT (LG/&3E8)*&3E8):LOCATE &,&1:PRINT N$; 2010 L=LG-HY:C=L*180/R/PI+B0:X=R*SIN C+Q:Y=O*R*(1-COS C)+P 2015 X3=X1-Y*SIN T+X*COS T:Y3=Y1+Y*COS T+X*SIN T 2020 LOCATE &10,&1:PRINT \&10,&2:PRINT \&1:WAIT:PRINT \& 2024 IF O>0 THEN K=T+A ELSE K=T-A 2025 IF K>=360 THEN K=K-360 2026 X1=X2+TH*COS K:Y1=Y2+TH*SIN K 2030 LG=INT (YH*&3E8)/&3E8:I=1 2040 IF I=1 THEN G=INT ((LG+W)/W)*W:LG=INT (HZ*&3E8)/&3E8:I=2 ELSE LG=INT (HZ*&3E8)/&3E8 2041 CLS:WAIT &:IF G>HZ THEN GOTO 3000 2045 N$=STR$ (INT (G/&3E8)):N$=N$+\(G-INT (G/&3E8)*& 以上是PC-E500S计算机程序，没有说明书，不知怎么用，我传上来，请各位大师帮个忙，告诉我怎么用，同时供广大测量工作者使用或参考，非常感激。 85 楼：hao728 咱向你们学习，因为你们都是高手！！！ 86 楼：hanniu_1973楼主 测量员： 如果多个程序使用扩展变量共存的系统，需注意：1、如果扩展变量只作中间变量的话，无所谓；2、如果扩展变量作为程序常量的话，应注意各个程序中的变量重名的问题，否则将出现变量相互覆盖的现象，因为CASIO系列的计算器系回 回 统，其扩展变量是采用不受保护的堆栈原理设计的。是的，在COMP状态下需一一输入扩展常量，不过只要不受初始化、Defm 0或Defm N或其它扩展变量的覆盖，其可永远使用。 87 楼：wang811736 请给我也发一份“Gauss-Legendre”到121767632@qq.com,谢谢！ 88 楼：radio_man 一个弯道两边缓和曲线不等长能算吗? 89 楼：测量员 yshy老师，您好！Gauss-Legendre公式能发给我一份吗，我的email:luozixiang@sina.com,谢谢！ 90 楼：wuheqi 回 回 回 回 评：+5 分 我认为你的主程序的 Defm 50的设置没必要,有点画蛇添足的味道.你在输入曲线要素的时候这个值的设置不是在程序里设置的 91 楼：wuheqi E=180/Π 还是N？ 92 楼：yshf radio_man：那个弯道不就是由三个线元组成吗，当然能算。 liuyouchun：你后面的那个问题，我没有相关的资料和程序。 wuheqi：“E=180/Π”中不是N，乃圆周率（3.14159267....)的符号也。 93 楼：陈秋 yshy老师： 请给我也发一份“Gauss-Legendre”到c_q_1_2@163.com好吗？ 谢谢啦！！ 94 楼：hanniu_1973楼主 wuheqi： 你说的程序中的“Defm 50”，的确在有些地方是可以取消的，但对于程序代码的可读性和安全性，有时是有必要的，具体情况应具体分析才行！ 95 楼：hanniu_1973楼主 各位网友： 因本人这边缺乏超过100个线元或25个平曲线以上的平、纵、 横面的线位数据进行程序调试，请问各位能提供这方面的资料吗？或者本愿意以《曲线任意里程中边桩坐标正反算(4850P计算器)程序(第五次修改版）》交换，不胜感谢！有意者可发至我的邮箱：hanniu_1973@163.com 回 回 回 评：+5 分 回 回

96 楼：陈秋 收到了谢谢！ 97 楼：wuheqi 今天对你的程序进行了调试,正算可以成功通过,但发算也成功通过! 回 回 现在开始研究你各个子程序,但在仔细阅读的过程中我碰到了一个问题,就是正算子程序中通过5点公式的推导出程序中的X 和Y值,对于这个问题我已经反复推导了好几次都无法正确的推导到程序中所写的公式,希望你能发一贴推导公式,或者能否把推导公式发一份到我的油箱wu_he_qi@163.com 98 楼：lizengrong 程序终于调试成功，我对你万分的感谢。 99 楼：hanniu_1973楼主 lizengrong: 恭喜你调试成功！但你有没有在我的第四版的基础上进行扩展呢？我本来是想通过网络共享的原则进行相互交流的，但我已发布一个月有余了，都没有见到一个像样的扩展程序，没有达到我预期的效果，我亦感到有点遗憾！其实，第四版中仍有许多需要改进的地方，我现提供给大家一些思路： 1、对于数据库：可分为扩展变量型和常量型。各有优缺点，扩展变量型：运行速度较快、维护方便，但数据安全性较差、占容量大；变量型：数据安全性较好、占容量相对较小、易于扩展，但其维护不方便、运行速度相对差一点等。 2、对程序功能：已有程序只有正反算坐标功能，为了更好的测放线提供数据，发必须对其进行扩展，比如：加入方向坐标法或称虚拟测站法、平行线解算和曲线求交功能，以任一角度任一地点解算坐标、里程和偏距。更高级别一就是整合平、纵、横面线位数据在CASIO系列（最好是4850P，因4800P容量太低，如线元多实现起来有点难度）上实现线路三维坐标动态控制系统 希望热爱CASIO系列编程的朋友继续努力为测工编制更多的便用的测量程序服务于现场工作！ 100 楼：wuheqi 回 回 回 评：+5 分 hanniu_1973 我对你的反算程序理解不通,能否具体的讲讲你的反算思路,!谢谢,如果有QQ,请可以加我49619996,加的时候请写清楚你的名字.希望能得到你的指点! 101 楼：wuheqi 有没有 Simpson公式 的介绍? 102 楼：lovewangziqiang 我们这里左右幅,路基隧道加桥梁,切多为斜交所以感觉不太实用, 哪位高手改一下,把扩充去掉,只保留核心部分(有左右边桩连线与线路中线的 右交角大于小于90度的边桩坐标正反算);我的邮箱lovewangziqiang@sina. 回 评：+5 分 回 com,lovewangziqiang@163.com有劳各位高手了! 103 楼：lovewangziqiang 急用,谢谢了 104 楼：hanniu_1973楼主 lovewangziqiang: 照你说的，你可参照我在本网站里的另一程序：《利用Simpson公式和垂距趋近正反算曲线程序》，其可正反算任意角度、支距的坐标点。 105 楼：lovewangziqiang 为什么坐标反算时,计算器没有反应了呢?还是速度太慢了. 106 楼：hanniu_1973楼主 lovewangziqiang: 反算时，的确有点慢，因为其一般是逐步趋近解算里程和偏距的，计算时间与距离待测点的长度成正比。但一般不超过30秒，如果超过的话，首先检查你的程序结构是否正确。如果能反算成功，但其速度太慢，需优化数据结构或算法。如果程序计算数据结构本身错误，需进一步检查、更正并不断调试，直到调试成功并满足运行速度和精度要求。 107 楼：ZHAOYU009 yshy老师： 请给我也发一份“Gauss-Legendre”到ZHAOYU009_002@163.COM好吗？ 非常感谢 108 楼：zhaoyu009 邮件已经收到,谢谢,非常感谢 109 楼：罗自翔 hanniu_1973 : 我发给你的东西你收到了吗??? 110 楼：kent 回 回 回 回 评：+5 分 回 回 评：+4 分 回 回 评：+5 分 请问各位大师，我正在用CASIO fx-4850P计算器施入此程序，发现程序S≤Z[8（N+1）+2]=>J=n 中的 J=n 小写 n 不知怎么样施入，请各位大师指教， 111 楼：kent 程序S≤Z[8（N+1）+2]=>J=n 中的 J=n 小写 n 不知怎么样施入，请各位大师指教， 回 112 楼：风雨屋 请问.........................怎么输入！！！！！ 113 楼：kent 回 回 请问yshy大师，我正在用CASIO fx-4850P计算器施入此 程序， 发现程序S≤Z[8（N+1）+2]=>J=n 中的 J=n 小写 n 不知 怎么样施入，请大师指教. 114 楼：wuheqi 你好，我现在正在考虑分离和匝道的条件判断，但总是觉得不怎么理想，希望能得到你的指点！ 115 楼：wuheqi Lbl 1：{SZ}：SZ S≤Z[18]=>J=1：Prog “DAT1”：Goto 3Δ S≤Z[26]=>J=2：Prog “DAT1”：Goto 3Δ S≤Z[34]=>J=3：Prog “DAT1”：Goto 3Δ S≤Z[42]=>J=4：Prog “DAT1”：Goto 3Δ …………… S≤Z[8（N+1）+2]=>J=n：Prog “DAT1”：Goto 3Δ 对于这句可以如下处理 Lbl 5：{SZ}：SZ:N=1:Lbl 1:Isz N:S≤Z[8(N+1)+2]=>J=N：Prog “DAT1”：Goto 3:：≠>Goto 1 数据库的ＤＡＴ１也可以这样处理 116 楼：KENT 请问各位大师，怎样才能调试成功，请指教。 117 楼：wuheqi 按步骤一步一步输入,如果遇到问题不要急着输入,赶紧查书看看格式有没有对, 118 楼：KENT 很难调试成功啊，有那位大师调试成功的传份上来参考啦，让大家一起学习吗 119 楼：G4WWWW yshf老师请你告诉我Lbl1:J=1=>U=Z[8]:V=Z[9]:怎么输入？ 是这样吗Lbl1:J=1=>U=Z[8]=1235.45:V=Z[9]=4534.455:谢谢。 120 楼：yshf Lbl1:J=1=>Z[8]=1235.45:U=Z[8]:Z[9]=4534.455:V=Z[9]:...... 回 回 评：+5 分 回 回 回 回 回 121 楼：wang811736 请给我也发一份“Gauss-Legendre”到121767632@qq.com,谢谢 122 楼：wuheqi 回 回 我现在通过原稿的修改改编很自己的程序,已经综合考虑了匝道,左右线分离的情况，正算速度很错,但反算速度还是差了一步,觉的有点慢!! 另外在反正程序中，我对起过程还是了解不是很清楚,能有有象Gauss-Legendre一样的资料说明,只有理解了， 才可以对程序进行跟好的改良! 123 楼：hanniu_1973楼主 G4WWWW: 你好！我在此纠正一下yshf朋友的答复如下：像“Lbl1:J=1=>U=Z[8]:V=Z[9]:”这样的语句中的Z[8]、Z[9]，其输入数据是在程序运行前在COMP状态下直接输入Z[8]=1234.56，其实我以前早就答复过了，请你仔细看看我上面的答述。good luck! 124 楼：胡跃籍 回 评：+5 分 回 好东西，我还没看太明白，测量大军里有你们这些慷慨的人无私奉献你们的知识，是大家的测量水平会更上一个台阶，衷心的谢谢你们。 125 楼：g4wwww hanniu_1973 谢了,我有不明白的地方再问你可以吗? 126 楼：g4wwww hanniu_1973 老师，我运行程序为什么不能反算！正算以，反算时好想只能算一段？ 不知道为什么计算下一段线元就不正确了。 127 楼：wuheqi 回 回 评：+5 分 回 如果这样的话,我判断你的程序问题不大,希望你能检查一下你的曲元要素有没有正确,另外反算是你的程序有没有 联接正确, 128 楼：dengyangya hunniu老师： 请给我也发一份“Gauss-Legendre”到dengyangyang@sina.COM好吗？ 非常感谢 129 楼：g4wwww 回 评：+5 分 回 hanniu_1973老师: 请您给我也发一份“Gauss-Legendre”到g4wwww@163.com好吗？ 谢谢 130 楼：zongze hanniu_1973老师: 请您给我也发一份“Gauss-Legendre”到www.zongje@126.com好吗？ 谢谢 131 楼：淡天 计算斜交时只要是Z=负值就会计算出错,这是什么原因! 是本程序的原因吗? 132 楼：xb65785355 hanniu_1973大师:请给我也发一份“Gauss-Legendre”到cexoyga068@.sona.COM好吗？ 非常感谢 133 楼：zz05 hanniu_1973老师: 请您给我也发一份“Gauss-Legendre”到zhangzong05@sohu.com好吗？ 谢谢 134 楼：淡天 hanniu_1973老师: 能帮我看一下吗?为什么计算左边桩斜交角度会出错啊! 135 楼：小沈 我编的方位角程序只显示Ｓ，不显示方位角Ｊ，请各位前辈指教． 程序如下（４８００）： ～～～ＸＸＸ～～～ ＦＩＸＭ：ＰＯＬ（Ａ－Ｃ，Ｂ－Ｄ）：Ｉ＂Ｓ＝＂▲← Ｊ≤０＝＞Ｊ＝Ｊ＋３６０：＂Ｊ＝＂▲ 姓名：小沈 05-12-13 17:6 还有就是４８００上的（ＤＭＳ）键再那里 两个都要 回 回 回 回 回 回 评：+5 分

136 楼：xb65785355 lizengrong:在上面看到你使用成功了,非常羡慕,能和你谈谈使用的体会吗?我的qq号是:549566755 137 楼：XIONGSR 公式实用性强 138 楼：XIONGSR 公式实用性强 139 楼：g4wwww hanniu_1973老师 我需要全线贯通的反算功能，但是上面这个程序我怎么都没能正确使用。你看 这个程序的反算功能，能改成象第3次修改那样吗？正算我都通过了。 我的一妹儿是：g4wwww@163.com 140 楼：yrc 请问小写的n字怎样输入，还有...........怎样输入,我是个菜鸟,请各位大虾帮助. 141 楼：yrc8888 各位大虾,如何才能调试成功啊? 142 楼：yrc8888 请问如何才能看到网友的全部留言? 143 楼：yrc8888 yshy老师，您好！请给我发一份“Gauss-Legendre”到yrc8888@163.com,谢谢！ 144 楼：yrc8888 回 回 回 回 评：+5 分 回 评：+5 分 回 回 回 回 评：+5 分 这个程序终于能运行了,正反算都可以,但就是运行结果对不到示例结果,郁闷啊!我也不知那里有问题了,不知谁能把输入示例线元要素和计算方式的程序完整的上传一遍那就太好了. 145 楼：yrc8888 回 请问主程序中的N和正算子程序中的N运行时会冲突吗?为什么我运行了一次程序后第二次运行时N=0.2844444444?另外A、B、N、K、L是固定值吗?请各位网友帮助!!! 146 楼：洛叶 hanniu_1973 请教一下， 回 评：+4 分 怎样才能在你的程序里只要输入一次已知数据 就可以计算对称曲线和非对称曲线的 中边桩坐标？如果能同时输入L01=?和L02=？好了 147 楼：lxwu hanniu_1973老师你好,我是菜鸟一只,用你的程序碰到断链,我该如何处理,能帮我解难吗?谢谢!! 148 楼：FDS_1 yshy 请给我也发一份“Gauss-Legendre”到.fengdisheng123@126.com中，谢谢 149 楼：bing507 回 回 回 太好了，我是新手，不能帮大家什么，但是我觉得这里真的是太好了。感谢无私奉贤的优秀测量工作者 150 楼：zbweixkj hanniu_1973老师： 我按你的程序一步步输入进去，可不知为何仅在第一二个曲线元是正确的再往后就不对了，能说明一下不，谢谢！！！ 151 楼：yzj37 yshy老师，您好！ 请给我发一份“Gauss-Legendre”到yzj37@163.com,谢谢！ 152 楼：wuheqi 回 回 评：+5 分 回 评：+5 分 我认为我对这个程序比较了解,还通过自己的想法对改程序进行了改进,如果有谁想交流的,可以 加的QQ:49619996. 谢谢,有好的建议都可以相互讨论, 加的时候请写名2加入原因,因为一般我不加人 153 楼：wuheqi 我的QQ:49619996 154 楼：zbweixkj hanniu_1973老师： 不知道怎么回事，我按你的程序一步步输入后，正算及反算均只在第一、二个线元内结果是正确的，超过后结果是错误的，能否说明一下不？谢谢。。 155 楼：wuheqi hanniu_1973老师 我们可否做一下深入的讨论,我的QQ:49619996 156 楼：wuheqi 大家对CASIO编程序感兴趣的 可以加这个群15514898 回 回 回 回 157 楼：zch720920 程序调试基本通过，但是在使用过程还出现如下几个问题，烦请hanniu_1973: 老师及各位大虾帮助解决。1、正算时程序运行过程中如果退出再进入的话，会出现 提示重新输入“X0”、“Y0”、“S0”……“Q0”等的问题，且有时会提示 子程序“DAT1”出错。2、正算中计算最后面两段线路中的点位时算出的结果精 度与您给出的算例中的结果差了不少（最多的差十几厘米），而且越往后差得 回 越多。3、正算中同一个位置两次算出的结果有时会不一样，这让我感到心里有点没底，不没底，不知道要用哪一个结果。3、反算时速度太慢。 158 楼：wangbin0798 回 评：+5 分 能不能给我发份完整的程序（可以全线放线的那种，坐标，标高可以连起来用的）我用的是4800。我的邮箱（yyyj1119@sohu.com) 谢谢了！！ 159 楼：lxwu hanniu_1973老师你好 这个程序真是太棒了,但在使用过程中发现一个问题: 当测点在第一线元起点的断面上,反算会出错. 160 楼：lxwu 我把Z[6]×Z[7]<0=>J=N改成Z[6]×Z[7]≤0=>J=N后 第一线元起点的中点的反算又是正确的,在其断面上的其它点反算还是出错. 能指点迷津吗? 161 楼：zch720920 怎么每一条回复内容中每行的最后几个字都被页面的最右边线挡住了呢？怎样才能 整的看到回复内容中的每行字呢？恳望各位大虾给予指点，谢谢 162 楼：wx2wxy 猛，我顶 163 楼：lxwu hanniu_1973老师你好 为什么最后一段线元上的点反算会出现Mem ERROR 而前面各段正反算均正确,能指点迷津吗? 回 回 回 回 回 评：+5 分 164 楼：Wenjian_Cao 回：zch720920 网页有BUG 现在弄好了！ 165 楼：wendi1685 发现问题: 1.正算精度还可以.有一个桩号780中.边桩坐标计算有误. 2.反算桩号和中距误差很大! 是不是还有什么问题? 166 楼：wuheqi 你把程序的Z[6]×Z[7]∠0改成Z[6]×Z[7]≤0 167 楼：forlu1314 hanniu_1973， yshy老师： 感谢你们两位的无私和慷慨 回 评：+5 分 回 评：+2 分 回 回 昨天我刚下了yshy老师的程序，也就是“曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序”今天早上主程序的子程序一可以正常运行，到晚上全部运行成功，先来报个喜。 真的是从内心里感谢你们 168 楼：forlu1314 hanniu_1973老师 请问你的程序可以放样的时候算出距离和方位角 如果可以的话就放样就非常方便了。 169 楼：唯能 hanniu_1973大师:请给我也发一份“Gauss-Legendre”到chenhongwu888@hotmail.com好吗？ 非常感谢 170 楼：wangbin0798 谁能告诉我？ 这个程序用在4800上，哪些地方要改？？？ 171 楼：wangbin0798 请教： 回 回 回 评：+5 分 回 ”FWI=”: F=F-M :F→DMS◢ 输入到4800里面因该改成什么？？？ 172 楼：wangbin0798 回 我4800里输入这程序，调试了好几遍还是不行呀，有哪位师傅可以把调试好的上传来学习下。发到我邮箱也行呀(yyyj1119@sohu.com) 173 楼：冬木 楼主 请问有没有4800计算高速公路上普通曲线的正算和反算程序。。。。。。谢谢。。。。有的话来个伊妹儿dongmujiyu@sohu.com 174 楼：ysh328328 缓和曲线的中点坐标算出，如何计算边桩坐标，有没有简单的算缓和曲线边桩坐标程序？ 175 楼：zbweixkj 非常感谢hanniu_1973大师，经过多次摸索，程序已基本调试通过，再次感谢。。 176 楼：sky 回 回 回 回 评：+5 分 非常感谢HANNIU 1973您的程序,我把范例计算一下正反算都可以很好算,为什么反算时919~999这段就不运行了?请帮我指点! 177 楼：JJL-XD 楼主，有没有适用CASIO FX-6300G的坐标计算程序呀？ 178 楼：hzg1237 yshy老师，您好！ 请给我发一份“Gauss-Legendre”到hzg1237@163.com,好吗？谢谢！ 179 楼：xiao hanniu 老师,你好! 请发一份\给我好吗? 我的QQ号是583279424 谢谢! 180 楼：放飞的鸟 哪4850呢？如何用扩张名呢？请高手指教~ 181 楼：谢了 回 评：+5 分 回 回 评：+5 分 回 回 hanniu 高手,如有方便,也请发一份\至本人油箱,不知妥否?如妥,本人将不胜感激!在此先谢了!本人油箱:1977_xxy@163.com

182 楼：lbs326 这个程序都说好用，我先试着用一下再说。 183 楼：Yangyb_001 yshy老师 请给我也发一份“Gauss-Legendre”到.yangyb_001@yahoo.com中，谢谢 初学者 184 楼：Yangyb_001 回 评：+3 分 回 回 4800里输入这程序，调试了好几遍还是不行呀，有哪位高手可以把调试好的上传来学习下。十分感谢,我邮箱yangyb_001@yahoo.com 185 楼：gcg315 请问楼主，能出示你的QQ号吗，以方便联系！ 186 楼：gcg315 yshy老师，您好！请给我发一份“Gauss-Legendre”到liuyouchun88@126.com,谢谢！ 187 楼：gcg315 yshy老师，您好！请给我发一份“Gauss-Legendre”到lightspring315@yahoo.com.cn,谢谢！ 188 楼：名字还没想?/font> 各位高手,谁有调试了的4800的程序,发一份给我.在下不胜感激. 我的邮箱88868498@sina.com 谢谢 189 楼：Yangyb_001 各位老师,谁有调试了的4800的程序,发一份给我.在下不胜感激.email:yangyb_001@yahoo.com 190 楼：Yangyb_001 各位老师,谁有调试了的4800的程序,发一份给我.在下不胜感激.email:yangyb_001@yahoo.com.cn 191 楼：测量新手 楼主：我有一套4850放样程序，但不完整，想请您帮忙改一下，可以吗？ 我还有两套4800程序想改成4850程序，但不知道怎么改，请您指点一下。 联系HDF13896667799@163.com QQ:310141891 回 回 回 回 评：+10 分 回 回 回 192 楼：hyz2008520 yshy老师: 回 我是新手,能给我发一个能全线放样的程序吗,还有不对称缓和曲线也行的,谢谢了.我的邮箱是fang570520@163.com 193 楼：hyz2008520 YSHY老师 U.W.G.Q.E.C.D在正算子程序运行时的符号代表什么啊,还有 反算的G.V,麻烦说明一下,我是初学者什么都不会,你可不要笑我啊 194 楼：ahui9225 yshy老师，hanniu_1973老师.您好！请给我发一份“Gauss-Legendre”到ahui9225@sina.com,谢谢！ 195 楼：chentaotsc 你好 hanniu_1973老师 我已用你的程序进行了输入运行，正算全部通过，但是反算无论如何就是不能正确算出，开始的直线部份能算出，但是后几段曲线反算桩号确实不能通过，我反复对了好几遍，我们这里的同事也一起进行输入过，就是不能进行反算。 还有请老师发一份Gauss-Legendre的公式到我的邮箱好吗？ 我的邮箱是yunfeiyangtsc@sina.com 如哪位知道的朋友能发一份感激不尽：） 如哪位反算在4800上通过的能介绍一下吗？ 不胜感激~！！！ 196 楼：溜洞 你好 hanniu_1973老师 请问在存在断链和坐标换代时怎么处理数据的选择啊 197 楼：测量新手 回 评：+8 分 回 回 回 回 第一次看到有这个论坛，感觉很兴奋、感兴趣，有大开眼界的感觉。我是个测量新手刚从学校毕业出来不久，对于4800编程我不懂。但我很想学它。 198 楼：小白 老师你好： 帮个忙吧，多谢了！ 这个导线平差只能计算10未知点以内的导线点数，我想计算10个以上的导线的点数。怎样来修改成序才行啊，我最初按程序说明去修改（本程序是按10个未知点来设计的，若未知点数超过10个时，按“4倍点数”来改回 变源程序中的“Defm 40）。可是不能正确计算，出现死循环啊。望个位老师来给于修改。多谢了!!!!!!!!!!!我的邮箱：baiyueming123@163.com

注：（摘自）广东交通学院精品课堂。

BI DAO XIAN（闭导线）

1．用途：该程序是“闭合导线平差”的计算程序。

2．程序清单：

Defm 40：X“X0”：Y“Y0”：F“A0”：N“UNKNOWN POINT”：N=N+1↙

I=0：E=0：C=0：Lbl 0 ↙

I=I+1：{L，D}：Z[2I-1]=L：Z[2I]=D：E=E+L：C=C+D：I

I+I+1◢ F=F+Z[2I-1]：F 180 F=F+180： F=F-180 Δ F 360 F+F-360 ΔZ[2N+2I-1]“A”=F+（I-1）B◢ I

I+I+1：Z[2I-1]=Z[2I]cosZ[2N+2I-1]：G=G+Z[2I-1]：Z[2N+2I]=Z[2I]sin Z[2N+2I-1]：H=H+Z[2N+2I]：I I=0：Lbl 3 ↙

I+I+1◢ Z[2I-1]=Z[2I-1]+U Z[2I]：Z[2N+2I]=Z[2N+2I]+V Z[2I]：X=X+ Z[2I-1] ◢ Y=Y+ Z[2N+2I] ◢ I

3．程序说明：

X0，Y0，A0——已知的起始点的纵、横坐标及起始边的坐标方位角。

UNKNOWN POINT——待求未知点的个数。

L，D——观测的导线各左转角及平距。 FW，R，FR——导线角度闭合差、角度闭合差的允许误差的系数及允许误差值。 I，A——各条导线边序号及其坐标方位角。 FX，FY——坐标增量ΔX、ΔY的闭合差。 F，K——导线全长闭合差、导线全长相对闭合差。 I，X，Y——各导线点的序号及其纵横坐标。 注： （1）本程序是按10个未知点来设计的，若未知点数超过10个时，按“4倍点数”来改变源程序中的“Defm 40”。 199 楼：长江不肖生 我是一新手,程序中第一行中的1,sz和2,xy这两点都是逗号吗? 谢谢 200 楼：长江不肖生 好 201 楼：feitiansos 老大，给我发一份 feitianyanting@163.com 小弟不胜感激！！！ 202 楼：wx2wxy 回 回 回 评：+5 分 回 这个我不知道会怎么样!能不能给我一份关于FX-4500PA的篇写程式正反算的啊!有给我的话请发wx2wxy@163.com .那样小弟我不胜感激啊! 203 楼：wanghuachao hanniu_1973 师傅： 你好！我是一个测量爱好者！看了你的程序！我十分佩服你对测量的研究那么的深奥！你是 我的 偶像！最后我有一回 个请求！能否把你的测量程序发给我一些呀！我的邮箱！wanghuachao@sohu.com.本人表示感谢！ 204 楼：lovethegod yshy 请给我也发一份“Gauss-Legendre”，在这里先谢谢了！ 邮箱地址：wween@163.com 205 楼：lovethegod 回 评：+10 分 回 我对编程可以说还不入门，我在输入这个程序的时候可能很多地方没有理解对作者的本意，不过我真的很想使用这个程序，哪们大哥有用4800P调试好过的吗？能否发一份给我，真的感谢了！ 邮箱地址：wween@163.com 206 楼：feitiansos 各位高手小弟想请教个简单问题请解答,先谢了哈!! 就上面的算列我把取Defm 50,正反算S0 X0 Y0 F0 LS R0 RN Q 500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0(正反算均通过)但在第二段和后面的曲线段都出现Ma......返回程序Lbl 2 ：U\：V\：O\：G\：H\：P\：R\：Q：C=1÷P：D=(P-R)÷(2HPR)：E=180÷π 在(2HPR)和E=180÷π之间的冒号上闪烁,我检查了程序都没问题,苦恼了几天,请解答!!!! 207 楼：feitiansos 请老师给我发一份该程序谢谢 feitiansos@163.com 208 楼：lovethegod 回 回 评：+5 分 回 请问一下：FW=”:F→DMS 在4800P里面怎么才能输入啊，我对4800和4850里的不一样的编程符号不知道怎么输，还有在4800里输入些程序还有哪些地方是需要改的啊，望好心人能告诉我，我急需这个程序啊！ wween@163.com QQ：28465609 209 楼：gui 请给我也发一份“Gauss-Legendre”，在这里先谢谢了！ 邮箱地址：gts456422@sohu.com 210 楼：gui hanniu_1973老师: 回 回 评：+5 分

本文来源：https://www.bwwdw.com/article/yxw6.html