第7章 MALAB符号计算.ppt

* * * * * * * * * 7.2.2 符号函数求导及其应用 diff函数用于对符号表达式求导数。 该函数的一般调用格式为： diff(f)：没有指定变量和导数阶数，则系统按findsym函数指示的默认变量对符号表达式f求一阶导数。 diff(f,v)：以v为自变量，对符号表达式f求一阶导数。 diff(f,n)：按findsym函数指示的默认变量对符号表达式f求n阶导数，n为正整数。 diff(f,v,n)：以v为自变量，对符号表达式f求n阶导数。 例7.3 求下列函数的导数。(1) y= ，求y (2) y= ，求y ，y  syms a b t x y z; f1=sqrt(1+exp(x)); diff(f1) %求(1) f2=x*cos(x); diff(f2,x,2) %求(2)中y diff(f2,x,3) %求(2)中y 例7.3 求下列函数的导数。(3) ,求yx ,yx (4) z= ，求zx ，zy g1=a*cos(t); g2=b*sin(t); diff(g2)/ diff(g1) %求(3)中y , %下行求y (diff(g1)*diff(g2,2)- diff(g1,2)*diff(g2))/ (diff(g1))^3 f4=x*exp(y)/y^2; diff(f4,x) %求(4) diff(f4,y) %求(4) 例7.3 求下列函数的导数。(5) z=f(x,y)由方程x2+y2+z2=a2定义 求zx ，zy f=x^2+y^2+z^2-a^2; zx=-diff(f,x)/diff(f,z) %求(5)。按隐函数求导公式求z对x的偏导数 zy=-diff(f,y)/diff(f,z) %求(5)。按隐函数求导公式求z对y的偏导数 例7.4 在曲线y=x3+3x-2上哪一点的切线与直线y=4x-1平行。 即求曲线哪一点的导数值为4 命令如下： x=sym(x); y=x^3+3*x-2; %定义曲线函数 f=diff(y); %对曲线求导数 g=f-4; solve(g) %求方程f-4=0的根，即求曲线何处的导数为4 实验 3 用符号方法求下列极限或导数。 7.3 符号积分 7.3.1 不定积分 符号积分由函数int来实现。调用格式为： int(f)：没有指定积分变量和积分阶数时，系统按findsym函数指示的默认变量对被积函数或符号表达式f求不定积分。 int(f,v)：以v为自变量，对被积函数或符号表达式f求不定积分。 例7.5 求不定积分。 命令如下： syms x t; f=(3-x^2)^3; int(f) %求不定积分 f=5*x*t/(1+x^2); int(f, t) 7.3.2 符号函数的定积分 定积分在实际工作中有广泛的应用。在MATLAB中，定积分的计算使用函数： int(f,x,a,b)：求定积分运算。 a,b分别表示定积分的下限和上限。该函数求被积函数在区间[a,b]上的定积分。a和b可以是两个具体的数，也可以是一个符号表达式，还可以是无穷(inf) 当函数f关于变量x在闭区间[a,b]上可积时，函数返回一个定积分结果。 当a,b中有一个是inf时，函数返回一个广义积分。 当a,b中有一个符号表达式时，函数返回一个符号函数。 例7.6 求定积分。 命令如下： x=sym(x); t=sym(t); int(abs(1-x),1,2) %求定积分(1) f=1/(1+x^2); int(f, -inf, inf) %求定积分(2) int(4*x/t, t, 2, sin(x)) %求定积分(3) f=x^3/(x-1)^10; I=int(f,2,3) %用符号积分的方法求定积分(4) double(I) %将上述符号结果转换为数值 例7.8 求空间曲线c从点(0，0，0)到点(3，3，2)的长度。求曲线c的长度是曲线一型 命令如下： syms t; x=3*t;y=3*t^2;z=2*t^3; f=diff([x,y,z],t) %求x,y,z对参数t的导数 g=sqrt(f*f) %计算一型积分公式中的根式部分 l=int(g,t,