matlab源程序代码

《matlab源程序代码》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、clearall;clc;tic;X=zeros(1000,4);Y=zeros(1000,4);U=zeros(1000,1);M=3;a=[2818747470726036121814907824546298363832];b=[4250346189850404207836205266014588854];Aeq=zeros(20,80);forj=1:20fork=1:20ifj==kAeq(j,k)=1;elseAeq(j,k)=0;endendfork=21:40ifk==j+20Aeq(j,k)=1;elseAeq(j,k)=0;endendfork=41:60ifk==j+40

2、Aeq(j,k)=1;elseAeq(j,k)=0;endendfork=61:80ifk==j+60Aeq(j,k)=1;elseAeq(j,k)=0;endendendAeq;A=zeros(4,80);fori=1:4fork=1:80ifi==1fork=1:20A(i,k)=1;endendifi==2fork=21:40A(i,k)=1;endendifi==3fork=41:60A(i,k)=1;endendifi==4fork=61:80A(i,k)=1;endendendendA;B=[40;50;60;70];%B=[s1s2s3s4]zm=zeros(80,1);for

3、kkk=1:Mx=unifrnd(2,100,1,4);X(kkk,:)=x;y=unifrnd(2,100,1,4);Y(kkk,:)=y;f=zeros(80,1);fori=1:4forj=1:20d(i,j)=sqrt((x(i)-a(j))^2+(y(i)-b(j))^2);endendd=[d(1:4,1:20)];f=reshape(d,4*20,1);kk=0;N=1000;u=0;while(kk

4、opt,key,c]=linprog(f,A,B,Aeq,Beq,zm,[]);u=u+f_opt;kk=kk+1;endU(kkk)=u/N;kkk=kkk+1;endU;toc;