道路中边桩坐标放样正反算程序ROAD--2.doc

道路中边桩坐标放样正反算程序ROAD-2（0908改进版）

一、前言
本次的“道路中边桩坐标放样正反算程序ROAD-2（0908改进版）”可理解为以下两个程序的升级版本：
1．HYPERLINK"http://user.qzone.qq.com/595077/blog/1251078433"\t"_blank"路线坐标放样计算程序ROAD-1（09-08改进版），09年8月24日发布
2．HYPERLINK"http://user.qzone.qq.com/595077/blog/1239980749%20­"\t"_blank"道路中边桩坐标放样正反算程序ROAD-2，09年4月17日发布
09年8月24日（也就是几天前）发布的“路线坐标放样计算程序ROAD-1（09-08改进版）”，是对《CASIOfx-5800P计算与道路坐标放样计算》书中的ROAD-1程序的改进，其改进的一些编程技巧和程序优化同样也被本程序所采纳。
09年4月17日发布的“道路中边桩坐标放样正反算程序ROAD-2”当时也是作为ROAD-1的升级替代程序出现，本次发布的ROAD-2较之前的版本又有较大改进，除了采用前面所述的一些编程技巧和程序优化，以及老版本的ROAD-2程序的特点之外，还有如下一些特点：
1．采用交点定位技术，程序开始执行输入一个定位桩号K0用于定位交点，后面的坐标正、反算均在该交点范围内进行，曲线要素计算只在程序开始时运行一次，避免了每次计算桩号都要重新定位交点数据、计算曲线要素的麻烦，提高了计算速度，虽然存在每次执行程序只能在一个交点计算范围内计算，要计算另一个交点范围的桩号，要重新执行程序的缺点，但比起由此带来的优势，这是微不足道的（交点定位技术还有其它优点，后面还会提到）；
2．交点定位技术的采用，结合交点计算桩号范围的确定，可圆满解决纯直线路段、断链、卵形曲线等公路中比较特殊的情况，具体解决方案会在后面的日志中详细解读；
3．数据库子程序采用矩阵变量的形式，使数据组织更加简洁，输入更加便捷；
4．程序考虑了调用数据库子程序和手工输入交点数据两种方法。

总而言之，ROAD-1程序功能要相对单一一些，可作为读者的道路计算入门程序进行学习和应用，而ROAD-2程序完全包含了ROAD-1的所有功能和特点，是可替代ROAD-1的更高级的版本，功能更加丰富，基本上覆盖了公路路线各种情况的坐标计算。

二、程序清单
道路坐标放样计算程序包括一个主程序和四个子程序（不含交点数据子程序）。
本程序中若有与ROAD-1相同的名称的子程序，则以本程序为准或对照修改。
1．主程序清单：ROAD－2

2．子程序1清单：ROAD－QXYS
功能：计算非对称型道路平曲线要素及主点桩号。

注意：这里删除了原ROAD-QXYS子程序第8行计算圆曲线长度的代码，请注意更改。

3．子程序2清单：ROAD－LXZB
功能：计算道路中线坐标。

注意：第16行笔误，改为：Z<0＝>Z+360→Z:Z>360＝>Z-360→Z。

4．子程序3清单：ROAD－XY
功能：计算缓和曲线上一点的局部坐标。

5．子程序4清单：ROAD－XS
功能：显示计算结果。

本子程序名称应为“ROAD-XS”，特此更正。

三、数据库子程序的编写
1．道路示例项目基本资料
基本资料同《CASIOfx-5800P计算与道路坐标放样计算》第6章HY高速公路第2合同段（合同段起止桩号：K4+800~K9+600）。这里摘取直线、曲线及转角表资料如下（若图片不清晰，请参见参见教材P161附录1）：

.

2．路线数据库子程序的编写与格式说明

注：路线数据库子程序ROAD－DATA1是根据计算示例项目的直曲表编写，大家使用时应按各自项目的直曲表改写或新建，本数据库子程序采用了简化坐标，根据项目的坐标范围，小数点前只取了四位数。
数据库采用给矩阵变量MatA赋值的形式，使数据组织更加简洁，为减少数据库子程序的输入量（因为要经常根据不同的项目编写数据库子程序），矩阵变量因子赋值给字母变量的程序代码写在主程序中（见ROAD-2程序第4行）。
矩阵变量简介：5800P计算器开始有矩阵计算功能，内部有七个矩阵存储器（MatA是其中之一），每个矩阵的行列数最大是10×10，有关矩阵具体内容见《CASIOfx-5800P计算器与道路坐标放样计算》34-38页。
数据库格式说明：

本交点计算起终点桩号的说明：
一般情况下，本交点的计算起点是指前一交点的HZ点或YZ点，本交点的计算终点是指后一交点的ZH点或ZY点。之所以要强调是一般情况下，是因为有些特殊情况下，要灵活处理，比如：整个路线就是一段直线、断链的处理、卵形曲线的处理等等，这些特殊情况后面我会以专题日