广告费用与销售价格调整程序设计

数学建模实践数学建模模课程设设计（程序设设计和论论文）题目 11对函数数进行麦麦克劳林林展开及及误差分分析 2 无变位位油罐中中油量确确定及误误差分析析 3评评卷成绩绩调整程程序设计计4 广告告费用与与销售价价格调整整程序设设计班级学号学生姓名名指导教师师60沈阳航空空航天大大学课程设计计任务书课程名称称数学建建模实践践院（系）理学院专业信息与计算科学班级学号号姓名课程设计计题目 1 对对函数进进行麦克克劳林展展开及误误差分析析 2 无无变位油油罐中油油量确定定及误差差分析 3 评评卷成绩绩调整程程序设计计 4 广广告费用用与销售售价格的的调整程程序设计计课程设计计时间: 220111 年年 6 月 277 日至至 20011 年 7 月 115 日日课程设计计的内容容及要求求：内容1（11）求函函数（2）编编写对任任意固定定的n计算多多项式函函数值的的函数MM文件（3）任任取n，在同同一平面面内画出出函数的的图形，并并进行比比较。2无变变位油罐罐中油量量确定设油罐中中油量VV与高度度h的关系系是其中，（1）编编写计算算体积VV(h)的函数数M文件fvv；（2）根根据“无无变位实实验采集集数据表表”中的的无变位位进油表表中的数数据计算算公式VV(h)与实验验数据之之间的误误差WCC(h)，并用用多项式式拟合确确定函数数WC(h)表达式式。（3）用用误差WWC(hh)调整V(h)，并并用“无无变位实实验采集集数据表表”中的的无变位位出油表表中的数数据检验验调整结结果。3评卷卷成绩调调整程序序设计设个专家家分别对对名学生生的试卷卷进行评评阅，设设表示教教师对学学生的试试卷所给给定的成成绩，这这样形成成成绩矩矩阵。由由于各专专家的评评分标准准不一致致，因此此需要对对成绩进进行一致致性调整整，具体体方法如如下：设分别表表示整体体成绩，是教师师j的平均均成绩和和标准差差，即第第j列数据据的平均均值和标标准差。调调整后的的成绩为为形成调整整后的成成绩矩阵阵，则的平平均值就就是第ii个学生生的最后后综合成成绩。而而是综合合成绩向向量，依依此确定定学生获获奖情况况。（1）编编写函数数M文件，收收入成绩绩矩阵，输输出是综综合成绩绩向量。（2）根根据下表表是成绩绩数据学生编号号专家1专家2专家3专家419090809729590809038585909048590809059580907069080907578090708588585708099570707510906060601175657060127570606513607060501465605050156560505016606050501760505050用上述方方法计算算综合成成绩向量量，并由由此确定定1个特等等奖，11个一等等奖，22个二等等奖；33个三等等奖。4.广告告模型某装饰材材料公司司欲以每每桶2元的价价钱购进进一批彩彩漆以供供日后销销售。为为了尽快快收回资资金并获获得较多多的赢利利，公司司经理李李先生打打算做广广告，于于是便找找到广告告公司的的王经理理进行咨咨询。李李经理认认为，随随彩漆售售价的提提高，预预期销售售量将减减少，并并对此进进行了估估算（见见表2）。他他问王经经理广告告有多大大的效应应。王经经理说“投投入一定定的广告告费后，销销售量将将有一个个增长，这这由销售售增长因因子来表表示。例例如，投投资3万元的的广告费费，销售售增长因因子为11.855，即销销售量将将是预期期销售量量的1.85倍倍。根据据经验，广广告费与与销售增增长因子子的关系系有表33。”售价2.0002.5003.0003.5004.0004.5005.0005.5006.000预期销售售量（千千桶）413834322928252220表2 售价与预期销售量广告费（元元）0100000200000300000400000500000600000700000销售增长长因子1.0001.4001.7001.8551.9552.0001.9551.800表3 广告费与销售增长因子问李经理理如何确确定彩漆漆的售价价和广告告费，才才能使公公司获得得的利润润最大？ 要求求1、学习习态度要要认真，要要积极参参与课程程设计，锻锻炼独立立思考能能力；2、严格格遵守上上机时间间安排；3、按照照MATTLABB编程训训练的任任务要求求来编写写程序；4、根据据任务来来完成数数学建模模论文；5、报告告书写格格式要求求按照沈沈阳航空空航天大大学“课课程设计计报告撰写写规范”；7、报告告上交时时间：课课程设计计结时上上交报告告。8、严谨谨抄袭行行为。指导教师师年月日负责教师师年月日学生签字字年月日沈阳航空空航天大大学课程设计计成绩评定单课程名称称数学建建模实践践院（系）理学院专业信息与计算科学课程设计计题目11 对函函数进行行麦克劳劳林展开开及误差差分析 2 无无变位油油罐中油油量确定定及误差差分析 3 评评卷成绩绩调整程程序设计计 4 广广告费用用与销售售价格的的调整程程序设计计学号 220099041140110022 姓名名郭婧指导教师师评语：课程设计计成绩指导教师师签字年月日目录目录V摘要1正文31 题目目一31.1 问题重重述31.2 问题求求解31.3 题目结结果42 题目目二52.1 问题重重述52.2 问题求求解62.3 题目结结果1003 题目目三1663.1 问题重重述1663.2 问题求求解1773.3 题目结结果2004 题目目四2114.1 问题重重述2114.2 问题求求解2224.3 题目结结果244参考文献献26源程序227摘要在本次课课程设计计中，我我的课程程设计题题目是四四道题。第第一道题题目里的的第一个个问号是是用Maatlaab编写写函数，根根据人为为设定的的n，函数数可以任任意展开开，并且且在Maatlaab运行行界面显显示的是是展开的的多项式式。第二二个问号号里要求求在任意意设定的的n阶下，带带入自变变量的值值，然后后求出的的麦克劳劳林展式式的函数数值。第第三个问问号里要要求我们们通过画画图对的的本来的的式子、麦麦克劳林林展开的的式子、以以及作比比较。在在题中已已经给定定画图区区间，在在这个区区间内画画出图形形，进行行比较。第二道题题目中，给给出了一一个Exxcell表格，里里面有无无变位进进油量表表和无变变量出油油量表。我我们知道道在一个个油罐中中罐中油油的高度度和体积积是有一一定的关关系的，题题中就把把这种关关系式给给了我们们，式式中的一一些参量量已经给给出，编编写这个个式子的的程序即即可。带带入不同同的高度度可以输输出不同同的体积积。接下下来根据据“无变变位实验验采集数数据表”中中的无变变位进油油表中的的数据，计计算公式式V(h)与实验验数据之之间的误误差WCC(h)。我们们可以先先把体积积数据保保存在MMatllab中中，然后后用表中中已经给给了的高高度带入入V(h)，这时时可以求求出一系系列的体体积，然然后与真真实值进进行作差差，得到到的数据据即为误误差。再再次用多多项式拟拟合确定定误差函函数WCC(h)表达式式。最后后，用误误差函数数WC(h)调整V(h)，并并用“无无变位实实验采集集数据表表”中的的无变位位出油表表中的数数据检验验调整结结果。第三道题题目是评评卷成绩绩调整程程序设计计，题目目中给出出了一些些学生的的由不同同专家给给出的阅阅卷成绩绩。要求求先求出出每一个个学生由由不同专专家给出出的成绩绩的平均均值，然然后求出出标准差差。再求求出第jj个专家家给出成成绩的平平均值，然然后求出出标准差差。调整整后的成成绩为。形形成调整整后的成成绩矩阵阵，则的平平均值就就是第ii个学生生的最后后综合成成绩。而而是综合合成绩向向量，依依此确定定学生获获奖情况况。由此此最后求求出的列列向量确确定1名特等等奖，11名一等等奖，22名二等等奖；33名三等等奖。第四题是是一个广广告模型型，某装装饰材料料公司欲欲以每桶桶2元的价价钱购进进一批彩彩漆以供供日后销销售。为为了可以以很快的的收益并并且收回回大量的的资金，公公司经理理李先生生打算做做广告，于于是便找找到广告告公司的的王经理理进行咨咨询。李李经理认认为，随随彩漆售售价的提提高，预预期销售售量将减减少，并并对此进进行了估估算，见见表格22，随着着销售价价格的增增加，销销售量下下降。他他问王经经理广告告有多大大的效应应。王经经理说：“投入入一定的的广告费费后，销销售量将将有一个个增长，这这由销售售增长因因子来表表示。通通过表格格3可知，随随着广告告费用的的升高，销销售因子子先上升升后下降降。问李李经理如如何确定定彩漆的的售价和和广告费费，才能能使公司司获得的的利润最最大。关键词：拟合函函数；误差分分析；调整矩矩阵；利益最最大化正文1 题目目一1.1 问题重重述第一个问问号是用用Mattlabb编写函函数，根根据人为为设定的的n，函数数可以任任意展开开，并且且在Maatlaab运行行界面显显示的是是展开的的多项式式。第二二个问号号里要求求在任意意设定的的n阶下，带带入自变变量的值值，然后后求出的的麦克劳劳林展式式的函数数值。第第三个问问号里要要求我们们通过画画图对的的本来的的式子、麦麦克劳林林展开的的式子、以以及作比比较。在在题中已已经给定定画图区区间，在在这个区区间内画画出图形形，进行行比较。1.2 问题求求解（1）根根据数学学分析课课程中学学到的麦麦克劳林林展开的的定义，可可知=（）（1）（）（22）（）（33）有了公式式（1）（2）（3）就可以以对编程程进行麦麦克劳林林展开，把把、用Mattlabb语言进进行编辑辑，然后后作差即即得得麦麦克劳林林展式。（2）再再上一个个问号中中已经把把麦克劳劳林展式式求出来来，在第第二个程程序中只只需把任任意自变变量值代代入求函函数值。把把Mattlabb中M文件的的函数名名由fuuncttionn Txx=myyfunn1(nn)改为为funnctiion Tx=myffun22(x,n)，输输入任意意的x和n就可以以求出任任意阶展展式的任任意函数数值。（3）第第三个问问号是画画图比较较，x的区间间已经给给出，用ploot命令令可以直直接画出出图形调调用格式式为a(k)=logg(11-x(k)/(11+x(k); ploot(xx,a,*)。画用用麦克劳劳林展开开的式子子调用格格式为TTn(kk)=mmyfuun2(x(kk),nn); ploot(xx,Tnn,*)。画图形形的调用用格式为为y(kk)=aa(k)-Tnn(k); pplott(x,y,*)。为了了更直观观的观察察图形之之间的关关系和差差距，最最后把三三个图形形画到同同一个图图形中，可可以用ssubpplott(m,n,pp)命令令把四个个画到同同一图中中，分为为四个小小子图，m，n为画几乘几的子图，p为第几个图。1.3 题目结结果（1）第第一个问问号的运运行结果果，n=10时时的的麦麦克劳林林展式。图1.33.1（2）第第一个问问号运行行结果，x=10，n=2时的麦克劳林展式的函数值。图1.33.2（3）第第三个问问号运行行结果，取取n=11时的、麦麦克劳林林展开式式、以及及同时在在一个图图形时的的图形。图1.33.32 题目目二2.1 问题重重述第二道题题目中，给给出了一一个Exxcell表格，里里面有无无变位进进油量表表和无变变量出油油量表。我我们知道道在一个个油罐中中罐中油油的高度度和体积积是有一一定的关关系的，题题中就把把这种VV与h的大致致关系式式给了我我们：，式式中的一一些参量量已经给给出，aa=177.8/2、b=112/22、L1=0.44、L2=2.005，编编写这个个式子的的程序即即可。带带入不同同的高度度可以输输出不同同的体积积。接下下来根据据“无变变位实验验采集数数据表”中中的无变变位进油油表中的的数据，计计算公式式V(h)与实验验数据之之间的误误差WCC(h)。我们们可以先先把体积积数据保保存在MMatllab中中，然后后用表中中已经给给了的高高度带入入V(h)，这时时可以求求出一系系列的体体积，然然后与真真实值进进行作差差，得到到的数据据即为误误差。再再次用多多项式拟拟合确定定误差函函数WCC(h)表达式式。最后后，用误误差函数数WC(h)调整V(h)，并并用“无无变位实实验采集集数据表表”中的的无变位位出油表表中的数数据检验验调整结结果。2.2 问题求求解（1）第第一个问问号中要要求编辑辑计算体体积的公公式的函函数M文件。根根据题目目中给出出的a、b、L1、L2的值值以及公公式V(h)进行编编辑。程程序为：funcctioon VVh=mmyfuun4(h)a=177.8/2;b=122/2;L1=00.4;L2=22.055;Vh=aa*b*(L11+L22)*asiin(h-bb)/bb)+(h-b)/b)*sqrrt(11-(hh-b)2/b22)+ppi/22*110;（2）要要求根据据“无变变位实验验采集数数据表”中中的无变变位进油油表中的的数据计计算公式式V(h)与实验验数据之之间的误误差WCC(h)，并用用多项式式进行拟拟合。误差值=|真实实值-公式求求解的函函数值| （4）所以我用用以下MM文件进进行求解解误差，此此M文件可可以求解解每一个个进油高高度所对对应的误误差。funcctioon mmyfuun5(X1,X2)for l=11:leengtth(XX1) a=117.88/2; b=112/22; L1=0.44; L2=2.005;Vh(ll)=aa*b*(L11+L22)*asiin(X1(l)-b)/b)+(XX1(ll)-bb)/bb)*ssqrtt(1-(X11(l)-b)2/b22)+ppi/22*110; a(ll)=aabs(X2(l)+2622-Vhh(l); dissp(aa(l)end根据高度度和误差差进行曲曲线拟合合，拟合合命令为为：x=X11;y=Y;plott(x,y,*)polyyfitt(x,y,nn)holdd onnfploot(fx,00,122,r)我分别对对曲线进进行了二二次，三三次，四四次，五五次拟合合，得到到以下拟拟合曲线线：二次拟合合曲线：图2.22.1三次拟合合曲线：图2.22.2四次拟合合曲线：图2.22.3五次拟合合曲线：图2.22.4由以上曲曲线拟合合可知：进行三三次和四四次多项项式拟合合的曲线线较好，更更贴合原原图。（3）用用误差WWC(hh)调整V(h)，并并用“无无变位实实验采集集数据表表”中的的无变位位出油表表中的数数据检验验调整结结果。把把编辑的的M文件的的求误差差的语句句中的绝绝对值去去掉后，误误差都为为负值，所所以应在在V(hh)后减减去WCC(h)，所以以，调整整后的VV(h)= VV(h)-WCC(h)。分别别用三次次多项和和四次多多项式进进行求解解误差。2.3 题目结结果（1）hh取10分米米时的结结果：图2.33.1（2）求求解的误误差值：图2.33.2图2.33.3三次拟合合出来的的曲线为为：WCC(h)= -0.0084* +11.50065* +55.82216* -11.71108四次拟合合出来的的曲线为为：WCC(h)= -0.000255*-00.01167* +00.88876*+8.08226*-4.338288（3）用用三次多多项式拟拟合出的的WC(h)调整的的结果：图2.33.4用三次多多项式拟拟合出的的WC(h)图形：图2.33.5用四次多多项式拟拟合出的的WC(h)调整的的结果：图2.33.6用四次多多项式拟拟合出的的WC(h)图形：图2.33.73 题目目三3.1 问题重重述第三道题题目是评评卷成绩绩调整程程序设计计，题目目中给出出了一些些学生的的由不同同专家给给出的阅阅卷成绩绩。要求求先求出出每一个个学生由由不同专专家给出出的成绩绩的平均均值，然然后求出出标准差差。再求求出第jj个专家家给出成成绩的平平均值，然然后求出出标准差差。调整整后的成成绩为。形形成调整整后的成成绩矩阵阵，则的平平均值就就是第ii个学生生的最后后综合成成绩。而而是综合合成绩向向量，依依此确定定学生获获奖情况况。由此此确定11个特等等奖，11个一等等奖，22个二等等奖；33个三等等奖。3.2 问题求求解这道题目目就是对对矩阵进进行变换换，根据据矩阵求求出要求求解出的的数值。如如每个学学生由不不同的专专家评卷卷得出的的成绩的的平均值值，每个学学生成绩绩的标准准差。每每个专家家对不同同学生评评卷成绩绩的平均均值，每每个专家家评卷成成绩的标标准差。求求出这些些必要的的数值后后，根据据题中所所给的公公式求出出Y矩阵，把把矩阵的的每一行行在作和和求平均均值，最最后得到到一个列列向量，列列向量的的每一行行就为最最后评定定学生成成绩的标标准，再再对这个个列向量量的每一一行的数数值进行行有大到到小的排排列，由由此最后后求出的的列向量量就可以以得出得得特等奖奖、一等等奖、二二等奖、三三等奖的的人。标准差的的求解公公式为：（5）平均值的的求解公公式为：（6）求解每个个学生成成绩的平平均值：B=A(:,11);for i=22:4 BB=B+A(:,i);endB=B/4;求解每个个学生成成绩的标标准差：for j=11:4 D(:,j)=(AA(:,j)-B).2;endF=0;for l=11:4 F=FF+D(:,ll);endF=sqqrt(F/44);求解每位位专家评评卷成绩绩的平均均值：C=A(1,:);for j=22:177 C=CC+A(j,:);endC=C/17;求解每位位专家评评卷成绩绩的标准准差：for k=11:177 E(kk,:)=(AA(k,:)-C).2./177;endG=0;for q=11:177 G=GG+E(q,:);endG=sqqrt(G/117);求解最后后的公式式：for p=11:4 H(:,p)=(AA(:,p)-C(pp)./G(p).*F+B;endI=0;for n=11:4 I=II+H(:,nn);endI=I./4;对最后求求出的列列向量进进行排列列并输出出学生号号：a=siize(I);a=a(1);for h=11:a J(hh,1)=h;endfor w=11:166for s=ww:177if II(w)I(s) t=I(ww); I(w)=I(ss); I(s)=t; z=J(ww); J(w)=J(ss); J(s)=z;endendend3.3 题目结结果（1）每每个学生生由不同同的专家家评卷得得出的成成绩的平平均值，每每个学生生成绩的的标准差差。每个个专家对对不同学学生评卷卷成绩的的平均值值，每个个专家评评卷成绩绩的标准准差。学生1234567889.22588.77587.5586.55083.77583.77581.22580.0009101112131415161777.55067.55067.55067.55060.00056.00056.22555.00052.550表3.33.1学生123456786.055705.444862.500004.144589.600146.499527.399516.122379101112131415161710.33078812.9990445.599025.599027.700116.499526.499525.000004.33301表3.33.2专家123479.44118873.88235568.88235569.882355专家12343.088143.122733.411813.95585表3.33.3（2）最最后求出出经过排排序的综综合成绩绩向量。学生15264738117.63332115.61118114.08337104.72001102.1622598.88903397.99249992.7710339111210131617141591.66423358.55699958.44798847.88612233.22867728.44566626.00512224.44039924.440399表3.33.4所以得到到特等奖奖的学生生为学生生1，得到到一等奖奖的学生生为学生生5，得到到二等奖奖的学生生为学生生2和学生生6，得到到三等奖奖的学生生为学生生4，学生生7，学生生3。4 题目目四4.1 问题重重述第四题是是一个广广告模型型，某装装饰材料料公司欲欲以每桶桶2元的价价钱购进进一批彩彩漆以供供日后销销售。为为了可以以很快的的收益并并且收回回大量的的资金，公公司经理理李先生生打算做做广告，于于是便找找到广告告公司的的王经理理进行咨咨询。李李经理认认为，随随彩漆售售价的提提高，预预期销售售量将随随之减少少，并对对此进行行了估算算，销售售价格与与销售量量之间的的关系见见表格22：售价2.0002.5003.0003.5004.0004.5005.0005.5006.000预期销售售量（千千桶）413834322928252220表2 售价与预期销售量通过表格格2知道，销销售量随随着销售售价格的的增加而而减少。他他问王经经理广告告有多大大的效应应。王经经理说：“投入入一定的的广告费费后，销销售量将将有一个个增长，这这有销售售因子决决定，如如投入的的广告费费用为2200000元时时，销售售量将为为原销售售量的11.7倍倍。广告告费用与与销售增增长因子子的关系系见表格格3：广告费（元元）0100000200000300000400000500000600000700000销售增长长因子1.0001.4001.7001.8551.9552.0001.9551.800表3 广告费与销售增长因子通过表格格3可知，随随着广告告费用的的升高，销销售因子子先上升升后下降降。问李李经理如如何确定定彩漆的的售价和和广告费费，才能能使公司司获得的的利润最最大。4.2 问题求求解根据题中中所给要要求可知知，此题题为一个个应用线线性规划划求解的的题目，所所以应使使用Liingoo软件进进行解题题。题中中问如何何确定彩彩漆的售售价和使使用的广广告费用用，才能能使公司司获得的的利益最最大，根根据这个个要求可可以写出出目标函函数。设设x为广告告费用，则则销售增增长因子子为，y为彩漆漆售价，则则为预期期销售量量。根据据表格中中所给的的数据可可以拟合合出和的曲线线。拟合合出来的的=-00.04426*x22+0.40992*xx+1.02888（这这里的xx以万元元为单位位），=-5.13333*xx+500.42222（这这里的以以千桶为为单位）。销售利润润=销售量量*销售因因子*销售价价格-广告费费用(77)目标函数数为：mmax=(-55.13333*y+550.442222)*110000*(-0.004266*x2+00.40096*x+00.01188)*y-x*1100000;因此只要要求出这这个目标标函数的的最大值值就是利利益最大大组合。下面是拟拟合的和和的图像像，根据据图像可可知，拟拟合出的的函数与与实际数数值有很很好的相相近关系系，所以以可以用用这两个个函数拟拟合表格格中的数数据。图4.22.1图4.22.2应用Liingoo求解：图4.22.3在此模型型求解中中，我假假设广告告费用xx可以在在0元到7000000元中任任意取值值，销售售价格可可以在22元到6元中任任意取值值，在由由线性规规划求解解后的值值向整数数靠近，此此时求得得的目标标函数的的值为线线性规划划中与现现实情况况比较接接近的解解。4.3 题目结结果由Linngo软软件解得得的结果果为：图4.33.1由线性规规划后的的结果可可知：在在广告费费用x=3.88595588万万元，销销售价格格y=44.91112885元时时，此时时的目标标函数值值最大：maxx=(-5.113333*y+50.42222)*10000*(-0.04226*xx2+0.440966*x+0.001888)*yy-x*100000=809901.71元元所以我们们取广告告费用xx=4万万元，销销售价格格y=55元，则则此时公公司获得得的利润润最大。最最大的利利润为8807558元。参考文献献1王王正东,数学软软件与数数学实验验.北京:科学出出版社,200042刘刘玉璉,傅沛仁仁等,数学分分析讲义义.北京.高等教教育出版版社,2200883吴吴建国,数学建建模案例例精编.北京:中国水水利水电电出版社社,20005源程序1 题目目一（1）ffuncctioon TTx=mmyfuun1(n)symss xfx=11;for k=11:n fx=fx+(-11)(2*kk-1)*xk/kk;%xx存在的的区间为为endgx=11;for l=11:n gx=gx+(-11)(l-11)*xxl/l; %x存存在的区区间为endTx=ffx-ggx;(2) funnctiion Tx=myffun22(x,n)fx=11;for k=11:n fx=fx+(-11)(2*kk-1)*xk/kk; %x存在在的区间间为endgx=11;for l=11:n gx=gx+(-11)(l-11)*xxl/l; %x存存在的区区间为endTx=ffx-ggx;(3) funnctiion myffun33(n)x=liinsppacee(-22/3,2/33,500);for k=11:leengtth(xx) Tn(k)=myffun22(x(k),n); a(kk)=llog(1-x(kk)/(1+x(kk); y(kk)=aa(k)-Tnn(k);endsubpplott(2,2,11);plott(x,a,*)titlle(ln(1-x)/(1+x)subpplott(2,2,22);plott(x,Tn,*) titlle(麦克劳劳林展式式)subpplott(2,2,33);plott(x,y,*)titlle(作差图图)subpplott(2,2,44);plott(x,Tn,r,x,a,b,x,yy,gg)2 题目目二（1）ffuncctioon VVh=mmyfuun4(h)a=177.8/2;b=122/2;L1=00.4;L2=22.055;Vh=aa*b*(L11+L22)*asiin(h-bb)/bb)+(h-b)/b)*sqrrt(11-(hh-b)2/b22)+ppi/22*110;（2）ffuncctioon mmyfuun5(X1,X2)for l=11:leengtth(XX1) a=117.88/2; b=112/22; L1=0.44; L2=2.005; Vh(l)=a*bb*(LL1+LL2)*assin(X11(l)-b)/b)+(X1(l)-b)/b)*sqrrt(11-(XX1(ll)-bb)22/b2)+pi/2*10; a(ll)=aabs(X2(l)+2622-Vhh(l); dissp(aa(l)end数据1：X1=1559.002176.14192.59208.50223.93238.97253.66268.04282.16296.03309.69323.15336.44349.57362.56375.42388.16400.79413.32425.76438.12450.40462.62474.78486.89498.95510.97522.95534.90546.82558.72570.61582.48594.35606.22618.09629.96641.85653.75665.67677.63678.54690.53690.82702.85714.91727.03739.19751.42763.70764.16776.53788.99801.54814.19826.95839.83852.84866.00879.32892.82892.84906.53920.45934.61949.05963.80978.91994.4310100.43310266.99910444.25510622.37710811.59911022.33311255.32211522.36611933.499./1100数据2：X2=5001001502002503003504004505005506006507007508008509009501000010500110001150012000125001300013500140001450015000155001600016500170001750018000185001900019500200002050020533.83321033.83321055.06621555.06622055.06622555.06623055.06623555.06624044.98824066.83324566.83325066.83325566.83326066.83326566.83327066.83327566.83328066.83328566.83329066.83329066.91129566.91130066.91130566.91131066.91131566.91132066.91132566.91133066.91133566.91134066.91134566.91135066.91135566.91136066.91136566.91137066.911数据3：Y=110.888266 112.663300 114.336499 116.113188 117.885199 119.660588 221.335200 223.008100 224.884688 226.557777 228.332888 330.005611 331.880233 333.554233 335.229711 337.005488 338.880811 440.555411 442.229399 444.003177 445.777499 447.449144 449.223544 550.997899 552.773766 554.448600 556.224255 557.998366 559.772977 661.445899 663.119299 664.995433 666.667900 668.443355 770.119788 771.995144 773.667411 775.443266 777.116277 778.888699 880.667033 880.777511 882.554100 882.556122 884.333533 886.005499 887.882388 889.553111 991.227755 993.003011 993.110988 994.883455 996.558933 998.333933 1000.008833 1001.883699 1003.558155 1005.331444 1007.006200 1008.880444 1110.555366 1110.554988 1112.226944 1114.002744 1115.776433 1117.551533 1119.225722 1221.001066 1222.775977 1224.550855 1226.222811 1227.998088 1229.773933 1331.447533 1333.222855 1334.995322 1336.770388 1338.445211（3）ffuncctioon mmyfuun7(X3,X4)for k=11:leengtth(XX3) a=117.88/2; b=112/22; L1=0.44; L2=2.005; Vh(k)=a*bb*(LL1+LL2)*assin(X33(k)-b)/b)+(X3(k)-b)/b)*sqrrt(11-(XX3(kk)-bb)22/b2)+pi/2*10; Wc(k)= -00.0884*XX3(kk)33+1.50665*XX3(kk)22+5.82116*XX3(kk)-11.71108; b(kk)=aabs(X4(k)-Vh(k)+Wc(k)+2622); dissp(bb(k);endfuncctioon mmyfuun9(X3,X4)for m=11:leengtth(XX3) a=117.88/2; b=112/22; L1=0.44; L2=2.005; Vh(m)=a*bb*(LL1+LL2)*assin(X33(m)-b)/b)+(X3(m)-b)/b)*sqrrt(11-(XX3(mm)-bb)22/b2)+pi/2*10; Wc(m)= -00.00025*X3(m)4-00.01167*X3(m)3+00.88876*X3(m)2+88.08826*X3(m)-4.338288; c(mm)=aabs(X4(m)-Vh(m)+Wc(m)+2622); dissp(cc(m);end数据1：X3=11550.77211233.99911011.15510800.51110611.36610433.29910266.08810099.544993.57978.08962.99948.26933.84919.69905.78892.10878.61865.30852.15839.14826.27813.52800.87788.33775.88763.51751.21738.98726.81714.70702.64690.61678.63666.68654.75642.84630.96619.08607.21595.35583.48571.61559.72547.82535.90523.95511.97499.96487.90475.80463.65451.43439.15426.80414.36401.84389.22376.49363.64350.67337.55324.27310.82297.18283.33269.24254.88240.21225.21209.81193.94177.54160.48142.62./1100数据2：X4=337066.911-552.772102.72152.72202.72252.72302.72352.72402.72452.72502.72552.72602.72652.72702.72752.72802.72852.72902.72952.7210022.72210522.72211022.72211522.72212022.72212522.72213022.72213522.72214022.72214522.72215022.72215522.72216022.72216522.72217022.72217522.72218022.72218522.72219022.72219522.72220022.72220522.72221022.72221522.72222022.72222522.72223022.72223522.72224022.72224522.72225022.72225522.72226022.72226522.72227022.72227522.72228022.72228522.72229022.72229522.72230022.72230522.72231022.72231522.72232022.72232522.72233022.72233522.72234022.72234522.72235022.72235522.72236022.72236522.72237022.7223 题目目三funcctioon I,JJ=mmyfuun8(A)B=A(:,11);for i=22:4 BB=B+A(:,i);endB=B/4;for j=11:4 D(:,j)=(AA(:,j)-B).2;endF=0;for l=11:4 F=FF+D(:,ll);endF=sqqrt(F/44);C=A(1,:);for j=22:177 C=CC+A(j,:);endC=C/17;for k=11:177 E(kk,:)=(AA(k,:)-C).2./177;endG=0;for q=11:177 G=GG+E(q,:);endG=sqqrt(G/117);for p=11:4 H(:,p)=(AA(:,p)-C(pp)./G(p).*F+B;endI=0;for n=11:4 I=II+H(:,nn);endI=I./4;a=siize(I);a=a(1);for h=11:a J(hh,1)=h;endfor w=11:166for s=ww:177if II(w)I(s) t=I(ww); I(w)=I(ss); I(s)=t; z=J(ww); J(w)=J(ss); J(s)=z;endendend数据：A=990 900 880 9795 90 800 99085 85 900 99085 90 800 99095 80 900 77090 80 900 77580 90 700 88585 85 700 88095 70 700 77590 60 600 66075 65 700 66075 70 600 66560 70 600 55065 60 500 55065 60 500 55060 60 500 55060 50 500 550