实验二-极限和导数.docVIP

实验二、极限与导数 【实验目的】 1．了解函数极限、导数的基本概念。 2．学习掌握MATLAB软件有关求极限、导数的命令。 【实验内容】 判断极限的存在性 验证极限 验证极限 求数的单调区间及极值 【实验准备】 极限、导数的基本概念 数列极限：如果对于，存在正整数，使得当时有。则称为数列的极限，或称收敛于，记为。直观上表示：趋于无穷大时，无限接近。 函数极限：如果当时有，则称为函数当时的极限。记为。若仅当且（或）时有，则称为函数当时的右极限（左极限），记为（或）。当且仅当时，当时的极限存在且等于这个值。 导数：函数在点的导数的定义为： 它反映了在点附近函数的变化率。当时，函数在点附近是上升的，反之时，函数在点附近是下降的，而当时，往往（但不一定）标志函数在点达到局部极大或局部极小。函数在点达到局部极大（或局部极小）的充分条件是且（或）。从几何意义上说，是函数在点切线的斜率。 2．求极限、导数的MATLAB命令 MATLAB中主要用limit,diff分别求函数的极限与导数。 limit(s,n,inf) 返回符号表达式当n趋于无穷大时表达式s的极限 limit(s,x,a) 返回符号表达式当x趋于a时表达式s的极限 limit(s,x,a,’left’) 返回符号表达式当x趋于a-0时表达式s的左极限 limit(s,x,a,’right’) 返回符号表达式当x趋于a-0时表达式s的右极限 diff(s,x,n) 返回符号表达式s对自变量x的n阶导数 可以用help limit, help diff查阅有关这些命令的详细信息 【实验方法与步骤】 练习1 判断极限的存在性 （1）作图观察 　x=-1:0.01:-0.01; y=cos(1./x); plot(x,y) （2）MATLAB程序直接求 clear syms x; limit(cos(1/x),x,0) ans = limit(cos(1/x), x = 0) 极限不存在。 同理可判断极限也不存在。 练习2判断极限的存在性 （1）作图观察 x=-1:0.0001:(-0.01); y=sin(x)./x; plot(x,y) （2）MATLAB程序直接求 clear syms x; limit(sin(x)/x,x,0) ans = 1 练习3 考虑时，数列的变化趋势，并求它们的极限。 （1）作图观察 clear x=1:0.1:10; y1=exp(x.*log(1+1./x)); y2=exp((x+1).*log(1+1./x)); y3=2.71828; plot(x,y1,-.,x,y2,:,x,y3) 或 x=1:0.1:10; y1=exp(x.*log(1+1./x)); y2=exp((x+1).*log(1+1./x)); y3=2.71828; plot(x,y1,r,x,y2,b,x,y3,k) （2）MATLAB程序直接求 求 clear syms n; limit((1+1/n)^n,n,inf ) %求 ans = exp(1) 求 clear syms n; limit((1+1/n)^(n+1),n,inf) %求 ans = exp(1) 练习4（极限的定义和判别） 用MATLAB语言来表达推理过程是比较困难的，它必须与实际的数值联系起来，比如无法用无穷小和高阶无穷小的概念，只能用等数职。极限的定义恰恰是用了和等数值的概念，因此不难用程序表述。 用函数极限的定义，对于函数，当任意给定一个正数时，有一个对应的正数存在，使得当 时，有， 则就是在时的极限，如果找不到这样的，就不是它的极限，只考虑左极限时，因必为正数，可去掉绝对值符号。检验左极限是否正确的程序为： （1）输入M文件 clear disp(A是否是f(xc)的左极限) eps=1.0e-10; A=input(A=); xc=input(xc=);%输入对应的自量变值; fxc=input(f(x)=); flag=1; delta=1; x=xc-delta; n=1; while flag==1 epsilon=input(任给一个小的数=) while abs(A-eval(fxc))epsilon delta=delta/2; x=xc-delta; if abs(delta)eps disp(找不到delta) n=0; break end end if n==0 disp(左极限不存