电器销售问题论文

1、PAGE PAGE 3 PAGE 3第1页，共37页2013年高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛

2、规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A我们的参赛报名号为（如果赛区设置报名号的话）：J1109所属学校（请填写完整的全名）：西安科技大学参赛队员 (打印并签名) ：1.闫旭斌 2.马鹏 3.郭卜源指导教师或指导教师组负责人 (打印并签名)：周彬（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修

3、改。如填写错误，论文可能被取消评奖资格。） 日期：2013年9月11日赛区评阅编号（由赛区组委会评阅前进行编号）：2013年高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：PAGE 69 PAGE 69电器销售问题摘要本论文针对大多数企业所面临的生产、销售问题做出了初步的假设分析，主要解决了公司的利润最大化、营销部的收入极大化等问题。在对具体问题的分析上，根据所求目标函数及约束的不同，列出不同的假设变量，结合相应的

4、条件，求解出了在不同要求下的最佳生产和销售方案。 本题中给出了公司的一些基本情况，包括公司生产能力，签约量（已签和意向），营销部的工作行为，发放的经费，客户的情况等。公司必须先满足已签约量，完成生产，由于生产能力的限制，对意向签约量须有选择的安排生产，以达到所得利润最大。营销部作为公司的销售主力，须努力提高签约量以尽最大可能的销售公司产品。首先，对表格中数据进行处理。通过对表3，表4等相应表格中的数据进行拟合，可以得到相关数据的函数关系，建立起数据之间的数学联系。对于第一，二，三，五问,通过对题目认真的分析,可以建立相应的非线性规划模型,利用模型求解即可得到其在相应要求下的最佳生产和销售方案。

5、对于问题四，为了兼顾公司、营销部的利益以及客户的需求，尽量做到均衡销售。首先定义客户的不满意度，并以其相反数作为满意度，同时乘以一个权重分子,使其可以进行不同量化指标之间的加和。从而使多目标决策转化为单目标决策,之后求解出了此条件下的最佳生产和销售方案。关键字：利润最大化非线性规划最优问题 期望值 满意度 一、问题重述本题所研究的是某电器公司根据企业的生产能力、已经签订的销售合同、意向签约量、计划外的销售量等制定相应的生产、销售方案。由于企业生产能力的制约，因此需要在满足已签约的销售合同量的基础上，完成生产，对意向签约量有选择的安排生产，使利润达到最大。营销部作为公司产品销售的主力军应尽可能的

6、增加签订量，提高产品的销售量为公司盈利。请分别按以下要求为该电器公司制定相应的生产、销售方案：（1）使公司的利润达到最大；（2）使营销部的总收入极大化；（3）兼顾公司和营销部二者的利益；（4）兼顾公司、营销部的利益以及客户的需求，尽量做到均衡销售；（5）公司为了调动营销部的积极性，为公司多创利润，规定对于计划外销售的产品，营销部可以自行定价。统计数据显示，自行定价后每单位产品的销售价格上涨或下降6，则三类小家电的可能销售量相应的减少或增加11。试确定使营销部总收入最大的定价、生产及销售方案。二、问题分析该题是企业在自身能力限制下，通过制定生产、销售方案使各方的利润都能够达到预期的目标，即一定约

7、束条件下的最优化问题，初步分析题意后可知所求目标函数是非线性的，是一个非线性优化问题。在对题目分析时，应理解企业一些运作情况。首先，公司是必须先完成已签约量之后，再对意向签约量有选择的安排生产。企业要达到最大利润，得选择盈利最多的产品进行生产并销售。除了之外，企业要多生产一些产品（计划外），让销售部额外销售。在产品的成本问题上，企业有自己的打算。一般企业单位产品成本会随着产量的增多而降低，故企业总希望尽量多的生产产品以减小生产成本。销售方面，单位产品的价格会随着订购量的增加而减小。在此题中，营销部的任务是尽量争取与意向签约的客户签订正式合约，还需努力多销售一些产品。对于所签约的销售合同和意向签

8、约量（计划内），企业根据销售量向营销部发放经费（包括工资以及宣传费用等）；对于计划外销售的产品，实行承包制，营销部向企业缴纳利润。这里要说明的是计划外的产品的销售额利润属于营销部，营销部上缴一定利润，要承担产品的宣传费用，但一般计划外卖出的产品利润较多，故营销部是赢利的，所以销售部希望在计划外能多卖出产品。对于意向签约订购合同，这里在问题一、二、三中是以最大的合同量计算求得利润然后再问题四中定义了不满意度，然后将不满意度和利润的单位去量纲化，通过一定的赋值找到最佳的权重，求得的解就是企业、营销部和客户的需求得到平衡的方案。问题五中营销部可以自行定价，然而产品的销量与产品单价有一定的关系，因此提

9、高或降低价格必然会减少或增加销量，根据题目中的约束建立模型。三、问题假设1、假设销售量、订购量、生产量相等；2、假设企业生产的产品都是合格的；3、假设企业内部的各种因素稳定，不会有资金链中断等情况；4、假设不满意度是由于生产量与意向订购量不匹配造成的；5、假设价格上涨或下降6%必然导致销售量相应的减少或增加11%；四、符号说明序号符号符号说明1产品的种类2第种产品订购量3第种产品产量4计划内第种产品销售量5计划外第种产品销售量6第种产品生产成本7第种产品生产成本8计划内产品宣传费9计划内企业发放的宣传费10计划外向企业缴纳的利润11计划外销售量12第种家电的定价13计划内的销售总额14第个客户

10、对第种产品的订购量期望15第个客户对第种可能出现的订购量16第个客户对第种产品的第种可能出现的订购量订购时的概率17不满意度18权重因子分别为1,10，20，30,60,10019第种意向签约量20第种计划外的销售量五、模型的建立与求解5.1数据的预处理5.1.1.1利用spss处理数据：通过用spss软件作散点图发现与存在线性关系，再利用多种函数拟合，比较分析，三次函数拟合的精度最高， 能够解释绝大部分数据。函数关系式：其中，为常数，、分别为的参数。1）、每种家电产量（千件）与成本（万元）的函数关系式的系数如表1所示：表1：产量与成本函数关系式系数系数家 电家电11.8547.79-0.07

11、40家电21.4078.709-0.0950家电32.0948.831-0.0830家电45.73220.381-0.2040家电58.0222.833-0.201-4.63E-05家电66.23624.174-0.2440家电7-3.4529.085-0.6370.006家电8-3.69731.163-0.6820.007家电9-3.92532.927-0.7210.007家电10-4.28233.797-0.7410.0082）、每种小家电的宣传费（万元）随销售量（千件）的函数关系式系数如表2所示：表2：宣传费随销售量的函数关系式系数系数家 电家电10.3740.659-0.0073.25

12、E-05家电20.360.737-0.0095.40E-05家电30.4290.746-0.0083.50E-05家电41.0661.723-0.0196.58E-05家电51.3491.93-0.0193.75E-05家电61.2162.043-0.0238.17E-05家电70.4252.453-0.0560.001家电80.452.63-0.060.001家电90.4822.777-0.0630.001家电100.4612.853-0.0650.0013）、每种小家电的销售额（万元）随订购量（千件）的函数关系系数如表3所示：表3：销售额随订购量的函数关系式系数系数家 电家电11.3581

13、1.443-0.2720家电21.21812.877-0.4140.008家电30.75613.941-0.560.017家电4-2.0836.01-1.9650.068家电5-0.59539.756-1.9020.058家电6-1.59441.966-2.240.078家电7-3.14444.648-2.7430.105家电8-2.68346.976-2.7240.099家电9-9.02955.79-4.3020.187家电10-2.7350.608-2.9050.1055.1.1.2利用excel处理数据通过用excel作出平滑线散点图，发现其具有一定的线性关系，再添加多种函数趋势线，分析

14、比较，三次函数的拟合优度超过99%，解释数据能力强。其函数关系式：其中，为常数，、分别为的参数。4）、企业向营销部发放计划内的销售量（千件）和经费（万元）函数关系系数如表4所示：表4：销售量和经费的函数关系式系数系数家 电家电11.83980.8422-0.01360.0001家电21.94510.9403-0.01630.0002家电32.07660.9539-0.01530.0001家电44.97822.1923-0.03630.0003家电55.85762.4607-0.03850.0003家电65.77522.6104-0.04370.0004家电74.95063.1029-0.084

15、40.001家电85.28473.3277-0.09060.0011家电95.59893.5112-0.09540.0011家电105.67613.611-0.09840.00125)、计划外营销部向企业缴纳利润(万元)和销售量（千件）函数关系系数如表5所示:表5：计划外缴纳利润和销售量函数关系系数系数家 电家电1-0.0142.3288-0.20240.007家电2-0.03462.6463-0.25720.01家电3-0.49153.2644-0.39860.0189家电4-0.86327.0365-0.73390.0293家电5-0.8678.2061-0.88430.0367家电6-0

16、.53518.1944-0.86310.0354家电7-0.9018.5683-0.92690.0387家电8-1.3159.4151-1.03730.0435家电9-1.42949.9706-1.0910.0446家电10-1.651110.439-1.17720.04985.1.1.3意向客户的需求量和计划外的需求量按最大处理表6：10种小家电的“已签约合同量”以及“意向签约量”表序 号12345678910家电名称1.5升热水壶1.8升热水壶2升热水壶0.9升豆浆机1.1升豆浆机1.3升豆浆机2升电饭煲25升电饭煲3升电饭煲3.5升电饭煲已签约量（万个）0.82.01.01.02.01.

17、51.11.71.00.5意向签约量（万个）0.30.30.20.50.50.20.41.00.30.4计划外最大签约量（万个）0.40.40.40.40.40.30.40.50.30.35.1.2模型一的建立及求解：根据题目分析企业的利润构成如下所示：企业利润=计划内销售额+计划外上缴利润-计划内经费-总成本由此，建立企业利润最大的非线性规划模型：其中，表示计划内销售总额，计划外上缴利润，表示生产总成本，计划内发放经费，热水壶、豆浆机、电饭煲的各自生产总量不能超过最大生产能力，计划内、计划外的销售量不能超过最大签约量，生产量、销售量、订购量均为整数。运用lingo软件编程求解得企业的最大利润

18、为：1129.233（万元），生产和销售方案如下表所示：表7：企业最大利润生产和销售方案 产量家电计划内已签约客户实际生产量（千件）意向签约客户实际生产量（千件）计划外客户实际生产量（千件）180022030310004100052050615207114081710091030105005.1.3模型二的建立及求解：根据题目分析企业的利润构成如下所示：销售部总收入=计划外的销售额-计划外上缴利润-计划外的宣传费由此，建立销售部总收入极大化规划模型：其中，表示计划外的销售额，是销售总额，是计划内销售额，计划外上缴利润，表示计划外宣传费，是总的宣传费，计划内的宣传费，热水壶、豆浆机、电饭煲的各自

19、生产总量不能超过最大生产能力，计划内、计划外的销售量不能超过最大签约量，生产量、销售量、订购量均为整数。运用lingo软件编程（附件二）求解得销售部极大利润为：953.2921（万元），生产和销售方案如下表所示：表8：营销部总输入极大化生产和销售方案 产量家电计划内已签约客户实际生产量（千件）意向签约客户实际生产量（千件）计划外客户实际生产量（千件）实际生产总量（千件）1804122200424310041441024165205429615232071100118171053291003131050165.1.4模型三的建立及求解：根据模型一、二中企业利润和销售利润，为了兼顾二者的利益分配，

20、建立二者利益平衡的规划模型：其中，表示销售总额，表示计划内发放经费，表示计划外宣传费，表示总的宣传费，表示计划内的宣传费，热水壶、豆浆机、电饭煲的各自生产总量不能超过最大生产能力，计划内、计划外的销售量不能超过最大签约量，生产量、销售量、订购量均为整数。运用lingo软件编程（附件三）求解得企业和销售部总和的最大利润为：1633.415（万元），生产和销售方案如下表所示：表9：兼顾企业和营销部二者利益的生产销售方案 产量家电计划内已签约客户实际生产量（千件）意向签约客户实际生产量（千件）计划外客户实际生产量（千件）实际生产总量（千件）180082203427310041441000105205

21、429615232071104158171053291000101050055.1.5模型四的建立及求解：为了兼顾公司、营销部的利益及客户的需求，尽量做到均衡销售。1）、首先将计划外可能订购合同作为意向客户6处理，如下表10所示：表10：意向客户6签订订购合同量的概率表产品合同量1234567891000.10.10.20.10.10.40.10.050.20.20.10.40.30.10.30.30.30.10.10.50.20.20.20.40.20.30.40.20.30.250.20.50.30.20.10.40.20.10.10.40.30.10.10.40.10.10.10.10.

22、10.10.20.50.12）、再次定义客户的不满意度为：其中，表示6个客户所能得到的第种产品的生产量之和，表示第个客户对第种产品的订购量期望，表示6个客户对第种产品订购量的期望之和，表示第个客户对第种可能出现的订购量，表示第个客户对第种产品的第种可能出现的订购量订购时的概率。3）、由此建立以下规划模型:其中，为不满意度的权重因子，表示销售总额，表示计划内发放经费，表示计划外宣传费，表示总的宣传费，表示计划内的宣传费，热水壶、豆浆机、电饭煲的各自生产总量不能超过最大生产能力，计划内、计划外的销售量不能超过最大签约量，生产量、销售量、订购量均为整数。运用lingo软件编程求解得到满足企业和销售部

23、总以及客户需求的均衡生产、销售方案。当时，生产和销售方案如下表11所示：表11：权重因子为1时的生产、销售方案 产量家电意向签约客户生产量（千件）计划外生产量（千件）实际生产总量（千件）108.80E-028.08798200.874499320.8745300.184750110.18475400105542962320700118105329131410005当分别赋予时，求解的生产和销售方案见附录五。通过对赋不同的值，来均衡公司、营销部的利益及客户的需求，通过观察分析比较不同的赋值对计算结果的影响，认为赋值为20时能较为恰当的均衡公司、营销部的利益及客户的需求。5.1.6模型五的建立及求

24、解：为了调动营销部门的积极性，且保证为公司多创造利润，根据每单位产品的销售价格上涨或下降6，则三类小家电的可能销售量相应的减少或增加11的规律，建立营销部确立商品价格与公司利益，营销部门利益的平衡规划模型：其中，为计划外的销售总额，则有，为第二问中求得的销售价格未改变时的销售量，计划外上缴利润，表示计划外宣传费，为该问所需求得的定价，为根据题意计算所得的销售量。运用LINGO软件编程（附件5），最大利润为：2122.570（万元）。生产和销售方案如下表12所示：表12：营销部最大定价、生产及销售方案 类型家电最大定价（万元）生产量（千件）销售量（千件）130.244332.0658612.06

25、58612141.03240.0000000.000000335.307712.0441842.044184486.464541.9738261.9738265246.00771.9876071.9876076147.76571.4744141.47441470.0000000.0000000.0000008309.92382.4900892.490089939.612501.3441381.3441381061.181650.44866020.4486602六、模型的评价与改进本文在对企业各种数据进行初步分析并假设后，建立此模型，运用lingo软件编程求解出符合实际的结果。6.1.1模型的优

26、点如下：1）、本题模型实用性很强，模型中应用文字以及表格表示解释了较为抽象的数学语言，使求解过程简单、易懂，容易推广应用到不同的公司的利润分配计算问题。2）、本题模型中数据分析和求解都经过充分的考虑，拟合函数的精度较高，从而提高结果的准确度；3)、本题模型思路清晰明了，逐步进行分析假设，有层次感，便于理解。6.1.2本模型的不足之处：本模型对实际问题的考虑略有不足，对市场中存在的不稳定因素没有考虑， 意向用户的概率以及价格变化对周围人群的影响变化不能予以更好的考虑，在实际问题中还有待改进，在函数拟合上本题采用三阶拟合以达到精度，实际问题中可能有更精确的拟合函数，还有待改进。七、参考文献1邓维斌

27、 唐兴艳 胡大权等,SPSS19(中文版)统计分析实用教程M，北京：电子工业出版社，2012.32庄君 蒋敏杰 李秀霞,EXCEL统计分析与应用M，北京：电子工业出版社，2013.13刘焕彬 库在强 廖小勇等.数学模型与实验M，北京：科学出版社，2008八、附录8.1.1模型一程序：model:sets:s1/1.4/;s2/1.10/;dd/1.10/:x;b(s2,s1):b11,b22,b33,b44,b55;endsetsdata:x=8 20101020151117105;b11=1.358 11.443-0.27201.218 12.877-0.4140.0080.756 13.9

28、41-0.56 0.017-2.08 36.01 -1.9650.068-0.59539.756-1.9020.058-1.59441.966-2.24 0.078-3.14444.648-2.7430.105-2.68346.976-2.7240.099-9.02955.79 -4.3020.187-2.73 50.608-2.9050.105;b22=-0.0142.3288-0.20240.007-0.03462.6463-0.25720.01-0.49153.2644-0.39860.0189-0.86327.0365-0.73390.0293-0.8678.2061-0.88430.

29、0367-0.53518.1944-0.86310.0354-0.9018.5683-0.92690.0387-1.3159.4151-1.03730.0435-1.42949.9706-1.0910.0446-1.651110.439-1.17720.0498;b33=1.854 7.79 -0.07401.407 8.709 -0.09502.094 8.831 -0.08305.732 20.381-0.20408.02 22.833-0.201-0.00004636.236 24.174-0.2440-3.45 29.085-0.6370.006-3.69731.163-0.6820.

30、007-3.92532.927-0.7210.007-4.28233.797-0.7410.008;b44=1.83980.8422-0.01360.00011.94510.9403-0.01630.00022.07660.9539-0.01530.00014.97822.1923-0.03630.00035.85762.4607-0.03850.00035.77522.6104-0.04370.00044.95063.1029-0.08440.0015.28473.3277-0.09060.00115.59893.5112-0.09540.00115.67613.611 -0.09840.0

31、012;enddatamax=(b11(1,1)+b11(1,2)*k1+b11(1,3)*k12+b11(1,4)*k13)+(b11(2,1)+b11(2,2)*k2+b11(2,3)*k22+b11(2,4)*k23)+(b11(3,1)+b11(3,2)*k3+b11(3,3)*k32+b11(3,4)*k33)+(b11(4,1)+b11(4,2)*k4+b11(4,3)*k42+b11(4,4)*k43)+(b11(5,1)+b11(5,2)*k5+b11(5,3)*k52+b11(5,4)*k53)+(b11(6,1)+b11(6,2)*k6+b11(6,3)*k62+b11(6

32、,4)*k63)+(b11(7,1)+b11(7,2)*k7+b11(7,3)*k72+b11(7,4)*k73)+(b11(8,1)+b11(8,2)*k8+b11(8,3)*k82+b11(8,4)*k83)+(b11(9,1)+b11(9,2)*k9+b11(9,3)*k92+b11(9,4)*k93)+(b11(10,1)+b11(10,2)*k10+b11(10,3)*k102+b11(10,4)*k103)+(b22(1,1)+b22(1,2)*c1+b22(1,3)*c12+b22(1,4)*c13)+(b22(2,1)+b22(2,2)*c2+b22(2,3)*c22+b22(

33、2,4)*c23)+(b22(3,1)+b22(3,2)*c3+b22(3,3)*c32+b22(3,4)*c33)+(b22(4,1)+b22(4,2)*c4+b22(4,3)*c42+b22(4,4)*c43)+(b22(5,1)+b22(5,2)*c5+b22(5,3)*c52+b22(5,4)*c53)+(b22(6,1)+b22(6,2)*c6+b22(6,3)*c62+b22(6,4)*c63)+(b22(7,1)+b22(7,2)*c7+b22(7,3)*c72+b22(7,4)*c73)+(b22(8,1)+b22(8,2)*c8+b22(8,3)*c82+b22(8,4)*c

34、83)+(b22(9,1)+b22(9,2)*c9+b22(9,3)*c92+b22(9,4)*c93)+(b22(10,1)+b22(10,2)*c10+b22(10,3)*c102+b22(10,4)*c103)-(b33(1,1)+b33(1,2)*n1+b33(1,3)*n12+b33(1,4)*n13)+(b33(2,1)+b33(2,2)*n2+b33(2,3)*n22+b33(2,4)*n23)+(b33(3,1)+b33(3,2)*n3+b33(3,3)*n32+b33(3,4)*n33)+(b33(4,1)+b33(4,2)*n4+b33(4,3)*n42+b33(4,4)*

35、n43)+(b33(5,1)+b33(5,2)*n5+b33(5,3)*n52+b33(5,4)*n53)+(b33(6,1)+b33(6,2)*n6+b33(6,3)*n62+b33(6,4)*n63)+(b33(7,1)+b33(7,2)*n7+b33(7,3)*n72+b33(7,4)*n73)+(b33(8,1)+b33(8,2)*n8+b33(8,3)*n82+b33(8,4)*n83)+(b33(9,1)+b33(9,2)*n9+b33(9,3)*n92+b33(9,4)*n93)+(b33(10,1)+b33(10,2)*n10+b33(10,3)*n102+b33(10,4)*

36、n103)-(b44(1,1)+b44(1,2)*k1+b44(1,3)*k12+b44(1,4)*k13)+(b44(2,1)+b44(2,2)*k2+b44(2,3)*k22+b44(2,4)*k23)+(b44(3,1)+b44(3,2)*k3+b44(3,3)*k32+b44(3,4)*k33)+(b44(4,1)+b44(4,2)*k4+b44(4,3)*k42+b44(4,4)*k43)+(b44(5,1)+b44(5,2)*k5+b44(5,3)*k52+b44(5,4)*k53)+(b44(6,1)+b44(6,2)*k6+b44(6,3)*k62+b44(6,4)*k63)+

37、(b44(7,1)+b44(7,2)*k7+b44(7,3)*k72+b44(7,4)*k73)+(b44(8,1)+b44(8,2)*k8+b44(8,3)*k82+b44(8,4)*k83)+(b44(9,1)+b44(9,2)*k9+b44(9,3)*k92+b44(9,4)*k93)+(b44(10,1)+b44(10,2)*k10+b44(10,3)*k102+b44(10,4)*k103);a13;a23;a32;a45;a55;a62;a74;a810;a93;a104;c14;c24;c34;c44;c54;c63;c74;c85;c93;c10=0;n1+n2+n3=0;n4

38、+n5+n6=0;n7+n8+n9+n10=62;end8.1.2模型二程序：model:sets:s1/1.4/;s2/1.10/;dd/1.10/:x;b(s2,s1):b11,b22,b33,b44,b55,b66;endsetsdata:x=8 20101020151117105;b11=1.358 11.443-0.27201.218 12.877-0.4140.0080.756 13.941-0.56 0.017-2.08 36.01 -1.9650.068-0.59539.756-1.9020.058-1.59441.966-2.24 0.078-3.14444.648-2.74

39、30.105-2.68346.976-2.7240.099-9.02955.79 -4.3020.187-2.73 50.608-2.9050.105;!t;b22=-0.0142.3288-0.20240.007-0.03462.6463-0.25720.01-0.49153.2644-0.39860.0189-0.86327.0365-0.73390.0293-0.8678.2061-0.88430.0367-0.53518.1944-0.86310.0354-0.9018.5683-0.92690.0387-1.3159.4151-1.03730.0435-1.42949.9706-1.

40、0910.0446-1.651110.439-1.17720.0498;!p;b33=1.854 7.79 -0.07401.407 8.709 -0.09502.094 8.831 -0.08305.732 20.381-0.20408.02 22.833-0.201-0.00004636.236 24.174-0.2440-3.45 29.085-0.6370.006-3.69731.163-0.6820.007-3.92532.927-0.7210.007-4.28233.797-0.7410.008;!f;b44=1.83980.8422-0.01360.00011.94510.940

41、3-0.01630.00022.07660.9539-0.01530.00014.97822.1923-0.03630.00035.85762.4607-0.03850.00035.77522.6104-0.04370.00044.95063.1029-0.08440.0015.28473.3277-0.09060.00115.59893.5112-0.09540.00115.67613.611 -0.09840.0012;!q;b55=-0.0142.3288-0.20240.007-0.03462.6463-0.25720.01-0.49153.2644-0.39860.0189-0.86

42、327.0365-0.73390.0293-0.8678.2061-0.88430.0367-0.53518.1944-0.86310.0354-0.9018.5683-0.92690.0387-1.3159.4151-1.03730.0435-1.42949.9706-1.0910.0446-1.651110.439-1.17720.0498;!p;b66=0.3740.659-0.0073.25E-050.360.737-0.0095.40E-050.4290.746-0.0083.50E-051.0661.723-0.0196.58E-051.3491.93-0.0193.75E-051

43、.2162.043-0.0238.17E-050.4252.453-0.0560.0010.452.63-0.060.0010.4822.777-0.0630.0010.4612.853-0.0650.001;!g;enddatamax=(b11(1,1)+b11(1,2)*n1+b11(1,3)*n12+b11(1,4)*n13)+(b11(2,1)+b11(2,2)*n2+b11(2,3)*n22+b11(2,4)*n23)+(b11(3,1)+b11(3,2)*n3+b11(3,3)*n32+b11(3,4)*n33)+(b11(4,1)+b11(4,2)*n4+b11(4,3)*n42

44、+b11(4,4)*n43)+(b11(5,1)+b11(5,2)*n5+b11(5,3)*n52+b11(5,4)*n53)+(b11(6,1)+b11(6,2)*n6+b11(6,3)*n62+b11(6,4)*n63)+(b11(7,1)+b11(7,2)*n7+b11(7,3)*n72+b11(7,4)*n73)+(b11(8,1)+b11(8,2)*n8+b11(8,3)*n82+b11(8,4)*n83)+(b11(9,1)+b11(9,2)*n9+b11(9,3)*n92+b11(9,4)*n93)+(b11(10,1)+b11(10,2)*n10+b11(10,3)*n102+

45、b11(10,4)*n103)-(b11(1,1)+b11(1,2)*k1+b11(1,3)*k12+b11(1,4)*k13)+(b11(2,1)+b11(2,2)*k2+b11(2,3)*k22+b11(2,4)*k23)+(b11(3,1)+b11(3,2)*k3+b11(3,3)*k32+b11(3,4)*k33)+(b11(4,1)+b11(4,2)*k4+b11(4,3)*k42+b11(4,4)*k43)+(b11(5,1)+b11(5,2)*k5+b11(5,3)*k52+b11(5,4)*k53)+(b11(6,1)+b11(6,2)*k6+b11(6,3)*k62+b11(

46、6,4)*k63)+(b11(7,1)+b11(7,2)*k7+b11(7,3)*k72+b11(7,4)*k73)+(b11(8,1)+b11(8,2)*k8+b11(8,3)*k82+b11(8,4)*k83)+(b11(9,1)+b11(9,2)*k9+b11(9,3)*k92+b11(9,4)*k93)+(b11(10,1)+b11(10,2)*k10+b11(10,3)*k102+b11(10,4)*k103)-(b55(1,1)+b55(1,2)*c1+b55(1,3)*c12+b55(1,4)*c13)+(b55(2,1)+b55(2,2)*c2+b55(2,3)*c22+b55

47、(2,4)*c23)+(b55(3,1)+b55(3,2)*c3+b55(3,3)*c32+b55(3,4)*c33)+(b55(4,1)+b55(4,2)*c4+b55(4,3)*c42+b55(4,4)*c43)+(b55(5,1)+b55(5,2)*c5+b55(5,3)*c52+b55(5,4)*c53)+(b55(6,1)+b55(6,2)*c6+b55(6,3)*c62+b55(6,4)*c63)+(b55(7,1)+b55(7,2)*c7+b55(7,3)*c72+b55(7,4)*c73)+(b55(8,1)+b55(8,2)*c8+b55(8,3)*c82+b55(8,4)*

48、c83)+(b55(9,1)+b55(9,2)*c9+b55(9,3)*c92+b55(9,4)*c93)+(b55(10,1)+b55(10,2)*c10+b55(10,3)*c102+b55(10,4)*c103)-(b66(1,1)+b66(1,2)*c1+b66(1,3)*c12+b66(1,4)*c13)+(b66(2,1)+b66(2,2)*c2+b66(2,3)*c22+b66(2,4)*c23)+(b66(3,1)+b66(3,2)*c3+b66(3,3)*c32+b66(3,4)*c33)+(b66(4,1)+b66(4,2)*c4+b66(4,3)*c42+b66(4,4)

49、*c43)+(b66(5,1)+b66(5,2)*c5+b66(5,3)*c52+b66(5,4)*c53)+(b66(6,1)+b66(6,2)*c6+b66(6,3)*c62+b66(6,4)*c63)+(b66(7,1)+b66(7,2)*c7+b66(7,3)*c72+b66(7,4)*c73)+(b66(8,1)+b66(8,2)*c8+b66(8,3)*c82+b66(8,4)*c83)+(b66(9,1)+b66(9,2)*c9+b66(9,3)*c92+b66(9,4)*c93)+(b66(10,1)+b66(10,2)*c10+b66(10,3)*c102+b66(10,4)

50、*c103)-(b66(1,1)+b66(1,2)*n1+b66(1,3)*n12+b66(1,4)*n13)+(b66(2,1)+b66(2,2)*n2+b66(2,3)*n22+b66(2,4)*n23)+(b66(3,1)+b66(3,2)*n3+b66(3,3)*n32+b66(3,4)*n33)+(b66(4,1)+b66(4,2)*n4+b66(4,3)*n42+b66(4,4)*n43)+(b66(5,1)+b66(5,2)*n5+b66(5,3)*n52+b66(5,4)*n53)+(b66(6,1)+b66(6,2)*n6+b66(6,3)*n62+b66(6,4)*n63)

51、+(b66(7,1)+b66(7,2)*n7+b66(7,3)*n72+b66(7,4)*n73)+(b66(8,1)+b66(8,2)*n8+b66(8,3)*n82+b66(8,4)*n83)+(b66(9,1)+b66(9,2)*n9+b66(9,3)*n92+b66(9,4)*n93)+(b66(10,1)+b66(10,2)*n10+b66(10,3)*n102+b66(10,4)*n103);a13;a23;a32;a45;a55;a62;a74;a810;a93;a104;c14;c24;c34;c44;c54;c63;c74;c85;c93;c10=0;n1+n2+n3=0;n

52、4+n5+n6=0;n7+n8+n9+n10=62;end8.1.3模型三程序：model:sets:s1/1.4/;s2/1.10/;dd/1.10/:x;b(s2,s1):b11,b22,b33,b44,b55,b66;endsetsdata:x=8 20101020151117105;b11=1.358 11.443-0.27201.218 12.877-0.4140.0080.756 13.941-0.56 0.017-2.08 36.01 -1.9650.068-0.59539.756-1.9020.058-1.59441.966-2.24 0.078-3.14444.648-2.7

53、430.105-2.68346.976-2.7240.099-9.02955.79 -4.3020.187-2.73 50.608-2.9050.105;!t;b22=-0.0142.3288-0.20240.007-0.03462.6463-0.25720.01-0.49153.2644-0.39860.0189-0.86327.0365-0.73390.0293-0.8678.2061-0.88430.0367-0.53518.1944-0.86310.0354-0.9018.5683-0.92690.0387-1.3159.4151-1.03730.0435-1.42949.9706-1

54、.0910.0446-1.651110.439-1.17720.0498;!p;b33=1.854 7.79 -0.07401.407 8.709 -0.09502.094 8.831 -0.08305.732 20.381-0.20408.02 22.833-0.201-0.00004636.236 24.174-0.2440-3.45 29.085-0.6370.006-3.69731.163-0.6820.007-3.92532.927-0.7210.007-4.28233.797-0.7410.008;!f;b44=1.83980.8422-0.01360.00011.94510.94

55、03-0.01630.00022.07660.9539-0.01530.00014.97822.1923-0.03630.00035.85762.4607-0.03850.00035.77522.6104-0.04370.00044.95063.1029-0.08440.0015.28473.3277-0.09060.00115.59893.5112-0.09540.00115.67613.611 -0.09840.0012;!q;b55=-0.0142.3288-0.20240.007-0.03462.6463-0.25720.01-0.49153.2644-0.39860.0189-0.8

56、6327.0365-0.73390.0293-0.8678.2061-0.88430.0367-0.53518.1944-0.86310.0354-0.9018.5683-0.92690.0387-1.3159.4151-1.03730.0435-1.42949.9706-1.0910.0446-1.651110.439-1.17720.0498;!p;b66=0.3740.659-0.0073.25E-050.360.737-0.0095.40E-050.4290.746-0.0083.50E-051.0661.723-0.0196.58E-051.3491.93-0.0193.75E-05

57、1.2162.043-0.0238.17E-050.4252.453-0.0560.0010.452.63-0.060.0010.4822.777-0.0630.0010.4612.853-0.0650.001;!g;enddatamax=(b11(1,1)+b11(1,2)*n1+b11(1,3)*n12+b11(1,4)*n13)+(b11(2,1)+b11(2,2)*n2+b11(2,3)*n22+b11(2,4)*n23)+(b11(3,1)+b11(3,2)*n3+b11(3,3)*n32+b11(3,4)*n33)+(b11(4,1)+b11(4,2)*n4+b11(4,3)*n4

58、2+b11(4,4)*n43)+(b11(5,1)+b11(5,2)*n5+b11(5,3)*n52+b11(5,4)*n53)+(b11(6,1)+b11(6,2)*n6+b11(6,3)*n62+b11(6,4)*n63)+(b11(7,1)+b11(7,2)*n7+b11(7,3)*n72+b11(7,4)*n73)+(b11(8,1)+b11(8,2)*n8+b11(8,3)*n82+b11(8,4)*n83)+(b11(9,1)+b11(9,2)*n9+b11(9,3)*n92+b11(9,4)*n93)+(b11(10,1)+b11(10,2)*n10+b11(10,3)*n102

59、+b11(10,4)*n103)-(b44(1,1)+b44(1,2)*k1+b44(1,3)*k12+b44(1,4)*k13)+(b44(2,1)+b44(2,2)*k2+b44(2,3)*k22+b44(2,4)*k23)+(b44(3,1)+b44(3,2)*k3+b44(3,3)*k32+b44(3,4)*k33)+(b44(4,1)+b44(4,2)*k4+b44(4,3)*k42+b44(4,4)*k43)+(b44(5,1)+b44(5,2)*k5+b44(5,3)*k52+b44(5,4)*k53)+(b44(6,1)+b44(6,2)*k6+b44(6,3)*k62+b44

60、(6,4)*k63)+(b44(7,1)+b44(7,2)*k7+b44(7,3)*k72+b44(7,4)*k73)+(b44(8,1)+b44(8,2)*k8+b44(8,3)*k82+b44(8,4)*k83)+(b44(9,1)+b44(9,2)*k9+b44(9,3)*k92+b44(9,4)*k93)+(b44(10,1)+b44(10,2)*k10+b44(10,3)*k102+b44(10,4)*k103)-(b66(1,1)+b66(1,2)*n1+b66(1,3)*n12+b66(1,4)*n13)+(b66(2,1)+b66(2,2)*n2+b66(2,3)*n22+b6