搜狗做题网matlab 函数 fmincon求解非线性规划,
来源：学生作业帮 编辑：搜狗做题网作业帮 分类：数学作业 时间：2024/04/28 15:50:19
matlab 函数 fmincon求解非线性规划,
求这个函数的最小值 function f = myfun(x)
f = (pi*(x(3)+2)*x(2)*x(1)^2)/4;
运行后返回（约束在最后面有写）
Warning:Trust-region-reflective method does not currently solve this type of problem,
using active-set (line search) instead.
> In fmincon at 439
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 300 (the default value).
x =
4.0499
14.0013
1.0013
FVAL =
541.3289
说是这个trust……解不了这个问题,那么该怎么改一下?help过没发现怎么改……
非线性约束 function [c,ceq] = mycon(x)
k=1.6*(x(1)/x(2))^0.14; %曲度系数
F=13; %最大力
f=20.2416; %气门最大升程
G=8100; %切变模量
w=1500; %凸轮轴转速
c(1) =8*k*F*x(2)/(pi*x(1)^3)-51; %切应力
c(2) =((x(3)+2-0.5)+1.2*f)/x(2)-5.3 ; %稳定性约束
c(3) =((G*x(1)^4*1.2*f/(8*x(2)^3*x(3)))*8*k*x(2))/(pi*x(1)^3); %极限切应力?
c(4) =10*w-3.56*10^5*x(1)/(x(2)^2*x(3)) ; %无共振
ceq = [];
主函数 [x,FVAL] = fmincon(@myfun,[2.5; 21; 7],[4 -1 0;-18 1 0],[0;0],[],[],[2.2;16;3],[3.5;28;Inf],@mycon(x))
求这个函数的最小值 function f = myfun(x)
f = (pi*(x(3)+2)*x(2)*x(1)^2)/4;
运行后返回（约束在最后面有写）
Warning:Trust-region-reflective method does not currently solve this type of problem,
using active-set (line search) instead.
> In fmincon at 439
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 300 (the default value).
x =
4.0499
14.0013
1.0013
FVAL =
541.3289
说是这个trust……解不了这个问题,那么该怎么改一下?help过没发现怎么改……
非线性约束 function [c,ceq] = mycon(x)
k=1.6*(x(1)/x(2))^0.14; %曲度系数
F=13; %最大力
f=20.2416; %气门最大升程
G=8100; %切变模量
w=1500; %凸轮轴转速
c(1) =8*k*F*x(2)/(pi*x(1)^3)-51; %切应力
c(2) =((x(3)+2-0.5)+1.2*f)/x(2)-5.3 ; %稳定性约束
c(3) =((G*x(1)^4*1.2*f/(8*x(2)^3*x(3)))*8*k*x(2))/(pi*x(1)^3); %极限切应力?
c(4) =10*w-3.56*10^5*x(1)/(x(2)^2*x(3)) ; %无共振
ceq = [];
主函数 [x,FVAL] = fmincon(@myfun,[2.5; 21; 7],[4 -1 0;-18 1 0],[0;0],[],[],[2.2;16;3],[3.5;28;Inf],@mycon(x))
主要问题：
1、关于Trust-region-reflective算法的警告
该类方法不适用于存在非线性约束的情况,所以用警告的方式提示你使用别的方法.如果不介意看到警告信息,直接忽略无妨；如果不想看到警告,主程序可以这样调用：
opt=optimset('Algorithm','active-set',);
x0=[2.5 21 7];
[x,FVAL] = fmincon(@myfun,x0,[4 -1 0;-18 1 0],[0;0],[],[],[2.2;16;3],[3.5;28;Inf],@mycon,opt)
2、程序现在最大的问题是,你给出的初始解不在可行域内,也就是说,不满足你的所有约束条件（很容易验证,不满足非线性不等式约束c(x)
再问： 多谢~~~ 呃……我发现是约束函数写错了…… 还想问一下opt=optimset('Algorithm', 'active-set',);这个里面的两个参数是怎么来的 See OPTIMSET for details. For a list of options accepted by FMINCON refer to the documentation. 我看有上面这句话，不过这个documentation在哪里找？
再答： 输入命令： doc optim/OPTIMSET 满意的话希望采纳。
再问： 嗯，多谢了~~我这个还得提升，方便留一下你的QQ或者邮箱么？以后有问题还想问你请教~我的QQ412757622.顺便问一下你对谢菲尔德遗传工具箱熟悉么？
再答： 我的邮箱是Q522319871，不常在线，如有问题探讨请发邮件，勿加好友。 官方的遗传算法工具箱我用过（从R2010a版开始，由原来的Genetic Algorithm and Direct Search Toolbox经扩充后改名为Global Optimization Toolbox），其它几个第三方遗传算法工具箱没用过。
1、关于Trust-region-reflective算法的警告
该类方法不适用于存在非线性约束的情况,所以用警告的方式提示你使用别的方法.如果不介意看到警告信息,直接忽略无妨；如果不想看到警告,主程序可以这样调用：
opt=optimset('Algorithm','active-set',);
x0=[2.5 21 7];
[x,FVAL] = fmincon(@myfun,x0,[4 -1 0;-18 1 0],[0;0],[],[],[2.2;16;3],[3.5;28;Inf],@mycon,opt)
2、程序现在最大的问题是,你给出的初始解不在可行域内,也就是说,不满足你的所有约束条件（很容易验证,不满足非线性不等式约束c(x)
再问： 多谢~~~ 呃……我发现是约束函数写错了…… 还想问一下opt=optimset('Algorithm', 'active-set',);这个里面的两个参数是怎么来的 See OPTIMSET for details. For a list of options accepted by FMINCON refer to the documentation. 我看有上面这句话，不过这个documentation在哪里找？
再答： 输入命令： doc optim/OPTIMSET 满意的话希望采纳。
再问： 嗯，多谢了~~我这个还得提升，方便留一下你的QQ或者邮箱么？以后有问题还想问你请教~我的QQ412757622.顺便问一下你对谢菲尔德遗传工具箱熟悉么？
再答： 我的邮箱是Q522319871，不常在线，如有问题探讨请发邮件，勿加好友。 官方的遗传算法工具箱我用过（从R2010a版开始，由原来的Genetic Algorithm and Direct Search Toolbox经扩充后改名为Global Optimization Toolbox），其它几个第三方遗传算法工具箱没用过。