帮忙看看哪里出现问题了？帮忙调试一下，谢谢了！
程序1 求初解clc,clearA1=eye(6,6);A2=eye(4,4);B1=A1(1:5,:);B2=A2(1,);X1=;X4=;%以下程序用于求解X2，X3中1-9个决策变量P1=X1;t=0;Whilet<2X=zeros(29,6);c11=zeros(1,6);c1=zeros(1,6);count=0;for m=1:9tag11=0;for n=1:6if P1(m,n)==0if c11(n)<2&tag11==0&count<3X(m,n)=1;tag11=1;c11(n)=c11(n)+1;c1(n)=c1(n)+1;if c11(n)==2count=count+1;endelseif c11(n)<1&tag11==0X(m,n)=1;tag11=1;c11(n)=c11(n)+1;c1(n)=c1(n)+1;endendendend%以下程序用于调整X3中前9行中不满足条件的决策变量ift==1a=sum(X(1:9,:)′);for m=1:9if a(m)==0b=m;break;endendY=X(1:9,:);for m=1:9kk=0;c=P1(b,:)+Y(m,:);for n=1:6if c(n)==2kk=1;breakendendif kk==0X(b,:)=Y(m,:);X(m,:)=Y(b,:);for n=1:6if P1(m,n)==0X(m,n)=1;breakendendbreakendendend%以下程序用于求解X2，X3中20-29个决策变量for m=10:29tag12=0;for n=1:6if P1(m,n)==0if c1(n)<5&tag12==0X(m,n)=1;tag12=1;c1(n)=c1(n)+1;elseX(m,n)=0;endelseX(m,n)=0;endendend%以下程序用于调整X3中20-29个不满足条件的决策变量if t==1Y=Xfor m=10:29if P1(m,6)==0X(m,:)=Y(29,:);X(29,:)=Y(m,:);X(m,6)=1;breakendendendif t==0X2=X;elseX3=X;endP1=P1+X;t=t+1;end%以下程序用以求得满足约束条件的X5，X6，X7；t=0;while t<3C2=randperm1(A2);B2= C2(1,:);Y=randperm1(A2); randperm1(A2); randperm1(A2);B2];if t==0X5=Y;elseif t==1X6=Y;elseX7=Y;endt=t+1;end%将初解保存在文本文件shuju中save shuju X1 X2 X3 X4 X5 X6 X7程序2 目标函数function f=F(A,B,C,D,E,F,G)Q1=A*A′-eye(29,29);Q2=B*B′-eye(29,29);Q3=C*C′-eye(29,29);Q4=D*D′-eye(29,29);Q5=E*E′-eye(29,29);Q6=F*F′-eye(29,29);Q7=G*G′-eye(29,29);Q_sum=Q1+Q2+Q3+Q4+Q5+Q6+Q7;Q_sum=Q_sum+eye(29,29);M=Q_sum1==zeros(29,29);count=sum(sum(M));g=norm(Q_sum);f=0.7*count+0.3*g;程序3 随即变换function=randperm1(A)n=size(A,1);y=randperm(n);for p=1:nX(p,:)=A(y(p),:);end程序4 退火过程%模拟退火算法求最优解tic;clc;clearload shuju;t0=100;tf=0.01k=1000;f=F(X1,X2,X3,X4,X5,X6,X7);t=t0;while t>tffor k=1:ky=randperm(9);for m=1:9X11(m,:)=X1(y(m),:);X22(m,:)=X2(y(m),:);X33(m,:)=X3(y(m),:);endz=randperm(20);z=z+9*ones(1,20);for m=10:29X11(m,:)=X1(z(m-9),:);X22(m,:)=X2(z(m-9),:);X33(m,:)=X3(z(m-9),:);endp=randperm(6);for m=1:6X11(:,m)=X1(:,p(m));X22(:,m)=X2(:,p(m));X33(:,m)=X3(:,p(m));endA2=eye(4,4);C4=randperm1(A2);B4=C4(1,:);C5=randperm1(A2);B5=C5(1,:);C6=randperm1(A2);B6=C6(1,:);C7=randperm1(A2);B7=C7(1,:);X44=;X55=;X66=;X77=;ff=F(X11,X22,X33,X44,X55,X66,X77);if ff<fX1=X11;X2=X22;X3=X33;X4=X44;X5=X55;X6=X66;X7=X77;elseif exp((f-ff)/t)>rand(1);X1=X11;X2=X22;X3=X33;X4=X44;X5=X55;X6=X66;X7=X77;endf=ffendt=0.87*t;endtoc;1=tocsave zhi X1 X2 X3 X4 X5 X6 X7
%将结果保留在zhi中程序5
计算见面次数load zhiQ1=X1*X1′-eye(29,29);Q2=X2*X2′-eye(29,29);Q3=X3*X3′-eye(29,29);Q4=X4*X4′-eye(29,29);Q5=X5*X5′-eye(29,29);Q6=X6*X6′-eye(29,29);Q7=X7*X7′-eye(29,29);Q_sum=Q1+Q2+Q3+Q4+Q5+Q6+Q7;Q_sum1=Q_sum+eye(29,29);M=Q_sum1==zeros(29,29);count=sum(sum(M));g=norm(Q_sum);f=0.7*count+0.3*g;f,count,gcount0=0;count1=0;count2=0;count3=0;count4=0;count5=0;count6=0;count7=0;for m=1:29for n=1:29if Q_sum(m,n)==0count0=count0+1;elseif Q_sum(m,n)==1count1=count1+1;elseif Q_sum(m,n)==2count2=count2+1;elseif Q_sum(m,n)==3count3=count3+1;elseif Q_sum(m,n)==4count4=count4+1;elseif Q_sum(m,n)==5count5=count5+1;elseif Q_sum(m,n)==6count6=count6+1;elseif Q_sum(m,n)==0count7=count7+1;endendend(count0-29)/2,count1/2,count2/2,count3/2,count4/2,count5/2,count6/2,count7/2 randperm1，运行时，说没有这个函数呢？帮忙调试一下，万分感谢！！！ randperm(n)作用是产生一个介于1~n的随机序列，比如randperm(6),可能得到序列。
但randperm1()这个函数在Matlab是没有的，你自己如果没有定义的话就肯定不能执行。你可以把randperm1改成randperm试一试，看结果有没有问题。 是生成n个在内的无重复的随机整数吗？是randperm，而不是randperm1... 回复 3# olh2008
randperm1()这个函数在Matlab是没有的，你自己如果没有定义的话就肯定不能执行。……那要怎么定义呢？ 回复 4# 数学者
%以下程序用以求得满足约束条件的X5，X6，X7；t=0;while t<3C2=randperm1(A2);B2= C2(1,:);Y=randperm1(A2); randperm1(A2); randperm1(A2);B2];if t==0X5=Y;elseif t==1X6=Y;elseX7=Y;endt=t+1;end 这一段有问题吗？因为这是别人的程序，不懂啊，帮忙看一下！ 回复 6# nashihuakaimm
把randperm1改成randperm试试看~
页:
[1]