求MATLAB程序,关于有约束优化设计的 主程序:clccleart s 分别是x1 x2A=[1-1;1 0];b=[-0.002;0.01];x0=[0.5;0.3];Aeq=[];beq=[];lb=[];ub=[];[x,fval]=fmincon(@fstfun,x0,A,b,Aeq,beq,lb,ub,@fstnonlcon)两个子程序:function f=fstfun(x)t s分别是x1 x2f=pi/4*(x(2).^2-x(1).^2);endfunction[c,ceq]=fstnonlcon(x)t s分别是x1 x2c=205*pi*(x(1).^4-x(2).^4)+0.016*x(2);ceq=[];end
matlab 非线性约束最优化问题,希望给出代码 先建立M文件fun.m定义目标函数:function f=fun(x)f=((x(1)+1)^2+4*(x(2)-1.5)^2)*((x(1)-1.2)^2+0.4*(x(2)-0.5)^2);再建立M文件mycon.m定义非线性约束:function[g,ceq]=mycon(x)g=[2*x(1)-x(1)*x(2)+5*x(2)-6;x(1)-x(2)+0.5;x(1)^2-4*x(2)^2+x(2)];ceq=0;主程序为:x0=[0;0];VLB=[0 0];VUB=[];[x,fval,exitflag,output]=fmincon('fun',x0,[],[],[],[],VLB,VUB,'mycon2')结果为:x=0.59251.0925fval=1.6306exitflag=1output=iterations:5funcCount:18lssteplength:1stepsize:8.6939e-007algorithm:'medium-scale:SQP,Quasi-Newton,line-search'firstorderopt:3.2838e-008constrviolation:1.6077e-011message:[1x144 char]
有约束最优化问题,用matlab求解 假设最优时候的a1不等于a2,那么取a1'=a2'=max{a1,a2}将是更优的解。因此,最优时候的a1与a2必定相等。给定角加速度a时,加速时间越长那么转过的角度越多。在加速度不大于0.5g的约束下,加速时间最多可以是:加速与减速过程所转过的角度是a*t(a)^2,是个随a递减的函数。假设最优时候的角加速度为a,加速时间t(a),那么可以增大a到某个值a',加速时间为t(a'),使得a'*t(a')^2=at^2。因此,最优时候的加速时间必取到最大值。综上,可得最终优化式子:代码如下:g=9.8;r=.056;t=(a)(g^2/(4*r^2*a^4)-1/a^2)^(1/4);f=(a)t(a)+22.2/a/t(a);a=fminsearch(f,1e-6);fprintf('a1=a2=f\\nt1=t3=f\\nt2=f\\n',a,t(a),22.2/a/t(a)-t(a))