第5章MATLAB符号运算PPTPPT.ppt

MATLAB 7.X程序设计 sym(1.5) sym(3.14) phi=(1+sqrt(5))/2 eval(phi) 5.4 微积分 1.符号极限 格式： limit(f,x,a)： limit(f) : limit(f,a) : limit(f,x,a,’left’) : limit(f,x,a,’right’) : sym x limit(x^x,x,0,’right’) ans=1 limit(sin(x)/x,x,0) ans= 1 limit(1/x,x,0,’ left’) ans=-Inf %负无穷大 例 计算 ， ，  函数极值 格式： fmin(‘fname’,x1,x2) fmin(‘fname’,x0) 求 在[0,5]内的最小值 2.微分 格式： diff(f):求f对默认自变量的一次微分值 diff(f,t):求f对自变量t的一次微分值 diff(f, n):求f对默认自变量的n次微分值 diff(f,t,n):求f对自变量t的n次微分值 例 求f(x)=ax2+bx+c的微分 f=sym（’a*x^2+b*x+c’）; diff(f) %对x求一次微分 diff(f,2) %对x求二次微分 diff(f,a) %对a求一次微分 diff(f,a,2) %对a求二次微分 diff(diff(f),a) %对x和a偏导 3.积分 符号积分格式： int(f):求f对默认自变量的不定积分 int(f,t):求f对t的不定积分 int(f,a,b):求f对默认自变量的定积分,积分区间[a,b] int(f,t,a,b):求f对t的定积分,积分区间[a,b] 例 求f(x)=ax2+bx+c的积分 syms a b c x f=sym(a*x^2+b*x+c) ; int(f) %f的不定积分，自变量是x int(f,0,2) %f在[0,2]的定积分，自变量是x int(f,a) %f的不定积分，自变量是a int(int(f,a),x) 数值积分 [I,n]=quad(‘fname’,a,b,tol,trace) [I,n]=quad8(‘fname’,a,b,tol,trace) dblquad(fun,a,b,c,d,tol) triplequad(fun,a,b,c,d,e,f,tol) 例 求定积分。 g=inline(exp(-0.5*x).*sin(x+pi/6)); [S,n]=quad(g,0,3*pi) [S,n]=quad8(g,0,3*pi) 例 求二重积分。 f=inline(exp(-x.^2/2).*sin(x.^2+y),x,y);I=dblquad(f,-2,2,-1,1) 4 级数 1） 级数符号求和 symsum(s,v,n,m) 其中s表示一个级数的通项，是一个符号表达式。v是求和变量，v省略时使用系统的默认变量。n和m是求和的开始项和末项。 例 求下列级数之和。 syms n; s=symsum(1/n^2,1,inf) s = 1/6*pi^2 eval(s) ans = 1.6449 2) 函数的泰勒级数 taylor(f,v,n,a) 该函数将函数f按变量v展开为泰勒级数，展开到第n项(即变量v的n-1次幂)为止，n的缺省值为6。v的缺省值与diff函数相同。参数a指定将函数f在自变量v=a处展开，a的缺省值是0。 例 将ln(x)在x=1处按5阶多项式展开。 syms x taylor(log(x),x,6,1) ans = x-1-1/2*(x-1)^2+1/3*(x-1)^3-1/4*(x-1)^4+1/5*(x-1)^5 5.5.1 多项式的四则运算 1．多项式的加减运算 2．多项式乘法运算 conv(P1,P2) P1、P2是两个多项式系数 向量。