MATLAB 应用 第三章 导数的求解[精品].ppt

MATLAB 高等数学实验 实验三 导数 实验目的 深入理解导数与微分的概念，导数的几何意义。掌握用MATLAB求导数与高阶导数的方法。深入理解和掌握求隐函数的导数及由参数方程定义的函数的导数的方法。 * 3.1 学习MATLAB命令 * 求导数命令是diff，常用格式为： syms x diff(f(x),x) diff(f,x)给出f关于x的导数，而将表达式f中的其他字母看作常量。因此，如果表达式是多元函数，则给出的是偏导数。 diff(f(x),x,n)给出f关于x的n阶导数或者偏导数。 * 3.2 实验内容 * 3.2.1 导数概念与导数的几何意义 【例1】用定义求 的导数。 输入： syms x diff(x^3-3*x^2+x+1) 执行以后得到导函数： ans= 3*x^2-6*x+1 * 再输入： x=-1:0.1:3; y1=x.^3-3*x.^2+x+1; y2=3*x.^2-6*x+1; plot(x,y1,b,x,y2,r:) 执行后便得到函数y1=g(x)和它的导数y2=g(x)的图形(见图3.1，图中虚线是曲线 g(x))。 * 图3-1 * 【例2】作函数 的图形和在x=-1处的切线。 输入： syms x hanshu=2*x^3+3*x^2-12*x+7; daoshu=diff(2*x^3+3*x^2-12*x+7); x=-1; hanshuzhi=eval(hanshu) daoshuzhi=eval(daoshu) * 执行后得到函数f(x)在x=-1处的函数值和导数值： hanshuzhi=20 daoshuzhi=-12 再执行： x=-4:0.1:3; y=2*x.^3+3*x.^2-12*x+7; y1=20-12*(x+1); plot(x,y,b,x,y1,r) 便在同一个坐标系内作出了函数f(x)的图形和它在x=-1处的切线(见图3.2，其中直线为切线)。 * 图3-2 * 3.2.2 求函数的高阶导数及函数  在某点的导数值 【例3】求函数 的一阶导数和二阶导数。 输入： syms x diff(x^n,1) diff(x^n,2) 执行后得一阶导数和二阶导数分别为： ans= n*x^(n - 1) ans= n*x^(n - 2)*(n - 1) * 【例4】求函数 的一阶导数， 并求 。 输入： syms x a b daoshu=diff(sin(a*x)*cos(b*x)) x=1/(a+b); daoshuzhi=eval(daoshu) 执行后分别得函数f(x)的一阶导数及其值： daoshu=a*cos(a*x)*cos(b*x) – b*sin(a*x)*sin(b*x) daoshuzhi=a*cos(a/(a + b))*cos(b/(a + b)) – b*sin(a/(a + b))*sin(b/(a + b)) * 【例5】求函数 的1阶到11阶导数。 为了将1阶到11阶导数一次都求出来，输入： syms x y=x^10+2*(x-10)^9; for n=1:11 diff(y,x,n) end * 输出为： ans =18*(x - 10)^8 + 10*x^9 ans =144*(x - 10)^7 + 90*x^8 ans =1008*(x - 10)^6 + 720*x^7 ans =6048*(x - 10)^5 + 5040*x^6 ans =30240*(x - 10)^4 + 30240*x^5 ans =120960*(x - 10)^3 + 151200*x^4 ans =362880*(x - 10)^2 + 604800*x^3 ans =1814400*x^2 + 725760*x - 7257600 ans =3628800*x + 725760 ans =3628800 ans =0 * 3.2.3 求隐函数的导数，由参数 方程定义的函数的导数 【例6】求方程 确定的隐函数的导数。 输入： syms x y z=2*x^2-2*x*y+y^2+x+2*y+1; daoshu=-diff(z,x)/diff(z,y) 执行后得到： daoshu= -(4*x-2*y