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

L*J+G*Z+Y*(P-F)
320 IF C>2LET X(C)=X:Y(C)=Y:GOTO 345-65*(Q=K)
325 X=X-X(1):Y=Y-Y(1):E=X(2)-X(1):F=Y(2)-Y(1)
330 W=ATN (Y/X)+90*(2-SGN X-SGN X*SGN Y)-ATN (F/(E-1E-9))-180*(E<=0)-360*(F<0AND E>0)
335 FOR J=1TO I-(C=2):F=K+J:X=X(F)-X(1):Y=Y(F)-Y(1):X(F)=X*COS W+Y*SIN W+X(1)
340 Y(F)=Y*COS W-X*SIN W+Y(1):NEXT J:Q=0
345 NEXT I:Z=0
375 FOR I=1TO S:C=VAL S$(I):D=VAL T$(I):IF C>KAND H(D-K)<D-CLET H(D-K)=D-C
380 NEXT I:IF T=0THEN 430
385 FOR I=1TO T:V=VAL D$(I):FOR J=1TO V:C=VAL P$(Z+J):Q=C-K:IF Q<1THEN 425
390 IF I<=KTHEN 405
395 IF C<IAND H(I-K)<I-CLET H(I-K)=I-C
400 IF C>IAND H(Q)<C-ILET H(Q)=C-I
405 FOR H=JTO V:D=VAL P$(Z+H):IF D<=KTHEN 420
410 IF D<CAND H(Q)<C-DLET H(Q)=C-D
415 IF D>CAND H(D-K)<D-CLET H(D-K)=D-C
420 NEXT H
425 NEXT J:Z=Z+V:NEXT I
430 W=0:H(0)=-1:H(U+1)=U:FOR I=1TO U+1:H(I)=H(I-1)+4*H(I)+3:NEXT I:V=0:IF X$="Y"LET F=1
435 P=180/PI :DIM S(1,S),N(INT ((H(U)+2*U+1)/10),9):U=2*U:FOR I=1TO S:GOSUB 6
440 S(U)=H-L(I):G=-COS E:H=SIN E:GOSUB 16:IF X$="Y"THEN 460
445 IF W$="N"GOSUB 29:GOTO 455
450 A=X(0):B=Y(0)
455 F=(L(0)/(A+L(I)*B/1E3+.003))^2
460 GOSUB 8:NEXT I:IF T=0THEN 485
465 FOR I=1TO T:C=VAL N$(I):L=1:V=VAL D$(I):GOSUB 19:Z=0
470 FOR J=1TO V:D=VAL P$(W+J):GOSUB 7:R=E-D(W+J):IF R<0LET R=R+360
475 IF LLET Z=Z+R/V:NEXT J:L=0:GOTO 470
480 S(U)=R-Z:GOSUB 15:GOSUB 8:NEXT J:F=-F/V:Q=1:W=W+V:GOSUB 8:NEXT I:V=0:W=0
485 IF OFOR I=1TO O:C=VAL #$(6+I):D=VAL #$(12+I):GOSUB 7:M=VAL #$(I)+.001:GOSUB 15:S(U)=E-A(I)
490 IF OGOSUB 19:GOSUB 8:NEXT I
500 FOR I=1TO U/2:S(2*I-2)=(H(I)+H(I-1)-1)/2:S(2*I-1)=H(I):NEXT I:B=H(I)-U-2:S(U)=B
505 A=INT (B/10):V=N(A,B-10*A):FOR H=0TO U-1:A=S(H)-S(H-SGN S(H))+(H=0)
510 FOR I=HTO U:Y=S(I)-I:X=S(I-SGN I)-Y-(I=0):IF H<=XTHEN 545
515 E=Y+H:F=INT (E/10):M=E-10*F:W=N(F,M):IF A=1THEN 535
520 FOR J=H-A+1TO H-1:IF J<=XTHEN 530
525 B=S(H)-H+J:C=INT (B/10):Q=Y+J:L=INT (Q/10):W=W-N(C,B-10*C)*N(L,Q-10*L)
530 NEXT J
535 IF I=HLET R=SQR W
540 W=W/R:N(F,M)=W:IF I=ULET V=V-W*W
545 NEXT I:NEXT H:Z=K+U/2:X$="Adjusted "
550 FOR I=U-1TO 0STEP -1:A=S(U-1)+I+1:B=INT (A/10):X=-N(B,A-10*B):IF I=U-1THEN 565
555 FOR J=I+1TO U-1:F=S(J)-J+I:IF F>S(J-1)LET E=INT (F/10):X=X-N(E,F-10*E)*S(1,J)
560 NEXT J
565 Q=INT (S(I)/10):X=X/N(Q,S(I)-10*Q):S(1,I)=X:IF INT (I/2)<I/2LET Y(Z)=Y(Z)+X:NEXT I
570 X(Z)=X(Z)+X:Z=Z-1:NEXT I:M=SQR (V/(S-U+N-T+O))*100:P=5E-4
575 USING U$:GOSUB 1:LPRINT :LPRINT X$;"X / Y":GOSUB 1:FOR I=1TO K+U/2:X=X(I):Y=Y(I):J=(I>K):IF J=0COLOR 3
580 LPRINT "P";STR$ I;TAB 4;X+P*SGN X*(X<1E6)*J:LPRINT TAB 4;Y+P*SGN Y*J:GOSUB 1:NEXT I
600 GOSUB 1:LPRINT " b ";X$:LPRINT TAB 9;"S":GOSUB 1:FOR I=1TO S:GOSUB 6:IF W$="N"USING :GOSUB 29:LPRINT B;
610 LPRINT TAB 2;H+P:NEXT I:GOSUB 1:GOSUB 1:LPRINT TAB 10;X$;TAB 13;"a":GOSUB 1:FOR I=1TO S:GOSUB 6:IF O=0THEN 630
620 FOR J=1TO O:IF C=VAL #$(6+J)AND D=VAL #$(12+J)COLOR 1:GOTO 630
625 NEXT J
630 G=DMS E:GOSUB 2:COLOR 0:NEXT I:GOSUB 1:IF T=0THEN 700
635 GOSUB 1:LPRINT TAB 10;X$;TAB 8;S$:GOSUB 1:W=0:FOR I=1TO T:LPRINT :C=VAL N$(I):V=VAL D$(I)
640 FOR J=1TO V:D=VAL P$(W+J):GOSUB 7:IF J=1LET O=E
645 E=E-O:IF E<0LET E=E+360
650 G=DMS E:GOSUB 2:NEXT J:GOSUB 1:W=W+V:NEXT I
700 GOSUB 1:LPRINT :LPRINT " Ms/S":GOSUB 1:H=0:COLOR 2
710 FOR O=1TO S:G=VAL S$(O):P=VAL T$(O):X=X(P)-X(G):Y=Y(P)-Y(G):Z=SQR (X*X+Y*Y)
715 FOR I=1TO U:S(I)=0:NEXT I:IF G>KLET A=2*(G-K)-2:S(A)=X/Z:S(A+1)=Y/Z
720 C=2*(P-K)-2:S(C)=-X/Z:S(C+1)=-Y/Z:IF G<=KLET A=C
725 X=0:GOSUB 20:LPRINT "1/";STR$ (INT (L(O)/M/SQR X+.5)*100):NEXT O:GOSUB 1:USING "####.#":COLOR 2:G=0
730 LPRINT "Mo=";M*10+.05;"mm":LPRINT :B=0:FOR O=1TO U/2:X=0:Y=0:Z=0:FOR H=0TO 1:A=2*O+H-2:GOSUB 20
735 NEXT H:A=ATN (2*Z/(X-Y))/2:A=A+45*(2-SGN A-SGN Z):LPRINT "P";STR$ (K+O);":";:GRAPH :ROTATE 1
740 LINE -(46,13),9:LPRINT "O":TEXT :LPRINT TAB 4;"le=";A+.05;"(DEG)"
745 LPRINT " E=";TAB 2;M*SQR (X+Z*TAN A)+.05;",F=";TAB 9;M*SQR (Y-Z*TAN A)+.05;"cm":LPRINT :NEXT O:USING :END 
10001 GRAPH :COLOR 0:LINE -(216,0):TEXT :LPRINT :RETURN 
5;",F=";TAB 9;M*SQR (Y-Z*TAN A)+.05;"cm":LPRINT :NEXT O:USING :END 
10001 GRAPH :COLOR 0:LINE -(216,0):TEXT :LPR

上一页  [1] [2] 

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

  • 下一篇文章:
  • 最新热点 最新推荐 相关文章
    匝道PCE500计算源程序
    任意三角网平差计算 [pc1500…
    PC1500单一附合导线平差计算
    PC1500任意高程网平差计算
    pc1500任意三角网平差计算
      网友评论:(只显示最新10条。评论内容只代表网友观点,与本站立场无关!)
    粤ICP备06026757号
    站长:白黑点击这里和我QQ聊天-进入[路桥吾爱]新版