2.2 一元函微积分学.docVIP

§2.2 一元函数微积分学 2.2.1 导数与微分 实验使用的函数与命令 1、极限函数p=limit(F,x,a)、绘图函数plot； 2、向量或矩阵求导指令：diff 1) Y = diff(X)：计算向量或矩阵的差分，若X是向量，则diff(X) 返回[X(2)-X(1) X(3)-X(2) ... X(n)-X(n-1)]；若X是矩阵，则diff(X) 返回X的各列差分矩阵[X(2:m,:)-X(1:m-1,:)]。 如：x = [1 2 3 4 5]; y = diff(x)其输出结果为y = 1 1 1 1 2) Y = diff(X,n)：计算n阶向量或矩阵的微分(或差分)； 如：x = [1 2 3 4 5]; z = diff(x,2) 其输出结果为z = 0 0 0 3、符号求导指令diff 1) diff(S,v) 或 diff(S,sym(v))：返回符号表达式S对自变量v的导数； 2) diff(S,n)：对于正整数n，求S的n阶导数； 3) diff(S,v,n) 与 diff(S,n,v)：返回符号表达式S对自变量v的n阶导数，缺省n时为求1阶导数、缺省v时默认变量x。 4、diff(y)./diff(x)是求近似导数。 实验指导 例1 求正弦函数的一阶导数与函数的二阶导数。 解：M文件程序处理如下： syms x; f=sin(5*x); g=exp(x)*cos(x); df=diff(f) dg=diff(g,2) 输出结果： df = 5*cos(5*x) dg = -2*exp(x)*sin(x) 例2 设函数，讨论在处的左、右导数与导数。 解：由若函数在点导数的定义和导数存在的充要条件，考察在在点的左右导数。M文件程序设计如下： syms x; diffleft=limit(abs(x)/x,x,0,left) diffright=limit(abs(x)/x,x,0,right) 输出结果： diffleft = -1 diffright = 1 由以上可知，函数在点左、右导数都存在，但，故函数在点导数不存在。 例3 已知抛射线的运动轨迹的参数方程为 求抛射体在时刻t的运动速度的大小和方向。 解：该问题的解决分为两步，先求在时刻t的运动速度的大小，再求轨道的切线方向，即(为切线的倾角)，由此M文件程序设计如下： syms t v1 v2 g; x=v1*t; y=v2*t-1/2*g*t^2; dx=diff(x,t); dy=diff(y,t); v=sqrt(dx^2+dy^2) %运动速度的大小 tan=dy/dx %轨道的切线方向 输出结果为： v = (v1^2+v2^2-2*v2*g*t+g^2*t^2)^(1/2) tan = (v2-g*t)/v1 对于该题，给定初始速度，即可画出抛射体在时刻t的运动轨迹。 命令diff(y)./diff(x)是求近似导数的一种方法。又由导数的几何意义：为函数在处的切线斜率，即，其中为切线的倾角。 例4 求正弦函数在处的导数，并利用导数几何意义，绘出正弦函数在处的切线。 解：对于正弦函数在处的导数，可以利用上述方法或定义形式来求得。M文件程序设计如下： x=0:0.001:0.001;%给出x在0点的向量 y=sin(x); diff(x);%计算向量x的差分 diff(y);%计算向量y的差分 dydx=diff(y)./diff(x)%求在x=0处近似导数 syms x; difflimit=limit(sin(x)/x,x,0)%利用定义求x=0处导数 x1=0:0.01:pi; y1=sin(x1); plot(x1,y1); x2=0:0.01:1; y2=dydx*x2+0; hold on plot(x2,y2,--r);%画出图像图2.1 输出数值结果为 dydx = 1.0000 difflimit = 1 注：由输出结果可以看到，近似计算dydx = difflimit ，因此diff(y)./diff(x)求近似导数有一定的实际意义。下面再看一例： 例 5 观察下列程序的结果，体会diff(y)./diff(x)用法。 x=0:0.1:5; y=cos(x); dydx=diff(y)./diff(x) hold on plot(x,y,--r); plot(x(1:50),dydx); 输出结果： 思考：1）在程序中最后语句plot(x(1:50),dydx);为什么要取x(1:50)而不是向量x ? 2）plot(x(1:50),dydx)所绘图形（图中实线）表示的意义是什么？ 例6 讨论函数在处的连续性和可导性