麻烦建模高手帮我看下下面程序哪错了

干一@者 · 发表于 2010-5-21 11:29:42

程序1 求初解
clc,clear
A1=eye(6,6);
A2=eye(4,4);
B1=A1(1:5,

;
B2=A2(1,

;
X1=[A1;A1;A1;A1;B1];
X4=[A2;A2;A2;A2;A2;A2;A2;B2];
%以下程序用于求解X2，X3中1-9个决策变量
P1=X1;t=0;
While  (t<2)
X=zeros(29,6);
c11=zeros(1,6);
c1=zeros(1,6);
count=0;
for  m=1:9
tag11=0;
for  n=1:6
if  P1(m,n)==0
if  c11(n)<2&tag11==0&count<3
X(m,n)=1;
tag11=1;
c11(n)=c11(n)+1;
c1(n)=c1(n)+1;
if  c11(n)==2
count=count+1;
end
elseif  c11(n)<1&tag11==0
X(m,n)=1;
tag11=1;
c11(n)=c11(n)+1;
c1(n)=c1(n)+1;
end
end
end
end
%以下程序用于调整X3中前9行中不满足条件的决策变量
If  t==1
a=sum(X(1:9,

′);
for  m=1:9
if  a(m)==0
b=m;
break;
end
end
Y=X(1:9,:);
for  m=1:9
kk=0;
c=P1(b,:)+Y(m,:);
for  n=1:6
if  c(n)==2
kk=1;
break
end
end
if  kk==0
X(b,:)=Y(m,:);
X(m,:)=Y(b,:);
for  n=1:6
if  P1(m,n)==0
X(m,n)=1;
break
end
end
break
end
end
end

%以下程序用于求解X2，X3中20-29个决策变量
For  m=10:29
tag12=0;
for  n=1:6
if  P1(m,n)==0
if  c1(n)<5&tag12==0
X(m,n)=1;
tag12=1;
c1(n)=c1(n)+1;
else
X(m,n)=0;
end
else
X(m,n)=0;
end
end
end
%以下程序用于调整X3中20-29个不满足条件的决策变量
if  t==1
Y=X
For  m=10:29
if  P1(m,6)==0
X(m,:)=Y(29,:);
X(29,:)=Y(m,:);
X(m,6)=1;
break
end
end
end

if  t==0
X2=X;
else
X3=X;
end
P1=P1+X;
t=t+1;
end
%以下程序用以求得满足约束条件的X5，X6，X7；
t=0;
while  (t<3)
C2=randperm1(A2);
B2= C2(1,:);
Y=[randperm1(A2); randperm1(A2); randperm1(A2); randperm1(A2);
randperm1(A2); randperm1(A2); randperm1(A2);B2];
if  t==0
X5=Y;
elseif  t==1
X6=Y;
else
X7=Y;
end
t=t+1;
end
%将初解保存在文本文件shuju中
Save  shuju X1 X2 X3 X4 X5 X6 X7

我的邮箱353079891@qq.com，懂的麻烦把结果发到我邮箱来，谢谢

		自动登录	找回密码
密码			注-册-帐-号

[问题征解] 麻烦建模高手帮我看下下面程序哪错了