第11章程序设计345 IF(N2,GE. LDGOTO 60 N2=NI-1 DO 80 I=LlN H=Xl(+1)-X1(D) AI(1,4)=Y1()!计算插值区间内插值节点处的x1,K2,K3,K (1,3)=A()H (1,2)=3.·(Y1(1+1)Y1())(2*A(I)H+A(+1)*H I(,1)=2·(Yl(I}Yl(+1)+(A(+A(+1)H 100 CONTINUE 0o200l=1,NSl-1 LI=KQ2) NI=KQ2(1+1) STEP=1/NI-Ll) Sx=X1(H+1)Xl(1) DO 200 J=LLNI 2(=X1()+S1·SX!计算滑裂面与土条侧面边界线交点的X,y值 Y2(J)=AI(,4)+S1*(AI(,3)HS1°(AI(,2H+SlAI(,1) 1=S1+STEP NW=KQ2(NSI WRITE(6, N=,NW WRITE(6, 105 M, X2(M), Y2(M), M=I, NW) 105 FORMAT(1X,1s,2F10.3) STOP END 3.光滑曲线滑裂面例题 (1)数据文件 Spline-ldat: 1-51522-35.000 6-46.314-28.235 计算成果 Spline-1 151 2-50480-33.466 -47.356-29.317 6-46.314-28.235 7-45.15727.280 8-44.000-26.545 9-42.842-25.975 10-41.685-25.511 l1-40.528-25.10 13-38.213-24.495 16-34.718-24.243 17-33.539-24.352 19-31.180-24.754 030.000-25.000第 11 章 程序设计 345 N2=N2-1 IF(N2.GE.L1)GOTO 60 N2=N1-1 DO 80 I=L1,N2 H=X1(I+1)-X1(I) AI(I,4)=Y1(I) ! 计算插值区间内插值节点处的κ1,κ2,κ3,κ4 AI(I,3)=A(I)*H AI(I,2)=3.*(Y1(I+1)-Y1(I))-(2.*A(I)+A(I+1))*H 80 AI(I,1)=2.*(Y1(I)-Y1(I+1))+(A(I)+A(I+1))*H 100 CONTINUE DO 200 I=1,NS1-1 L1=KQ2(I) N1=KQ2(I+1) STEP=1./(N1-L1) SX=X1(I+1)-X1(I) S1=0. DO 200 J=L1,N1 X2(J)=X1(I)+S1*SX ! 计算滑裂面与土条侧面边界线交点的 X Y 值 Y2(J)=AI(I,4)+S1*(AI(I,3)+S1*(AI(I,2)+S1*AI(I,1))) 200 S1=S1+STEP NW=KQ2(NS1) WRITE(6,*)'N= ',NW WRITE(6,105)(M,X2(M),Y2(M),M=1,NW) 105 FORMAT(1X,I5,2F10.3) STOP END 3. 光滑曲线滑裂面例题 (1) 数据文件 Spline-1.dat 5 1 -51.522 -35.000 6 -46.314 -28.235 10 -41.685 -25.511 15 -35.898 -24.224 20 -30.000 -25.000 0 计算成果 Spline-1.c N= 20 1 -51.522 -35.000 2 -50.480 -33.466 3 -49.439 -31.978 4 -48.397 -30.580 5 -47.356 -29.317 6 -46.314 -28.235 7 -45.157 -27.280 8 -44.000 -26.545 9 -42.842 -25.975 10 -41.685 -25.511 11 -40.528 -25.107 12 -39.370 -24.765 13 -38.213 -24.495 14 -37.055 -24.311 15 -35.898 -24.224 16 -34.718 -24.243 17 -33.539 -24.352 18 -32.359 -24.530 19 -31.180 -24.754 20 -30.000 -25.000