资源描述：
《机械优化设计,惩罚函数法,多变量,约束优化》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、C==============PROGRAMCOMPLEC==============DIMENSIONX(25),GX(50),XCOM(1250)COMMON/ONE/ITE,IXE,ILI,NPE,NFX,NGRREAD(*,*)N,KG,KWRITE(*,10001)N,KG,K10001FORMAT(25X,'==========PRIMARYDATA=========='//5X,1'N=',I4,5X,'KG=',I4,5X,'K=',I4)CALLMAISUB(N,K,KG,X,GX,XCOM)STOPENDC================================
2、===SUBROUTINEMAISUB(N,K,KG,X,GX,XCOM)C===================================DIMENSIONX(N),GX(KG),XCOM(N,K),FXK(50),XR(25)DIMENSIONXO(25),XH(25),XL(25),BL(25),BU(25)COMMON/ONE/ITE,IXE,ILI,NPE,NFX,NGRCOMMON/TWO/ISEREAD(*,*)(X(I),I=1,N)READ(*,*)EPSREAD(*,*)KWR,ISEREAD(*,*)(BL(I),I=1,N),(BU(I),I=1,N)WRIT
3、E(*,1010)(BL(I),I=1,N)WRITE(*,1015)(BU(I),I=1,N)WRITE(*,1020)EPS1010FORMAT('BL:'/(5X,5E15.6))1015FORMAT('BU:'/(5X,5E15.6))1020FORMAT(5X,'EPS=',E10.2)ITE=0NFX=0IXE=0RM=2657863.01025CALLPRICOM(N,K,KG,X,GX,XCOM,FXK,BL,BU,RM)IF(KWR.LT.0)GOTO1041WRITE(*,1030)1030FORMAT(/25X,'==========PRIMARYCOMPLEX===
4、======='/)WRITE(*,1080)KDO1031L=1,KWRITE(*,1085)L,(XCOM(I,L),I=1,N)1031CONTINUEWRITE(*,1035)(FXK(I),I=1,K)1035FORMAT(4X,'FXK:'/(5X,5E15.6))1041WRITE(*,1042)1042FORMAT(/25X,'==========ITEATIONCOMPUTE=========='/)1045ITE=ITE+1CALLFXSEGU(N,K,XCOM,FXK)DO1050I=1,N1050XL(I)=XCOM(I,K)FXL=FXK(K)SDX=0.0DO1
5、055I=1,K-11055SDX=SDX+(FXL-FXK(I))**2SDX=SQRT(SDX/FLOAT(K-1))IF(SDX.LE.EPS)GOTO1210IF(KWR.GT.0)GOTO1056IF(ITE/10*10.NE.ITE)GOTO10901056WRITE(*,1060)ITE,FXLIF(KWR.LT.0)GOTO1090WRITE(*,1065)(XL(I),I=1,N)WRITE(*,1070)FXLWRITE(*,1075)(GX(I),I=1,KG)WRITE(*,1035)(FXK(I),I=1,K)1060FORMAT(/1X,'*****ITE=',
6、I4,5X,'FXL=',E15.7)1065FORMAT('X:'/(5X,5E15.6))1070FORMAT('FX:'/(5X,5E15.6))1075FORMAT('GX:'/(5X,5E15.6))1080FORMAT('XCOM:(K=',I3,')')1085FORMAT(2X,I2/(5X,5E15.6))1090LH=11095DO1100I=1,N1100XH(I)=XCOM(I,LH)FXH=FXK(LH)CALLXCENTE(N,K,K,LH,XO,XCOM)CALLFFX(N,XO,FXO)CALLGGX(N,KG,XO,GX)DO1105J=1,KGIF(GX
7、(J).GE.0.0)GOTO11701105CONTINUE1140PHI=1.31145DO1150I=1,N1150XR(I)=XO(I)+PHI*(XO(I)-XH(I))CALLFFX(N,XR,FXR)CALLGGX(N,KG,XR,GX)DO1151J=1,KGIF(GX(J).GE.0.0)GOTO11521151CONTINUEGOTO11551152PHI=0.5*PHIGOTO11451155IF(
