[MATLAB工具箱] 利用函数fmincon求解非线性方程时遇到计算提前终止的问题

[alafazhen]

我想利用fmincon求解非线性方程，a1=(0.09-(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019).*x+x.*x./4);
a2=(0.5*exp(-x).*(1-exp(x/2))+(0.25+0.019)*exp(-x).*(1-exp(-x))-0.019*x+0.09);
a3=log(sqrt(8*(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019)./x+4+0.37./x./x)+2*(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019)./x+2);
a4=log(2*(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019)./x+sqrt(3)./x);
待求解方程：
sqrt(a1./a2).*(a3-a4)-2.98=0;
程序如下：
主函数main:
clc
clear
close all
%带非线性约束的最小值@fun是优化的目标函数 见fun.m ；@mycon是非线性约束见mycon.m
options = optimoptions('fmincon','MaxIter',10000);
x = fmincon(@fun,[1],[],[],[],[],[],[],@mycon);
ymin=fun(x);
x
ymin
函数文件fun:
function y=fun(x)

%约束条件 就是原公式的最小值
a1=(0.09-(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019).*x+x.*x./4);
a2=(0.5*exp(-x).*(1-exp(x/2))+(0.25+0.019)*exp(-x).*(1-exp(-x))-0.019*x+0.09);
a3=log(sqrt(8*(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019)./x+4+0.37./x./x)+2*(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019)./x+2);
a4=log(2*(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019)./x+sqrt(3)./x);


y=(sqrt(a1./a2).*(a3-a4)-2.98);
约束条件文件mycon:
function [c,ceq] = mycon(x)

%非线性约束，要保证根号内的大于0，ln内的也大于0 ，否则就会求出复数
c(1)=(0.09-(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019).*x+x.*x./4)./(0.5*exp(-x).*(1-exp(x/2))+(0.25+0.019)*exp(-x).*(1-exp(-x))-0.019*x+0.09); %不等式约束1
c(2)=8*(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019)./x+4+0.37./x./x; %不等式约束2
c(3)=sqrt(8*(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019)./x+4+0.37./x./x)+2*(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019)./x+2;
c(4)=(2*(0.25*exp(-x)-(0.25+0.019)*exp(-2*x)+0.019)./x+sqrt(3)./x);
ceq = [];%不存在等式约束

运行程序后，出现Solver stopped prematurely.

fmincon stopped because it exceeded the iteration limit,
options.MaxIter = 1000 (the default value).
我把迭代次数改为1次，运行后还是出现上述报错，不知道哪里出了问题，论坛里有牛人能帮忙看一下吗？感激不尽了，我刚接触这块，不是很精通！