matlab求多元非线性，求和里面的m也是变量吗

Question

目标函数：
function y=myfun(k)
y=x(1)*x(2)
约束条件：

function [c,ceq]=cons(x)
eq=symsum((400^2+800*(m-1)*y+(m-1)^2*y^2)/(800*(m-1)*y+(m-1)^2*y^2),m,2,x);
c=-double(eq)-6;
ceq=[];

主程序：

>> syms x y m
ceq=[];
A=[];
b=[];
Aeq=[];
beq=[];
x(1)=6,x(2)=8;
k0=[x(1),x(2)];
lb=[6,8];
ub=[12,16];
[x,fval,exitflag,output]=fmincon(@fun,k0,A,b,[],[],[],[],@cons);
disp('x(1)、x(2)、y值')
 
x =
 
6
 
错误使用 fmincon (line 224)
FMINCON requires the following inputs to be of data type double: 'X0'.
 
![img](https://img-mid.csdnimg.cn/release/static/image/mid/ask/19219581353612.png "=600 #left")

Answer 1

同学，你把你的优化的目标函数和约束条件都说一下，目前来看你的这些代码比较凌乱，看不出具体在优化什么，以及约束什么，还有sym变量，非常忌讳的东西。
全局找最小，不用优化：

% 这个是全局找最小
[X,Y] = meshgrid(6:12,8:16);
F = X.*Y;
C = true(size(X));
for i = 1:numel(X)
    n = X(i);
    y = Y(i);
    c = 0;
    for k = 1:1:n
        c = c+ (184900+860*(k-1)*y^2+(k-1)^2*y^2)/...
            (860*(k-1)*y+(k-1)^2*y^2+249600);
    end
    C(i) = 24-c<=0;
end
minf = min(F(C));
I = find(minf==F&C);
x = X(I)
y = Y(I)
fmin = x*y

最后结果：

x =

    12


y =

    10


fmin =

   120

优化的代码在评论区，只是fmincon只能找到局部最优

Answer 2

老师您好，我现在需要对这个进行优化，麻烦老师看一下，因为我是机械专业在做课题，小白入门，麻烦最好写下代码

Answer 3

你好，根据fmincon得到的matlab函数如下：

function main
x0 = [12;16];
A=[];
b=[];
Aeq = [];
beq = [];
lb = [6;8];
ub = [12; 16];
count = 1;
fvalmin = 1e10;
while(count<100)
    x0(1) = randi([6,12],1);
    x0(2) = randi([8,16],1);
    [x,fval,exitflag,output]=fmincon(@myfun,x0,A,b,Aeq,beq,lb,ub,@cons);
    if(fvalmin>fval && exitflag>0)
        xmin = round(x);
        fvalmin =fval;
    end
    count = count + 1;
end
xmin
fvalmin
end

function y = myfun(x)
y = floor(x(1))*floor(x(2));
end

function [c,ceq]=cons(x)
n = floor(x(1));
y = floor(x(2));
c = 0;
for k = 1:1:n
    c = c+ (184900+860*(k-1)*y^2+(k-1)^2*y^2)/...
        (860*(k-1)*y+(k-1)^2*y^2+249600);
end
c = 24-c;
ceq = [];
end

但是这个fmincon并不能得到最小值，只能得到局部域里面的最小，只有循环若干次之后才有全局最小（其实也不怎么算）：

xmin =

    11
    11


fvalmin =

   121