最近我们小组做了建模训练，选择的题目是21年的国赛C题，我负责软件部分，这题我刚好会做，第一次独自完成了代码，有点激动，和大家分享一下，共同学习，不过因为我个人比较懒，讲解的不怎么多。。。不过会分享一下完整的代码

其实我的代码还有点问题，转运方案没有按照它的要求，当时没注意到这个问题，所以输出的结果和题目给的表格不太一样。

我路子可能比较野，写的代码有点杂乱无章，感觉算小屎堆了（悲）

因为事后没有整理，所以代码会有一些多余的部分，是当时用过后发，之后又不用的。。。

我也是个建模新手，还有很多不足的地方，欢迎大家指出我的错误，一起交流

大家会看到一个Excel，叫‘问题1评分排名.xlsx’，当时鬼迷心窍弄成这样，现在一看实属脱裤子放屁。这个表格大家复制就能用了

第一问

第一问选取了

几个指标，利用熵权topsis法进行排名，取前五十名

不过数据处理的代码里还有一些计算其他的，写了但是没用上。。。

熵权topsis代码非原创，是借鉴了某位大佬的，不过我忘记是哪位的了。。。。

如果大佬不满意的话请联系我，我把它删掉

第一问代码

%%第一问代码
%处理数据
clc;clear
[num,txt,raw]=xlsread('附件1 近5年402家供应商的相关数据.xlsx',1);
[num2,txt2,raw2]=xlsread('附件1 近5年402家供应商的相关数据.xlsx',2);
%计算供给达标周数
sum_hege=zeros(402,1);
[m,n]=size(num);
for j=1:m
    for i=1:n
        if num(j,i)~=0&&num2(j,i)>=num(j,i)
            sum_hege(j,1)=sum_hege(j,1)+1;%计算供给达标的周数
        end
    end
end
num_gongji=sum(num~=0,2);%计算企业要求供应商供给的周数
num2_sum=sum(num2,2);%求出供应商一共提供了多少材料；
num2_mean=mean(num2,2);%得出平均每个星期提供的材料数量；
num2_gongji=sum(num2~=0,2);%计算供应商有供货的周数；
num2_wugongji=sum(num2==0,2);%计算没有供货的周数；
num2_month=num2_mean*4;
dabiaolv=sum_hege./num_gongji;
%将材料量化
cailiao=ones(402,1);
xunzhaoA = strcmp(txt2(2:end,2),'A');%找出A，标记为1
mA=find(xunzhaoA==1);%找出字符串所在的索引
xunzhaoB = strcmp(txt2(2:end,2),'B');%找出A，标记为1
mB=find(xunzhaoB==1);%找出字符串所在的索引
cailiao(mA)=1.2;
cailiao(mB)=1.1;
data=[cailiao num2_sum  num2_gongji dabiaolv];
bingtu=[sum(num2_sum<100),sum(num2_sum>=100&num2_sum<1000),sum(num2_sum>=1000&num2_sum<10000),sum(num2_sum>=10000)];
pie3(bingtu)
legend('供货量小于100','供货量100-1000','供货量1000-10000','供货量大于10000')
%%
%熵权topsis
x=data;%导入数据
[n,m]=size(x);
% 数据的归一化处理
[X,ps]=mapminmax(x',0,1);
 ps.ymin=0.002; % 归一化后的最小值
ps.ymax=0.996; % 归一化后的最大值
ps.yrange=ps.ymax-ps.ymin; % 归一化后的极差
X=mapminmax(x',ps);
X=X';  % X为归一化后的数据
 % 计算比重
 for i=1:n
     for j=1:m
         p(i,j)=X(i,j)/sum(X(:,j));
     end
 end
 %计算熵值
 k=1/log(n);
 for j=1:m
     e(j)=-k*sum(p(:,j).*log(p(:,j)));
 end
 d=ones(1,m)-e;  % 计算信息熵冗余度
w=d./sum(d);    % 求权值w
%TOPSIS综合评价
A=data;%评价矩阵
W=w%权重
 [ma,na]=size(A);
 for i=1:na
     B(:,i)=A(:,i)*W(i); %得到加权标准化矩阵
end
V1=zeros(1,na);            %初始化正理想解和负理想解
V2=zeros(1,na);
BMAX=max(B);               %取加权标准化矩阵每列的最大值和最小值
BMIN=min(B);               
for i=1:na
     V1(i)=BMAX(i);
     V2(i)=BMIN(i);
end
for i=1:ma   
     C1=B(i,:)-V1;
     S1(i)=norm(C1);   %S1,S2分别为离正理想点和负理想点的距离
     C2=B(i,:)-V2;
     S2(i)=norm(C2);
     T(i)=S2(i)/(S1(i)+S2(i));
end
output=T'%得出评分
xlswrite('评分.xlsx',output)

第二问

为了省事，我直接选择只从以上五十家购买材料，我发现，好像，应该，可能只购买材料A或者只购买材料C，最便宜，不过我当时在储存费用方面考虑欠佳。然后限制条件吧啦吧啦，懒得码字了，直接贴代码了。。。大家有不明白的，或者发现我问题的，欢迎在评论区和我交流

第二问代码

%若出现错误，可能是内存问题，重新运行即可
clc;clear;
[num,txt,raw]=xlsread('问题1评分排名.xlsx');
data=xlsread('附件1 近5年402家供应商的相关数据.xlsx',2);
index=1:402;
num2=[index' num];
num3=sortrows(num2,7);%对排名进行升序
wushi=num3(1:50,:);
indexC=find(wushi(:,2)==1);
C=wushi(indexC,:);
for i=1:20
Max(i)=max(data(C(i,1),:));
end
T=28200*0.72;
changku=0;%表示出库材料存量
buy=zeros(20,24);
A=0;
f=100000000;%初始化
X=1000;%模拟次数
count=1;
for Xunhuan=1:X
    
for i=1:20
buy(i,1)=min(6000,randi(round([Max(i)*0.3,Max(i)*0.6])));
end
if sum(buy(:,1))>=T
    changku=sum(buy(:,1))-T;
    Changku(1)=changku;
for j=2:24%星期
    for i=1:20%供应商
        fanwei=min(6000,round([Max(i)*0.3,Max(i)*0.6]));
        if sum(buy(:,j))<=T||changku<=2*T
            buy(i,j)=randi(fanwei);
            if sum(buy(:,j))>=T
                changku=changku+sum(buy(:,j))-T;
                if changku>=2*T
                    continue
                end
            end
        end
    end
    if sum(buy(:,j))<=T
        changku=changku-(T-sum(buy(:,j)));
    end
    Changku(j)=changku;%记录每周仓库的存储数量
end
end
%数据合并
Buy=zeros(402,24)*nan;
Buy(C(:,1),:)=buy;
%%
%模拟转运损耗
into=xlsread('附件2 近5年8家转运商的相关数据.xlsx','B2:IG9');
yunshu=zeros(5,24)
for j=1:8%第j个转运商
hang=into(j,:);
%sunhao=unique(hang);
%[m,n]=size(sunhao);
for i=1:24
xunzhe=randi(240)
yunshu(j,i)=hang(xunzhe)
end
end
%%
%实施转运
p=1:8;;
a=1;%第a列
transport=zeros(20,24);
for i=1:24%24周
    a=1;%第a行
    R=0;%运输量
yunshu2=[p' yunshu(:,i)];
yunshu2(find(yunshu2(:,2)==0),:)=[];%剔除不运输的转运商
[m,n]=size(yunshu2);
yunshu3=sortrows(yunshu2,2)%进行降序
for j=1:20%20家供应商
indexT=yunshu3(a,1);
if R<=6000
R=R+buy(j,i);
transport(j,i)=indexT;
else
    a=a+1;
     if a>m
        a=m;
    end
    transport(j,i)=yunshu3(a,1);
    R=0;
    R=R+buy(j,i);
end
end
end
Transport=zeros(402,24)*nan;
Transport(C(:,1),:)=transport;
%%
%计算损失
Sunshi=zeros(20,24)*nan;
for j=1:24;%周
    for i=1:20%供应商
        Shunshi(i,j)=buy(i,j)*yunshu(transport(i,j),j)
    end
end
Shunshi2=zeros(402,24);
Shunshi2(C(:,1),:)=Shunshi;
if length(find(Changku>=0))==24 && changku>2*T
    F(count)=sum(sum(buy))*1.5+sum(Changku)*0.5;
    count=count+1
    if F(count-1) 
第三问代码
 
我的第三问其实和第二问差不多，还是选第一问得出的五十家供应商，只不过变成三种材料都采购了，并且给采购ABC的数量加了个权重，大概就是随机出三个大于0的数，它们的和为1，之后排序，储存在r2里，通过乘这个权重，达到多买A，少买C
 
第三问代码
 
这个代码感觉还有非常大大大的优化空间
 
因为有的时候会模拟出一些离谱的订单方案，比如说连仓库都欠账了之类的所以在最后加了一个判定，通过判定了才是有效模拟，直接count+1，有效模拟次数是count-1。
 
其实这个判定和 有效模拟次数 这玩意，第二问原本是没的，但是我在第三问写了这玩意之后，就有一种感觉“不行，第二问也要写，多水点图片”，所以每个问都有了
 
%第三问代码
%若出现某些错误，导致终止程序，可能是内存问题，重新运行即可
clc;clear;
[num,txt,raw]=xlsread('问题1评分排名.xlsx');
data=xlsread('附件1 近5年402家供应商的相关数据.xlsx',2);
index=1:402;
num2=[index' num];
num3=sortrows(num2,7);%对排名进行升序
wushi=num3(1:50,:);
%找出卖A的
indexA=find(wushi(:,2)==1.2);
A=wushi(indexA,:);
%找出卖B的
indexB=find(wushi(:,2)==1.1);
B=wushi(indexB,:);
%找出卖C的
indexC=find(wushi(:,2)==1);
C=wushi(indexC,:);
for i=1:16
MaxA(i)=max(data(A(i,1),:));
end
for i=1:14
MaxB(i)=max(data(B(i,1),:));
end
for i=1:20
MaxC(i)=max(data(C(i,1),:));
end
F=[]
Buy2=[];
R2=[];
Transport=[];
Shunshi2=[];
f=1000000000;
count=1;
X=100%模拟次数
for Xunhuan=1:X
x=rand(1,3);y=sum(x);r=x/y;
r2=sort(r,2)
%TA=28200*0.6*r2(3);
%TB=28200*0.66*r2(2);
%TC=28200*0.72*r2(1);
if TC<1
    TC=1
end
changkuA=0;
changkuB=0;
changkuC=0;%表示出库材料存量
buyA=zeros(16,24);
buyB=zeros(14,24);
buyC=zeros(20,24);
%%
%购买材料A
for i=1:16
buyA(i,1)=min(round(TC*0.5),randi(round([MaxA(i)*0.5,MaxA(i)])));
end
if sum(buyA(:,1))>=TA
    changkuA=sum(buyA(:,1))-TA;
    ChangkuA(1)=changkuA;
for j=2:24%星期
    for i=1:16%供应商
        fanwei=min(round(TC*0.5),randi(round([MaxA(i)*0.5,MaxA(i)])));
        if fanwei==0;
            fanwei=1
        end
        if sum(buyA(:,j))<=TA||changkuA<=2*TA
            buyA(i,j)=randi(fanwei);
            if sum(buyA(:,j))>=TA
                changkuA=changkuA+sum(buyA(:,j))-TA;
                if changkuA>=2*TA
                    continue
                end
            end
        end
    end
    if sum(buyA(:,j))<=TA
        changkuA=changkuA-(TA-sum(buyA(:,j)));
    end
    ChangkuA(j)=changkuA;%记录每周仓库的存储数量
end
end
%%
%购买材料B
for i=1:14
buyB(i,1)=min(round(TB*0.5),randi(round([MaxB(i)*0.5,MaxB(i)])));
end
if sum(buyB(:,1))>=TB
    changkuB=sum(buyB(:,1))-TB;
    ChangkuB(1)=changkuB;
for j=2:24%星期
    for i=1:14%供应商
        fanwei=min(round(TC*0.5),randi(round([MaxB(i)*0.5,MaxB(i)])));
        if fanwei==0;
            fanwei=1;
        end
        if sum(buyB(:,j))<=TB||changkuB<=2*TB
            buyB(i,j)=randi(fanwei);
            if sum(buyB(:,j))>=TB
                changkuB=changkuB+sum(buyB(:,j))-TB;
                if changkuB>=2*TB
                    continue
                end
            end
        end
    end
    if sum(buyB(:,j))<=TB
        changkuB=changkuB-(TB-sum(buyB(:,j)));
    end
    ChangkuB(j)=changkuB;%记录每周仓库的存储数量
end
end
%%
%购买材料C
for i=1:20
buyC(i,1)=min(round(TC*0.5),randi(round([MaxC(i)*0.5,MaxC(i)])));
end
if sum(buyC(:,1))>=TC
    changkuC=sum(buyC(:,1))-TC;
    ChangkuC(1)=changkuC;
for j=2:24%星期
    for i=1:20%供应商
        fanwei=min(round(TC*0.5),randi(round([MaxC(i)*0.5,MaxC(i)])));
        if fanwei==0
            fanwei=1
        end
        if sum(buyC(:,j))<=TC||changkuC<=2*TC
            buyC(i,j)=randi(fanwei);
            if sum(buyC(:,j))>=TC
                changkuC=changkuC+sum(buyC(:,j))-TC;
                if changkuC>=2*TC
                    continue
                end
            end
        end
    end
    if sum(buyC(:,j))<=TC
        changkuC=changkuC-(TC-sum(buyC(:,j)));
    end
    ChangkuC(j)=changkuC;%记录每周仓库的存储数量
end
end
Buy=zeros(50,24);
Buy(A(:,7),:)=buyA;
Buy(B(:,7),:)=buyB;
Buy(C(:,7),:)=buyC;
%数据合并
buy=zeros(402,24)*nan;
buy(A(:,1),:)=buyA;
buy(B(:,1),:)=buyB;
buy(C(:,1),:)=buyC;
%xlswrite('订购方案2',buy);
%%
%模拟转运损耗
fujian2=xlsread('附件2 近5年8家转运商的相关数据.xlsx','B2:IG9');
yunshu=zeros(5,24);
for j=1:8%第j个转运商
hang=fujian2(j,:);
%sunhao=unique(hang);
%[m,n]=size(sunhao);
for i=1:24
xunzhe=randi(240);
yunshu(j,i)=hang(xunzhe);
end
end
%%
%实施转运
p=1:8;%运输排名
transport=zeros(50,24);
for i=1:24%24周
    a=1;%第a行
    R=0;%运输量
yunshu2=[p' yunshu(:,i)];
yunshu2(find(yunshu2(:,2)==0),:)=[];%剔除不运输的转运商
[m,n]=size(yunshu2);
yunshu3=sortrows(yunshu2,2);%进行降序
for j=1:50%50家供应商
indexT=yunshu3(a,1);
if Buy(j,i)==0
    continue
end
if R<=6000
R=R+Buy(j,i);
transport(j,i)=indexT;
else
    a=a+1;
    if a>m
        a=m;
    end
    transport(j,i)=yunshu3(a,1);
    R=0;
    R=R+Buy(j,i);
end
end
end
Transport2=ones(402,24)*nan;
Transport2(wushi(:,1),:)=transport;
%%
%计算损失
Sunshi=zeros(50,24);
for j=1:24%周
    for i=1:50%供应商
        if transport(i,j)==0
            continue
        end
        Shunshi(i,j)=Buy(i,j)*yunshu(transport(i,j),j);
    end
end
Shunshi2=ones(402,24)*nan;
Shunshi2(wushi(:,1),:)=Shunshi;
if length(find(ChangkuA>=0))==24&&length(find(ChangkuB>=0))==24&&length(find(ChangkuC>=0))==24&&sum(sum(Buy))>406080
    F(count)=sum(sum(buyA))*1.2+sum(sum(buyB))*1.1+sum(sum(buyC))+0.5*(sum(ChangkuA)+sum(ChangkuB)+sum(ChangkuC))+sum(sum(Buy))*0.5;
    count=count+1;
    if F(count-1) 
第四问
 
第四问和第三问类似，受到另外一篇文章启发，画图发现购买材料C越来越多，B相对平稳，A逐步下降，所以那个权重就变成了多买C，少买A，然后为了产能增加，我经过乱七八糟的计算，我直接取了个产量最多的年份的平均周产量，其实就是算年的之后除48（小声bb），然后取了40625，我就不负责的假设他为了有这么多生产能力
 
关于生产力提高这方面，我个人认为也能用预测模型预测，未来能达到的产量，预测好之后也能直接放到我这个代码，这个改起来很简单
 
第四问代码
 
杂项
 
这个代码被我删删改改，丢了一部分。。。。
 
%%
%前五十供应商画图
clc;clear
[num2,txt2,~]=xlsread('问题1评分排名.xlsx');
index=1:402;
num3=[index' num2];
num4=sortrows(num3,7);%对排名进行升序
num5=num4(1:50,:);
wushi=sortrows(num5,1);
indexA=find(wushi(:,2)==1.2);
A=wushi(indexA,:);
indexB=find(wushi(:,2)==1.1);
B=wushi(indexB,:);
indexC=find(wushi(:,2)==1);
C=wushi(indexC,:);
%平均每周
a1=A(:,3)./420
b1=B(:,3)./420
c1=C(:,3)./420
%平均每月
a2=A(:,3)./60
b2=B(:,3)./60
c2=C(:,3)./60
%每月画图
plot(A(:,1),a2,'-k*','linewidth',1)
title('\fontname{宋体}材料A每月平均供应量');
xlabel('\fontname{宋体}供应商\fontname{Times New Roman}ID');
ylabel('\fontname{宋体}材料供应量');
set(gca,'FontName','Times New Roman','fontsize',10.5)
box on
grid on
xlim([1,240])
figure
plot(B(:,1),b2,'-k*','linewidth',1)
title('\fontname{宋体}材料B每月平均供应量');
xlabel('\fontname{宋体}供应商\fontname{Times New Roman}ID');;
ylabel('\fontname{宋体}材料供应量');
set(gca,'FontName','Times New Roman','fontsize',10.5)
box on
grid on
xlim([1,240])
figure
plot(C(:,1),c2,'-k*','linewidth',1)
title('\fontname{宋体}材料C每月平均供应量');
xlabel('\fontname{宋体}供应商\fontname{Times New Roman}ID');
ylabel('\fontname{宋体}材料供应量');
set(gca,'FontName','Times New Roman','fontsize',10.5)
box on
grid on
xlim([1,240])
%%
%分离附件1订单材料ABC
clc;clear
[num,txt,~]=xlsread('附件1 近5年402家供应商的相关数据.xlsx',1);
cailiao=ones(402,1);
xunzhaoA = strcmp(txt(2:end,2),'A');%找出A，标记为1
mA=find(xunzhaoA==1);%找出字符串所在的索引
xunzhaoB = strcmp(txt(2:end,2),'B');%找出A，标记为1
mB=find(xunzhaoB==1);%找出字符串所在的索引
xunzhaoC = strcmp(txt(2:end,2),'C');%找出A，标记为1
mC=find(xunzhaoC==1);%找出字符串所在的索引
cailiao(mA)=1.2;
cailiao(mB)=1.1;
cailiao(mC)=1;
indexA=find(cailiao(:,1)==1.2);
A=num(indexA,:);
indexB=find(cailiao(:,1)==1.1);
B=num(indexB,:);
indexC=find(cailiao(:,1)==1);
C=num(indexC,:);
shengchanli=sum(reshape(sum(A),48,5))./0.6+sum(reshape(sum(B),48,5))./0.66+sum(reshape(sum(C),48,5))./0.72;
plot(shengchanli)
 
模拟代码
 
%第四问模拟代码
%若出现某些错误，终止程序，可能是内存问题，重新运行即可
clc;clear;
[num,txt,raw]=xlsread('问题1评分排名.xlsx');
data=xlsread('附件1 近5年402家供应商的相关数据.xlsx',2);
index=1:402;
num2=[index' num];
num3=sortrows(num2,7);%对排名进行升序
wushi=num3(1:50,:);
%找出卖A的
indexA=find(wushi(:,2)==1.2);
A=wushi(indexA,:);
%找出卖B的
indexB=find(wushi(:,2)==1.1);
B=wushi(indexB,:);
%找出卖C的
indexC=find(wushi(:,2)==1);
C=wushi(indexC,:);
for i=1:16
MaxA(i)=max(data(A(i,1),:));
end
for i=1:14
MaxB(i)=max(data(B(i,1),:));
end
for i=1:20
MaxC(i)=max(data(C(i,1),:));
end
F=[]
Buy2=[];
R2=[];
Transport=[];
Shunshi2=[];
f=1000000000;
count=1;
X=1000%模拟次数
for Xunhuan=1:X
x=rand(1,3);
y=sum(x);
r=x/y;
r2=sort(r,2)
TA=40625*0.6*r2(1);
TB=40625*0.66*r2(2);
TC=40625*0.72*r2(3);
if TA<1
    TA=1
end
changkuA=0;
changkuB=0;
changkuC=0;%表示出库材料存量
buyA=zeros(16,24);
buyB=zeros(14,24);
buyC=zeros(20,24);
%%
%购买材料A
for i=1:16
buyA(i,1)=min(round(TC*0.5),randi(round([MaxA(i)*0.5,MaxA(i)])));
end
if sum(buyA(:,1))>=TA
    changkuA=sum(buyA(:,1))-TA;
    ChangkuA(1)=changkuA;
for j=2:24%星期
    for i=1:16%供应商
        fanwei=min(round(TC*0.5),randi(round([MaxA(i)*0.5,MaxA(i)])));
        if fanwei==0;
            fanwei=1
        end
        if sum(buyA(:,j))<=TA||changkuA<=2*TA
            buyA(i,j)=randi(fanwei);
            if sum(buyA(:,j))>=TA
                changkuA=changkuA+sum(buyA(:,j))-TA;
                if changkuA>=2*TA
                    continue
                end
            end
        end
    end
    if sum(buyA(:,j))<=TA
        changkuA=changkuA-(TA-sum(buyA(:,j)));
    end
    ChangkuA(j)=changkuA;%记录每周仓库的存储数量
end
end
%%
%购买材料B
for i=1:14
buyB(i,1)=min(round(TB*0.5),randi(round([MaxB(i)*0.5,MaxB(i)])));
end
if sum(buyB(:,1))>=TB
    changkuB=sum(buyB(:,1))-TB;
    ChangkuB(1)=changkuB;
for j=2:24%星期
    for i=1:14%供应商
        fanwei=min(round(TC*0.5),randi(round([MaxB(i)*0.5,MaxB(i)])));
        if fanwei==0;
            fanwei=1;
        end
        if sum(buyB(:,j))<=TB||changkuB<=2*TB
            buyB(i,j)=randi(fanwei);
            if sum(buyB(:,j))>=TB
                changkuB=changkuB+sum(buyB(:,j))-TB;
                if changkuB>=2*TB
                    continue
                end
            end
        end
    end
    if sum(buyB(:,j))<=TB
        changkuB=changkuB-(TB-sum(buyB(:,j)));
    end
    ChangkuB(j)=changkuB;%记录每周仓库的存储数量
end
end
%%
%购买材料C
for i=1:20
buyC(i,1)=min(round(TC*0.5),randi(round([MaxC(i)*0.5,MaxC(i)])));
end
if sum(buyC(:,1))>=TC
    changkuC=sum(buyC(:,1))-TC;
    ChangkuC(1)=changkuC;
for j=2:24%星期
    for i=1:20%供应商
        fanwei=min(round(TC*0.5),randi(round([MaxC(i)*0.5,MaxC(i)])));
        if fanwei==0
            fanwei=1
        end
        if sum(buyC(:,j))<=TC||changkuC<=2*TC
            buyC(i,j)=randi(fanwei);
            if sum(buyC(:,j))>=TC
                changkuC=changkuC+sum(buyC(:,j))-TC;
                if changkuC>=2*TC
                    continue
                end
            end
        end
    end
    if sum(buyC(:,j))<=TC
        changkuC=changkuC-(TC-sum(buyC(:,j)));
    end
    ChangkuC(j)=changkuC;%记录每周仓库的存储数量
end
end
Buy=zeros(50,24);
Buy(A(:,7),:)=buyA;
Buy(B(:,7),:)=buyB;
Buy(C(:,7),:)=buyC;
%数据合并
buy=zeros(402,24)*nan;
buy(A(:,1),:)=buyA;
buy(B(:,1),:)=buyB;
buy(C(:,1),:)=buyC;
%xlswrite('订购方案2',buy);
%%
%模拟转运损耗
fujian2=xlsread('附件2 近5年8家转运商的相关数据.xlsx','B2:IG9');
yunshu=zeros(5,24);
for j=1:8%第j个转运商
hang=fujian2(j,:);
%sunhao=unique(hang);
%[m,n]=size(sunhao);
for i=1:24
xunzhe=randi(240);
yunshu(j,i)=hang(xunzhe);
end
end
%%
%实施转运
p=1:8;%运输排名
transport=zeros(50,24);
for i=1:24%24周
    a=1;%第a行
    R=0;%运输量
yunshu2=[p' yunshu(:,i)];
yunshu2(find(yunshu2(:,2)==0),:)=[];%剔除不运输的转运商
[m,n]=size(yunshu2);
yunshu3=sortrows(yunshu2,2);%进行降序
for j=1:50%50家供应商
indexT=yunshu3(a,1);
if Buy(j,i)==0
    continue
end
if R<=6000
R=R+Buy(j,i);
transport(j,i)=indexT;
else
    a=a+1;
    if a>m
        a=m;
    end
    transport(j,i)=yunshu3(a,1);
    R=0;
    R=R+Buy(j,i);
end
end
end
Transport2=ones(402,24)*nan;
Transport2(wushi(:,1),:)=transport;
%%
%计算损失
Sunshi=zeros(50,24);
for j=1:24%周
    for i=1:50%供应商
        if transport(i,j)==0
            continue
        end
        Shunshi(i,j)=Buy(i,j)*yunshu(transport(i,j),j);
    end
end
Shunshi2=ones(402,24)*nan;
Shunshi2(wushi(:,1),:)=Shunshi;
if length(find(ChangkuA>=0))==24&&length(find(ChangkuB>=0))==24&&length(find(ChangkuC>=0))==24&&sum(sum(Buy))>702000
    F(count)=sum(sum(buyA))*1.2+sum(sum(buyB))*1.1+sum(sum(buyC))+0.5*(sum(ChangkuA)+sum(ChangkuB)+sum(ChangkuC))+sum(sum(Buy))*0.5;
    count=count+1;
    if F(count-1) 
其实还有一个模型检验的代码，不过太敷衍了，我就不放出来了
 
关于转运方案，改一改代码就能换成对应的表格了，我当时没注意它的表格格式，弄成了和订单一样的样子。。。
 
我这个用是蒙特卡洛模拟的思想
 
最后，欢迎大家和我交流
 
 品牌网店加盟  小程序开发教程官网  建外贸网站公司  怎么做小程序  php企业网站  创建自己的个人网站