数学实验报告

1、高等数学数学实验报告实验人员：院（系）_ _ 学号 _姓名 _ _实验地点： 计算机中心机房实验一空间曲线与曲面的绘制一、实验题目：（实验习题 1-2 ）利用参数方程作图，做出由下列曲面所围成的立体图形：(1)z1 x 2y2 , x2y 2x 及xOy平面；(2)zxy , xy 1 0 及 z0.二、实验目的和意义1、利用数学软件 Mathematica 绘制三维图形来观察空间曲线和空间曲面图形的特点，以加强几何的直观性。2、学会用 Mathematica 绘制空间立体图形。三、程序设计空间曲面的绘制xx(u, v)yy(u,v), uumin , max, v vmin ,vmax 作参

2、数方程zz(u,v)所确定的曲面图形的Mathematica 命令为：ParametricPlot3Dxu,v,yu,v,zu,v,u,umin,umax,v,vmin,vmax,选项 (1)(2)四、程序运行结果(1)（ 2）五、结果的讨论和分析1、通过参数方程的方法做出的图形，可以比较完整的显示出空间中的曲面和立体图形。2、可以通过 mathematica 软件作出多重积分的积分区域，使积分能够较直观的被观察。3、从 (1) 中的实验结果可以看出，所围成的立体图形是球面和圆柱面所围成的立体空间。4、从 (2) 中的实验结果可以看出围成的立体图形的上面曲面的方程是z xy ，下底面的方程是

3、z=0，右边的平面是 x y1 0 。实验一 空间曲线与曲面的绘制一、实验题目：（实验习题1-3 ）观察二次曲面族 zx 2y 2kxy 的图形。特别注意确定k 的这样一些值，当 k经过这些值时，曲面从一种类型变成了另一种类型。二、实验目的和意义1.学会利用Mathematica 软件绘制三维图形来观察空间曲线和空间曲线图形的特点。2. 学会通过表达式辨别不同类型的曲线。三、程序设计这里为了更好地分辨出曲线的类型，我们采用题目中曲线的参数方程来画图，即 zr 2kr 2 cos t sin t输入代码：ParametricPlot3Dr*Cost，r*Sint,r2+ k*r2*Cost*Si

4、nt,t, 0, 2*Pi, r, 0, 1，PlotPoints -> 30式中k 选择不同的值：-4到 4 的整数带入。四、程序运行结果k=4：k=3：k=2：k=1:k=0:k=-1:k=-2:k=-3:k=-4:五、结果的讨论和分析k 取不同值，得到不同的图形。我们发现，当|k|<2时，曲面为椭圆抛物面；当|k|=2 时，曲面为抛物柱面；当 |k|>2 时，曲面为双曲抛物面。实验二 无穷级数与函数逼近一、实验题目：（实验习题 2-2 ）改变例 2 中 m及 x0 的数值来求函数的幂级数及观察其幂级数逼近函数的情况。二、实验目的和意义1. 利用 Mathematica

5、显示级数部分和的变化趋势。2.学会如何利用幂级数的部分和对函数进行逼近以及函数值的近似计算。三、程序设计若函数 f (x)(1x) m 能展开成 x- x0 的幂级数（这里不验证），则根据函数展开为幂级数的展开公式，其展开式为 f ( x)f ( n) ( x0 ) (xx0 ) n 。因此首先定义 f ( x) 的 n 阶导数的函n 0n!数 g(n,x0 ) ，最后再构成和式即得f (x) 的幂级数展开式。用Mathematica 观察幂级数部分和逼近函数的情况。m=2， x0 =2 时输入如下命令：m=- 2 ;f x _ : =( 1+ x ) m; x 0= 2 ;g n _ , x

6、 0 _ : = D f x , x , n /.xx 0 ;s n _ , x _ : = Sum gk , x0k !* ( x - x 0 ) k , k , 0 , n ;t =Ta bl e s n , x , n , 20 ;p 1= Pl o t Ev a l u at e t , x , - 1 / 2 , 1/ 2 ;p 2= Pl o t ( 1 + x ) m, x , - 1/ 2,1 /2 , Pl ot St y l e RGBCo l o r 0, 0, 1 ;Sho w p 1,p2四、程序运行结果从输出的图形观察f ( x) 展开的幂级数的部分和逼近函数f (

7、x) 的情况：五、结果的讨论和分析从图中可以看到，当n 越大时，幂级数越逼近函数。实验二无穷级数与函数逼近一、实验题目：（实验习题 2-3 ）观察函数 f ( x )x ,x 0 展成的傅里叶级数的部分和逼近f ( x ) 的情况。1,0x二、实验目的和意义1. 利用 Mathematica 显示级数部分和的变化趋势。2. 学会展示傅里叶级数对周期函数的逼近情况。三、计算公式f ( x ) 可以展开成傅里叶级数： a0( an cos nx bnsin nx ) ，其中2n 1ak11kxdx( k0,1,2, )f ( x ) cos kxdx ( k0,1,2, ) ， bkf ( x )

8、 sin四、程序设计输入代码：fx_ := Which-Pi <= x < 0, -x, 0 <= x < Pi, 1; an_ := Integrate-x*Cosn*x, x, -Pi, 0/Pi +IntegrateCosn*x, x, 0, Pi/Pi;bn_ := Integrate-x*Sinn*x, x, -Pi, 0/Pi + IntegrateSinn*x, x, 0, Pi/Pi;sx_, n_ :=a0/2+Sumak*Cosk*x + bk*Sink*x, k, 1, n; g1 = Plotfx, x, -2Pi, 2Pi, PlotStyle

9、 -> RGBColor0, 0, 1,DisplayFunction -> Identity; m = 18;Fori = 1, i <= m, i += 2,g2 = PlotEvaluatesx, i, x, -Pi, Pi, DisplayFunction -> Identity;Showg1, g2, DisplayFunction -> $DisplayFunction五、程序运行结果六、结果的讨论和分析从图表可以看出， n 越大逼近函数的效果越好， 还可以注意到傅里叶级数的逼近是整体性的。实验三最小二乘法一、实验题目： ( 实验习题 3-2)一种合金

10、在某种添加剂的不同浓度下进行实验，得到如下数据：浓度 x10.015.020.025.030.0抗压强度 y27.026.826.526.326.1已知函数 y 与 x 的关系适合模型： ya bx cx 2 ，试用最小二乘法确定系数a，b，c，并求出拟合曲线。二、实验目的和意义1. 学会利用最小二乘法求拟合曲线。2. 学会画数据点的散点图及拟合函数的图形，并将两个图画在同一坐标下。三、计算公式n( a bxi cx i2 ) y i 2 取最小值，令此函数对根据最小二乘法， 要求 Q( a, b, c )i1各个参数的偏导等于0，解 n+1 元的方程组便可求得这些参数的最小二乘解。四、程序设

11、计输入代码：x = Table10.0 + 5.0*i, i, 0, 4;y = 27.0, 26.8, 26.5, 26.3, 26.1;xy = Tablexi, yi, i, 1, 5;qa_, b_, c_ := Sum(a + b*xi + c*xi2 - yi)2, i, 1, 5 NSolveDqa, b, c, a = 0, Dqa, b, c, b = 0, Dqa, b, c, c = 0, a, b, ct1 = ListPlotxy, PlotStyle -> PointSize0.02,DisplayFunction -> Identity;fx_ :=

12、27.56 + -0.0574286*x + 0.000285714*x2; t2 = Plotfx, x, 5, 35, AxesOrigin -> 5, 25, DisplayFunction -> Identity;Showt1, t2, DisplayFunction -> $DisplayFunction五、程序运行结果首先得到 a，b，c 三个值：a -> 27.56, b -> -0.0574286, c -> 0.000285714然后得到同一坐标系下的数据点散点图及拟合函数的图形：六、结果的讨论和分析观察 a，b，c 的值以及图像可以发现，

13、二次方项的系数非常小，而所得的图像也非常接近于直线。实验三最小二乘法一、实验题目： ( 实验习题 3-3)在研究化学反应速度时，得到下列数据：369121518212457.641.931.022.716.612.28.96.5其中 xi 表示实验中作记录的时间，yi 表示在相应时刻反应混合物中物质的量，试根据这些数据建立经验公式。二、实验目的和意义1. 学会利用最小二乘法求拟合曲线。2. 学会由实际经验或相关的学科理论，能够提供拟合函数的可取类型，通过适当的变量代换将拟合函数线性化，建立经验公式。三、计算公式在许多场合下，拟合函数不具有线性形式，但是由实际经验或相关的学科理论，能够提供拟合函

14、数的可取类型，而且可以通过适当的变量代换将拟合函数线性化，同样可以建立经验公式。模型 yaebx 可以用变量替换Yln y, Xx 将函数化为线性函数：Yln abX 。四、程序设计输入代码：（1）生成数据并作图观察t1=3,6,9,12,15,18,21,24;y1=57.6,41.9,31.0,22.7,16.6,12.2,8.9,6.5;data1=Transposet1,y1;d2=ListPlotdata1,PlotStyle->RGBColor0,0,1,PointSize0.02;（2）确定回归函数的类型logy=Logy1;data2=Transposet1,logy;d3=ListPlotdata2,PlotStyle-&g