N/,单位长度允许挠度/ =1/400.Р(3)求解方法和结果用内点惩罚函数法来求解这个问题.其惩罚函数按式Р Р初始点取,是一个可行点.取惩罚因子初始值,降低系数C=0.7,收敛精度.其计算结果如下:Р Р 3Р1.0000Р30.0000Р157.5689Р 2.1Р0.6519Р33.3545Р127.1697Р 1.47Р0.6460Р32.1997Р122.1886Р 1.029Р0.6391Р30.2894Р115.0050Р Р Р 0.0415Р0.6349Р26.1520Р103.5664Р 0.17x10^-11Р0.6366Р24.9685Р 101.3706Р0.06827x10^-12Р0.6366Р 24.9685Р101.3706Р 最优解为=0.6366 =24.9685 Р(4) 内点法的流程图:РMATLAB程序:Р利用文件编辑器编写目标函数myfun.m函数文件Рfunction f=myfun(x)Рf=120*x(1)+x(2)Р编写约束条件mycon.m函数文件Рfunction[c,ceq]=mycon(x)Рc(1)=1-(7/45)*x(1)*x(2);Рc(2)=1-(7/45)*x(1)^3*x(2);Рc(3)=1-(1/321)*x(1)*x(2)^2;Рceq=[ ]Р在窗口中调用优化程序РA=[-1 0Р 0 -1Р 0 -0.25];Рb=[0;0;-1];Рx0=[1;30];РIb=zeros(2,1);Р[x,fval,exitflag,output]=fmincon(@myfun,x0,A,b,[],[],Ib,[],@mycon)Р得到的结果为:Рx =Р 0.6332Р 25.3264Рfval =Р 101.3056Р 经过笔算与matlab程序结果相比较.笔算有一定的误差,可见笔算精确度不高.
全文预览
