主脚本写:fun = @(X)X(1)^3+X(2)^4;fmincon(fun,[0 0],[],[],[],[],[],[],@con) 然后写一个con.m(别嫌麻烦,非线性约束必须这么写,没办法)function [c,ceq]=con(X)c = -1;ceq = X(1)^2+X(2)^2-1;