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

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

我想利用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次，运行后还是出现上述报错，不知道哪里出了问题，论坛里有牛人能帮忙看一下吗？感激不尽了，我刚接触这块，不是很精通！