坐标法面积、周长计算程序
- 格式:doc
- 大小:22.50 KB
- 文档页数:1
坐标法面积、周长计算程序(命名为ZBMJ)
第1行:FreqOn:ClrStat:ClrMemory
第2行:“N=1”: 1→N: “X=”?→C: “Y=”?→D: C→A:D→B:A→List X[1] :B→List Y[1]
第3行:Lbl 0:“N=”:N+1→N◢“X=”?→E: “Y=”?→F:E→List X[N] :F→List Y[N]
第4行:Pol(E-C,F-D+10^(-45)): I→G:J→H:L+G→L
第5行:Pol(E-A,F-B+10^(-45)): E→C:F→D:N=2=>J►DMS◢
第6行:“PMT=”: N=2=>0→I: L+I◢
第7行:“ARE=”:S+GIsin(J-H)÷2→S:Abs(S)→List Freq[N]◢
第8行:Goto 0
N——取样点的顺序号
PMT——多边形的周长
ARE——多边形的面积
X——对应点号N的取样点X坐标
Y——对应点号N的取样点Y坐标
以上为5800程序,在9860中程序如下
第1行:ClrList:0→L:0→S:0→List 3[1]:0→List 4[1]
第2行:“N=1”: 1→N: “X=”?→C: “Y=”?→D: C→A:D→B:A→List 1[1] :B→List 2[1]
第3行:Lbl 0:“N=”:N+1→N◢“X=”?→E: “Y=”?→F:E→List 1[N] :F→List 2[N]
第4行:Pol(E-C,F-D+10^(-45)) : List Ans[1]→G:List Ans[2]→H:L+G→L
第5行:Pol(E-A,F-B+10^(-45)) : List Ans[1]→I:List Ans[2]→J: E→C:F→D:N=2=>J►DMS◢
第6行:“PMT=”: N=2=>0→I: L+I→List 3[N]◢
第7行:“ARE=”:S+GIsin(J-H)÷2→S:Abs S→List 4[N]◢
第8行:Goto 0
2.计算实例:
N=1 X = 940.5392 Y = 898.4861
N=2 X = 1114.8811 Y = 996.4907 PMT=200.000 ARE=0 (29°20′31.54″)
N=3 X = 1065.8788 Y = 1083.6617 PMT=523.607 ARE=10000
N=4 X = 1338.2253 Y = 1007.3243 PMT=995.153 ARE=20000
N=5 X = 1261.8880 Y = 734.9778 PMT=1226.241 ARE=70000
N=6 X = 1212.8857 Y = 822.1488 PMT=1248.528 ARE=60000
N=7 X = 1038.5438 Y = 724.1442 PMT=1365.686 ARE=80000
3. 说明:
(1)动态显示结果:按照一定顺序(逆时针或顺时针均可)依次取样输入到计算器,可动态显示周长和面积,并可以随时停止,最后一次显示结果为最终结果。
(2)程序转换功能:当取样点为两点时,此时本程序功能变为坐标反算距离和方位角的程序,当取样点多于两点时为计算周长和面积程序
(3)纠错功能:如果当前取样点是错误点且已输入到计算器中时,可以输入该点的前一点坐标或在其前一点的位置附近重新取样再次输入到计算器内,基本不影响面积结果,精度取决于您再次放置位置与上次放置位置的偏移量,但周长需人工去除边长数据。
(4)取样点存储功能:每次取样点存储到X、Y两个序列中,对应面积显示到Freq序列中。
4. 注意事项:
(1)是否需要闭合数据:取点完成后,无论是否回归到第1点,均能正确计算结果。
(2)若不想存储可以将程序中含有→List的语句去掉(5800中有5处,9860中有8处)。