Question

matlab 关于fmincon的问题

function [c,ceq]=count2(x);
c=[-(x(5)+sqrt(40^2-(x(2)/2)^2)-sqrt(40^2-(x(4)/2)^2));(x(5)+sqrt(40^2-(x(2)/2)^2))-(x(1)/2); sqrt(((sqrt(40^2-(x(2)/2)^2)+x(5)*sqrt(-(70-x(3))^2+(x(1)/2-sqrt(40^2-(x(2)/2)^2))^2)/(x(1)/2-sqrt(40^2-(x(2)/2)^2))-sqrt(40^2-(0.5)^2*x(4)^2))^2)+(x(5)*(70-x(3))/(x(1)/2-sqrt(40^2-(x(2)/2)^2)))^2)+sqrt(40^2-(x(4)/2)^2)-(x(1)/2)];
ceg=[];
end

f=@(x)x(1)*x(2)*x(3);
A=[180 60 1 1 1 ;0 0.5 0 1 0 ];
b=[40;70];
% c=[-(x(5)+sqrt(40^2-(x(2)/2)^2)-sqrt(40^2-(x(4)/2)^2));(x(5)+sqrt(40^2-(x(2)/2)^2))-(x(1)/2); sqrt(((sqrt(40^2-(x(2)/2)^2)+x(5)*sqrt(-(70-x(3))^2+(x(1)/2-sqrt(40^2-(x(2)/2)^2))^2)/(x(1)/2-sqrt(40^2-(x(2)/2)^2))-sqrt(40^2-(0.5)^2*x(4)^2))^2)+(x(5)*(70-x(3))/(x(1)/2-sqrt(40^2-(x(2)/2)^2)))^2)+sqrt(40^2-(x(4)/2)^2)-(x(1)/2)];
% ceg=[];
x0=[180 60 1 1 1];
[x,y]=fmincon(f,x0,A,b,[],[],[],[],@count2)
x=fmincon('fun22',x0,l,a,l2,b,[],[],[],c,ceg)

为什么会出错呢 求大神
不对 最后一句没有

Answer 1

一定要使用matlab吗？我使用mathematica计算了一下，由于mathematica中x(1)+4500x(2)+3375x(3)=4570条件报错，所以我修改了一下条件，近似采用x(1)+4500x(2)+3375x(3)介于4569.9与4570.1之间来约束，结果如下：

希望对你能有所帮助。

追问

对不起哈 我看不到结果