机械matlab程序,机械优化设计matlab程序.doc

机械优化设计matlab程序.doc

机械优化设计作业1用二次插值法求函数极小值，精度E001。21TT在MATLAB的M文件编辑器中编写的M文件，如下FINLINE T1T22 , T A0B3EPSILON001T1AF1FT1T3BF3FT3T205T1T3F2FT2C1F3F1/T3T1C2F2F1/T2T1C1/T2T3T405T1T3C1/C2F4FT4K0WHILEABST4T2EPSILONIFT2F4F1F2T1T2T2T4F2F4ELSEF3F4T3T4ENDELSEIFF2F4F3F2T3T2T2T4F2F4ELSEF1F4T2T4ENDENDC1F3F1/T3T1C2F2F1/T2T1C1/T2T3T405T1T3C1/C2F4FT4KK1END输出最优解IFF2F4TT4FFT4ELSETT2FFT2ENDFPRINTF1, 迭代计算K30F\N ,KFPRINTF1, 极小点坐标T30F\N ,TFPRINTF1, 函数值F34F\N ,F运行结果如下迭代计算K7极小点坐标T2函数值F000012用黄金分割法求函数的极小值，精度E001。321TT在MATLAB的M文件编辑器中编写的M文件，如下FINLINE T2/3T211/3 , T A0B3EPSILON001T1B0618BAF1FT1T2A0618BAF2FT2K1WHILEABSBAEPSILONIFF1EPSONPG1\G1X0X0PG1SUBSDF,{X1,X2},{X01,1,X02,1}G1SUBSG,{X1,X2},{X01,1,X02,1}KK1ENDXX0FXSUBSF,{X1,X2},{X1,1,X2,1}运行结果如下X,FX,KNIUDUNFA11X19999554476059523381489991377897099997772380297616907449956889483FX00000000000000000039398907941382470301534502947647K232用阻尼牛顿法在MATLAB的M文件编辑器中编写的M文件，如下FUNCTIONX,FX,KZUNINIUDUNFAX0阻尼牛顿法SYMSX1X2FX124X12X22FX0VX1,X2DFJACOBIANF,VDFDF GJACOBIANDF,VEPSON1E12停机原则G1SUBSDF,{X1,X2},{X01,1,X02,1}G1SUBSG,{X1,X2},{X01,1,X02,1}K0迭代次数PG1\G1A0P G1/P G1PX0X0A0PWHILENORMA0PEPSONPG1\G1A0P G1/P G1PX0X0A0PG1SUBSDF,{X1,X2},{X01,1,X02,1}G1SUBSG,{X1,X2},{X01,1,X02,1}KK1ENDXX0FXSUBSF,{X1,X2},{X01,1,X02,1}运行结果如下X,FX,KZUNINIUDUNFA11X19999554476059523381489991377897099997772380297616907449956889483FX00000000000000000039398907941382470301534502947647K233用变尺度法在MATLAB的M文件编辑器中编写的M文件，如下4用共轭梯度法求函数的极小点122113,XXXF(1)用共轭梯度法在MATLAB的M文件编辑器中编写的M文件，如下FUNCTIONY,X,KCGA,B,C,X0共轭梯度法解MINF(X)05X AXB XCEPS1E6迭代停机原则FX05X0 AX0B X0CR0AX0BIFNORMR0EPSBETAR1 R1/R0 R0P1R1BETAP0ALPHAR1 P1/P1 AP1X1X1ALPHAP1R2AX1BP0P1R0R1R1R2KK1ENDXX1Y05X AXB XC运行结果如下Y,X,KCG3111,20,0,21Y1X1000010000K12用变尺度法在MATLAB的M文件编辑器中编写的M文件，如下FUNCTIONX,FX,KBIANCHIDUFAA,B,C,X0用变尺度法求FX05X AXB XCEPSON1E12G0AX0BG0AH0EYE2K0D0H0G0A0D0 G0/D0 G0D0S0A0D0XK1XKY0AA0D0GK1GKX1X0A0D0WHILENORMS0EPSONSWITCHKCASE{10}X0X1G0AX0BH0EYE2K0D0H0G0A0D0 G0/D0 G0D0S0A0D0X1X0A0D0BREAKOTHERWISEG1AX1BY0AA0D0S0A0D0H1H0S0S0 /S0 Y0H0Y0Y0 H0/Y0 H0Y0H1H01Y0 H0Y0/S0 Y0S0S0 H0Y0S0 S0Y0 H0/S0 Y0KK1D1H1G1A1D1 G1/D1 G0D1A0A1D0D1H0H1S0A0D0X1X1A0D0BREAKENDENDXX1FX05X1 AX1B X1C运行结果如下X,FX,KBIANCHIDUFA3111,20,0,21H104031025780257808945FX1X1000010000FX1K1故函数极小点是点(1,1)5用鲍威尔法求函数的极小点。2121214,XXXF用鲍威尔法在MATLAB的M文件编辑器中编写的M文件，如下FUNCTIONX,FX,KBOWELLA,B,C,X0鲍威尔法D0110D0201X0200ESP1E12停机原则K0迭代次数WHILENORMX0X02ESPKK1G01AX0BA01D01 G01/D01 AD01X01X0A01D01G02AX01BA02D02 G02/D02 AD02X02X01A02D02D10X02X0G10AX02BA10D10 G10/D10 AD10X10X0A01D01D01D02D02D10X0X10ENDXX0FX05X AXB XC运行结果如下X,FX,KBOWELL2224,40,0,21FX8X42FX8K36用单纯形法求线性规划问题4,321054MIN1321JXXTSXXFJ用单纯形法在MATLAB的M文件编辑器中编写的M文件，如下单纯形法MATLAB程序DANCHUNXINGFA求解标准型线性规划MAXCXSTAXBX0本函数中的A是单纯初始表，包括最后一行是初始的检验数，最后一列是资源向量BN是初始的基变量的下标输出变量SOL是最优解,其中松弛变量(或剩余变量)可能不为0输出变量VAL是最优目标值，KK是迭代次数FUNCTIONSOL,VAL,KKDANCHUNXINGFAA,NMA,NASIZEAKK0迭代次数FLAG1WHILEFLAGKKKK1IFAMA,0INBI进基变量的下标ENDENDSITAZEROS1,MA1FORI1MA1IFAI,INB0SITAIAI,NA/AI,INBENDENDTEMPINFFORI1MA1IFSITAI0122535112233440N34SOL,VAL,KKDANCHUNXINGFAA,NSOL004000016667VAL73333K