住房的合理定价-20088387郭诚.doc

西南交通大学峨眉校区 数学建模 姓名： 学号： 联系电话计算机与通信工程系 2010年4月日  当冰雹由高空落下时，它受到地球引力和空气阻力的作用，阻力的大小与冰雹的形状和速度有关，一般可以对阻力作两种假设： （1）阻力大小与下落的速度成正比； （2）阻力大小与速度的平方成正比； 请根据两种不同的假设，建立速度满足的微分方程，并计算冰雹下落的极限速度（已知初速度v(0)=0、冰雹质量m、重力加速度g、正比例系数k0）。 解1： 由物理学可知，冰雹受到地球引力和空气阻力的影响，其加速度应该为： (1.1) 而显然 (1.2) 联立上述两式子得到： (1.3) 显然3式就是所要求解的微分方程模型。 在（3）两边作Laplace 变换并且代入v(0)=0，得到： (1.4) 解代数方程，得 (1.5) 最后对做Laplace 反变换，最终求得： (1.6) 从式中不难看出，当时，v(t)的值为。 解二： 分析方法和第一题很类似，只不过将上述第一式改为： (2.1) 这个时候所要求解的模型变成了 (2.2) 由于这是一个非线性的微分方程，以我的数学水平很难给出它的精确解，但是，我们可以来讨论一下这个冰雹的下落情况，在一开始速度比较小时（阻力小于重力），冰雹的速度总是增加的。当速度达到一定时（阻力等于重力），速度不再增加，显然这个时候就是冰雹速度的最大值。 所以，若要冰雹的速度达到最大，必有： (2.3) 即： (2.4) 即： (2.5) 但是只如此，我们便知知道冰雹的最终速度，而对其运动过程了解不多，所以可以采用数值计算的办法来求解这个非线性微分方程。此时，取m=0.1kg,k=0.1,g=9.8。(k是随意取的) 数值解如图表 1所示： 图表 1 模型的改进： 对于阻力，我们所能够知道的是，它一定是与速度成正比的。但是具体是怎么样一个“正比”法呢？ 请想一想，当您在散步的时候，能感觉到空气的阻力吗？显然是不能的！而当您参加运动会百米赛跑，就会比较明显的感受到阻力。而当您骑着摩托车在公路上飙车，更是能亲身体会到空气阻力的“强大”。 所以不妨先将阻力设成速度的函数，即令： （3.1） 然后我们利用Taylor级数将其理论上展开成为如下的形式： (3.2） 取到v的多少次方取决于速度的大小，就像我们刚才讨论的那样，走路的时候就一项都不取，跑步时取到,飙车时取到，要是要讨论一架飞机或者一个UFO，就一定要取到了。所以对这个冰雹的模型，既然题目第二题让我们讨论到了速度的平方项，故不妨将阻力设为： ； （3.3） 则我们的微分方程变为： （3.4） 不幸的是，它看起来更难求解了，所以我们还是只能用数值计算的方法求出数值解，取m=0.1kg,k1=0.1,k2=0.05. 数值解的结果如图表 2所示 图表 2 附录 附录1：数值计算的两个matlab程序： function xdot = diffq(t,x) /*第二个微分方程*/ xdot=9.8-0.05*x.^2; function xdot = difeq(t,x) /*改进后的微分方程*/ xdot=9.8-0.05*x.^2-0.1*x; [t,x]=ode23(diffq,[t0,tf],x0); /用低阶法求微分方程二的数值解/ plot(t,x) [t,x]=ode23(difeq,[t0,tf],x0); /用低阶法求微分方程（改进）的数值解/ plot(t,x)  住房的合理定价 摘要 本文一开始根据历年住房的平均价格建立了房价与时间的函数关系，并根据这个关系预测了2010的房价。之后根据各组数据分别建立了两个模型并讨论了房价与GDP的关系。最后根据题目的要求，在基本假设及多个条件的约束下建立了2010房价的最优化模型。 问题一中