VB常用算法源代码.doc

《VB常用算法源代码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、1.高斯函数GAUSSJSubGAUSSJ(A(),N,B())DimIPIV(50),INDXR(50),INDXC(50)ForJ=1ToNIPIV(J)=0NextJForI=1ToNBIG=0#ForJ=1ToNIfIPIV(J)<>1ThenForK=1ToNIfIPIV(K)=0ThenIfAbs(A(J,K))>=BIGThenBIG=Abs(A(J,K))IROW=JICOL=KEndIfElseIfIPIV(K)>1ThenPrint"Singularmatrix"EndIfNextK

2、EndIfNextJIPIV(ICOL)=IPIV(ICOL)+1IfIROW<>ICOLThenForL=1ToNDUM=A(IROW,L)A(IROW,L)=A(ICOL,L)A(ICOL,L)=DUMNextLDUM=B(IROW)B(IROW)=B(ICOL)B(ICOL)=DUMEndIfINDXR(I)=IROWINDXC(I)=ICOLIfA(ICOL,ICOL)=0#ThenPrint"Singularmatrix."PIVINV=1#/A(ICOL,ICOL)A(ICOL,ICOL)=

3、1#ForL=1ToNA(ICOL,L)=A(ICOL,L)*PIVINVNextLB(ICOL)=B(ICOL)*PIVINVForLL=1ToNIfLL<>ICOLThenDUM=A(LL,ICOL)A(LL,ICOL)=0#ForL=1ToNA(LL,L)=A(LL,L)-A(ICOL,L)*DUMNextLB(LL)=B(LL)-B(ICOL)*DUMEndIfNextLLNextIForL=NTo1Step-1IfINDXR(L)<>INDXC(L)ThenForK=1ToNDUM=A(K,I

4、NDXR(L))A(K,INDXR(L))=A(K,INDXC(L))A(K,INDXC(L))=DUMNextKEndIfNextLEndSub2.松弛迭代SubSSOR(A(),N,B(),X(),EPS,OM,II)IMAX=200ForI=1ToNR=1#/A(I,I)B(I)=B(I)*RForJ=1ToNA(I,J)=A(I,J)*RNextJNextIForII=1ToIMAXRX=0#ForI=1ToNR=B(I)ForJ=1ToNR=R-A(I,J)*X(J)NextJIfAbs(R)

5、>RXThenRX=Abs(R)X(I)=X(I)+OM*RNextIIfOM*RX<=EPSThenExitSubNextIIPrint"Toomanyiterations"EndSub3.SVDCMPSubSVDCMP(A(),M,N,W(),V())DimRV1(100)IfM

6、#IfI<=MThenForK=IToMSCALE1=SCALE1+Abs(A(K,I))NextKIfSCALE1<>0#ThenForK=IToMA(K,I)=A(K,I)/SCALE1S=S+A(K,I)*A(K,I)NextKF=A(I,I)G=-Sqr(S)*Sgn(F)H=F*G-SA(I,I)=F-GIfI<>NThenForJ=LToNS=0#ForK=IToMS=S+A(K,I)*A(K,J)NextKF=S/HForK=IToMA(K,J)=A(K,J)+F*A(K,I)NextKN

7、extJEndIfForK=IToMA(K,I)=SCALE1*A(K,I)NextKEndIfEndIfW(I)=SCALE1*GG=0#S=0#SCALE1=0#IfI<=MAndI<>NThenForK=LToNSCALE1=SCALE1+Abs(A(I,K))NextKIfSCALE1<>0#ThenForK=LToNA(I,K)=A(I,K)/SCALE1S=S+A(I,K)*A(I,K)NextKF=A(I,L)G=-Sqr(S)*Sgn(F)H=F*G-SA(I,L)=F-GForK=LT

8、oNRV1(K)=A(I,K)/HNextKIfI<>MThenForJ=LToMS=0#ForK=LToNS=S+A(J,K)*A(I,K)NextKForK=LToNA(J,K)=A(J,K)+S*RV1(K)NextKNextJEndIfForK=LToNA(I,K)=SCALE1*A(I,K)NextKEndIfEndIfIfANORM>Abs(W(I))+Abs(RV1(I))ThenANORM=ANORMElseANORM=Ab