主程序
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));
|