实验一 Mathematica系统使用入门 【实验目的】 熟悉 Mathematica的一些基本命令的格式与作用并能熟练使用这些命令 【基本命令】 1帮助命令 命令格式 页码 作用注意事项 ?大写字母* 显示以该字母开头的所有命令 该命令的基本格式 ??命令 该命令的详细介绍 2系统内部常数E→e,Pi→, infinity→,1→1 3代数运算命令 匚命令格式 作用 注意事项 Na, P3求数值a的近似值(保留n位 有效数字) Nal 「求数值a的保留6位有效数字也可写为/N 的近似值 P3将数值a赋予 使用后注意及时清除 Clearly 清除对ⅹ的赋值 Clear x,y,z…]可一次 清除多个变量的赋值 表达式{x>a}P3将表达式中的x换成数值a后1{xay>b…}可代 的结果 入多个变量的值 Factor[表达式]p对表达式进行因式分解 也可写为/ Factor Expandl表达式]|P4展开表达式 也可写为/ xpand Smp表达式P4化简表达式 也可写为/ Simplify 例1求代数式x2+3x-10在x=2处的值 解:按序号依次输入下列10个语句观察每个语句执行后的结果,认真体会赋值语句与代入 规则的不同 命令 结果 命令 结果 1.x^2+3x-10 x^2+3x-10 3.x^2+3x-10 4. 5.x*y 6.x=或 Clearly 7.x^2+3x-10x^2+3x-10 8.x^2+3x-10/,x->2 0 9.x^2+3x-10 x^2+3x-10 x^3-1 x^3-1 X*y 12.x*y. 6 3解方程与方程组命令 命令格式 注意事项 Solve[方程,x P4求解以x为自变量的代数方程,只能求解代数方 也可用于求解方程组 程 SOlve[方程x]|P4求以x为自变量的代数方程解的只能求解代数方 近似值,也可用于求解方程组程 Reduce[方程x|P4求解以x为自变量的代数方程,只能求解代数方 也可用于求解方程组 Eliminata[{方程|P5在给定的方程组中消去x 不能用这一命令 求方程组的解 4求和命令
实验一 Mathemetica 系统使用入门 【实验目的】 熟悉 Mathemetica 的一些基本命令的格式与作用并能熟练使用这些命令。 【基本命令】 1 帮助命令 命令格式 页码 作用 注意事项 ？大写字母* 显示以该字母开头的所有命令 ？命令 该命令的基本格式 ？？命令 该命令的详细介绍 2.系统内部常数 E →e , Pi → , Inifinity → , I → i 3.代数运算命令 命令格式 页码 作用 注意事项 N[a,n] P3 求数值 a 的近似值（保留 n 位 有效数字） N[a] 求数值a的保留6位有效数字 的近似值 也可写为//N x=a P3 将数值 a 赋予 x 使用后注意及时清除 Clear[x] x =. P3 清除对 x 的赋值 Clear[x,y,z,…]可一次 清除多个变量的赋值 表达式/.{x->a} P3 将表达式中的x换成数值 a后 的结果 /.{x->a,y->b,…}可代 入多个变量的值 Factor[表达式] P4 对表达式进行因式分解 也可写为//Factor Expand[表达式] P4 展开表达式 也可写为//Expand Simplify[表达式] P4 化简表达式 也可写为//Simplify 例 1 求代数式 3 10 2 x + x − 在 x = 2 处的值. 解: 按序号依次输入下列 10 个语句,观察每个语句执行后的结果,认真体会赋值语句与代入 规则的不同. 命令 结果 命令 结果 1. x^2+3x-10 x^2+3x-10 2. x=2 2 3. x^2+3x-10 0 4. x^3-1 7 5. x*y 2y 6. x=.或 Clear[x] 7. x^2+3x-10 x^2+3x-10 8. x^2+3x-10/.x->2 0 9. x^2+3x-10 x^2+3x-10 10. x^3-1 x^3-1 11. x*y x*y 12. x*y/.{x->2,y->3} 6 3 解方程与方程组命令 命令格式 页码 作用 注意事项 Solve[方程,x] P4 求解以 x 为自变量的代数方程， 也可用于求解方程组 只能求解代数方 程 NSolve[方程,x] P4 求以 x 为自变量的代数方程解的 近似值，也可用于求解方程组 只能求解代数方 程 Reduce[方程,x] P4 求解以 x 为自变量的代数方程， 也可用于求解方程组 只能求解代数方 程 Eliminata[{方程 组},x] P5 在给定的方程组中消去 x. 不能用这一命令 求方程组的解 4 求和命令
命令格式 页码 作用 注意事项 Sum["n, in, 1,n2, di, 求通项为n,n从到当1=1,d=1时 可省略,即 步长为d的和的精 确值 ,{n,h2 P5求上式的近似值 n1,n2 5作图命令 命令格式 码 注意事项 定义函数fx) Plot[f[x],(x, a, b1 P7 画出函数f(x)在区间 ab]上的图形 Parametric Plot[(x[t], y[]),P8 画出参数方程 h1,2 x=x(1) =y()在[,12l上 的图形 ListPlot[list] P9|画出以所给表lst为坐该命令中的选项 标的点的散点图 Plotstyle->True可 画出折线图 Showel,g2,g3,…] P9将g1,g2;g3…等图形组 6关于数表的命令 命令格式 码 作用 注意事项 P10 Table[un, (n, 1, n2, di)1 生成以为通项的数表B=di=1时 可省略 P10生成一个有n+1个元素的 NestList[f(x], 0, n 数表,表中第个元素为以 xo为初值,以x为迭代 函数,迭代n-1次的值 list[n) P10取出表lst中的第n个元 素 list[[1, m2n1 取出二层表lst中第h个 元素中的第个元素 Take[list, n P11n>0时,取出表list中前n 个元素,n<0时,取出表
命令格式 页码 作用 注意事项 Sum[ n a ,{n, 1 2 n , n ,di}] P5 求通项为 n a ，n 从 1 n 到 2 n ，步长为 di 的和的精 确值 当 1 n =1,di=1 时 可省略，即 Sum[ n a ,{n, 2 n }] NSum[ n a ,{n, 1 2 n , n ,di}] ] P5 求上式的近似值 5 作图命令 命令格式 页码 作用 注意事项 f[x_]:= P7 定义函数 f(x) Plot[f[x],{x, a, b}] P7 画出函数 f(x)在区间 [a,b]上的图形 ParametricPlot[{x[t],y[t]}, {t, 1 2 t ,t }] P8 画出参数方程 = = ( ) ( ) y y t x x t 在 [ , ] 1 2 t t 上 的图形 ListPlot[list] P9 画出以所给表 list 为坐 标的点的散点图 该命令中的选项 PlotStyle->True 可 画出折线图 Show[g1,g2,g3, …] P9 将 g1,g2,g3…等图形组 合显示在一张图中 6 关于数表的命令 命令格式 页码 作用 注意事项 Table[ n a ,{n, 1 2 n ,n ,di}] P10 生成以 n a 为通项的数表 n1 = 1 ,di=1 时 可省略 NestList[f[x], 0 x ,n] P10 生成一个有n+1个元素的 数表,表中第个元素为以 0 x 为初值，以 f(x)为迭代 函数，迭代 n-1 次的值 list[[n]] P10 取出表 list 中的第 n 个元 素 list[[ 1 2 n ,n ]] P11 取出二层表 list 中第 1 n 个 元素中的第 2 n 个元素 Take[list,n] P11 n>0 时，取出表 list 中前 n 个元素，n<0 时，取出表
中后n个元素 ength list P11求表ls的长度(表中元 素的个数) P将a加在表的最前面 Pll 加在表list的最后面 rt[list, a, n] P1将a加在表ls的第n个 位置 Transpose[ list] Pl6将二层表的元素进行转 例2例2设数列求由该数列的前10项构成的数表(表中元素保留4位有效数字) n: tt=Table(N[1/(2n-1),4 n, 1, 10, 1) 1,0.30.2,0.1429,0.111009091,0.07692,0.0670.05882,0.05263} 命令 结果 ttl4ll 0.1429 Length tt 分别输入以下语句,并观察结果 Prepend tt, oI Append[tt, x 2+y 21 Insert tt, ** 3 tt=Table[1/(2n-1), 4,n, 10fI 例3求函数f(x)=e-h1+x)+x在x=0,x=02,x=04…x=1.8x=2处的近似 值将它们与x的对应值作成一个二层表并画出相应的图形 解:nx1:=E^x-Log+x]+x ff-TableRx, fx //N, x, 0, 2, 0.231 0,1.},{0.2,1.23908},{0.4,1.55535},{06,1,95212},{0.8,2.43775},{1,3.02513},{1.2, 373166},{1.4,4.57973},{1.6,5.59752},{1.8,6.82003},{2.,8.29044} ListPlot[ff, PlotJoined->Truel gg=Transpose ff 0,0.2,0.4,0.6,0.8,1,1.2,14,1.6,1,8,2.},{1,1.23908,1.55535,1.95212,2.43775 3.02513,3.73166457973,5.59752,6.82003,829044} ggll2ll 1,1.23908,1.55535,1.95212,2.43775,3.02513,3.73166,4.57973,5.59752,682003, 829044} Transposelggl 10,1.},{0.2,1,23908},{0.4,1.55535},{0.6,1,95212},{0.8,243775},{1,3.02513},{1 373166},{14,457973},{1.6,5.59752},{1.8,6.82003},{2.,8.29044} 7条件语句与循环语句 命令格式 页码 作用 注意事项 条件,语句1,语句2]P12条件成立时执行语句 1,不成立时执行语句2 Whch条件1语句1,条P2条件1成立时执行语 件2,语句2…条件n,语句 句1条件2成立时执行语 句2…条件n成立时执行 语句n P3循环语句 Fo=5,i<=n=i+32,循 i=+S2,可记为 环体语句] 计=52,。特别的 i=i+1也可记为
中后 n 个元素 Length[list] P11 求表 list 的长度（表中元 素的个数） Prepend[list,a] P11 将 a 加在表 list 的最前面 Append[list,a] P11 将 a 加在表 list 的最后面 Insert[list,a,n] P11 将 a 加在表 list 的第 n 个 位置 Transpose[list] P16 将二层表的元素进行转 换 例2 例 2 设数列,求由该数列的前 10 项构成的数表(表中元素保留 4 位有效数字) 解: tt=Table[N[1/(2n-1),4],{n,1,10,1}] {1., 0.3333, 0.2, 0.1429, 0.1111, 0.09091, 0.07692, 0.06667, 0.05882, 0.05263} 命令 结果 tt[[4]] 0.1429 Length[tt] 10 分别输入以下语句,并观察结果. Prepend[tt,0] Append[tt,x^2+y^2] Insert[tt,**,3] tt=Table[N[1/(2n-1),4],{n,10}] 例 3 求函数 f x e x x x ( ) = − ln(1+ ) + 在 x = 0, x = 0.2, x = 0.4, x = 1.8, x = 2 处的近似 值,将它们与 x 的对应值作成一个二层表,并画出相应的图形 解: f[x_]:=E^x-Log[1+x]+x; ff=Table[{x,f[x]//N},{x,0,2,0.2}] {{0, 1.}, {0.2, 1.23908}, {0.4, 1.55535}, {0.6, 1.95212}, {0.8, 2.43775}, {1,3.02513}, {1.2, 3.73166}, {1.4, 4.57973}, {1.6, 5.59752}, {1.8, 6.82003}, {2., 8.29044}} ListPlot[ff,PlotJoined->True] gg=Transpose[ff] {{0, 0.2, 0.4, 0.6, 0.8, 1., 1.2, 1.4, 1.6, 1.8, 2.}, {1., 1.23908, 1.55535, 1.95212, 2.43775, 3.02513, 3.73166, 4.57973, 5.59752, 6.82003, 8.29044}} gg[[2]] {1., 1.23908, 1.55535, 1.95212, 2.43775, 3.02513, 3.73166, 4.57973, 5.59752, 6.82003, 8.29044} Transpose[gg] {{0, 1.}, {0.2, 1.23908}, {0.4, 1.55535}, {0.6, 1.95212}, {0.8, 2.43775}, {1,3.02513}, {1.2, 3.73166}, {1.4, 4.57973}, {1.6, 5.59752}, {1.8, 6.82003}, {2., 8.29044}} 7 条件语句与循环语句 命令格式 页码 作用 注意事项 If[条件，语句 1，语句 2] P12 条件成立时执行语句 1，不成立时执行语句 2 Which[条件 1,语句 1,条 件 2,语句 2,…条件 n,语句 n] P12 条件 1 成立时执行语 句1,条件2成立时执行语 句2,…条件n成立时执行 语句 n For[i= 1 s ,i<=n,i=i+ 2 s ,循 环体语句] P13 循环语句 i=i+ 2 s ,可记为 i+= 2 s ,。特别的 i=i+1 也可记为
P13 循环语句 Do[循环体语句n,, H=1di=1时 可省 P13 循环语句 划线部分不可 n=h; While1
i++ Do[循环体语句,{n, 1 n , 2 n , di}] P13 循环语句 n1 = 1 ,di=1 时 可省略 n= 1 n ；While[n< 2 n ,循环 体语句；n=n+di] P13 循环语句 划线部分不可 省略，否则进入死 循环。 例 4 分别用 For,Do,While 语句打印 1—10 解： For[i=1,i<=10,i=i+1,Print[i]] Do[Print[i],{i,1,10,1}] i=1;While[i<=10,Print[i];i=i+1] 例 5 用 For 语句求 100 1 3 1 2 1 1+ + ++ 的近似值。 解： a=0;For[i=1,i<=100,i++,a=a+1/i];Print[a//N] 5.18738 【初学者容易出现的一些错误】 1. 该大写的地方不大写。 2. 括号用错：在 Mathemetica 中，大、中、小括号是有不同的用处的。大括号用于数表，中 括号用于命令或函数，小括号用于在表达式中改变运算顺序。 3. 缺少乘号*,两个字母相乘中间必须加 * 或空格，否则在运算中这两个字母作 为一 个整体处理。 4. 用 Solve,Nsolve,Reduce 等命令求解超越方程。 5. 没有及时清除变量的赋值。 6. 系统内部的函数名拼写或函数格式错误。 7. 注意标点符号，特别是逗号“，”与分号“；”的用法（参考教材 P13）。 【练习题】 1. 求下列各式的精确值与近似值(分别保留 4 位,6 位及 15 位有效数字) (1) 2 1 2 3 ) 7 3 ) 27 (3 2 1 − 2 ( + + − (2) arctan 3.2 ln 4 lg 2 1 ) 4 3 (sin 2 + + − − 2. 用三种不同的方法求有理分式 12 3 6 9 3 5 − + − x x x 在 x = 2 时的值,并体会这三种方法的区别。 (三种方法为 (i)赋值语句, (ii)代入规则, (iii)自定义函数) 3. 用两个不同的命令求解方程 3 2 13 3 2 4 3 3 2 x − x − x + x − =0, 并对结果进行比较。 4. 求 ( ) 3 1 5 f x = x − x + 的所有零点。 5. 已知数列的一般项为 n n an + + = 2 1 2 (1) 试求出由该数列的前 10 项构成的数表(表中的元素保留 3 位有效数字) (2) 取出该表中的第 3 个元素和第 7 个元素 (3) 分别把 0,1,-1 放在该表的第 1 位, 最后一位及第 8 位 (4) 求出最后一个表的长度 6. 分别画出函数 f (x) = sin 3x + 4cos x , − = 1 0 , 0, ( ) 2 e x x x g x x , + − = 1 , 1 1, 0 1, 1 , 0, ( ) x x x x x h x
的图形并用两种不同的方法将三个图形组合在一个图中(给每条曲线设计出不同的颜色)。 7.画出函数1+7在2上的图形 8.画图求曲线y=x,y-e-x=0的所有交点。 ∫(x) 9.用描点法画出函数 1+x的图形,并与用Pot语句画出的图形进行比较 (提示:(1)在某个范围内例如3,3]选取x的值,(2)生成以(x(x)为元素的二层表() 画出该表的折线图(4)将其与f(x)的图形组合进行比较) 10.用For语句打印 (1)1-10的平方 (2)自己的姓名(用拼音)10次 11.用三种不同的循环语句计算下列各题 11、+…+100(精确到0) 2 1+ 的前10;20;100项的和 实验二与极限、连续有关的问题 【实验目的】 1.能通过极限的动画演示理解极限的概念,特别是两个重要极限 2.熟练运用 Mathematica软件命令求极限(左、右极限),与理论知识相结合能够运用 极限存在的充要条件判断极限的存在性 3.能够用作图命令和极限命令结合理论知识求函数的间断点及判断间断点的类型 4.了解二分法的基本思想,会利用 Mathematica软件用二分法求方程的近似根; 5.熟练利用 Mathematica软件解决与极限、连续有关的一些简单实际问题 【实验的预备内容】 复习的数学内容 (1)极限的有关概念、数列极限与函数极限的关系、极限存在的充要条件; (2)一元函数连续的有关概念、间断点及其类型:一元函数的零点概念及零点定理 复习 Mathematica语句 Table、 ListPlot、/N与N[%]、f[x]:=、ShoW、Plot、For、Do、 While、If、 Which Print等语句的格式和用法 【 Mathematica新命令】 命令格式 功能说明 Limit [f[x],x->xo, Direction ->-1 lim f(x 求右极限x→x Limit f[], x->xo, Direction->+l Im f(x) 求左极限x→而 Point[ix, yI] Line[(,y).(r2,22)J1 画连接点(x,y)与(x2y)的直线段。 Pointsize[r],其中0≤r≤1 表示点的大小
的图形,并用两种不同的方法将三个图形组合在一个图中(给每条曲线设计出不同的颜色)。 7. 画出函数 + = + = 3 2 3 1 3 1 3 t t y t t x 在[-2,2]上的图形。 8. 画图求曲线 , 0 2 y = x y − e − x = x 的所有交点。 9. 用描点法画出函数 2 1 ( ) x x f x + = 的图形,并与用 Plot 语句画出的图形进行比较。 (提示: (1) 在某个范围内例如[-3，3]选取 x 的值, (2) 生成以 (x, f (x)) 为元素的二层表 (3) 画出该表的折线图 (4) 将其与 f (x) 的图形组合,进行比较) 10. 用 For 语句打印 (1) 1—10 的平方 (2) 自己的姓名(用拼音)10 次 11. 用三种不同的循环语句计算下列各题 (1) 2 2 2 100 1 3 1 2 1 1+ + ++ (精确到 5 10− ) (2) n 1 3 1 2 1 1+ + ++ 的前 10；20；100 项的和 实验二 与极限、连续有关的问题 【实验目的】 1．能通过极限的动画演示理解极限的概念，特别是两个重要极限； 2．熟练运用 Mathematica 软件命令求极限（左、右极限），与理论知识相结合能够运用 极限存在的充要条件判断极限的存在性； 3．能够用作图命令和极限命令结合理论知识求函数的间断点及判断间断点的类型； 4．了解二分法的基本思想，会利用 Mathematica 软件用二分法求方程的近似根； 5．熟练利用 Mathematica 软件解决与极限、连续有关的一些简单实际问题。 【实验的预备内容】 复习的数学内容 （1）极限的有关概念、数列极限与函数极限的关系、极限存在的充要条件； （2）一元函数连续的有关概念、间断点及其类型；一元函数的零点概念及零点定理。 复习 Mathematica 语句 Table、ListPlot、//N 与 N[%]、f[x]:=、Show、Plot、For、Do、While、If、Which、 Print 等语句的格式和用法。 【Mathematica 新命令】 命令格式 功能说明 [ [ ], , 1] Limit f x x− x0 Direction− − 求右极限 lim ( ) 0 f x x x → + [ [ ], , 1] Limit f x x− x0 Direction− + 求左极限 lim ( ) 0 f x x x → − Point[{x,y}] 画一点。 Line[{ ( ) 1 1 x , y , ( ) 2 2 x , y }] 画连接点 ( ) 1 1 x , y 与 ( ) 2 2 x , y 的直线段。 PointSize[r]，其中 0 r 1 表示点的大小
GRBColor[h27],其中05≤1i=123表示一种颜色。 Mathematica新命令】 通过动画演示理解极限概念 A,则。如何理解这句话呢?首先我们以数列248…、1一个确定的常数 数列极限概念的通俗说法是:若当n充分大时,如m能充分接近与唯 l11 为例来做实验。 例1观察数列当n→>∞时的变化趋势 解用 Table命令先选取数列的前10项进行观察 aa= Table[1/n^2,{n,1,10}]//N 结果:{1.,0.25,0.111110.0625,0.04,0.027778,0.0204082,0.015625, 0.0123457,0.01} 用 Listplot画出其散点图,执行后观察数列的变化趋势。 ListPlot/aa, PlotRange->110, 118/-0.2, 18/ 如图 10 从上图可以看出,当n增大时,数列(2”逐渐趋近于0 再选取数列的前20项进行观察 bb=Table(1/n 2,n, 1, 201//N 结果:{1.,0.25,0.111100625,0.04,0.0277780.02040820015625,0.0123457, 001,0.00826446,0.00694444,0.00591716,0.00510204,0.0044440.00390625, 000346021,0.00308642,0.00277008,0.0025} 再用 ListPlot画出其散点图,执行后观察数列的变化趋势 ListPlot/bb, PlotRange-> 0, 201, -0.2, B1 如图: 0.8 2.557.51012.51517.520
GRBColor[ 1 2 3 r ,r ,r ]，其中 0 ri 1,i =1,2,3 表示一种颜色。 【Mathematica 新命令】 一、通过动画演示理解极限概念 数列极限概念的通俗说法是：若当 n 充分大时， n a 能充分接近与唯一一个确定的常数 A ，则 an A n = → lim 。如何理解这句话呢？首先我们以数列 n 2 1 , 8 1 , 4 1 , 2 1 为例来做实验。 例 1 观察数列 2 1 n an = 当 n → 时的变化趋势。 解 用 Table 命令先选取数列的前 10 项进行观察。 用 ListPlot 画出其散点图，执行后观察数列的变化趋势。 从上图可以看出，当 n 增大时，数列 n 2 1 逐渐趋近于 0。 再选取数列的前 20 项进行观察。 再用 ListPlot 画出其散点图，执行后观察数列的变化趋势。 aa=Table[1/n^2,{n,1,10}]//N 结果：{1., 0.25, 0.111111, 0.0625, 0.04, 0.0277778, 0.0204082, 0.015625, 0.0123457, 0.01} ListPlot[aa,PlotRange->{{0,11},{-0.2,1}}] 如图： bb=Table[1/n^2,{n,1,20}]//N 结果：{1., 0.25, 0.111111, 0.0625, 0.04, 0.0277778, 0.0204082, 0.015625, 0.0123457, 0.01, 0.00826446, 0.00694444, 0.00591716, 0.00510204, 0.00444444, 0.00390625, 0.00346021, 0.00308642, 0.00277008, 0.0025} ListPlot[bb,PlotRange->{{0,20},{-0.2,1}}] 如图：
从上图再次可以看出,当乃增大时,数列(2确实逐渐趋近于 再逐次作数列前30、40、50等项的散点图,注意纵横上刻度的变化,可以更加确信数 列趋于0。因此可得n20 可通过下面的动画程序理解数列极限的概念 Forli=100, i>=10, i=i-10, aa=Table[ 1/n 2, n, 1, i//N ListPlot]aa, PlotRange->((0, 1001, (-0.2, 13311 运行程序后,得到一系列图形,选中所有图形(用鼠标单击最右边括号,使之变色),再单 击胶卷图标,即可动画演示(想终止动画,单击空白处) f∫(x)= 例2分析函数 当x→>∞时的变化趋势。 解先取一个较小的范围,如区间[-10,10],作出函数在这一区间上的图形 Plo1x^2=Sin{x,x,-10.,10 如图 手住区间 函数的图 Plot(1/2*Sin/x/,/x,40, 408/ 结果 f(x) sin x 从上面两个图形可以看出当增大时,函数 逐渐趋于0。 可通过下面的动画程序理解函数极限的概念 fx: =1/x 2*Sin x|; ForlF10, if(-30, 303,1-0.2,0.233II 运行程序后,得到一系列图形,选中所有图形(用鼠标单击最右边括号,使之变色),再单
从上图再次可以看出，当 P3 增大时，数列 n 2 1 确实逐渐趋近于 0。 再逐次作数列前 30、40、50 等项的散点图，注意纵横上刻度的变化，可以更加确信数 列趋于 0。因此可得 0 2 1 lim = → n n 。 可通过下面的动画程序理解数列极限的概念。 运行程序后，得到一系列图形，选中所有图形（用鼠标单击最右边括号，使之变色），再单 击胶卷图标，即可动画演示（想终止动画，单击空白处）。 例 2 分析函数 x x f x sin 1 ( ) 2 = 当 x → 时的变化趋势。 解 先取一个较小的范围，如区间[-10，10]，作出函数在这一区间上的图形。 再在区间[-20，20]内作出这一函数的图形 从上面两个图形可以看出当 x 增大时，函数 x x f x sin 1 ( ) 2 = 逐渐趋于 0。 可通过下面的动画程序理解函数极限的概念。 运行程序后，得到一系列图形，选中所有图形（用鼠标单击最右边括号，使之变色），再单 For[i=100,i>=10,i=i-10,aa=Table[1/n^2,{n,1,i}]//N; ListPlot[aa,PlotRange->{{0,100},{-0.2,1}}]] Plot[1/x^2*Sin[x],{x,-10,10}] 如图： Plot[1/x^2*Sin[x],{x,-40,40}] 结果： f[x_]:=1/x^2*Sin[x];For[i=10,i{{-30,30},{-0.2,0.2}}]]
f(x=-sin x 击胶卷图标,即可动画演示(想终止动画,单击空白处)。可得函数 xX→ 时的极限为0。 练习1借助于上面两个动画演示程序考察下列函数(或数列)在给定极限过程中的变 化趋势 sin x f(x)=(1+-) f(x)=xsin an=(1+-) (3) x,x→>0(4) n→0 练习2已知数列x=2,xm=2+x ,画图观察数列的极限是否存在? 二、极限的计算 例3判极限x→0的存在性,若存在,求出其极限。注意左右极限的用法。 解先注意下列结果是否正确 Limit[Exp-1/x],x->0结果:0 下面区分极限方向再判断 Limit Expl-1/xl,x->0, Direction->-11结果:0 Limit(Expl-1/x1, x->0, Direction->+1| t: Infinity lim e I 从上面的结果可以判断x∽0不存在。很显然第一次计算的结果是错误的 例4计算下列极限 xsin x+cosx-1 m (1)Limi(3a,n> Infinit, Direction-x1结果:E Limit(+)~(3m),n-> Infinity, Direction-)+结果:E 由此可得 (2)Limit((x*Sin(x)+ Cos x] -1)/x 2, x->0, Direction->-1| 4R: 2 lim -sin x+ cosx-l 1 由此可得
击胶卷图标，即可动画演示（想终止动画，单击空白处）。可得函数 x x f x sin 1 ( ) 2 = 当 x → 时的极限为 0。 练习 1 借助于上面两个动画演示程序考察下列函数（或数列）在给定极限过程中的变 化趋势。 （1） x x f x sin ( ) = ， x →0 （2） x x f x ) 1 ( ) = (1+ ， x → （3） x f x x 1 ( ) = sin ， x →0 （4） n n n a ) 1 = (1+ ，n → 练习 2 已知数列 x1 = 2 ， n n x = + x +1 2 ，画图观察数列的极限是否存在？ 二、极限的计算 例 3 判极限 x x e 1 0 lim − → 的存在性，若存在，求出其极限。注意左右极限的用法。 解 先注意下列结果是否正确？ 下面区分极限方向再判断： 从上面的结果可以判断 x x e 1 0 lim − → 不存在。很显然第一次计算的结果是错误的。 例 4 计算下列极限 （1） n n n 3 1 lim 1 + → （2） 2 0 sin cos 1 lim x x x x x + − → （3） x x 1 lim sin →0 由此可得 3 3 1 lim 1 e n n n = + → 。 由此可得 2 sin cos 1 1 lim 2 0 = + − → x x x x x 。 Limit[Exp[-1/x],x->0] 结果：0 Limit[Exp[-1/x],x->0,Direction->-1] 结果：0 Limit[Exp[-1/x],x->0,Direction->+1] 结果：Infinity （1）Limit[(1+1/n)^(3n),n->Infinity,Direction->-1] 结果： 3 E Limit[(1+1/n)^(3n),n->Infinity,Direction->+1] 结果： 3 E （2）Limit[(x*Sin[x]+Cos[x]-1)/x^2,x->0,Direction->-1] 结果： 2 1 Limit[(x*Sin[x]+Cos[x]-1)/x^2,x->0,Direction->+1] 结果： 2 1
(3)Limit(Sin(1/x], x->0, Direction->-1| tA: Interval[( -1, 1)1 Limit Sin 1/xl,x->0. Direction+11结果: Interval-1,1(在(-11)之间振荡) Plot(Sin1/x,x,0.5,0.5}如下图 直接计算不能凑效 (1) Limit(-1)^(2n,n-> Infinity, Direction->-1结果: Indeterminate(不定式) Limit(1)^(2n),n-> Infinity, Direction->+1结果: Indeterminate(不定式) (2) Limit[Logix/x,x->0, Direction->-l不能直接计算出结果。 Limit[ Log|x]/x,x->0, Direction->+1]不能直接计算出结果 改变函数形式后在计算,问题就可得到解决 (1)Limit(((-1)2)n, n->Infinity, Direction->-1| A:I Limit(((-1)/2)n, n->Infinity, Direction->+1| #R:1 im(-1)=1 由此可得n→ (2)Limit[ Log[x (1/x)l, x->0, Direction->-1 R:-Infinity Limit LogIx"(1/x)1, x->0, Direction->+1 R: Infinity lim nr 由此可得x0x (不存在)。 例6求下列极限的近似值 Iim +3x3-5 (2)n→n”(精度为10 im1+-+-+…+ (1)法一用 Limit命令先求左、右极限 Limit((x*3-1)/(x4+3x* 2+2x-1)x->1, Direction->-1 LR: 5
例 5 计算下列极限。 （1） n n 2 lim (−1) → （2） x x x ln lim →0 直接计算不能凑效。 改变函数形式后在计算，问题就可得到解决。 由此可得 lim ( 1) 1 2 − = → n n 。 由此可得 = → x x x ln lim 0 （不存在）。 例 6 求下列极限的近似值。 （1） 3 5 2 1 1 lim 4 3 2 3 1 + − + − − → x x x x x x （2） n n n n! lim → （精度为 6 10 − ） （3） + + + + n→ n 1 3 1 2 1 lim 1 （1）法一 用 Limit 命令先求左、右极限。 （1）Limit[(-1)^(2n),n->Infinity,Direction->-1] 结果：Indeterminate（不定式） Limit[(-1)^(2n),n->Infinity,Direction->+1] 结果：Indeterminate（不定式） （2）Limit[Log[x]/x,x->0,Direction->-1] 不能直接计算出结果。 Limit[Log[x]/x,x->0,Direction->+1] 不能直接计算出结果。 （1）Limit[((-1)^2)^n,n->Infinity,Direction->-1] 结果：1 Limit[((-1)^2)^n,n->Infinity,Direction->+1] 结果：1 （2）Limit[Log[x^(1/x)],x->0,Direction->-1] 结果：-Infinity Limit[Log[x^(1/x)],x->0,Direction->+1] 结果：Infinity Limit[(x^3-1)/(x^4+3x^3-5x^2+2x-1),x->1,Direction->-1] 结果： 5 3 Limit[(x^3-1)/(x^4+3x^3-5x^2+2x-1),x->1,Direction->+1] 结果： 5 3 （3）Limit[Sin[1/x],x->0,Direction->-1] 结果：Interval[{-1, 1}] Limit[Sin[1/x],x->0,Direction->+1] 结果：Interval[{-1, 1}]（在 (−1,1) 之间振荡） Plot[Sin[1/x],{x,-0.5,0.5}] 如下图：
执行后结果均为5,从而可知左、右极限均存在且相等,故 m x小1x++3x3-5x2+2x-15 法二在x=1附近,用Plot命令画出函数的图形求极限 nx|:=x^3-1)(x^4+3x^3-5x^2+2x-1); Plot fx{x,0,2 从图形看不出趋向情况,缩小画图区间: 印x]=(x^3-1)/(x~4+3x^3-5x^2+2x-1)Plot[x],{x0.5,1.5} 从图知纵轴上所代表的函数值从0.54变化到0.66之间,再缩小区间画图 印x]=(x^3-1)/(x~4+3x^-3-5x2+2x-1) Plot[fx]{x0.9991001 0,599 限为0.6 (2)法一用 Limit命令先求左、右极限。 Limit]n! /n n, n->Infinity, Direction->-1I 不能直接计算出结果。 Limit[n! /nn, n->Infinity, Direction-+ll 不能直接计算出结果。 x|n|:=n!/n^n; Forli=l, i<=20, i++, IfN(x(i+1l-xillP10(6), mi+1; Return 1; PrintI"m=, m Printl"xIm]=", x|m///N 结果:m=16xm]=113423×106
执行后结果均为 5 3 ，从而可知左、右极限均存在且相等，故 5 3 3 5 2 1 1 lim 4 3 2 3 1 = + − + − − → x x x x x x 。 法二 在 x=1 附近，用 Plot 命令画出函数的图形求极限。 从图形看不出趋向情况，缩小画图区间： 从图知纵轴上所代表的函数值从 0.54 变化到 0.66 之间，再缩小区间画图： 函数值从 0.59975 变化到 0.60025，如此可再作下去，知函数值逐渐缩小至 0.6，故所求极 限为 0.6 （2）法一 用 Limit 命令先求左、右极限。 法二 用循环语句求极限的近似值。精度要求为 6 10 − ，即往往要求 6 | 1 | 10− xn+ − xn . f[x_]:=(x^3-1)/(x^4+3x^3-5x^2+2x-1);Plot[f[x],{x,0,2}] 如图： f[x_]:=(x^3-1)/(x^4+3x^3-5x^2+2x-1);Plot[f[x],{x,0.5,1.5}] 如图： f[x_]:=(x^3-1)/(x^4+3x^3-5x^2+2x-1);Plot[f[x],{x,0.999,1.001}] 如图： Limit[n!/n^n,n->Infinity,Direction->-1] 不能直接计算出结果。 Limit[n!/n^n,n->Infinity,Direction->+1] 不能直接计算出结果。 x[n_]:=n!/n^n; For[i=1,i10^(-6),m=i+1;Return];]; Print["m=",m] Print["x[m]=",x[m]//N] 结果：m=16 x[m]=1.13423×10-6