PARO,K)要2R3 Q=AI(I.K)+BICK) S=(ARCK)+AKLK))(BR(K)+B(K) 眼D= ()-P+Q B0)=B)-S+P+Q 0 CONTINUE 100 CONTINUE D-AK(N.N):AR(N.N)+AICN.N)+AION.N) n F(W41.a0.1.)TEN L-0 WRTTE(+.20> RETURN END IH P=ARN,N)m级CN) Q--AT(N.N)BION) S-(AR(N .N- AIN.N).(BRON)+BION) BR(N)-(P-Q)/D BI(N)-(S-P-Q)/D D02001-N-11,-1 D0150J=+1,N P=AR1,J》+R Q=AIC.J)+BI) +AJ(IJ))(R ()-B]d)S+P+Q 150 CONTINUE 200 CONTINUE SON)N D0110K=N,1,1 BR(K) BRC)R(JS(K) BRUS(K))=P P=BI(K) BK)BIOS(K) (S))-P 120 NTINUE RETURN END 六、例 求解下列复系数方程组 10
AX-8 其中A=AR+iAI,B=BR+iB 132 f3-216 R 21 9153-2 A1= -2758 9-3151 L-2-2115 -2-276d BR=2,1,3,9,B=(1,2,-,3 主程序（文件名：ACGAS0,FOR)为 DIMENSION AR(4.4),AI(4.4).BR(4).BI(4).JS(4) DATA AR/1.0,.D,-20,0,20,15,0,-么0 20.1.0.3.0,11-0.13.0,-20,-20.5.0/ DATA A1/30.-2.,9.0,-20,-2.0,.0,-30,-20, 1.05,0,15.0,7.0,6.0,80,1.0,i.c/ DATA BR/.20,7,0.a0,9.0/ DATA81.0,20.-a.0,1./ CALL ACGAS(AR.AL.4.BR.BI.LJS) LNE)HN WRITE(10)(BR(D)I14) RTE(.,10(O,l-1,4) END IF 100RMAT(1X,4DI4.6) 运行结果为 .67R233D-01-.182341D+o0 ,59852D+00.245455D+00 .707823D-01-.75194D+00 -.437131D+00 .11395D+60 1.4复系数方程组的全选主元高斯-约当消去法 一、功 用全选主元高所约当(Gauss-】orn)消去法求鳏系数矩阵相同顶具有多组右鞘常 数向量的复系数线性代数方程组AX=B。 二、方法说明 基本方齿同1.？节。 为了计算两个复数的秦积 e+f=a+)+0 本子程序来用如下算法： p-ac:q=bd,s=(a+6)(c+d) =p-gf-s-p-g 11
为了计算两个复数的 e+if-(c+id)/(a+b) 本子程序采用如下算法： P-ac.q--bd.-(a-b)(c+d)ww=d+B e=(p-q)fwf=(s-p-q)/w 三，子程序语句 SUBROUTLNE ACJDN(AR.AI.N.BR.Bl.M.L.JS> 四、形参说明 AR,AI一均为双精度实型二维教组，体积为NXN,输入参数。分别存放方程组系 数矩阵的实部与虚部，返回时将被坏。 N一整型变量，输入参数。方程组阶数。 BR,B一均为双精度实型二堆数组，体积为NXM,输入兼行出参数，润月时分别 存放M组右装数向的实部与虚都：返回时分则存放M组解向量的实部与崖部。 M一整型变量，静入参数。方程组右给常数向量的组数。 L一整型变量，输出鑫数。若返回L=0,说明方程组的系数矩阵奇异，求解失败，若 L≠0，表示正常返向。 S一型一维数组，长度为N。本子程序的工作数组。 五、芋程序（文件名，ACDN,FOR) SUBROUTINE ACIDN(AR,AI,N,RR,图，M,L.IS) DIMENSTON AR(N.N).A1(N.N).BR(N.M),BIN.M).JS(N) DOUBLE FRECISION AR.A .Bl.D.P.Q.S D-0.0 DO 101-K.N DO 10J-K.N P-AR(L.J)ARO0J》+M0D￥AIJ) IF OP.GT.D)THEN KK)-I END F 10 CONTINUE W=D IF (W+1.0.EQ.1.0)THEN WRITE(20) 1=0 RETURN FORMAT() 12
DO 3I-K.N P=ARK, RKJ AR(S.J) R(S.J》 AI(K.I)AI(S.] AI(IS,D=P 30 CONTINUE D0 35J-1.M PR(K BRCK.D-BRUIS.D RQS.J)-F P-BIK,J S》一上 35 CONTINUE D 50I-1.N P=ARQ,K) AR1,K)=A.SK》 AROIS(K=P P=AICK) AI(I.K)-AI(I.JS(K)) K-P 50 DO60J-K-IN P=AR(K,J】ARK,K Q-A(KJ)AI(K,K) S-(AR(K,K)-AI(K,K))+(AR(K.J)+AICK.J)) AR(K.J)-(PQ)/ A0K.J)=5-P-Q)石 CONTINUE D065J=1M Q=-BK,J·0K,K) S=(ARK,K》一A1(K,K)e(BR(X,J》风K,I) BR(K=(P-0)/D B团(K,J》=sP-Q/D CONTINUE D090I=1,N 压L.NE.)THE 080J=K+1. P-ARKL.K)+AR(K.D)
Q-AIO,K)+AIKJ) S-(AR(L.K)+AIC.K))(AR(K.D)+AI(K.J)) ARIJ》=RL,n-P+Q AldD=AIO.J)-S+P+Q CONTINUE 0851=1,M Q-AIC S-(AR(,K)+ IO.K))*(BR(K.J)+BICK.J> RA.J》=RI,D-P+Q BI(1,])=队0，D-s+P+Q CONTINUE END IF CONTINUE Do110K=N.i.-1 D01oJ P-RR(K.J) )0), 欧K,-灯 P-B[区，JD BI(KJ)-BI(JS(K).J) BIUSK)D=P 110 CONTINUE RETURN 六.到 求解复系数方程组X=日，其中 A=AR十jAI,B=BR+BI 2 3 -21可 AR= 21-2 -2758 915a-2 ,A= 9-3151 -2-15J L-2-27 211 「123 72 27 BR- 3-2 BI- -2 3 93 L3 主程序（文件名：ACJDNO..FOR)为 D0 MENSION AR(4,4),AI(,4),B眼(4,2》BI4,2》,S(4) DOUBLE PRECISION AR.AI.BR.BI DATA AR/1,0.7.0,1.0,-2.0,3.0,2.015.01-2.0 1