数学建模论文-企业的营销管理问题.doc

企业的营销管理问题 一：摘要摘要：本题为典型的求生产“最优”问题，这就用到数学建模中的规划模型, 求得在公司或者营销部门的最大利润。销售的产品分为计划外与计划内两类，计划内又分已签约量与意向签约量，销售部门可从计划外的产品销售中获得利润，再在这样的基础上运用目标规划的方法来解决题中所给的五个问题。针对问题（1）（2）分别建立以求公司的最大利润为目标的规划模型和营销部的总收入极大化为目标的规划模型，得到目标函数为：公司的利润=计划内的销售额+计划外销售部缴纳的利润-所有的成本-发放的经费与营销部的总收入=公司发放的经费+计划外卖出的销售额-计划外销售部缴纳的利润-宣传费， 由lingo软件求解,得到公司的最大的利润为:2520.580（万元）和营销部的最大利润为：931.9888（万元）。针对问题（3）是建立以兼顾公司和营销部二者的利益为目标的双目标规划模型，即是将公司利润和营销部门利润之和达到最大，建立目标函数。则我们将双目标规划问题转化为单目标问题求解。由lingo软件求解，兼顾公司和营销部的总利润最大：2689.400（万元）针对问题（4）：是建立在问题三的基础上的，兼顾公司、营销部的利益以及客户的需求，尽量做到均衡销售。我们转化为求公司、营销部、客户的最大满意度建立模型，用lingo软件求解当满意度为最大时的利润为：2512.910（万元）针对问题（5）：是建立在问题（2）的基础上的，由于价格波动影响销售量，所以引入变量变量来表示价格的上涨率，则为销售价格的变化率，建立模型并用lingo软件求解得到使营销部总收入最大的定价对应值为0.2198512，对应定价分别为：7.31、7.92、8.23、19.3、22.2、23.7、23.1、23.9、27.2、28.0（十元）。关键词：拟合函数 利润最大化 非线性规划 最优问题 满意度二：问题重述任何一个企业都面临着将产品推向市场，销售出去。一般的，企业下设营销部，企业生产的产品要通过营销部进行销售。由于企业生产能力的制约，因此需要在满足已签约的销售合同量的基础上，对意向签约量有选择的安排生产。一方面，企业会组织安排生产，完成已签约的销售合同；另一方面，企业希望营销部门尽力争取与意向签约的客户签订正式销售合同；除此之外，企业还希望销售部门努力再多销售一些产品（计划外）。对于所签约的销售合同和意向签约量（计划内），企业根据销售量向营销部发放经费（包括工资以及宣传费用等）；对于计划外销售的产品，实行承包制，营销部向企业缴纳利润。产品的生产费用由企业承担，与销售有关的费用（如产品的宣传费用等）由营销部承担。宏宇电器公司2011年计划生产三类10种小家电，其中包括：热水壶（1.5升，1.8升，2升）、豆浆机（0.9升，1.1升，1.3升）、电饭煲（2升，2.5升，3升，3.5升）。三类小家电的年最大生产能力分别为：热水壶：5万个；豆浆机：6.5万个；电饭煲6.2万个目前，10种小家电的“已签约合同量”以及“意向签约量”见附件的表1。公司的生产成本详见附件的表2。每种小家电的宣传费用随着销售量而变化，详见附件的表3。每种产品的销售额随着每单合同的订购量而变化，详见附件中的表4。公司现已签约订单八个，每单的签约量见附件的表5。公司现有意向签约客户五家，每家的意向签约量见附件的表6。对于每个有签约意向的客户，可能签订订购合同量的概率见附件中的表7-表11。对于营销部销售计划外的产品，可能签订订购合同量（万个）的概率见附件中的表12。对于计划内销售的产品，企业向营销部发放的经费（万元）见附件的表13。对于计划外销售的部分，营销部向企业缴纳的利润（万元）见附件的表14。请分别按以下要求为宏宇电器公司制定相应的生产、销售方案：（1）使公司的利润达到最大；（2）使营销部的总收入极大化；（3）兼顾公司和营销部二者的利益；（4）兼顾公司、营销部的利益以及客户的需求，尽量做到均衡销售；（5）公司为了调动营销部的积极性，为公司多创利润，规定对于计划外销售的产品，营销部可以自行定价。统计数据显示，自行定价后每单位产品的销售价格上涨或下降5，则三类小家电的可能销售量相应的减少或增加10。试确定使营销部总收入最大的定价、生产及销售方案。 三：模型假设 1:经费只包括题目中的工资和宣传费，没有其他的费用2:生产出的产品都是合格的3:在探讨期间内公司持续经营4: 在一定时间内，市场售价保持不变(物价稳定) 四：符号说明 1. 表示由表4拟合的函数,即第种小家电在相应订购量为时的 销售额（单位：万元）2. 表示由表2拟合的函数,即第种小家电在相应产量为时的总成 本（单位：万元）3. 表示由表13拟合的函数,即企业向营销部发放计划内销售第种 小家电在相应销售量为时的经费（万元）4. 表示第个订单已签约第种小家电的签约量5. 表示第个客户对第种小家电的意向签约量6. 表示计划外对第种小家电的签约量7. 表示第种小家电的总的销售量的最大生产值8. 表示第个客户对第种小家电的意向签约量的最大值9. 表示第12个表中第个小家电的最大合同量10. 表示由表3拟合的函数，即第种小家电在相应销售量为时的宣传费用（单位：万元）11. 表示销售价格的变化率12. 表示每一题所求利润 五：问题分析先观察附件中表里的数据，分析表里各变量的关系，用matlab将其各表的变量拟合成相关的连续函数。对于第（1）（2）问都是求最优解，则我们用线性规划的模型来建立。问题（1）使公司的利润达到最大，则我们牵涉到公司计划内的销售量，计划外销售部缴纳的利润，所有的成本，发放的经费问题，在这样的基础上我们建立线性规划的模型。问题（2）使营销部的总收入极大化，则也是考虑得到的利润最大的问题，则这主要牵涉到公司发放的经费，计划外卖出的销售额，计划外销售部缴纳的利润，宣传费问题，在这样的基础上我们建立线性规划的模型。对于第（3）问是要求兼顾公司和营销部二者的利益，则也和（1）（2）问题一样求最大的利润，只是第（3）问是要求公司和营销部的总利润的最大。 对于第（4）问要求兼顾公司、营销部的利益以及客户的需求，尽量做到均衡销售。我们转化为求公司、营销部、客户的最大满意度。公司、营销部的满意度可以用是否达到最大利润来刻画，公司、营销部的满意度=。客户的满意度用客户是否达到最大的签约量来刻画，客户的满意度= 。均衡销售用各类家电销售量的方差来刻画，刻画方程=其中均值为:，。第（5）问，是对计划外销售的产品，营销部可以自行定价。由销售部自行定产品的价格，由数据统计显示价格波动影响了销售量，并可以知道价格量的变化会引起相应的销售量两倍的影响 。每单位产品的销售价格=，当每单位产品的销售价格上涨与减少后，得到的销售价格=，销售量=。所以价格变化后的销售额=计划外销售部缴纳的利润、宣传费只与家电的销售量有关，所以计划外销售部缴纳的利润=宣传费= 六：模型建立及求解（一）附件中各表数据的处理 针对表2、3、4、13、14，分别对相应成本费用、宣传费用、销售额、经费、上缴利润建立以各产品产量为自变量的函数关系，并用matlab对函数进行拟合，通过对比，找出最佳拟合函数为二次函数。具体程序见附件1-5。 各表对应拟合函数见附件11（二）问题（1）目标函数：公司的利润=计划内的销售额+计划外销售部缴纳的利润-所有的成本-发放的经费计划内的销售额：表示由各小家电已签约量产生的销售额和意向签约量产生的销售额组成所有成本：表示计划外核计划内的所有小家电的成本发放的经费：表示各小家电的意向签约量的经费和已签约量的经费即： 约束条件： 运用lingo软件对模型求解，得到公司的最大的利润为:2520.580（单位：万元）对应方案为：序号12345678910对应已签约生产量15.015.917.116.817.217.913.914.114.714.6意向签约生产量0.670.690.513.954.604.420.541.300.831.79计划外生产量0000000000总产量15.616.617.620.721.822.314.415.415.516.4注：上表中各生产量单位：千件；具体求解程序见附件6。（三）问题（2）目标函数：营销部的总收入=公司发放的经费+计划外卖出的销售额-计划外销售部缴纳的利润-宣传费宣传费：表示各小家电已签约和意向签约量所花费的宣传费用约束条件： 运用lingo软件对模型求解，得营销部的最大利润为：931.9888（单位：万元）对应方案为：序号12345678910对应已签约生产量1.967.359.685.0210.211.72.835.417.168.08对应意向签约产量1.967.359.685.0210.211.72.835.417.168.08对应计划外生产量4.04.04.04.04.03.04.05.03.03.0总产量7.9218.723.314.024.426.49.6615.817.319.2注：上表中各生产量单位：千件；具体求解程序见附件7。（四）问题（3）目标函数： 约束条件： 运用lingo软件对模型求解，兼顾公司和营销部的总利润最大：2689.400（单位：万元）对应方案为：序号12345678910对应已签约生产量14.815.716.915.516.116.713.613.814.414.3对应意向签约产量0.450.490.312.243.152.900.100.910.431.44对应计划外生产量0.450.490.312.243.152.900.100.910.431.44总产量15.716.717.520.022.422.513.815.615.317.2注：上表中各生产量单位：千件；具体求解程序见附件8。（五）：问题（4）目标函数：其中，约束条件： 运用lingo软件对模型求解，均衡销售的最大利润：2512.910 （单位：万元）对应方案为：序号12345678910对应已签约生产量000015.00022.100对应意向签约产量046.00045.90034.800对应计划外生产量04.0004.0005.000总产量0.0050.00.0000.064.90.000.0061.70.000.00注：上表中各生产量单位：千件；具体求解程序见附件9。（六）：问题（5）建立在问题（2）的基础上建立模型目标函数约束条件： 运用lingo软件对模型求解，均衡销售的最大利润：1031.033 （单位：万元）,对应值：0.2198512对应方案为：序号12345678910对应已签约生产量1.967.359.685.0210.211.72.835.417.168.08对应意向签约产量1.967.359.685.0210.211.72.835.417.168.08对应计划外生产量3.993.993.993.993.992.993.994.992.992.99总产量7.9118.623.414.024.426.49.6515.817.319.1对应定价(十元)7.317.928.2319.322.223.723.123.927.228.0注：上表中各生产量单位：千件；具体求解程序见附件10。七：模型的评价及改进模型评价优点：1、该模型有效的刻画公司的最大利润，营销部的最大利润，公司和营销部二者的最大利益，公司、营销部的利益以及客户的需求。2、模型具有通用性缺点：1、 拟合的函数存在误差，使得得到的结果不是最精确的，只能是较精确的解。模型的改进 由于数据较少、较分散，使得拟合函数存在一定的误差。所以应将数据进行细化，得到更好的拟合函数。附件1：clearclfclcx=1 3 5 7 10 15 20 25 30 35 40;a=zeros(10,11);a=xlsread(e:数模b2.xls);%1绿色 plot(x,a(1,:),g) hold on ; a=polyfit(x,a(1,:),2); aa=poly2str(a,x) %显示拟合的表达式1 %2蓝色 plot(x,a(2,:),b) hold on ; b=polyfit(x,a(2,:),2); bb=poly2str(b,x) %显示拟合的表达式2%3红色 plot(x,a(3,:),r) hold on ; c=polyfit(x,a(3,:),2); cc=poly2str(c,x) %显示拟合的表达式3%4青色 plot(x,a(4,:),c) hold on ; d=polyfit(x,a(4,:),2);dd=poly2str(d,x) %显示拟合的表达式4 %5黑色 plot(x,a(5,:),k) hold on ; e=polyfit(x,a(5,:),2); ee=poly2str(e,x) %显示拟合的表达式5 %6黄色 plot(x,a(6,:),y) hold on ; f=polyfit(x,a(6,:),2); ff=poly2str(f,x) %显示拟合的表达式6 %7叉号 plot(x,a(7,:),x) hold on ; g=polyfit(x,a(7,:),2); gg=poly2str(g,x) %显示拟合的表达式7 %8* plot(x,a(8,:),*) hold on ; h=polyfit(x,a(8,:),2); hh=poly2str(h,x) %显示拟合的表达式8 %9 圆圈 plot(x,a(9,:),o) hold on ; i=polyfit(x,a(9,:),2); ii=poly2str(i,x) %显示拟合的表达式9 %10十字号 plot(x,a(10,:),+) hold on ; j=polyfit(x,a(10,:),2); jj=poly2str(j,x) %显示拟合的表达式10附件2：clearclfclcx=1 3 5 7 10 15 20 25 30 35 40;a=zeros(10,11);a=xlsread(e:数模b3.xls);%1绿色 plot(x,a(1,:),g) hold on ; a=polyfit(x,a(1,:),2); aa=poly2str(a,x) %显示拟合的表达式1 %2蓝色 plot(x,a(2,:),b) hold on ; b=polyfit(x,a(2,:),2); bb=poly2str(b,x) %显示拟合的表达式2%3红色 plot(x,a(3,:),r) hold on ; c=polyfit(x,a(3,:),2); cc=poly2str(c,x) %显示拟合的表达式3%4青色 plot(x,a(4,:),c) hold on ; d=polyfit(x,a(4,:),2);dd=poly2str(d,x) %显示拟合的表达式4 %5黑色 plot(x,a(5,:),k) hold on ; e=polyfit(x,a(5,:),2); ee=poly2str(e,x) %显示拟合的表达式5 %6黄色 plot(x,a(6,:),y) hold on ; f=polyfit(x,a(6,:),2); ff=poly2str(f,x) %显示拟合的表达式6 %7叉号 plot(x,a(7,:),x) hold on ; g=polyfit(x,a(7,:),2); gg=poly2str(g,x) %显示拟合的表达式7 %8* plot(x,a(8,:),*) hold on ; h=polyfit(x,a(8,:),2); hh=poly2str(h,x) %显示拟合的表达式8 %9 圆圈 plot(x,a(9,:),o) hold on ; i=polyfit(x,a(9,:),2); ii=poly2str(i,x) %显示拟合的表达式9 %10十字号 plot(x,a(10,:),+) hold on ; j=polyfit(x,a(10,:),2); jj=poly2str(j,x) %显示拟合的表达式10附件3：clearclfclcx=1 3 5 7 10 15 20 25 30 35 40;a=zeros(10,11);a=xlsread(e:数模b4.xls);%1绿色 plot(x,a(1,:),g) hold on ; a=polyfit(x,a(1,:),2); aa=poly2str(a,x) %显示拟合的表达式1 %2蓝色 plot(x,a(2,:),b) hold on ; b=polyfit(x,a(2,:),2); bb=poly2str(b,x) %显示拟合的表达式2%3红色 plot(x,a(3,:),r) hold on ; c=polyfit(x,a(3,:),2); cc=poly2str(c,x) %显示拟合的表达式3%4青色 plot(x,a(4,:),c) hold on ; d=polyfit(x,a(4,:),2);dd=poly2str(d,x) %显示拟合的表达式4 %5黑色 plot(x,a(5,:),k) hold on ; e=polyfit(x,a(5,:),2); ee=poly2str(e,x) %显示拟合的表达式5 %6黄色 plot(x,a(6,:),y) hold on ; f=polyfit(x,a(6,:),2); ff=poly2str(f,x) %显示拟合的表达式6 %7叉号 plot(x,a(7,:),x) hold on ; g=polyfit(x,a(7,:),2); gg=poly2str(g,x) %显示拟合的表达式7 %8* plot(x,a(8,:),*) hold on ; h=polyfit(x,a(8,:),2); hh=poly2str(h,x) %显示拟合的表达式8 %9 圆圈 plot(x,a(9,:),o) hold on ; i=polyfit(x,a(9,:),2); ii=poly2str(i,x) %显示拟合的表达式9 %10十字号 plot(x,a(10,:),+) hold on ; j=polyfit(x,a(10,:),2); jj=poly2str(j,x) %显示拟合的表达式10附件4：clearclfclcx=1 3 5 7 10 15 20 25 30 35 40;a=zeros(10,11);a=xlsread(e:数模b13.xls);%1绿色 plot(x,a(1,:),g) hold on ; a=polyfit(x,a(1,:),2); aa=poly2str(a,x) %显示拟合的表达式1 %2蓝色 plot(x,a(2,:),b) hold on ; b=polyfit(x,a(2,:),2); bb=poly2str(b,x) %显示拟合的表达式2%3红色 plot(x,a(3,:),r) hold on ; c=polyfit(x,a(3,:),2); cc=poly2str(c,x) %显示拟合的表达式3%4青色 plot(x,a(4,:),c) hold on ; d=polyfit(x,a(4,:),2);dd=poly2str(d,x) %显示拟合的表达式4 %5黑色 plot(x,a(5,:),k) hold on ; e=polyfit(x,a(5,:),2); ee=poly2str(e,x) %显示拟合的表达式5 %6黄色 plot(x,a(6,:),y) hold on ; f=polyfit(x,a(6,:),2); ff=poly2str(f,x) %显示拟合的表达式6 %7叉号 plot(x,a(7,:),x) hold on ; g=polyfit(x,a(7,:),2); gg=poly2str(g,x) %显示拟合的表达式7 %8* plot(x,a(8,:),*) hold on ; h=polyfit(x,a(8,:),2); hh=poly2str(h,x) %显示拟合的表达式8 %9 圆圈 plot(x,a(9,:),o) hold on ; i=polyfit(x,a(9,:),2); ii=poly2str(i,x) %显示拟合的表达式9 %10十字号 plot(x,a(10,:),+) hold on ; j=polyfit(x,a(10,:),2); jj=poly2str(j,x) %显示拟合的表达式10附件5clearclfclcx=1 2 3 4 5 6 7 8 9 10 11;a=zeros(10,11);a=xlsread(e:数模b14.xls);%1绿色 plot(x,a(1,:),g) hold on ; a=polyfit(x,a(1,:),2); aa=poly2str(a,x) %显示拟合的表达式1 %2蓝色 plot(x,a(2,:),b) hold on ; b=polyfit(x,a(2,:),2); bb=poly2str(b,x) %显示拟合的表达式2%3红色 plot(x,a(3,:),r) hold on ; c=polyfit(x,a(3,:),2); cc=poly2str(c,x) %显示拟合的表达式3%4青色 plot(x,a(4,:),c) hold on ; d=polyfit(x,a(4,:),2);dd=poly2str(d,x) %显示拟合的表达式4 %5黑色 plot(x,a(5,:),k) hold on ; e=polyfit(x,a(5,:),2); ee=poly2str(e,x) %显示拟合的表达式5 %6黄色 plot(x,a(6,:),y) hold on ; f=polyfit(x,a(6,:),2); ff=poly2str(f,x) %显示拟合的表达式6 %7叉号 plot(x,a(7,:),x) hold on ; g=polyfit(x,a(7,:),2); gg=poly2str(g,x) %显示拟合的表达式7 %8* plot(x,a(8,:),*) hold on ; h=polyfit(x,a(8,:),2); hh=poly2str(h,x) %显示拟合的表达式8 %9 圆圈 plot(x,a(9,:),o) hold on ; i=polyfit(x,a(9,:),2); ii=poly2str(i,x) %显示拟合的表达式9 %10十字号 plot(x,a(10,:),+) hold on ; j=polyfit(x,a(10,:),2); jj=poly2str(j,x) %显示拟合的表达式10附件6：model:sets:cp/1.10/;dg/1.8/;yx/1.5/;link(dg,yx);endsetsmax=a1+a2-b1-b2-b3-c1-c2+(-0.075559*z12 + 1.6933*z1+ 0.75536-0.076346 *z22 + 1.7398*z2+ 1.0628 -0.059056 *z32 + 1.5629*z3 + 1.5684-0.20613*z42 + 4.3918*z4 + 2.3385-0.22296*z52 + 4.892*z5 + 3.1453-0.22629*z62 + 5.0031*z6+ 3.3284 -0.23018*z72 + 5.0769*z7 + 3.3259 -0.25493*z82 + 5.4943*z8 + 3.4316-0.28777*z92 + 5.9454*z9+ 3.4438-0.28056*z102 + 5.9462*z10+ 3.7885);!计划内已签约销售额a1;for(dg(j):a1=-0.26713 *xj12 + 11.4183 *xj1 + 1.3873-0.27436* xj22 + 12.1778* xj2+ 2.0648-0.26305 *xj32 + 12.4521* xj3+ 2.5588-0.73776*xj42 + 29.8586*xj4 + 5.3673-0.85117 *xj52 + 34.4913*xj5 + 5.788 -0.83007*xj62 + 34.9027* xj6+ 6.9564-0.86107*xj72 + 35.2201*xj7 + 8.2703-0.94988*xj82 + 38.0877 *xj8 + 8.0776-0.93555*xj92 + 38.9202*xj9 + 11.3939-1.0184*xj102 + 41.1537*xj10+ 8.7158); !计划内意向签约销售额a2;for(yx(k):a2=-0.26713 *yk12 + 11.4183 *yk1 + 1.3873-0.27436* yk22 + 12.1778* yk2+ 2.0648-0.26305 *yk32 + 12.4521* yk3 + 2.5588-0.73776*yk42 + 29.8586*yk4 + 5.3673-0.85117*yk52+34.4913*yk5+5.7788 -0.83007*yk62 + 34.9027* yk6 + 6.9564-0.86107*yk72 + 35.2201*yk7 + 8.2703-0.94988*yk82+ 38.0877 *yk8 + 8.0776-0.93555*yk92 + 38.9202*yk9 + 11.3939-1.0184*yk102 + 41.1537*yk10+ 8.7158);!计划内已签约对应成本b1;for(dg(j):b1=-0.0603*xj12 + 7.5807*xj1+ 2.4845-0.067029*xj22 + 8.2675*xj2+ 2.7346-0.068437* xj32 + 8.5944*xj3+ 2.8063-0.18187 *xj42 + 20.0371*xj4+ 6.7665-0.20381*xj52 + 22.8783*xj5+ 7.8851-0.21575*xj62 + 23.7317*xj6+ 7.5642-0.2396*xj72 + 22.7827*xj7+ 15.4914-0.25671*xj82 + 24.41*xj8+ 16.598-0.27125 *xj92 + 25.7934*xj9+ 17.5138-0.2784 *xj102 + 26.4602*xj10+ 17.7686);!计划内意向签约对应成本b2;for(yx(k):b2=-0.0603*yk12 + 7.5807*yk1+ 2.4845-0.067029*yk22 + 8.2675*yk2+ 2.7346-0.068437*yk32 + 8.5944*yk3+ 2.8063-0.18187 *yk42 + 20.0371*yk4+ 6.7665-0.20381*yk52 + 22.8783*yk5+ 7.8851-0.21575*yk62 + 23.7317*yk6+ 7.5642-0.2396*yk72 + 22.7827*yk7+ 15.4914-0.25671*yk82 + 24.41*yk8+ 16.598-0.27125 *yk92 + 25.7934*yk9+ 17.5138-0.2784 *yk102 + 26.4602*yk10+ 17.7686);!计划外成本b3;b3=-0.0603*z12 + 7.5807*z1+ 2.4845-0.067029*z22 + 8.2675*z2+ 2.7346-0.068437*z32 + 8.5944*z3+ 2.8063-0.18187 *z42 + 20.0371*z4+ 6.7665-0.20381*z52 + 22.8783*z5+ 7.8851-0.21575*z62 + 23.7317*z6+ 7.5642-0.2396*z72 + 22.7827*z7+ 15.4914-0.25671*z82 + 24.41*z8+ 16.598-0.27125 *z92 + 25.7934*z9+ 17.5138-0.2784 *z102 + 26.4602*z10+ 17.7686;!计划内已签约对应经费c1;for(dg(j):c1=-0.0057999*xj12 + 0.71783*xj1+ 2.2137-0.0065679*xj22 + 0.78584*xj2+ 2.4092-0.0066856*xj32+0.81785*xj3+2.4855-0.017717*xj42+1.8979*xj4+5.8632-0.01979*xj52+2.1646*xj5+6.7475-0.021084*xj62 +2.251*xj6+ 6.8553-0.023093*xj72 + 2.1299*xj7+ 7.8747-0.024779*xj82+2.2839*xj8+8.4219-0.02616*xj92+2.4125*xj9+8.9007-0.026904*xj102 + 2.4774*xj10+ 9.0828);!计划内意向签约对应经费c2;for(yx(k):c2=-0.0057999*yk12+0.71783*yk1 + 2.2137-0.0065679*yk22 + 0.78584*yk2 +2.4092-0.0066856*yk32+0.81785*yk3+2.4855-0.017717*yk42+1.8979*yk4+5.8632-0.01979*yk52+2.1646*yk5+6.7475-0.021084*yk62+2.251*yk6+6.8553-0.023093*yk72+2.1299*yk7+7.8747-0.024779*yk82+2.2839*yk8+8.4219-0.02616*yk92+2.4125*yk9+8.9007-0.026904*yk102 + 2.4774*yk10+ 9.0828);x11=4.5;x31=3.5;x12=8;x22=5;x32=4.5;x42=2.5;x13=1;x23=3.5;x33=2;x43=3.5;x44=2;x54=3.5;x64=4.5;x65=4;x75=6.5;x85=9.5;x56=4;x66=2.5;x76=3.5;x86=5;x17=3.5;x37=4.5;x57=3;x18=5;x38=4.5;x58=3.5;x68=4;y111.5;y211.5;y121;y222;y131;y231;y343;y442;y351.5;y453.5;y361;y461;y171;y271;y572;y182.5;y284;y583.5;y192;y291;y2102;y5102;z14;z24;z34;z44;z54;z63;z74;z85;z93;z103;for(link(j,k):xj1+xj2+xj3+yk1+yk2+yk3+z1+z2+z350);for(link(j,k):xj4+xj5+xj6+yk4+yk5+yk6+z4+z5+z665);for(link(j,k):xj7+xj8+xj9+xj10+yk7+yk8+yk9+yk10+z7+z8+z9+z1062;end附件7：model:sets:cp/1.10/;dg/1.8/;yx/1.5/;link(dg,yx);endsetsmax=a1+a2+b-c1-c2-c3-(-0.075559*z12+1.6933*z1+ 0.75536-0.076346 *z22 + 1.7398*z2+ 1.0628 -0.059056 *z32+1.5629*z3 + 1.5684-0.20613*z42 + 4.3918*z4 + 2.3385-0.22296*z52 + 4.892*z5 + 3.1453-0.22629*z62 + 5.0031*z6+ 3.3284 -0.23018*z72 + 5.0769*z7 + 3.3259 -0.25493*z82 + 5.4943*z8 + 3.4316-0.28777*z92 + 5.9454*z9+ 3.4438-0.28056*z102 + 5.9462*z10+ 3.7885);!计划内已签约经费a1;for(dg(j):a1=-0.0057999*xj12 + 0.71783*xj1+ 2.2137-0.0065679*xj22 + 0.78584*xj2+ 2.4092-0.0066856*xj32+0.81785*xj3+2.4855-0.017717*xj42+1.8979*xj4+5.8632-0.01979*xj52 + 2.1646*xj5+ 6.7475 -0.021084*xj62 + 2.251*xj6+ 6.8553-0.023093*xj72 + 2.1299*xj7+ 7.8747-0.024779*xj82+2.2839*xj8+8.4219-0.02616*xj92+2.4125*xj9+8.9007-0.026904*xj102 + 2.4774*xj10+ 9.0828);!计划内意向签约经费a2;for(yx(k):a2=-0.0057999*yk12+0.71783*yk1+2.2137-0.0065679*yk22 + 0.78584*yk2+ 2.4092-0.0066856*yk32+0.81785*yk3+2.4855-0.017717*yk42+1.8979*yk4+5.8632-0.01979*yk52 + 2.1646*yk5+ 6.7475 -0.021084*yk62 + 2.251*yk6+ 6.8553-0.023093*yk72 + 2.1299*yk7+7.8747-0.024779*yk82+2.2839*yk8+ 8.4219-0.02616*yk92 + 2.4125*yk9+ 8.9007-0.026904*yk102 + 2.4774*yk10+ 9.0828);!计划外销售额b;b=-0.26713 *z12 + 11.4183 *z1 + 1.3873-0.27436* z22 + 12.1778* z2+ 2.0648-0.26305 *z32 + 12.4521* z3+ 2.5588-0.73776*z42 + 29.8586*z4 + 5.3673-0.85117 *z52 + 34.4913*z5 + 5.788 -0.83007*z62 + 34.9027* z6+ 6.9564-0.86107*z72 + 35.2201*z7 + 8.2703-0.94988*z82 + 38.0877 *z8 + 8.0776-0.93555*z92 + 38.9202*z9 + 11.3939-1.0184*z102 + 41.1537*z10+ 8.7158;!计划外宣传费c1;c1=-0.0050148*z12 + 0.62783*z1+ 0.46871-0.0055698 *z22 + 0.68425*z2 + 0.51762-0.0056906* z32 + 0.71167 *z3 + 0.53081-0.015133 *z42 + 1.6588*z4+ 1.2577 -0.016938 *z52 + 1.8931*z5+ 1.4585-0.017923*z62 + 1.9636 *z6+ 1.4544-0.019882 *z72 + 1.8817*z7+ 2.1416-0.021309 *z82 + 2.0163*z8+ 2.2941-0.022505 *z92 + 2.1301*z9+ 2.4257-0.02312 *z102 + 2.1862*z10+ 2.4633;!计划内已签约宣传费c2;for(dg(j):c2=-0.0050148*xj12 + 0.62783*xj1+ 0.46871-0.0055698 *xj22 + 0.68425*xj2 + 0.51762-0.0056906* xj32 + 0.71167 *xj3 + 0.53081-0.015133 *xj42 + 1.6588*xj4+ 1.2577 -0.016938 *xj52 + 1.8931*xj5+ 1.4585-0.017923*xj62 + 1.9636 *xj6+ 1.4544-0.019882 *xj72 + 1.8817*xj7+ 2.1416-0.021309 *xj82 + 2.0163*xj8+ 2.2941-0.022505 *xj92 + 2.1301*xj9+ 2.4257-0.02312 *xj102 + 2.1862*xj10+ 2.4633);!计划内意向签约宣传费c3;for(yx(k):c3=-0.0050148*yk12 + 0.62783*yk1+ 0.46871-0.0055698 *yk22 + 0.68425*yk2 + 0.51762-0.0056906* yk32 + 0.71167 *yk3 + 0.53081-0.015133 *yk42 + 1.6588*yk4+ 1.2577-0.016938 *yk52 + 1.8931*yk5+ 1.4585-0.017923*yk62 + 1.9636 *yk6+ 1.4544-0.019882 *yk72 + 1