基于数据分析的建模方法问题* 在建立数学模型的过程中,经常需要建立变量之间的关系. *由于对研究对象的内部机理不甚了解,不能通过合理的假设,或根据物理定律、原理, 经过机理分析法而得到.解决思路1* 借助于由实验或测量得到的一批离散数据.*通过对数据充分观察和分析, 获得数据所含信息;*揭示变量间的内在联系;*选择适当的数学式对变量间的关系进行拟合.yox2两类变量关系两类数据确定性关系确定的函数关系相关关系存在相依关系,但未达到相互确定的程度.已知规律(函数)的测试数据(在特定时间点或距离上的数据)呈现随机性的数据,可看成具有某种概率分布的随机样本值.3针对两种不同类型的数据, 有不同的建立模型方法:1. 数据拟合法(适用于第一类数据)基本思想已知函数y= f(x) 的一组测试数据(xi ,yi),(i=1,2,…,n),寻求一个函数ψ(x),使ψ(x)对上述测试数据的误差较小,即ψ(xi)≈yi,于是可以用ψ(x)来近似替代f (x).常用的数据拟合方法:一般插值法、最小二乘法、样条函数光顺法等.4插值法的基本思想寻找f(x)的近似替代函数φ(x), 在插值节点xi 上满足φ( xi )=yi,(i=1,2,…,n),其余点用φ(x)近似替代f (x), 称φ(x)为f (x)的插值函数.yixif(x)x5最小二乘法基本思想寻找f(x)的近似替代函数φ(x), 使n2min.(f(xi)(xi))2. 随机分析方法对于随机数据进行拟合,可用统计学中的回归分析方法或时间序列分析方法. 二.经验模型的建立以上两种建模方法都是建立在对数据进行充分分析的基础上.6i1寻找或选择适当的函数拟合变量之间的关系(函数关系或回归关系)是重要的环节.一般步骤1)绘制数据散布图;见p1562)分析数据散布图;3)选择函数关系形式.1)通过分析数据散布图可以获得对变量间关系的感性认识, 形成初步的看法, 以便于对问题做进一步的分析.7氮施肥量N-土豆产量数据散布图8磷施肥量-土豆产量数据散布图92)分析数据散布图;对数据散布图进行分析,可以分析出变量的关系是:1)线性的还是非线性的?2)有无周期性?3)呈现何种变化趋势?变化率如何?…,等等有用的初步结论.10例6.2.1 建立一个简洁的函数关系式来描述某个地区人的身高和体重的对应关系, 数据见表7.4(p156). 曲线特征是体重W随身高H 的增长而单调增长,但可以观察到是非线性增长.身高-体重数据散布图11练习试分析以下问题1. 氮施肥量N、磷施肥量P 关于土豆产量的数据散布图(P153例7.1.1).2. 海浪潮高度x随时间t的数据散布图.123)选择函数关系形式原1. 形式尽可能简洁, 尽可能线性化;则2. 依据实际问题的精度要求,合乎实际规律.续例6.2.1高体重关系.a选择幂函数W= , cH描述身优点此函数可以线性化.两边取对数, 有lnWalnHlnc13令ylnW,xlnH,blncyaxb.2变换为线性函数例6.2.2可选二次函数yb0b1xb2x注:其中b0= y(0) = 15.18.描述氮肥施肥量与土豆产量间的变量关系.关于磷肥施肥量和土豆产量的变量关系可选择威布尔模型:14yABeKx,x0合理性?有1.当x,yA;2.y(0)AB;3.y 是单调升函数.也可以选择S函数:y1abex,x015S函数也满足:11.当x,y;a12.y(0);ab3. y 是单调升函数;哪个模型更好?分析S 模型所含参数更少, 另外若令1xy,xe,y可得线性模型yabx.16重要定理(维尔斯脱拉斯) 若函数f(x)在有限闭区间上连续, 则存在一个多项式序列{Pn(x)}在有限闭区间[a , b]上一致收敛于f(x).称f(x) 在[a,b] 上可由多项式函数逼近.nlimPn(x)f(x),x[a,b]例6.2.3 估计供水塔的水流量试用以下数据估计任意时刻(包括水泵正在输水的时间内)从水塔流出的流量f(t),并估计一天的总用水量.17时间(秒) 水位(0.01英尺) 时间(秒) 水位(英尺)0 3175 46636 33503316 3110 49953 32606635 3054 53936 316710619 2994 57254 308713937 2947 60574 301217921 2892 64554 292721240 2850 68535 284225223 2795 71854 276728543 2752 75021 269732284 2697 79254 35932 水泵开动39332 水泵开动39435 3550 89953 339743318 3445 92370 3340水泵开动82649 水泵开动85968 347518某小镇某天水塔水位散布图19假设水位高度(或水塔的水容量)是连续变化的.可以选择n 次多项式Pn(x)2nPn(x)a0a1xa2xanx来近似描述水位随时间的变化规律.问题归结为选择足够大的n及估计各个系数值.思考为什么考虑用多项式函数?有什么优点?20参见电子科技大学《概率统计》p228 “非线性交调的频率设计”问题.21