Matlab高手请进——请教多目标求解问题！

lionheartyin 发表于 2004-12-19 19:10:42

多目标函数求解问题，问题如下：
 
主程序M文件如下：(dqtp.m文件)
A=[0.10452.1 0.19
0.9593 1.48 0.129
2.8041 4.4546 0.6578
0.8729 1.1668 0.68
1.858 2.2566 0.0247]; 
B=[1.439 15.85 0.384
 3.806 14.82 0.267
 8.632 19.79 3.982
 4.954 5.993 2.324
 6.716 12.38 0.239]; 
C=[22.811.9 3.2
 22.8 11.9 3.2
 22.8 11.9 3.2
 22.8 11.9 3.2
 22.8 11.9 3.2];
Goals1=; 
Goals2=; 
Num1=; 
Num2=; 
Xc=; lb=; 
ub=; 
MFE=;%最大函数评价值。
GEA=;%目标达到个数。
MI=;%最大循环测算次数。
=size(A);
=size(B);
=size(C);
=size(Goals1);
for i=1:Num1
 A1(:,:,i)=;
end 
for i=1:Num1
 A2(:,:,i)=;
end
for i=1:Num1
 A3(:,:,1)=;
 A3(:,:,i+1)=A3(:,:,i)+A2(:,:,i); 
end 
A0=A3(:,:,Num1+1);
for i=1:Num1
 B1(:,:,i)=;
end 
for i=1:Num1
 B2(:,:,i)=;
end
for i=1:Num1
 B3(:,:,1)=;
 B3(:,:,i+1)=B3(:,:,i)+B2(:,:,i); 
end 
B0=B3(:,:,Num1+1);
for i=1:Num1
 C1(:,:,i)=;
end 
for i=1:Num1
 C2(:,:,i)=;
end
for i=1:Num1
 C3(:,:,1)=;
 C3(:,:,i+1)=C3(:,:,i)+C2(:,:,i); 
end 
C0=C3(:,:,Num1+1);
Aeq1=;%zeros(1,na) ones(1,ma) zeros(1,mb+mc)];
Beq1=;%zeros(1,nb+ma) ones(1,mb) zeros(1,mc)];
Ceq1=;%zeros(1,nc+ma+mb) ones(1,mc)];
ABCeq0=;
Deq0=;
=size(A0);
=size(B0);
=size(C0);
A0qk=A0(:,1:Num1*Num2);
for i=1:Num1
 A1qk(:,:,i)=A0qk(i,:);
end
for i=1:Num1
 A2qk(:,:,1)=zeros(1,Num1*Num2);
 A2qk(:,:,i+1)=A2qk(:,:,i)+A1qk(:,:,i);
end
A4=-A2qk(:,:,Num1+1);
B0qk=B0(:,1:Num1*Num2);
for i=1:Num1
 B1qk(:,:,i)=B0qk(i,:);
end
for i=1:Num1
 B2qk(:,:,1)=zeros(1,Num1*Num2);
 B2qk(:,:,i+1)=B2qk(:,:,i)+B1qk(:,:,i);
end
B4=B2qk(:,:,Num1+1);
C0qk=C0(:,1:Num1*Num2);
for i=1:Num1
 C1qk(:,:,i)=C0qk(i,:);
end
for i=1:Num1
 C2qk(:,:,1)=zeros(1,Num1*Num2);
 C2qk(:,:,i+1)=C2qk(:,:,i)+C1qk(:,:,i);
end
C4=C2qk(:,:,Num1+1);
=size(A4);
ABC=[A0 zeros(ma,ma+mb+mc);
 A4 zeros(1,ma+mb+mc);
 B4 zeros(1,ma+mb+mc);
 C4 zeros(1,ma+mb+mc);
 zeros(1,a0n) -ones(1,ma) zeros(1,mb+mc);
 zeros(1,b0n) zeros(1,ma) ones(1,mb) zeros(1,mc);
 zeros(1,c0n) zeros(1,ma+mb) ones(1,mc)
 -eye(na4+ma+mb+mc,na4+ma+mb+mc);
 eye(na4+ma+mb+mc,na4+ma+mb+mc)];
for i=1:nd
 goals1(:,:,i)=;
end
goals1(:,:,1)=-goals1(:,:,1);
for i=1:nd
 goals2(:,:,1)=zeros((md*nd),1);
 goals2(:,:,i+1)=goals2(:,:,i)+goals1(:,:,i);
end
goals0=goals2(:,:,nd+1);
Deq1=;
goal=goals0;
weight=goals0;
options=optimset('MaxFunEvals',MFE,'GoalsExactAchieve',GEA,'MaxIter',MI);
=fgoalattain(@dqtaim,Xc,goal,weight,ABC,Deq1,ABCeq0,Deq0,lb,ub,[],options)

lionheartyin 发表于 2004-12-19 19:15:02

目标函数程序如下所示：（dqtaim.m文件）
function f=myfun(x)
A=[0.10452.1 0.19
0.9593 1.48 0.129
2.8041 4.4546 0.6578
0.8729 1.1668 0.68
1.858 2.2566 0.0247]; 
B=[1.439 15.85 0.384
 3.806 14.82 0.267
 8.632 19.79 3.982
 4.954 5.993 2.324
 6.716 12.38 0.239];
C=[22.811.9 3.2
 22.8 11.9 3.2
 22.8 11.9 3.2
 22.8 11.9 3.2
 22.8 11.9 3.2];
 Num1=; 
Num2=; 
=size(A);
=size(B);
=size(C);
for i=1:Num1
 A1(:,:,i)=;
end 
for i=1:Num1
 A2(:,:,i)=;
end
for i=1:Num1
 A3(:,:,1)=;
 A3(:,:,i+1)=A3(:,:,i)+A2(:,:,i); 
end 
A0=A3(:,:,Num1+1);
for i=1:Num1
 B1(:,:,i)=;
end 
for i=1:Num1
 B2(:,:,i)=;
end
for i=1:Num1
 B3(:,:,1)=;
 B3(:,:,i+1)=B3(:,:,i)+B2(:,:,i); 
end 
B0=B3(:,:,Num1+1);
for i=1:Num1
 C1(:,:,i)=;
end 
for i=1:Num1
 C2(:,:,i)=;
end
for i=1:Num1
 C3(:,:,1)=;
 C3(:,:,i+1)=C3(:,:,i)+C2(:,:,i); 
end 
C0=C3(:,:,Num1+1);
=size(A0);
=size(B0);
=size(C0);
A0qk=A0(:,1:Num1*Num2);
for i=1:Num1
 A1qk(:,:,i)=A0qk(i,:);
end
for i=1:Num1
 A2qk(:,:,1)=zeros(1,Num1*Num2);
 A2qk(:,:,i+1)=A2qk(:,:,i)+A1qk(:,:,i);
end
A4=A2qk(:,:,Num1+1);
B0qk=B0(:,1:Num1*Num2);
for i=1:Num1
 B1qk(:,:,i)=B0qk(i,:);
end
for i=1:Num1
 B2qk(:,:,1)=zeros(1,Num1*Num2);
 B2qk(:,:,i+1)=B2qk(:,:,i)+B1qk(:,:,i);
end
B4=B2qk(:,:,Num1+1);
C0qk=C0(:,1:Num1*Num2);
for i=1:Num1
 C1qk(:,:,i)=C0qk(i,:);
end
for i=1:Num1
 C2qk(:,:,1)=zeros(1,Num1*Num2);
 C2qk(:,:,i+1)=C2qk(:,:,i)+C1qk(:,:,i);
end
C4=C2qk(:,:,Num1+1);
D=;
=size(D);
for i=1:nd
 X1(:,:,i)=;
end
for i=1:nd
 X2(:,:,1)=zeros(nd,1);
 X2(:,:,i+1)=X2(:,:,i)+X1(:,:,i);
end
X=X2(:,:,nd+1);
for i=1:md
 f(i)=D(i,:)*abs(X);end

lionheartyin 发表于 2004-12-19 20:58:27

运行结果如下所是：
increase OPTIONS.MaxFunEvals
 
x =
 
1.0e+003 *
 
Columns 1 through 7 
 
 -0.0323 0.0631 0.0518 0.0391 0.0398 0.0518 -0.0075
 
Columns 8 through 14 
 
 0.1198 -0.0371 -0.0194 0.0881 0.0601 0.0760 -0.0054
 
Columns 15 through 21 
 
 -0.0439 0.1390 0.2921 0.3686 0.1346 0.3009 0.9453
 
Columns 22 through 28 
 
 0.7514 2.1586 1.0900 0.1202 0.0555 1.4783 0.2722
 
Columns 29 through 30 
 
 0.8253 1.3574
 
 
fval =
 
1.0e+003 *
 
 1.1428 5.7530 8.5181
 
 
attainfactor =
 
1.5920e+003
 
 
exitflag =
 
0
 
问题如下：
求解的x数值严重越界，造成目标函数无法达到预期目标；我反复调整了MaxFunEvals，MaxIter参数仍然出现越界问题！请高手教我！多谢了！ 
油箱：<a href="mailtyinhang19791013@yahoo.com.cn" target="_blank" >yinhang19791013@yahoo.com.cn</A> oicq:24079022

standing 发表于 2004-12-23 20:46:11

自己摸摸索好了

页: [1]

数模论坛's Archiver

Matlab高手请进——请教多目标求解问题！