实验库1：一元函数微分学2

1、实验一 一元函数微积分学实验2 极限与连续（基础实验）实验目的 通过计算与作图, 从直观上揭示极限的本质,加深对极限概念的理解. 掌握用Mathematica画散点图, 以及计算极限的方法. 深入理解函数连续的概念,熟悉几种间断点的图形特征,理解闭区间上连续函数的几个重要性质.1、作散点图例2.1 (教材 例2.1) 分别画出坐标为的散点图, 并画出折线图.分别输入命令t1=Tablei2,i,10; g1=ListPlott1,PlotStyle->PointSize0.02;g2=ListPlott1,PlotJoined->True;Showg1,g2;t2=Tablei2,

2、4i2+i3,i,10;g1=ListPlott2,PlotStyle->PointSize0.02;g2=ListPlott2,PlotJoined->True;Showg1,g2;则分别输出所求图形.例2.2 画出前25个素数的散点图.输入命令TablePrimen,n,25;ListPlotTablePrimen,n,25,PlotStyle->PointSize0.015;则分别输出所求图形.2、数列极限的概念例2.3 观察数列的前100项变化趋势.输入命令t=NTablen(1/n),n,1,100;ListPlott,PlotStyle->PointSize

3、0.015;则分别输出所求图形. 从图中可看出, 这个数列似乎收敛于1. 下面我们以数值的方式来说明这一变化趋势. 输入以下语句, 并观察其数值结果.m=2;xn=0;Fori=1,i<=1000,i+=50,IfAbsxn-1>10(-m),xn=Nn(1/n),20;Printi, " ",xn;设该数列收敛于不妨取下面考察与A的接近程度. 输入以下Mathematica语句.u = 109(-2); A = 1 + u; m = 5; n = 3; an = Sqrt3;WhileAbsA-an >= 10(-m), n+; an = Nn(1/n)

4、;Print" n=", n, " an=", an, "|A-an|=", AbsA - an;结果表明: 当时, 与的距离小于例2.4 观察Fibonacci数列的变化趋势.Fibonacci数列具有递推关系令.输入命令fn1=1;fn2=1;rn=1;Fori=3,i<=14,i+,Fn=fn2+fn1;fn2=fn1;fn1=fn;rn=Nfn2/fn1,20;dn=rn-rn1;rn1=rn;Printi, " ", fn1, " ", rn, " ",dn;

5、其中第二列给出了Fibonacci数列的前14项, 第3列给出了的值, 由第4列可以看出, 我们也可以用散点图来观察Fibonacci数列的变化趋势如图所示, 输入命令Clearf; fn_:= fn-1+fn-2;f0=1;f1=1;fab20=Tablefi,i,0,20;ListPlotfab20,PlotStyle ->PointSize0.02;Infab20=Logfab20;ListPlotInfab20,PlotStyle-> PointSize0.02;则输出所求散点图.为了更好地观察数列的变化趋势, 我们可以利用Mathematica的动画功能来进一步观察数列随

6、着n的增大的变化趋势.例2.5 通过动画观察当时数列的变化趋势.输入Cleartt;tt=1,1/22,1/32;Dott=Appendtt,N1/i2;ListPlottt,PlotRange->0,1,PlotStyle->PointSize0.02,i,4,20则输出所求图形动画. 从图中可以看出所画出的点逐渐接近于x轴.例2.6 研究极限输入Printn, " ", Ai, " ",0.4-Ai;Fori=1, i<=15, i+,Aii=N(2 i3+1)/(5 i3+1),10;Bii=0.4-Aii; Printi, &q

7、uot; ", Aii, " ", Bii则输出nAi0.4-Ai10.50.120.4146340.014634130.4044120.0044117640.4018690.0018691650.4009580.00095846660.4005550.00055504270.400350.0003496580.4002340.00023428390.4001650.000164564100.400120.000119976110.400090.0000901442120.4000690.0000694364130.4000550.000054615140.400

8、0440.0000437286150.4000360.0000355534观察所得数表. 第一列是下标n. 第二列是数列的第n项它与0.4越来越接近. 第三列是数列的极限0.4与数列的项的差, 逐渐接近0.再输入fn=Table(2 n3+1)/(5 n3+1),n,15;ListPlotfn,PlotStyle->PointSize0.02则输出散点图. 观察所得散点图, 可见表示数列的点逐渐接近于直线注:命令For的格式见项目二中实验1的基本命令.3、函数的极限例2.7 在区间上作出函数的图形, 并研究 和 输入命令Clearf;fx_=(x3-9x)/(x3-x);Plotfx,x

9、,-4,4;则输出的图形. 从图可猜测 不存在.例2.8 观察函数当时的变化趋势.取一个较小的区间1, 10, 输入命令fx_=Sinx/x2;Plotfx,x,1,20;则输出在这一区间上的图形. 从图中可以看出图形逐渐趋于0. 事实上, 逐次取更大的区间, 可以更有力地说明当时, 作动画: 分别取区间画出函数的图形, 输入以下命令:i=3;Whilei<=20,Plotfx,x,10,5*i,PlotRange->10,100,-0.008,0.004;i+则输出17幅图, 点黑右边的线框, 并选择从前向后的播放方式播放这些图形, 可得函数当时变化趋势的动画, 从而可以更好地理

10、解此时函数的变化趋势.例2.9 考虑函数 输入PlotArcTanx,x,-50,50则输出该函数的图形. 观察当时, 函数值的变化趋势. 分别输入LimitArcTanx,x->Infinity,Direction->+1LimitArcTanx,x->Infinity,Direction->-1输出分别为与考虑函数分别输入LimitSignx,x->0,Direction->+1LimitSignx,x->0,Direction->-1输出分别为-1与1.4、两个重要极限例2.10 考虑第一个重要极限输入PlotSinx/x,x,-Pi,Pi

11、则输出函数的图形. 观察图中当时, 函数值的变化趋势. 输入LimitSinx/x,x->0输出为1, 结论与图形一致.例2.11 研究第二个重要极限 输入Limit(1+1/n)n,n->Infinity输出为e. 再输入Plot(1+1/x)x,x,1,100则输出函数的图形. 观察图中函数的单调性. 理解第二个重要极限4、无穷大例2.12 考虑无穷大. 分别输入Plot(1+2 x)/(1-x),x,-3,4Plotx3-x,x,-20,20则分别输出两个给定函数的图形. 在第一个函数的图形中,时函数的绝对值无限增大,在第二个函数的图形中,时函数的绝对值在无限增大. 输入Li

12、mit(1+2x)/(1-x),x->1Mathematica输出的是. 这个结果应该是右极限.例2.13 考虑单侧无穷大. 分别输入PlotE(1/x),x,-20,20,PlotRange->-1,4LimitE(1/x),x->0,Direction->+1LimitE(1/x),x->0,Direction->-1则输出所给函数的图形、左极限0和右极限值. 再输入LimitE(1/x),x->0Mathematica的输出仍然为.这又是右极限（同上例）. 因此在没有指明是左右极限时, 命令Limit给出的是右极限.例2.14 输入Plotx+4

13、*Sinx,x,0,20 Pi则输出所给函数的图形. 观察函数值的变化趋势. 当时, 这个函数是无穷大. 但是, 它并不是单调增加. 于是, 无穷大并不要求函数单调.例2.15 输入Plotx*Sinx,x,0,20 Pi则输出所给函数的图形.观察图中函数的变化趋势. 这个函数无界, 但是, 当时, 这个函数不是无穷大. 即趋向于无穷大的函数当然无界, 而无界函数并不一定是无穷大.5、连续与间断例2.16 考察函数在处的连续性.选取几个考察当时, 的变化趋势, 依次取当时, 他们的极限均为5.输入命令g1 = ListPlotTableSin5 + 1/n, n, 1, 1000, 5,Plo

14、tStyle -> RGBColor1, 0, 0;g2 = ListPlotTableSin5 + (-1)n/Sqrtn, n, 1, 1000, 5,PlotStyle -> RGBColor0, 1, 0;g3 = ListPlotTableSin5*n*Log(1 + 1/n), n, 1, 1000, 5,PlotStyle -> RGBColor0, 0, 1;g = Showg1, g2, g3;则输出相应的的散点图. 由图可看出它们趋于同一极限值.例2.17 观察可去间断. 分别输入PlotTanx/x,x,-1,1Plot(Sinx-x)/x2,x,-Pi

15、,Pi则输出所给函数的图形. 从图可见，是所给函数的可去间断点.例2.18 观察跳跃间断. 分别输入PlotSignx,x,-2,2Plot(E(1/x)-1)/(E(1/x)+1),x,-2,2则分别输出所给函数的图形. 从图可见，是所给函数的跳跃间断点.例2.19 观察无穷间断. 分别输入Plot1/(1-x2),x,-3,3则输出所给函数的图形. 从图可见，是所给函数的跳跃间断点.例2.20 观察振荡间断. 分别输入PlotCos1/x,x,-Pi,Pi则输出所给函数的图形. 从图可见，是所给函数的跳跃间断点. 再输入LimitSin1/x,x->0Mathematica4.0输出为Interval-1,1. 读者可猜测这是什么意思.例2.21 有界量乘以无穷小. 分别输入Plotx*Sin1/x,x,-Pi,PiLimi