机械优化设计MATLAB程序

机械优化设计作业1.用二次插值法求函数极小值，精度e=0.01。在MATLAB的M文件编辑器中编写的M文件，如下：f=内嵌(1)*(t-2)2,t)a=0 .b=3；=0.01；t1=a .f1=f(t1)；T3=b。F3=f(T3)；T2=0.5 *(t1 T3)；F2=f(T2)；C1=(F3-f1)/(T3-t1)；C2=(F2-f1)/(T2-t1)-C1)/(T2-T3)；T4=0.5 *(t1 T3-C1/C2)；F4=f(T4)；k=0；而(澳大利亚统计局(T4-T2)=)如果t2f4f1=f2t1=t2t2=t4f2=f4其他f3=f4t3=t4目标其他if f2f4f3=f2t3=t2t2=t4f2=f4其他f1=f4t2=t4目标目标C1=(F3-f1)/(T3-t1)；C2=(F2-f1)/(T2-t1)-C1)/(T2-T3)；T4=0.5 *(t1 T3-C1/C2)；F4=f(T4)；k=k1目标%输出最优解if f2f4t=t4f=f(T4)；其他t=t2f=f(T2)；目标fprintf(1，迭代计算k=%3.0fn，k)fprintf(1，极小点坐标t=%3.0fn，t)fprintf(1，函数值f=%3.4fn，f)运行结果如下：迭代计算k=7极小点坐标t=2函数值f=0.00012.用黄金分割法求函数的极小值，精度e=0.01。在MATLAB的M文件编辑器中编写的M文件，如下：f=inline(t(2/3)-(t2 1)(1/3),t)；a=0 .b=3；=0.01；t1=b-0.618 *(b-a)；f1=f(t1)；T2=a 0.618 *(b-a)；F2=f(T2)；k=1；而防抱死制动系统(b-a)=if f1epson)G1 G1；x0=x0 p；g1=subs(df，x1，x2，x0(1，1)，x0(2，1)；G1=subs(G，x1，x2，x0(1，1)，x0(2，1)；k=k1目标x=x0fx=subs(f，x1，x2，x(1，1)，x(2，1)；运行结果如下：x，fx，k=纽敦法(1；1)x=1.99995544760595233814899913778970.999977238029761690749956889483外汇=0.0000000000000000039398907941382470301534502947647k=23(2)用阻尼牛顿法在MATLAB的M文件编辑器中编写的M文件，如下：函数x，fx，k=zuniudunf(x0)%阻尼牛顿法syms x1 x2f=(x1-2)4(x1-2*x2)2；FX=0；v=x1，x2；df=雅可比(f，v)；df=df .g=雅可比(df，v)；爱普生=1e-12；%停机原则g1=subs(df，x1，x2，x0(1，1)，x0(2，1)；G1=subs(G，x1，x2，x0(1，1)，x0(2，1)；k=0；%迭代次数G1 G1；A0=-p * G1/(p * G1 * p)；x0=x0 A0 * p；同时（标准(a0*p)爱普生)G1 G1；A0=-p * G1/(p * G1 * p)；x0=x0 A0 * p；g1=subs(df，x1，x2，x0(1，1)，x0(2，1)；G1=subs(G，x1，x2，x0(1，1)，x0(2，1)；k=k1目标x=x0fx=subs(f，x1，x2，x0(1，1)，x0(2，1)；运行结果如下：x，fx，k=祖纽邓发(1；1)x=1.99995544760595233814899913778970.999977238029761690749956889483外汇=0.0000000000000000039398907941382470301534502947647k=23(3)用变尺度法在MATLAB的M文件编辑器中编写的M文件，如下：4.用共轭梯度法求函数的极小点(1)用共轭梯度法在MATLAB的M文件编辑器中编写的M文件，如下：函数y，x，k=重心(A，b，c，x0)%共轭梯度法解最小值(x)=0.5*X*A*X bx cEPS=1e-6；%迭代停机原则%fx=0.5*x0 .* a . * x0 b . * x0 c；r0=A * x0 b；if norm(r0)=epsx=x0y=0.5 * x * A * b * c；k=0；目标P0=-r0；=-r0 * P0/(P0 * P0)；x1=x0 a * p0R1=A * x1 b；k=0；而norm(r1)eps=(R1 * R1)/(r0 * r0)；P1=-R1* P0；=-(R1 * P1)/(P1 * P1)；x1=x1R2=A * x1 b；p0=p1r0=r1 .r1=r2k=k1目标x=x1y=0.5 * x * A * b * c；运行结果如下：y，x，k=CG(3-1；-11)，-2；0，0，2；1)y=-1x=1.00001.0000k=1(2)用变尺度法在MATLAB的M文件编辑器中编写的M文件，如下：函数x，fx，k=bianchidufa(A，b，c，x0)%用变尺度法求FX=0.5 * x * A * b * c；爱普生=1e-12；G0=A * x0 b；G0=A；H0=眼睛(2)；k=0；d0=-H0 * G0；A0=-d0 * G0/(d0 * G0 * d0)；s0=A0 * d0 % x(k1)-x(k)；y0=A * A0 * d0 % g(k1)-g(k)；x1=x0 a0 * d0而(norm(s0)=epson)开关k案例10x0=x1G0=A * x0 b；H0=眼睛(2)；k=0；d0=-H0 * G0；A0=-d0 * G0/(d0 * G0 * d0)；s0=a0 * d0x1=x0 a0 * d0破裂否则G1=A * x1 b；y0=A * a0 * d0s0=a0 * d0% H1=H0 s0 * s0/(s0 * y0)-H0 * y0 * y0 * H0/(y0 * H0 * y0)；H1=H0(1 y0 * H0 * y0/(s0 * y0)* s0 * s0-H0 * y0 * s0-s0 * y0 * H0)/(s0 * y0)；k=k1D1=-H1 * G1；a1=-D1 * G1/(D1 * G0 * D1)；a0=a1d0=d1 .H0=H1；s0=a0 * d0x1=x1 a0 * d0破裂目标目标x=x1FX=0.5 * x1 * A * x1 b * x1 c；运行结果如下：x，fx，k=比安奇都法(3-1；-11)，-2；0，0，2；1)H1=0.40310.25780.25780.8945fx=-1x=1.00001.0000fx=-1k=1因此，函数的最小点是点(1，1)5.用鲍威尔法求函数的最小点。用鲍威尔法在MATLAB的M文件编辑器中编写的M文件如下：函数x，fx，k=鲍威尔(a，b，c，x0)%鲍威尔方法d01=1；0；d02=0；1；x02=0；0；esp=1e-12；%关机原则k=0；%迭代而norm(x0-x02)=espk=k1；g01=A * x0 b；a01=-d01 * g01/(d01 * A * d01)；x01=x0 a01 * d01g02=A * x01b；A02=-d02 * g02/(d02 * A * d02)；x02=x01 a02 * d02D10=x02-x0；G10=A * x02b；a10=-D10 * G10/(D10 * D10)；x10=x0 a01 * d01d01=d02d02=D10；x0=x10目标x=x0FX=0.5 * x * A * b * c；运行结果如下：x，fx，k=鲍尔(2-2；-24，-4；0，0，2；1)fx=-8x=42fx=-8k=36.用单纯形法求解线性规划问题在MATLAB的M文件编辑器中，用单纯形法编写的M文件如下：%单纯形法matlab程序%求解标准线性规划： max c * x；标准时间A * x=b；x=0这个函数中的% A是一个简单的初始表，包括：最后一行是初始检查号，最后一列是资源向量b% N是初始基变量的