设施位置（misole03）
/*****************************************************************//* *//*S A S A M P L E L I B R A R Y公司*//* *//*名称：misole03*//*标题：设施位置（misole03）*//*产品：或*//*系统：所有*//*关键字：或*//*进程：OPTMODEL，GPLOT*//*数据：*//* *//*支持：更新：*//*参考：*//*MISC：混合整数线性规划示例3*//*数学规划的求解章节*//* *//*****************************************************************/标题“设施位置问题”；%让NumCustomers=50；%设站点数=10；%设SiteCapacity=35；%设MaxDemand=10；%设xmax=200；%设ymax=100；%设种子=938；/*生成随机客户位置*/数据cdata（drop=i）；长度名称$8；do i=1到&NumCustomers；name=压缩（'C'||put（i，最好））；x=兰尼（&seed）*&xmax；y=兰尼（&seed）*&ymax；需求=拉努尼（&seed）*&MaxDemand；产出；结束；运行；/*生成随机站点位置和固定费用*/数据sdata（drop=i）；长度名称$8；i=1到&NumSites；name=压缩（'SITE'||put（i，best.））；x=兰尼（&seed）*&xmax；y=兰尼（&seed）*&ymax；固定电荷=30*（绝对值（&xmax/2-x）+绝对值（&ymax/2-y））；产出；结束；运行；过程操作模型；设置<str>客户；设置<str>SITES初始化{}；/*客户和场地的x和y坐标*/num x{客户联盟站点}；num y{客户联盟站点}；num需求{客户}；num固定电荷{SITES}；/*客户i到站点j的距离*/客户中的num dist{i，站点中的j}=平方（（x[i]-x[j]）^2+（y[i]-y[j]，^2）；将数据cdata读入CUSTOMERS=[name]x y demand；将数据sdata读入SITES=[name]x y fixed_charge；var分配{CUSTOMERS，SITES}二进制；var Build{SITES}二进制；最小成本无固定费用=总和{i在CUSTOMERS中，j在SITES中}dist[i，j]*分配[i，j]；最低成本固定费用=成本NoFixedCharge+现场费用总额{j}固定费用[j]*建造费用[j]；/*每个客户只分配到一个站点*/conassign_def{i在客户}中：站点中的和{j}赋值[i，j]=1；/*如果客户i分配到站点j，则设施必须在j建造*/con链接{i在客户中，j在站点}中：分配[i，j]<=构建[j]；/*每个站点最多可以处理SiteCapacity需求*/容量{j现场}：客户需求[i]*分配[i，j]中的总和{i<=&站点容量*构建[j]；/*无固定费用解决MILP*/以milp/logfreq=500求解对象成本NoFixedCharge；/*清理溶液*/对于用户中的{i，站点中的j分配[i，j]=round（分配[i、j]）；对于{j in SITES}Build[j]=圆形（Build[j]）；调用信号（“varcostNo”，put（CostNoFixedCharge，6.1））；/*创建供GPLOT使用的数据集*/从创建数据CostNoFixedCharge_data[客户站点]={i在客户中，j在站点中：分配[i，j]=1}xi=x[i]yi=y[i]xj=x[j]yj=y[j]；/*解决MILP问题，使用热启动的固定电荷*/用milp/primalin logfreq=500求解对象CostFixedCharge；/*清理溶液*/对于用户中的{i，站点中的j分配[i，j]=round（分配[i、j]）；对于{j in SITES}Build[j]=圆形（Build[j]）；num varcost=sum{i in CUSTOMERS，j in SITES}dist[i，j]*Assign[i，j].sol；num fixcost=sum{j in SITES}fixed_charge[j]*构建[j].sol；调用信号（“varcost”，put（varcost，6.1））；调用symput（“固定成本”，put（固定成本，5.1））；呼叫信号（“总成本”，put（CostFixedCharge，6.1））；/*创建供GPLOT使用的数据集*/从创建数据CostFixedCharge_data[客户站点]={i在客户中，j在站点中：分配[i，j]=1}xi=x[i]yi=y[i]xj=x[j]yj=y[j]；退出；title 1 h=1.5“设施位置问题”；title 2“总成本=&varcostNo（变量=&varcostNo，固定=0）”；数据csdata；设置cdata（重命名=（y=cy））sdata（rename=（y=sy））；运行；/*创建注释数据集以在客户和指定站点之间绘制线*/%annomac；数据anno（drop=xi yi xj yj）；%系统（2，2，2）；设置CostNoFixedCharge_Data（keep=xi yi xj yj）；%线（xi，yi，xj，yj，*，1，1）；运行；proc gplot data=csdata anno=anno；axis1标签=无顺序=（0到&xmax乘以10）；轴2标签=无顺序=（0到&ymax乘以10）；symbol1值=点interpol=无pointlabel=（“#name”nodropcollisions height=1）cv=黑色；符号2值=菱形插值=无pointlabel=（“#name”nodropcollisions color=蓝色高度=1）cv=蓝色；绘制cy*xsy*x/覆盖haxis=axis1 vaxis=axis2；运行；退出；title1“设施选址问题”；title2“总成本=&TotalCost（变量=&varcost，固定=&fixcost）”；/*创建注释数据集以在客户和指定站点之间绘制线*/数据anno（drop=xi yi xj yj）；%系统（2，2，2）；设置CostFixedCharge_Data（keep=xi yi xj yj）；%线（xi，yi，xj，yj，*，1，1）；运行；proc gplot data=csdata anno=anno；axis1标签=无顺序=（0到&xmax乘以10）；轴2标签=无顺序=（0到&ymax乘以10）；symbol1值=点interpol=无pointlabel=（“#name”nodropcollisions height=1）cv=黑色；符号2值=菱形插值=无pointlabel=（“#name”nodropcollisions color=蓝色高度=1）cv=蓝色；绘制cy*xsy*x/覆盖haxis=axis1 vaxis=axis2；运行；退出；%让NumCustomers=45；%设站点数=8；%设SiteCapacity=35；%设MaxDemand=10；%设xmax=200；%设ymax=100；%设种子=2345；/*生成随机客户位置*/数据cdata（drop=i）；长度名称$8；do i=1到&NumCustomers；name=压缩（'C'||put（i，最好））；x=兰尼（&seed）*&xmax；y=兰尼（&seed）*&ymax；需求=拉努尼（&seed）*&MaxDemand；产出；结束；运行；/*生成随机站点位置和固定费用*/数据sdata（drop=i）；长度名称$8；i=1到&NumSites；name=压缩（'SITE'||put（i，best.））；x=兰尼（&seed）*&xmax；y=兰尼（&seed）*&ymax；固定电荷=（abs（&xmax/2-x）+abs（-ymax/2-y））/2；产出；结束；运行；过程操作模型；设置<str>客户；设置<str>SITES初始化{}；/*客户和场地的x和y坐标*/num x{客户联盟站点}；num y{客户联盟站点}；num需求{客户}；num固定电荷{SITES}；/*客户i到站点j的距离*/客户中的num dist{i，站点中的j}=平方（（x[i]-x[j]）^2+（y[i]-y[j]，^2）；将数据cdata读入CUSTOMERS=[name]x y demand；将数据sdata读入SITES=[name]x y fixed_charge；var分配{CUSTOMERS，SITES}二进制；var Build{SITES}二进制；最低成本固定费用=客户中的sum{i，站点中的j}dist[i，j]*赋值[i，j]+sum{j in SITES}fixed_charge[j]*构建[j]；/*每个客户只分配到一个站点*/conassign_def{i在客户}中：站点中的和{j}赋值[i，j]=1；/*如果客户i分配到站点j，则设施必须在j建造*/con链接{i在客户中，j在站点}中：分配[i，j]<=构建[j]；/*每个站点最多可以处理SiteCapacity需求*/容量{j现场}：客户需求[i]*分配[i，j]<=&SiteCapacity*Build[j]；/*为构建变量分配优先级（y）*/对于{j in SITES}Build[j].priority=10；/*使用分支优先级，用固定费用求解MILP*/以milp/logfreq=1000求解对象CostFixedCharge；退出；