搜搜做题作业网请问你可以帮我编一下下面的函数吗,
来源：学生作业帮 编辑：搜搜做题作业网作业帮 分类：综合作业 时间：2024/05/08 02:34:09
请问你可以帮我编一下下面的函数吗,
function y=optm(Xij)
options = optimset('Algorithm','interior-point');
X=Xij;
LB=zeros(1,31);
UB=960000*ones(1,31);
X0=10000*rand(1,31);
[y,FVAL,EXITFLAG] = FMINCON(@(x)obfun(x,X),X0,[],[],[],[],LB,UB,@(x)con(x,X),options)
function GG=obfun(x,X)
W=[0.1440 0.1380 0.1440 0.1440 0.1440 0.1420 0.1440];
Sx=sum(x);
for i=1:31
Y(i)=x(i)/Sx;
end
for j=1:7
SG=0;
for i=2:31
SG=SG+(X(j,i)-X(j,i-1))*(Y(i)+Y(i-1));
end
G(j)=1-SG;
end
GG=sum(W.*G);
function [c,ceq]=con(x,X)
Gu=[0.3160 0.5960 0.5310 0.3060 0.4850 0.4550 0.4020];
Sx=sum(x);
for i=1:31
Y(i)=x(i)/Sx;
end
for j=1:7
SG=0;
for i=2:31
SG=SG+(X(j,i)-X(j,i-1))*(Y(i)+Y(i-1));
end
G(j)=1-SG;
end
for i=1:7
c(i)=G(i)-Gu(i);
end
for i=8:14
c(i)=-G(i-7);
end
c(15)=-Sx+120000;
c(16)=Sx-960000;
ceq=[];
options = optimset('Algorithm','interior-point');
X=Xij;
LB=zeros(1,31);
UB=960000*ones(1,31);
X0=10000*rand(1,31);
[y,FVAL,EXITFLAG] = FMINCON(@(x)obfun(x,X),X0,[],[],[],[],LB,UB,@(x)con(x,X),options)
function GG=obfun(x,X)
W=[0.1440 0.1380 0.1440 0.1440 0.1440 0.1420 0.1440];
Sx=sum(x);
for i=1:31
Y(i)=x(i)/Sx;
end
for j=1:7
SG=0;
for i=2:31
SG=SG+(X(j,i)-X(j,i-1))*(Y(i)+Y(i-1));
end
G(j)=1-SG;
end
GG=sum(W.*G);
function [c,ceq]=con(x,X)
Gu=[0.3160 0.5960 0.5310 0.3060 0.4850 0.4550 0.4020];
Sx=sum(x);
for i=1:31
Y(i)=x(i)/Sx;
end
for j=1:7
SG=0;
for i=2:31
SG=SG+(X(j,i)-X(j,i-1))*(Y(i)+Y(i-1));
end
G(j)=1-SG;
end
for i=1:7
c(i)=G(i)-Gu(i);
end
for i=8:14
c(i)=-G(i-7);
end
c(15)=-Sx+120000;
c(16)=Sx-960000;
ceq=[];
美国人的姓氏和名字有哪些.如果可以帮我编一下名字
帮我把下面的中文,翻译成英文.你可以改变一下位置和写法等.（60字左右）
帮我给下面的作文分类取上三四个小标题,要有诗意,好听点.重编一下目录.
可以帮我翻译一下下面的英语句子吗?谢谢!
英语好的帮我翻译一下下面这句话可以吗,谢谢了啊
英语翻译请帮我翻译一下下面这段文字,请问你现在有时间吗?我想跟你学英语,请先不要拒决我,因为我需要你的帮助,每天只需要花
我的名字是钰婕你可以帮我找一下他的英文名吗?
你可以帮我普及一下奢侈品知识吗?
可以帮我画一下它的大致函数图像吗?
请问一下可以帮我估下2010年12月的四级分数吗?
请问,你可以教我一下‘雨中情’的作文.
英语翻译可不可以帮我翻一下下面的广东话?尹家嘅咁多位,鉴于我手机唔见咗,喺我有新手机之前都唔可以join你哋嘅group