数学建模作业姓名：王士彬学院：计算机科学与技术班级：2014级计科2班学号：2014001300701.在区域x∈[-2,2],y∈[-2,3]内绘制函数z=exp^(-x2-y2)曲面图及等值线图。

解：曲面图如下：>> x=-2:0.5:2;>> y=-2:0.5:3;>> [X,Y]=meshgrid(x,y);>> Z=exp(-X.^2-``Y.^2);>> mesh(X,Y,Z)>>等值线图如下：>> x=-2:0.5:2;>> y=-2:0.5:3;>> [X,Y]=meshgrid(x,y);>> Z=exp(-X.^2-Y.^2);>> mesh(X,Y,Z)>> surf(X,Y,Z)>> surf(X,Y,Z) >> contour(X,Y,Z) >>2.已知一组观测数据，如表1所示.(1)试用差值方法绘制出x ∈[-2,4.9]区间内的光滑曲线，并比较各种差值算法的优劣.(2)试用最小二乘多项式拟合的方法拟合表中的数据，选择一个能较好拟合数据点的多项式的阶次，给出相应多项式的系数和偏差平方和.(3)若表中数据满足正态分布函数222/)(21)(σμπσ--=x e x y .试用最小二乘非线性拟合的方法求出分布参数σμ,值，并利用锁求参数值绘制拟合曲线，观察拟合效果.解：(1)分别用最领近插值，分段线性插值（缺省值），分段三次样条插值，保形分段三次插值方法绘制在x ∈[-2,4.9]的光滑曲线，图形如下：样条插值效果最好，其次线性插值，最近点插值效果最差，在这里效果好像不太明显。

最近点插值优点就是速度快，线性插值速度稍微慢一点，但效果好不少。

所以线性插值是个不错的折中方法。

样条插值，它的目的是试图让插值的曲线显得更平滑，为了这个目的，它们不得不利用到周围若干范围内的点，不过计算显然要比前两种大许多。

MATLAB 文件如下：>> x0=-2:0.3:4.9;>> y0=[0.10289 0.11741 0.13158 0.14483 0.15656 0.16622 0.17332 0.17750 0.17853 ... 0.17635 0.17109 0.16302 0.15255 0.1402 0.12655 0.11219 0.09768 0.08353 ... 0.07015 0.05876 0.04687 0.03729 0.02914 0.02236]; >> cx=-2:0.3:4.9;>> y1=interp1(cx,y0,cx,'nearest'); >> y2=interp1(cx,y0,cx,'linear'); >> y3=interp1(cx,y0,cx,'spline'); >> y4=interp1(cx,y0,cx,'cubic');>> subplot(2,2,1),plot(cx,y0,'o',cx,y1,'-r'),title('Nearest Interpolant'); >> subplot(2,2,2),plot(cx,y0,'o',cx,y1,'-k'),title('Linear Interpolant'); >> subplot(2,2,3),plot(cx,y0,'o',cx,y1,'-b'),title('Spline Interpolant'); >> subplot(2,2,4),plot(cx,y0,'o',cx,y1,'-k'),title('Cubic Interpolant'); >> subplot(2,2,1),plot(cx,y0,'o',cx,y1,'-r'),title('Nearest Interpolant');(2),从图形可以看出曲线函数遵从幂函数的形式，设幂函数形式为：βαx y =可化为.ln ln ln x y βα+=即把非线性函数转化为线性函数，原线性函数形式为01)(a x a x p +=由此我们可以得出p(x)等价于lny;x 等价于lnx;1a =β，0ln a =α 我们可以先求出0,1a a 。

求一个线性多项式01)(a x a x p +=使之在最小二乘准则下拟合这些观测值，问题即化为求1,0a a 使E(1,0a a)=21011,0)]([min ∑=+-m i i i a a a x a y 利用多元函数极值原理可知，若目标函数E(1,0a a)的极小值存在，一定有,00=a E ϑϑ01=a Eϑϑ，用MA TLAB 工具我们可以求得最后的结果。

>> log(x0); >> log(y0); >> x0=log(x0); >> y0=log(y0); >> n=length(x0); >> a=sum(x0); >> b=sum(y0); >> c=sum(x0.*y0); >> d=sum(x0.^2);>> a0=(d*b-c*a)*(n*d-a^2); >> a1=(n*c-a*b)/(n*d-a^2); >> a0,a1 a0 =-2.5891e+05 - 1.7515e+06ia1 =0.1045 - 0.3558i即系数a0为 -2.5891e+05 - 1.7515e+06i ，a1为0.1045 - 0.3558i其相应多项式的系数和偏差平方和.我们可以求出E= -7.2019e+13 + 2.1767e+13i 其MA TLAB 文件如下： >> Y=a1*x0+a0; >> e=Y-y0; >> E=sum(e.^2) E =-7.2019e+13 + 2.1767e+13i即其相应多项式的系数和偏差平方和. 为 -7.2019e+13 + 2.1767e+13i (3)？3.将某物体放置在空气中，在t=0时刻测得其温度u 0=150度，10min 后测得温度u 1=87度，假设空气的温度为24度。

试建立数学模型给出物体的温度u 与时间t 的关系，并计算20min 后物体的温度。

解：为了解决上述问题，我们首先需要了解有关热力学的一些基本规律：比如：热量总是从温度高的物体向温度低的物体传导的；在一定的温度范围（其中包括了上述问题的温度在内），一个物体的温度与这物体的温度和其所在介质的温度的差值成正比例。

这是已为实验证明了的牛顿冷却定律。

设空气的温度为u a ,物体在时刻t 的温度为)(t u u =，则温度的变化速度为dtdu。

注意热量总是从温度高的物体向温度低的物体传导的，因而初始温度大于空气温度，即（u 0>u a ），所以温差u-u a 恒正；又因为物体的温度将随时间而逐渐冷却，故温度变化速度dtdu恒负。

因此，由牛顿冷却定律得到)(ua u K dtdu--=............(1) 这里的K>0是比例常数。

此（1）方程就是冷却过程的数学模型。

为了确定温度u 与时间t 的关系，我们需要从上面（1）的方程中解出u 。

又因为ua 是常数，并且u-ua>0,所以我们可以将上述式子改写成Kdt uau ua u d -=--)( 将此式积分可得到如下式子1)ln(c Kt ua u +-=-)(^)1(^Kt ce c Kt e ua u -=+-=-即u=u a +ce^(-Kt)根据初始条件：t=0时，u=u 0代入上式得 c=u 0-u a于是u=u 0+(u 0-u a )e^(-Kt)又根据条件，当t=10时，u=u 1代入上式得u 1=u a +(u 0-u a )e^(-10K) ⇒ln 101=K [(u 0-u a )/(u 1-u a )]根据题意我们可知u 0=150，u 1=87,u a =24,代入得到K=248724150ln 101--=2ln 101=0.069 从而u=24+126e^(-0.069t)这就是物体冷却时温度u 随着时间t 的变化规律。

用t=20代入得u=55.7度4.假设在某商场中，某种商品在t 时刻的价格为P(t)，若假定其变化率与商品的需求量D 和供给量S 之差成正比(比例系数为k),若dP c S bP a D +-=-=,其中d c b a ,,,均为正常数，若已知初始价格为Po,求任意时刻t 时该商品的价格。

解：一般情况下，某种商品的价格主要服从市场供求关系，由题意我们可知商品需求量D 是价格P 的单调递减函数，商品供给量S 是价格P 的单调递增函数，即dP c S bP a D +-=-=,------------------------------------------------------------------(1)其中d c b a ,,,均为常数，且b>0,d>0.当需求量与供给量相等时，由(1)可得供求平衡时的价格P e =db ca ++，并称P e为均衡价格。

由题意得：)]()([p S p D k dtdp-= 其中比例系数k>0，用来反应价格的调整进度。

将(1)式代入方程可得其中常数=k(b+d)λ>0，所以此方程的通解为P(t)=P e +Ce^(-λt)由于初始价格P(0)=P 0代入上式，得C=P 0-P e 于是我们可以求出任意时刻价格P 与时刻t 之间的函数为：P(t)=P e +(P 0-P e )^(-λt)，并且我们可以得出，因为λ>0知，+∞→t 时P(t)→P e ，说明随着时间的不断推延，实际价格P(t)将逐渐趋近均衡价格P e 。

5.农场种植计划问题某农场根据土地的肥沃程度，把耕地分为I II III 三等，相应的耕地面积分别为100、300和200km 2,计划种植水稻、大豆和玉米.要求三种作物的最低收获量分别为190、130和350吨(t).I 、 II 、III 等耕地种植三种作物的单产如表所示.若三种作物的售价分别为水稻1.2元/kg,大豆1.50元/kg ，玉米0.80元/kg.那么(1)如何制订种植计划，才能使总产量最大？ (2)如何制订种植计划，才能使总产值最大？解：(1):•问题分析：确定种植最佳土地分配，即每种等级耕地分别种植水稻、大豆、玉米的面积•模型建立： 1，决策变量：令x1,x2,x3分别为I II III 三等耕地上种植的水稻面积，令x4,x5,x6分别为I II III 三等耕地上种植的大豆面积，令x7,x8,x9分别为I II III 三等耕地上种植的玉米面积。