如何用matlab解决这样的非线性优化问题！ - MATLAB爱好者论坛-LabFans.com

yfchenth2 · 2008-07-26, 10:25

主程序
clc
clear
p0=[160.06 64.94 147.5];
a=[];
b=[];
aeq=[];
beq=[];
xl=[60 64 64];
xu=[180 180 180];
i=1;
global h
for h=150:50:3500
[p,fmin]=fmincon(@obcon,p0,a,b,aeq,beq,xl,xu,@constraint);
c(i,1:3)=p;
p0=p;
i=i+1;
end
c

子程序
function [g u]=constraint(x)
global h
g=[];
u=762.5*cos(x(1)+x(2)+x(3))-1260*cos(x(1))-1182.5*cos(x(1)+x(2))+511.5-h+1260*sin(x(1))-1182.5*sin(x(1)+x(2))+762.5*sin(x(1)+x(2)+x(3))+1400;
end

function y=obcon(x)
y=-36962.04*cos(x(1))+15868.3*cos(x(1)+x(2))-3385*cos(x(1)+x(2)+x(3));

20052110 · 2008-07-28, 10:37

先谢谢这位哥们了，很高兴你能回帖
不过，为什么运行的时候前面会出现很多警告，而且用两个非线性等式约束为什么要加在一起，我看书上好像是用分号隔开的！

相似的主题
主题	主题作者	版面	回复	最后发表
Inside MATLAB Objects in R2008a	TechnicalArticles	MATLAB技术文章	0	2008-10-22 02:59
人脸识别分类的问题	feiyan2000000	MATLAB论坛	0	2008-06-06 11:43
【求助】新手请教，谢谢	zhugeping	MATLAB论坛	0	2007-07-15 16:01
【求助】走遍了所有的Matlab社区都没人能够解答的一个小问题	TD——SCDMA	MATLAB论坛	2	2007-05-19 12:56

2008-07-26, 10:25	#1
yfchenth2 初级会员注册日期: 2008-07-03 年龄: 48 帖子: 16 声望力: 18	主程序 clc clear p0=[160.06 64.94 147.5]; a=[]; b=[]; aeq=[]; beq=[]; xl=[60 64 64]; xu=[180 180 180]; i=1; global h for h=150:50:3500 [p,fmin]=fmincon(@obcon,p0,a,b,aeq,beq,xl,xu,@constraint); c(i,1:3)=p; p0=p; i=i+1; end c 子程序 function [g u]=constraint(x) global h g=[]; u=762.5cos(x(1)+x(2)+x(3))-1260cos(x(1))-1182.5cos(x(1)+x(2))+511.5-h+1260sin(x(1))-1182.5sin(x(1)+x(2))+762.5sin(x(1)+x(2)+x(3))+1400; end function y=obcon(x) y=-36962.04cos(x(1))+15868.3cos(x(1)+x(2))-3385*cos(x(1)+x(2)+x(3));

2008-07-28, 10:37	#2
20052110 初级会员注册日期: 2008-07-24 年龄: 39 帖子: 2 声望力: 0	先谢谢这位哥们了，很高兴你能回帖不过，为什么运行的时候前面会出现很多警告，而且用两个非线性等式约束为什么要加在一起，我看书上好像是用分号隔开的！