第五章 MATLAB符号运算3.pptVIP

第五章 MATLAB符号运算 靖辉 5.1 数值运算与符号运算 数值运算：运算前必先对变量赋值，才能运算 符号运算：对字符串进行数学分析，运算结果 以标准符号形式表述。 5.2 符号变量和符号表达式 格式： (1)符号变量 syms 符号变量1 符号变量2…符号变量n sym(‘符号变量’) 说明： syms可同时定义多个符号变量，各变量之间用空格分开 sym只定义一个符号变量 (2)符号表达式 f=sym(‘表达式’) 例1创建圆的半径符号函数r和面积符号函数area syms x y r=sqrt(x^2+y^2); area=pi*r^2; r r = (x^2+y^2)^(1/2) area area = pi*(x^2+y^2) 5.3 微积分 1.极限 格式： limit(f) : limit(f,a) : limit(f,x,a,’left’) : limit(f,x,a,’right’) : 例3 计算 、 、 syms x limit(1/x) ans=NaN %不定量 limit(1/x,x,0,’left’) ans= - Inf %负无穷大 limit(1/x,x,0,’ right’) ans=Inf %无穷大 2.微分 格式： diff(f):求f对默认自变量的一次微分值 diff(f,t):求f对自变量t的一次微分值 diff(f, n):求f对默认自变量的n次微分值 diff(f,t,n):求f对自变量t的n次微分值 例4 求f(x)=ax2+bx+c的微分 f=sym（’a*x^2+b*x+c’）; diff(f) %对x求一次微分 ans=2*a*x+b diff(f,2) %对x求二次微分 ans=2*a diff(f,a) %对a求一次微分 ans=x^2 diff(f,a,2) %对a求二次微分 ans=0 diff(diff(f),a) %对x和a偏导 ans=2*x 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] 例5 求f(x)=ax2+bx+c的积分 clear syms a b c x f=sym(a*x^2+b*x+c) ; int(f) %f的不定积分，自变量是x ans =1/3*a*x^3+1/2*b*x^2+c*x int(f,0,2) %f在[0,2]的定积分，自变量是x ans = 8/3*a+2*b+2*c int(f,a) %f的不定积分，自变量是a ans =1/2*a^2*x^2+b*x*a+c*a int(int(f,a),x) ans =1/6*a^2*x^3+1/2*b*x^2*a+c*a*x 5.4 符号表达式的运算 1.算术运算 与数值运算一样，注意”\”,”/”的区别 例6 表达式的算术运算 f1=sym(1/(a-b)) ; f2=sym(a/(a+b)); f1+f2 f1-f2 f1*f2 f1/f2 ans = 1/(a-b)/a*(a+b) f1\f2 ans = a/(a+b)*(a-b) f1^2 ans = 1/(a-b)^2 2.关系运算 运算符：== ~= 关系成立为1，不成立为0 f1==f2 ans=0 f1~=f2 ans=1 3.函数运算 (1)合并 collect 格式：collect(f) 对默认变量同幂项系数合并 collect(f,v) 对变量v同幂项系数合并 例7 合并 syms x y; R1=collect([(x+1)*(y+1),x+y]) R1 =[ (y+1)*x+y+1, x+y] R2=collect( (x+y)*(x^2+y^2+1),y) R2 =y^3+x*y^2+(x^2+1)*y+x*(x^2+1) (2)展开expand 格式： expand(f) 将表达式f展开 例