如何用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));

主题工具
显示可打印版本邮寄本页给好友
显示模式
切换到平板模式切换到混合模式树形模式

相似的主题
主题	主题作者	版面	回复	最后发表
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	#2
yfchenth2 初级会员注册日期: 2008-07-03 年龄: 47 帖子: 16 声望力: 17	主程序 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));