数学建模第三篇数学分支中的相关数学模型

1、第三篇第三篇 数学分支中的相关数学模型数学分支中的相关数学模型1 1 高等数学相关模型高等数学相关模型 1.11.1卫星轨道长度卫星轨道长度 1.21.2射击命中概率射击命中概率 1.31.3人口增长率人口增长率 2 2 线性代数相关模型线性代数相关模型 2.12.1投入产出综合平衡分析投入产出综合平衡分析 2.22.2输电网络输电网络3 3 概率统计相关模型概率统计相关模型 3.13.1合金强度与碳含量合金强度与碳含量 3.23.2年龄与运动能力年龄与运动能力 3.33.3商品销售量与价格商品销售量与价格 1 1 高等数学相关模型高等数学相关模型问题问题1.1 1.1 卫星轨道长度卫星轨道长

2、度人造地球卫星轨道可视为平面上的椭圆人造地球卫星轨道可视为平面上的椭圆.分析分析卫星轨道椭圆的参数方程卫星轨道椭圆的参数方程 )20(sin,costtbytax椭圆长度椭圆长度 分别是长、短半轴分别是长、短半轴 ba,dttbtadlL2/0212222)cossin(4椭圆积分椭圆积分无法解析计算无法解析计算 近地点距地球表面近地点距地球表面439km.远地点距地球表面远地点距地球表面2384km. .地球半径地球半径6371km. 求该卫星的轨道长度求该卫星的轨道长度. . 输出输出MATLAB程序程序 function y=x5(t)a = 8 7 5 5 ; b = 6 8 1 0

3、; y=sqrt(a2*sin(t).2+b2*cos(t).2); cleart=0:pi/10:pi/2;y1=x5(t);L1=4*trapz(t,y1) L2=4*quad(x5,0,pi/2,1e-6)L1=4.908996526785276e+004 L2=4.908996531830460e+004 输出输出求解求解梯形公式梯形公式 辛普森公式辛普森公式 68104396371,875523846371ba评注评注问题问题1.2 1.2 射击命中概率射击命中概率射击目标为正椭圆形区域射击目标为正椭圆形区域, ,弹着点与中心有随机偏差弹着点与中心有随机偏差. 分析分析设目标中心设目

4、标中心x=0,=0,y=0, =0, )(21222221),(yxyxyxeyxp无法解析计算无法解析计算 弹着点围绕中心成二维正态分布弹着点围绕中心成二维正态分布, ,偏差在偏差在X、Y方向独立方向独立.求炮弹落在椭圆形区域内的概率求炮弹落在椭圆形区域内的概率. 则弹着点则弹着点( (x, ,y) )概率密度函数概率密度函数 myx1001:,),(2222byaxdxdyyxpP炮弹命中椭圆形区域的概率炮弹命中椭圆形区域的概率 80,120ba椭圆在椭圆在X方向半轴长方向半轴长120m,Y方向半轴长方向半轴长80m.设弹着点偏差的均方差在设弹着点偏差的均方差在X和和Y方向均为方向均为10

5、0m.求解：求解：蒙特卡罗方法蒙特卡罗方法作变换作变换 ,bvyaux以以100(m)100(m)为为1 1单位，则单位，则 8 . 0, 2 . 1, 1bayxdudvvupabP),(1:,21),(22)(212222vuevupvbuaMATLAB程序程序 输出输出clear, n=100000;a=1.2;b=0.8;m=0;z=0;for i=1:n x=rand(1,2); y=0; if x(1)2+x(2)2=1 y=exp(- 0.5*(a2*P=0.3752, m=78552 x(1)2+b2*x(2)2); z=z+y;m=m+1; endendp=4*a*b*z/2

6、/pi/n,m 评注评注问题问题1 11.3 1.3 人口增长率人口增长率20世纪美国人口数据世纪美国人口数据(10(106 6 ), ), 年份年份 1900 1910 1920 1930 人口人口 76.0 92.0 106.5 123.2 1940 1950 1960 1970 1980 1990 131.7 150.7 179.3 204.0 226.5 251.4 计算各年份人口增长率计算各年份人口增长率. 记时刻记时刻t t人口为人口为x( (t t) )，则人口相对增长率为，则人口相对增长率为分析分析)(/)(txdtdxtr记记19001900年为年为k=0 =0 求解：求解：

7、数值微分三点公式数值微分三点公式 8 , 2 , 1,2011kxxxrkkkk99879021002034,2043xxxxrxxxxr 年增长率年增长率 2.20 1.66 1.46 1.02 1.04 1.58 1.49 1.16 1.05 1.04 评注评注计算计算程序程序 问题问题2 2已知某地区已知某地区20世纪世纪70年代的人口增长率年代的人口增长率, ,且且19701970年人口为年人口为210210（百万），（百万）， 年份年份 1970 1972 1974 年增长率（年增长率（%） 0.87 0.85 0.89 1976 1978 1980 0.91 0.95 1.10 试

8、估计试估计1980年的人口年的人口. 记时刻记时刻t t人口为人口为x( (t t) )，则人口增长满足微分方程，则人口增长满足微分方程 分析分析)()(txtrdtdx记记19701970年为年为k=0 =0 求解求解tduurextx0)(0)(评注评注0)0(xx1980年该地区人口为年该地区人口为230.2（百万）（百万） 数值积分梯形公式数值积分梯形公式 计算计算程序程序 clearx=76 92 106.5 123.2 131.7 150.7 179.3 204 226.5 251.4;r(1)=(-3*x(1)+4*x(2)-x(3)/20/x(1);for i=2:9 r(i)

9、=(x(i+1)-x(i-1)/20/x(i);endr(10)=(x(8)-4*x(9)+3*x(10)/20/x(10); r=100*r数值微分计算数值微分计算程序程序：people_model.m数值积分计算数值积分计算程序程序：people_model.m clear,x0=210;t=1970:2:1980;r=0.87 0.85 0.89 0.91 0.95 1.10;s=trapz(t,r);x1980=x0*exp(s/100) 为算出瑞士的国土面积，首先对瑞士地图作如下测量：为算出瑞士的国土面积，首先对瑞士地图作如下测量：以由西向东方向为以由西向东方向为x轴，由南到北方向为

10、轴，由南到北方向为y轴，选择方便轴，选择方便的原点，并将从最西边界点到最东边界点在的原点，并将从最西边界点到最东边界点在x轴上的区轴上的区间适当地划分为若干段，在每个分点的间适当地划分为若干段，在每个分点的y方向测出南边方向测出南边界点和北边界点的坐标，得到数据（界点和北边界点的坐标，得到数据（mm）见下页）见下页. 习题：习题： 国土面积问题国土面积问题 根据地图比例根据地图比例, ,18mm相当于相当于40km，试由测量数据计，试由测量数据计算瑞士国土的近似面积，与它的精确值算瑞士国土的近似面积，与它的精确值41288km比较比较. 计算计算程序（数据见下页程序（数据见下页, ,省略）省略

11、）clear,x=7.0 158;y1=44 68;y2=44 68;s=trapz(x,y2)-trapz(x,y1);s=s*402/182, error=s-412887.0 10.5 13.0 17.5 34.0 40.5 44.5 48 56 61 68.5 76.5 80.5 9144 45 47 50 50 3 8 30 30 34 36 34 41 45 4644 59 70 72 93 100 110 110 110 117 118 116 118 11896 101 104 106.5 111.5 118 123.5 136.5 142 146 150 157 158 43

12、 37 33 28 32 65 55 54 52 50 66 66 68121 124 121 121 121 122 116 83 81 82 86 85 68x1y2yx1y2y2 2 线性代数相关模型线性代数相关模型背景背景2.1 2.1 投入产出综合平衡分析投入产出综合平衡分析国民经济各个部门之间存在着相互依存的关系国民经济各个部门之间存在着相互依存的关系. .投入产出综合平衡模型投入产出综合平衡模型: :根据各部门间的投入根据各部门间的投入产出关产出关系，确定各部门的产出水平，以满足社会的需求系，确定各部门的产出水平，以满足社会的需求. .设国民经济仅由农业、制造业、和服务业三个设国

13、民经济仅由农业、制造业、和服务业三个部门构成，已知某年它们之间的投入产出关系、部门构成，已知某年它们之间的投入产出关系、外部需求、初始投入等如表（产值单位为亿元）外部需求、初始投入等如表（产值单位为亿元） 简化简化问题问题每个部门在运转中将其他部门的产品或半成品经每个部门在运转中将其他部门的产品或半成品经过加工（投入）变为自己的产品（产出）过加工（投入）变为自己的产品（产出）. . 产出产出 投入投入农业农业制造业制造业服务业服务业外部外部需求需求总产出总产出农业农造业制造业301045115200服务业服务业2060/70150外部外部 需求需求3511075总产出

14、总产出100200150说明说明假定每个部门的产出与投入是成正比的，由上表假定每个部门的产出与投入是成正比的，由上表能够确定这三个部门的投入产出表能够确定这三个部门的投入产出表 产产出出投入投入农业农业制造业制造业服务业服务业农业农业0.150.100.20制造业制造业0.300.050.30服务业服务业0.200.300说明说明表中数字称为表中数字称为投入系数投入系数或或消耗系数消耗系数 假设系数是常数假设系数是常数 设有设有n个部门，已知投入系数，给定外部需求，建个部门，已知投入系数，给定外部需求，建 立求解各部门总产出的模型立求解各部门总产出的模型. . 如果今年对农业、制造业、服务业的

15、外部需求分如果今年对农业、制造业、服务业的外部需求分别为别为5050，150150，100100亿元，三个部门总产出？亿元，三个部门总产出？ 模型可行模型可行：对于任意给定的、非负的外部需求：对于任意给定的、非负的外部需求, ,都都能得到非负的总产出能得到非负的总产出. .为使模型可行为使模型可行, ,投入系数满足？投入系数满足？ 如果三个部门的外部需求分别增加如果三个部门的外部需求分别增加1 1个单位，他们个单位，他们的总产出应分别增加多少？的总产出应分别增加多少？ 分析分析投入产出综合平衡分析投入产出综合平衡分析 若有若有n个部门，记一定时期内第个部门，记一定时期内第i个部门的总产出为个部

16、门的总产出为xi，其中对第其中对第j个部门的投入为个部门的投入为xij，满足的外部需求为，满足的外部需求为di，则，则nidxxinjiji, 2 , 1,1投入产出表每一行都满足投入产出表每一行都满足 记第记第j个部门的单位产出需要第个部门的单位产出需要第i个部门的投入为个部门的投入为aij ，在每个部门的产出与投入成正比的假定下，有在每个部门的产出与投入成正比的假定下，有njixxajijij, 2 , 1,nidxaxinjjiji, 2 , 1,1记记投入系数矩阵投入系数矩阵nnijaA)(产出向量产出向量 Tnxxx),(1需求向量需求向量 Tnddd),(1dAxxdxAI)(则则

17、 或或 若若I-A可逆，则可逆，则 dAIx1)( 各部门总产出各部门总产出 MATLAB程序程序 a=0.15 0.1 0.2;0.3 0.05 0.3;0.2 0.3 0;a=0.15 0.1 0.2;0.3 0.05 0.3;0.2 0.3 0; d=50;150;100;d=50;150;100; b=eye(3)-a;b=eye(3)-a; x=bdx=bd , c=inv(b)c=inv(b) 三部门总产出三部门总产出:139.2801,267.6056,208.1377:139.2801,267.6056,208.1377亿元亿元 外部需求分别增加外部需求分别增加1 1个单位时，

18、总产出分别增加个单位时，总产出分别增加C=1.3459 0.2504 0.3443C=1.3459 0.2504 0.3443 0.5634 1.2676 0.4930 0.5634 1.2676 0.4930 0.4382 0.4304 1.2167 0.4382 0.4304 1.2167 部门关联系数部门关联系数 当对农业的需求增加当对农业的需求增加1 1个单位时个单位时, ,农业、制造业、和服农业、制造业、和服务业的总产出分别增加务业的总产出分别增加1.34591.3459，0.56340.5634，0.43820.4382单位单位 dAIx1)( 模型可行模型可行njaniij, 2

19、 , 1, 11若若问题问题2.2 2.2 输电网络输电网络一种大型输电网络可简化为电路一种大型输电网络可简化为电路 负载电阻负载电阻 nRRR,21线路内阻线路内阻 nrrr,21电源电压电源电压V TnIII),(1负载电流负载电流 列出各负载上电流的方程列出各负载上电流的方程 设设 讨论情况讨论情况 18, 6, 1,11VRrrrrRRRnnn=10,=10,求求 nIII,21及总电流及总电流 0In2r1rnr2R1RnR2InI1IV2i1ini分析分析记记 nrr,1上的电流为上的电流为 niii,21根据电路中电流、电压关系根据电路中电流、电压关系, ,列出列出 111122

20、221111nnnnnnIRIRirIRIRirVIRi rnnnnniIiiIiiIiiI112321210)(0)()(11222211121111nnnnnnnIrRIRIrIrRIRVIrIrIrR消消 niii,21和和 求电流方程求电流方程 ERI 求电流方程求电流方程 其中其中 nnnrRRrrrRRrrrrRR12222111111TTnVEIIII0 , 0 ,21MATLAB计算电流计算电流程序程序 r=1;R=6;v=18;n=10; b1=sparse(1,1,v,n,1);b=full(b1);a1=triu(r*ones(n,n);a2=diag(R*ones(1,

21、n);a3=-tril(R*ones(n,n),-1)+ tril(R*ones(n,n),-2);a=a1+a2+a3; I=ab;I0=sum(I) k0 01 12 23 34 45 56 67 78 89 910105.99705.99706.00006.00002.00052.00052.00002.00001.33441.33441.33331.33330.89070.89070.88890.88890.59550.59550.59260.59260.39950.39950.39510.39510.27020.27020.26340.26340.18580.18580.17560.

22、17560.13240.13240.11710.11710.1010.1011 10.0780.0780 00.08670.08670.05200.0520k11111212131314141515161617171818191920200.03470.03470.02310.02310.01540.01540.01030.01030.00690.00690.00470.00470.00320.00320.00230.00230.0010.0018 80.00150.0015)20()10(nInIkk)20( nIk说明说明从从n=10=10到到n=20=20，I0几乎不变，几乎不变，I1-

23、 -I5变化也很小变化也很小 Ik+1差不多是差不多是Ik的的2/32/3倍倍如果如果n增加到增加到50,100?50,100? 可以得到类似的结论可以得到类似的结论 证明证明 习题：种群的繁殖与稳定收获习题：种群的繁殖与稳定收获 种群的数量因繁殖而增加，因自然死亡而减少，种群的数量因繁殖而增加，因自然死亡而减少，对于人工饲养的种群（比如家畜）而言，为了对于人工饲养的种群（比如家畜）而言，为了保证稳定的收获，各个年龄的种群数量应维持保证稳定的收获，各个年龄的种群数量应维持不变不变.种群因雌性个体的繁殖而改变，为方便起种群因雌性个体的繁殖而改变，为方便起见以下种群数量均指其中的雌性见以下种群数量

24、均指其中的雌性. 种群年龄记作种群年龄记作 当年年龄当年年龄 的种群数量的种群数量记作记作 ，繁殖率记作，繁殖率记作 （每个雌性个体一年繁（每个雌性个体一年繁殖的数量），自然存活率记作殖的数量），自然存活率记作 为为一年的死亡率），收获量记作一年的死亡率），收获量记作 ，则来年年龄，则来年年龄 的种群数量的种群数量 应为应为, 2 , 1nkkkxkbkkkdds,1(khkkx) 1, 2 , 1(,11111nkhxsxhxbxkkkknkkk （1）若）若 已知，给定收获量已知，给定收获量 ，建立求各年龄的，建立求各年龄的稳定种群数量稳定种群数量 的模型（用矩阵、向量表示）。的模型（用矩

25、阵、向量表示）。 （2）设）设 如要求如要求 为为500，400，200，100，100,求求 。 （3）使）使 均为均为500，如何达到，如何达到?kksb ,khkx, 6 . 0, 4 . 0, 3, 5, 0, 5324143521ssssbbbbbn51 hh51 xx51 hh要求各个年龄种群数要求各个年龄种群数量每年维持不变，即量每年维持不变，即) 1, 2 , 1(nkxxkk543215432143215432154321hhhhhxxxxxssssbbbbbxxxxx提示提示 hAxIxhIAx1)(3 3 概率统计相关模型概率统计相关模型问题问题3.13.1合金强度与碳含

26、量合金强度与碳含量合金的强度合金的强度y( (kg/mmkg/mm) )与其中的碳含量与其中的碳含量x( ( %)%)有比有比较密切的关系，从生产中收集一批数据较密切的关系，从生产中收集一批数据. .求拟合函数求拟合函数y(x) ，再用回归分析进行检验，再用回归分析进行检验. .x0.10 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.20 0.21 0.23 0.10 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.20 0.21 0.23 y42.0 41.5 45.0 45.5 45.0 47.5 49.0 55.

27、0 50.0 55.0 55.5 60.542.0 41.5 45.0 45.5 45.0 47.5 49.0 55.0 50.0 55.0 55.5 60.5描点作图描点作图 分析分析0.10.120.140.160.180.20.220.24404550556065y与与x近似近似为线性为线性 拟合拟合 y=ax+b MATLAB程序程序 x=0.1:0.01:0.23;x=0.1:0.01:0.23;x=x(1:9),x(11:12),x(14);x=x(1:9),x(11:12),x(14);y=42,41.5,45,45.5,45,47.5,y=42,41.5,45,45.5,45,

28、47.5,49,55,50,55,55.5,60.5;49,55,50,55,55.5,60.5;pp=pp=polyfitpolyfit(x,y,1);(x,y,1);xxxx=0.08:0.01:0.25;=0.08:0.01:0.25;yy=polyval(pp,xxyy=polyval(pp,xx););plot(x,y,rplot(x,y,r* *,xx,yy,xx,yy) ) 0.080.10.120.140.160.180.20.220.240.2635404550556065拟合拟合 y=ax+ba=140.6194,b=27.0269 评注评注是否线性显著是否线性显著 有无异

29、常点有无异常点预测预测 MATLAB统计工具箱统计工具箱 多元线性回归多元线性回归语法语法 b=regress(Y,X)b=regress(Y,X) b,bint,r,rintb,bint,r,rint,stats=regress(Y,X,alpha) ,stats=regress(Y,X,alpha) Y,X为按列排列的数据为按列排列的数据说明说明b,bintb,bint为回归系数估计值及其置信区间为回归系数估计值及其置信区间 m10,alphaalpha为显著性水平（缺省时设定为为显著性水平（缺省时设定为0.050.05）statsstats包括包括R R2 2，F F，概率，概率p pr

30、,rintr,rint为残差及置信区间，可用为残差及置信区间，可用rcoplot(r,rintrcoplot(r,rint) )画图画图 合金强度与碳含量问题合金强度与碳含量问题回归模型回归模型 xy10 回归模型与统计检验回归模型与统计检验 MATLAB程序程序 clear,x1=0.1:0.01:0.18;clear,x1=0.1:0.01:0.18;x=x1 0.2 0.21 0.23;y=42 41.5 45 45.5 45 47.5 49 55 50 55 55.5 60.5;x=ones(12,1) x;x=ones(12,1) x;b,bint,r,rintb,bint,r,ri

31、nt,stats=,stats=regressregress(y,x);(y,x);b,bint,stats,rcoplot(r,rintb,bint,stats,rcoplot(r,rint) ) b= 27.0269 140.6194b= 27.0269 140.6194bintbint=22.3226 31.7313=22.3226 31.7313 111.7842 169.4546 111.7842 169.4546stats=0.9219 118.0670 0.0000stats=0.9219 118.0670 0.0000 6194.140,0269.2710 y=27.0269+

32、140.6194x线性显著线性显著模型成立模型成立 有无异常点有无异常点画残差分布图画残差分布图24681012-6-4-20246Residual Case Order PlotResidualsCase Number除第除第8个数据外其余残差个数据外其余残差的置信区间均包含零点的置信区间均包含零点第第8个点应视为异常点，个点应视为异常点，剔除后重新计算，可得剔除后重新计算，可得 b= 26.8968 139.9043 bint=24.1330 29.6606 122.7939 157.0148stats=0.9744 342.1259 0.0000 评注评注 预测预测 clear,x1=0

33、.1:0.01:0.18;clear,x1=0.1:0.01:0.18;x=x1 0.2 0.21 0.23;y=42 41.5 45 45.5 45 47.5 49 55 50 55 55.5 60.5;x(8)=;y(8)=;x(8)=;y(8)=;x=ones(11,1) x;x=ones(11,1) x;b,bint,r,rintb,bint,r,rint,stats=,stats=regressregress(y,x);(y,x);b,bint,stats,b,bint,stats,rcoplot(r,rintrcoplot(r,rint) ) 问题问题3.2 3.2 年龄与运动能力

34、年龄与运动能力 将将1717至至2929岁的运动员每两岁一组分为岁的运动员每两岁一组分为7 7组，组，求年龄对这种运动能力的影响关系求年龄对这种运动能力的影响关系. .多项式回归多项式回归 年龄年龄 17 19 21 23 25 27 2917 19 21 23 25 27 29第一人第一人 20.48 25.13 26.15 30.0 26.1 20.3 19.3520.48 25.13 26.15 30.0 26.1 20.3 19.35第二人第二人 24.35 28.11 26.3 31.4 26.92 25.7 21.324.35 28.11 26.3 31.4 26.92 25.7

35、21.3分析分析MATLAB散点图程序散点图程序 每组两人测量其旋转定向能力每组两人测量其旋转定向能力. .clear,x=17:2:29;clear,x=17:2:29;y1=20.48,25.13,26.15,30.0,y1=20.48,25.13,26.15,30.0,26.1,20.3,19.35;26.1,20.3,19.35;y2=24.35,28.11,26.3,31.4,y2=24.35,28.11,26.3,31.4,26.92,25.7,21.3;26.92,25.7,21.3;plot(x,y1,+,x,y2,+)plot(x,y1,+,x,y2,+)axis(15 30

36、 15 35)axis(15 30 15 35) 152025301520253035应拟合一条二次曲线应拟合一条二次曲线 可利用可利用polyfit polyfit 一元多项式回归一元多项式回归年龄与运动能力的二次模型年龄与运动能力的二次模型 3221axaxayMATLAB程序程序 x1=17:2:29;x=x1,x1;x1=17:2:29;x=x1,x1;y=20.48 25.13 26.15 30.0 26.1 20.3 19.35 y=20.48 25.13 26.15 30.0 26.1 20.3 19.35 24.35 28.11 26.3 31.4 26.92 25.7 21.

37、3;24.35 28.11 26.3 31.4 26.92 25.7 21.3;p,s=polyfitp,s=polyfit(x,y,2); p(x,y,2); p p=-0.2003 8.9782 -72.2150 p=-0.2003 8.9782 -72.2150 a1=-0.2003 a2=8.9782 a3=-72.2150 a1=-0.2003 a2=8.9782 a3=-72.2150 S S是一个数据结构是一个数据结构 Y,delta=polyconfY,delta=polyconf(p,x,s);Y (p,x,s);Y 得到得到x与与y的拟合效果的拟合效果求解求解Y=22.52

38、43 26.0582 Y=22.5243 26.0582 统计检验统计检验y1=mean(y); resquare=sum(Y-y1).2)./sum(y-y1).2), s=sqrt(sum(y-Y).2)./12), rsquare=0.6980 s=2.0831 衡量拟合优劣的指标衡量拟合优劣的指标 x=17:2:29;y1=20.48,25.13,26.15,30.0,26.1,20.3,19.35;y2=24.35,28.11,26.3,31.4,26.92,25.7,21.3;plot(x,y1,+,x,y2,*),hold onx=x x;y=y1 y2;p=polyfit(x,

39、y,2); xx=17:0.2:29;yy=polyval(p,xx);plot(xx,yy),hold off1618202224262830161820222426283032交互界面交互界面: : polytoolpolytool(x,y,2)(x,y,2) 问题问题3.3 3.3 商品销售量与价商品销售量与价格格 某厂生产电器的销售量某厂生产电器的销售量y与竞争对手的价格与竞争对手的价格x1 1和本厂的价格和本厂的价格x2 2有关有关. .在在1010个城市的销售记录个城市的销售记录 建立建立y与与x1 1和和x2 2的关系式的关系式. . x1 1元元120 140 190 130

40、155 175 125 145 180 150120 140 190 130 155 175 125 145 180 150 x2 2元元100 110 90 150 210 150 250 270 300 250100 110 90 150 210 150 250 270 300 250y个个102 100 120 77 46 93 26 69 65 85102 100 120 77 46 93 26 69 65 85分析分析 对模型和系数进行检验对模型和系数进行检验. . 若本厂售价若本厂售价160160元，对手售价元，对手售价170170元，预测销售量元，预测销售量. . 画画散点图散点

41、图 1001502002030405060708090100110120130(x1,y)01002003002030405060708090100110120130(x2,y)y与与x2 2有较明显的线性关系，而有较明显的线性关系，而y与与x1 1之间的关系则难之间的关系则难以确定，作几种尝试，用统计分析决定优劣以确定，作几种尝试，用统计分析决定优劣. . 设回归模型设回归模型22110 xxyMATLAB程序程序 clear,x1=120 140 190 130 155 175 125 145 180 150; x2=100 110 90 150 210 150 250 270 300 250; x2=100 110 90 150 210 150 250 270 300 250; y=102 100 120 77 46 93 26 69 65 85; x=ones(1