§4.常见的数学建模方法(1)---数据拟合(曲线拟合)法.ppt

4.常见的数学建模方法(1)-数据拟合（曲线拟合）法,在建立数学模型时,实际问题有时仅给出一组数据.处理这类问题的较简单易行的方法是通过数据拟合法求得“最佳”的近似函数式-经验公式.从几何上看就是找一条“最佳”的曲线,使之和给定的数据点靠得最近,即进行曲线拟合.根据一组数据来确定其经验公式,一般可分为三步进行:,（1）决定经验公式的形式.根据所描绘的系统固有的特点,参照已知数据的图形和特点或者它应服从的规律来决定经验公式的形式.这一步是关键的一步.,（2）决定经验公式中的待定参数.一般可用线性情况下的最小二乘法.它误差较小,适用于测定数据比较精确的情况.在使用最小二乘法时,如遇到数学模型是非线性经验公式时其中参数的待定,通常是尝试能否经适当的变量替换,将之化为线性模型来计算.,（3）进行模型检验.求得确定的经验公式后,将实际测定值与用公式算出的理论值进行比较.,线性模型下的最小二乘法法则是：如果一组数据为：（xi,yi）,(i=0,n)，它服从线性函数y=kx+b模型,则,在决定经验公式的形式时,大致思路是:a)利用所研究系统的有关问题在理论上已有的结论,来确定经验公式的形式.b)在无现成理论情况下,最简单的处理手段是用描图的方法,将数据点连成光滑曲线,把它与已知函数曲线进行比较,找出与之比较接近的曲线.c)如要考虑所建立的模型必要的逻辑性与理论价值,可利用合适的数学方法,对所研究系统的有关问题进行定量化的机理分析,导出较为严密的数学公式.,实例1.找出基于下列数据的铜棒长度l与温度t之间关系的经验公式.,建模过程:利用已有的物理学固体热胀冷缩定律:l=l0(1+at)作为该组数据应服从的数学模型，如记l-1000=l,l01000=b,al0=k,则有l=b+kt.可以算得：,根据最小二乘法公式，,可得：l=999.804（1+0.0000212t）.,最后检验该模型(经验公式):,残差的平方和为:v2=0.00194,这个结果应该说是较好的.,实例2.找出基于下列数据的油的粘度y与温度x之间关系的经验公式.,建模过程:无现成机理明确的公式,使用描点比较法：,可以认为该光滑曲线相似于一条双曲线,故设其数学模型为y=axb(b0).为了将它化为线性模型,两边取对数，再作变量替换：Y=lny,X=lnx,即得线性模型：Y=A+bX,其中A=lna,而(X,Y)的数据为:（lnxi,lnyi）,(i=1,8).,用最小二乘法算得：a=17.2463,b=-0.6048.由此最后可得到油的粘度y与温度x之间依赖关系的数学模型为:,检验该模型(经验公式):y=17.2463x-0.6048.,残差的平方和为:v2=0.0146,这个结果应该说也是较好的.,说明:该例中的变量替换方法运用,使得线性模型的最小二乘法公式应用范围大大扩大.常见的非线性模型的变换方式如下表所列:,实例.找出基于下列数据的美国马萨诸塞州生产量、劳动力和投资之间变化的经济增长模型（道格拉斯Douglas生产函数模型）,实例3.某研究所为了研究三种肥料氮,磷,钾对于土豆和生菜的作用,分别对每种作物进行了三组试验.实验数据如下列表格所示,其中ha表示公顷,t表示吨,kg表示千克.试建立反映施肥量与产量关系的数学模型.,氮施肥量(公斤/公顷)与土豆产量（吨/公顷）关系的实验数据,磷施肥量(公斤/公顷)与土豆产量（吨/公顷）关系的实验数据,钾施肥量(公斤/公顷)与土豆产量（吨/公顷）关系的实验数据,氮施肥量(公斤/公顷)与生菜产量（吨/公顷）关系的实验数据,磷施肥量(公斤/公顷)与生菜产量（吨/公顷）关系的实验数据,钾施肥量(公斤/公顷)与生菜产量（吨/公顷）关系的实验数据,1.磷施肥量x关于土豆产量y的情况.,描点图为:,可选择,作为经验公式.,为了运用线性模型的最小二乘法公式,令,最终的数学模型是:,根据这个模型,可以得到土豆的最高极限产量是43吨.这个结论从定性角度看,与农业资料的结论是一致的,即在一定的范围内磷施肥量可以使土豆产量增长,但过多地施磷肥对土豆产量不起作用.在这一点上,该模型是经得起实际检验的.,得y=a+bx.由此可算得:a=0.0232,b=0.0073.,2.磷施肥量x关于生菜产量y的情况.,描点图为:,由描点图可知,在模型建立中应注意以下两个因素:,当磷肥施肥量为零时,生菜产量并非为零,这说明土壤中原来就含有一定量的磷肥成分;,(2)实验数据说明,磷肥施肥量再多不会引起产量明显下降,而使生菜产量趋于一个渐近值,即极限产量.,考虑到上面一些分析,可采用双曲线模型:,这里a为生菜极限产量数.,为了利用线性模型的最小二乘法,令X=1/(1+x),Y=y,k=y0a,化为线性函数模型:Y=a+kX.,根据最小二乘法计算公式和统计数据,先算得a和k,然后再算出y0.,相应的统计数据为:,在建立曲线拟合法的数学模型时,如果能尽量做一些定量化的机理分析,然后运用数学手段推导出合理的数学模型,则建模的效果会更好一些.,实例4.利用例3的资料,建立土豆产量y和生菜产量y依赖于氮施肥量x的数学模型,并由此求出氮和磷的最佳施肥数量.(氮肥价格350元/吨,土豆价格0.8元/公斤,生菜价格0.2元/公斤),建模假设:根据所给的数据和实际经验,当施肥量合适时,土豆产量和生菜产量随氮施肥量的增加而增加,但是氮施肥量的过量会造成土豆产量和生菜产量的下降.如设xm为达到土豆或生菜最高产量时的施肥量,现在假定边际产量yx与xmx成正比(Nicklas和Miller理论).,模型建立:,模型求解:y=ax2+bx+c其中c=y(0)=15.18(或11.02).,令x=x,y=(y-c)/x,可以化为线性模型：y=ax+b.,根据所给数据,运用线性模型的最小二乘法公式,得土豆产量y依赖于氮施肥量x的数学模型:y=-0.00034x2+0.197x+15.18;生菜产量y依赖于氮施肥量x的数学模型:y=-0.00024x2+0.101x+11.02.,模型分析与模型决策:当下列关系y(x0)=Tx/Ty(Tx,Ty分别为氮肥和生产作物的价格)成立时,投入一吨肥料得到的效益最大,此时的施肥量即为最佳施肥量;这是因为：,利润S(x)=土豆收益-氮肥支出=y(x)Ty-xTx,当S(x)=0时，利润最大,即：y(x)Ty-Tx=0，解得y(x0)=Tx/Ty。,由土豆产量y依赖于氮施肥量x的数学模型0.1970.00068x0=350/0.8,对土豆的最佳氮施肥量x0=290.57kg/ha.,由生菜产量y依赖于氮施肥量x的数学模型0.1010.00048x=350/0.2,对生菜的最佳氮施肥量x0=203.57kg/ha.,实例5.利用例3的资料,建立生菜产量y依赖于磷施肥量x的数学模型.,建模假设:经过查阅农业资料,有关文献中的米采利希理论认为,只增加某种养分时,引起产量的增加与该种养分供应充足时达到的最高产量和现在产量之差相关.现在,假定边际产量dy/dx与ymy成正比(ym:极限产量),模型建立:,模型求解:,这是一个非线性模型,无论使用怎样的变换,都不能化为线性模型.这样的模型,称为本质非线性模型.在实例7之后，我们再来看如何用已知数据来算得ym和k.,实例6.建立赛艇比赛成绩的模型(McMahon模型),赛艇是一种靠桨手划桨前进的小船,分单人艇,双人艇,四人艇,八人艇四种.T.A.McMahon比较了各种赛艇19641970年四次2000米比赛的最好成绩数据表后,认为比赛成绩与桨手数量之间存在着某种关系,并建立了一个数学模型来解释这种关系.,建模目的:寻求桨手数n与比赛成绩t之间的函数关系t=t(n).,建模假定:（1）艇的几何形状相同，艇长:艇宽=l:b=constant,艇重w0与桨手数n成正比；（2）艇速v=constant,阻力f与sv2成正比，s为艇浸没部分面积；（3）桨手体重w相同，划桨功率p不变，且p与w成正比。,建模过程：总功率与阻力和速度的乘积成正比：。,假设（2）,由此可得：,浸没部分面积应与艇的特征尺寸的平方成正比：,另一方面，艇的排水体积应与艇的特征尺寸的立方成正比：,故可得：,艇和桨手的总重w=w0+nw,根据阿基米德定律：,上面的,由于比赛成绩t（时间）与v成反比，所以,根据以上机理分析，可以认为有经验公式：,令X=lnn,Y=lnt,A=lna,得Y=A+bX.,由实际数据代入最小二乘法公式可得：,假设（3）,实例7.人口Logistic模型,根据下面一组统计数据（ti,xi）(i=1,2,n)：,寻求人口数x随时间t变化的具体数学模型。,建模过程：,模型求解：,建模假定（1）人口x(t)关于时间t的增长率与人类生存空间容纳度成正比。（2）人口基数x0是已知的。,这是一个本质非线性模型,无论使用怎样的变换,都不能化为线性模型.为了用所给的统计数字通过线性模型的最小二乘法公式