网站公告列表 [路桥吾爱]全新改版,正在为打造一个真正属于路桥人自己并完全免费的网络空间而努力。新站启用了多线路镜像空间,让您无论使用哪种线路都能体会到不错的速度。本站继续坚持“我为人人,人人为我”的精神,鼓励大家发布并免费共享自己收藏的资料。现在加入让我们与[路桥吾爱]共同成长,有了您的加入我们会更加优秀!
 
收藏本站
会员中心
用户注册
您现在的位置: 路桥吾爱 >> 文章中心 >> 其他相关 >> 测绘 >> 文章正文
  任意三角网平差计算 [pc1500]          【字体:
任意三角网平差计算 [pc1500]

0
440 V=F-W:GOSUB 9:H=1:GOSUB 25:NEXT J:M=1:H=-1/N:GOSUB 25:Z=Z+N
445 NEXT I:M=0:GOSUB 3:R=R-3+D:IF T=0THEN 460
450 FOR J=1TO T:GOSUB 13:V=C-DEG A(J):GOSUB 9:H=1E4:Z=PEEK (A+2)+PEEK (A+3)/100:IF ZLET H=VAL M$^2/Z/Z
455 GOSUB 25:NEXT J
460 W=G-1:R=R+4*T:IF L=0THEN 485
465 FOR J=1TO L:GOSUB 13:IF I>KLET C(Q-1)=-COS C:C(Q)=-SIN C
470 S=SQR (X*X+Y*Y):C(G)=S-S(J):IF O>KLET F=2*(O-K)-1:C(F-1)=COS C:C(F)=SIN C
475 H=10:Z=PEEK (A+3)/1E3:IF ZLET H=(VAL M$/(PEEK (A+2)+S(J)*Z))^2/.1296
480 GOSUB 25:NEXT J
485 FOR H=0TO W:Z=VAL I$(H):Q=INT (Z/10):A=Z-VAL I$(H-SGN Z)+(H=0):FOR I=HTO W:Y=VAL I$(I)-I
490 X=VAL I$(I-SGN I)-Y-(I=0):IF H<=XTHEN 525
495 E=Y+H:F=INT (E/10):M=E-10*F:IF A=1THEN 515
500 FOR J=H-A+1TO H-1:IF J<=XTHEN 510
505 B=Z-H+J:C=INT (B/10):V=Y+J:S=INT (V/10):N(F,M)=N(F,M)-N(C,B-10*C)*N(S,V-10*S)
510 NEXT J
515 IF I=HLET O=SQR N(Q,Z-10*Q)
520 N(F,M)=N(F,M)/O
525 NEXT I:IF A>1FOR J=H-A+1TO H-1:B=Z-H+J:F=INT (B/10):L(H)=L(H)-N(F,B-10*F)*L(J):NEXT J
530 L(H)=L(H)/O:L(G)=L(G)-L(H)^2:NEXT H:J=K+U
535 FOR H=WTO 0STEP -1:X=-L(H):FOR I=H+1TO W:B=VAL I$(H)
540 IF I<GLET F=VAL I$(I)-I+H:IF F>VAL I$(I-1)LET E=INT (F/10):X=X-N(E,F-10*E)*L(I)
545 NEXT I:Q=INT (B/10):X=X/N(Q,B-10*Q):L(H)=X:IF INT (H/2)<H/2LET Y(J)=Y(J)+X:NEXT H
550 X(J)=X(J)+X:J=J-1:NEXT H:USING U$:M=SQR (L(G)/(D-X(0)-G+T+L))*100
555 GOSUB 4:LPRINT TAB 6;"X / Y":LPRINT :GOSUB 1:Z=1:FOR I=1TO K+U:X=X(I):Y=Y(I):J=1:IF I<=KCOLOR 3:J=0
560 LPRINT X+5E-4*SGN X*(X<1E6)*J:LPRINT Y+5E-4*SGN Y*J:GOSUB 5:IF NLF N-2
565 GOSUB 1:NEXT I:GOSUB 4:LPRINT TAB 6;"a":LPRINT :GOSUB 1:F=T:GOSUB 3:W=R-3+D
570 FOR I=1TO K+U:GOSUB 5:IF N=0LF 2:GOTO 605
575 FOR H=1TO N:R=Z:GOSUB 24:O=P:GOSUB 7
580 G=DMS (ATN (Y/(X-1E-5))+180*(X<=0)+360*(X>0AND Y<0)):IF I<=KAND O<=KCOLOR 3:GOTO 590
585 IF TGOSUB 31
590 GOSUB 16:COLOR 0:IF H=1LET E=G:NEXT H
595 C=DEG G-DEG E:IF C<0LET C=C+360
600 D(Z+H)=C:NEXT H:Z=Z+N
605 GOSUB 1:NEXT I:GOSUB 4:LPRINT TAB 8;"S":LPRINT :GOSUB 1:F=L
610 W=W+4*T:FOR I=1TO K+U:GOSUB 5:IF N=0LF 2:GOTO 630
615 FOR H=1TO N:R=Z:GOSUB 24:O=P:GOSUB 7:IF I<=KAND O<=KCOLOR 3:GOTO 625
620 IF LGOSUB 31
625 LPRINT TAB 1;SQR (X*X+Y*Y)+5E-4:COLOR 0:NEXT H:Z=Z+N
630 GOSUB 1:NEXT I:GOSUB 4:LPRINT TAB 2;D$:LPRINT :GOSUB 1:USING Y$:FOR I=1TO K+U:GOSUB 5:IF 
N=0LF 2:GOTO 645
640 LPRINT " 0 00 00.00":FOR J=2TO N:G=DMS D(Z+J):GOSUB 16:NEXT J:Z=Z+N
645 GOSUB 1:NEXT I:Z=0:GOSUB 1:LF 2:LPRINT " Ms/S":LPRINT :GOSUB 1:V=0:USING :GOSUB 3:G=R-3+D+4*T
650 FOR O=1TO K+U:IF PEEK (196+O)=0LF 2:GOTO 715
655 FOR H=1TO PEEK (196+O):R=Z:GOSUB 24:IF P<=KAND O<=KLPRINT :GOTO 710
660 IF O<POR PEEK (196+P)=0THEN 675
665 W=P:GOSUB 23:FOR I=1TO PEEK (196+P):IF PEEK (196+2*(K+U)+R+I)=OLPRINT :GOTO 710
670 NEXT I
675 IF L=0THEN 690
680 FOR J=1TO L:Q=G+4*J:A=PEEK Q:B=PEEK (Q+1):IF O=AAND P=BAND PEEK (Q+2)+PEEK (Q+3)=0
LPRINT :GOTO 710
685 NEXT J
690 X=X(P)-X(O):Y=Y(P)-Y(O):S=SQR (X*X+Y*Y):FOR I=0TO 2*U:C(I)=0:NEXT I
695 IF P>KLET A=2*(P-K)-2:C(A)=X/S:C(A+1)=Y/S
700 IF O>KLET C=2*(O-K)-2:C(C)=-X/S:C(C+1)=-Y/S:IF O<PLET A=C
705 X=0:GOSUB 43:LPRINT "1/";INT (S/M/SQR X+.5)*100
710 NEXT H:Z=Z+PEEK (196+O)
715 GOSUB 1:NEXT O:USING Y$:COLOR 2:LPRINT "Mo=";M*36+.005:H=0:B=0:FOR O=1TO U:X=0:Y=0:Z=0:LPRINT 
720 FOR V=0TO 1:A=2*O-2+V:GOSUB 43:NEXT V
725 A=ATN (2*Z/(X-Y))/2:A=A+45*(2-SGN A-SGN Z):LPRINT "P";STR$ (K+O);":";:GRAPH :ROTATE 1
730 LINE -(46,13),9:LPRINT "O":TEXT :USING "####.#":LPRINT TAB 4;"le=";A+.05;"(DEG)"
735 LPRINT " E=";TAB 2;M*SQR (X+Z*TAN A)+.05;",F=";TAB 10;M*SQR (Y-Z*TAN A)+.05;"cm":NEXT O:USING :END 
10001 GRAPH :COLOR 0:LINE -(216,0),(Z>0):TEXT :LPRINT :RETURN 

上一页  [1] [2] 

文章录入:admin    责任编辑:admin 
  • 上一篇文章:

  • 下一篇文章:
  • 最新热点 最新推荐 相关文章
    曲线任意里程中边桩坐标正反…
    曲线任意里程中边桩坐标正反…
    CASIO fx4800P曲线任意里程中…
    任意里程中边桩坐标正反算(…
    全站仪对向观测法三角高程测…
    全站仪三角高程测量在公路测…
    PC1500任意高程网平差计算
    pc1500任意三角网平差计算
      网友评论:(只显示最新10条。评论内容只代表网友观点,与本站立场无关!)
    粤ICP备06026757号
    站长:白黑点击这里和我QQ聊天-进入[路桥吾爱]新版