第九章matlab符号计算-席2015详解.pptVIP

级数求和 求无穷级数的和，利用符号表达式求和函数symsum，其调用格式为： symsum(s,v,a,b) 其中s表示一个级数的通项，是一个符号表达式。v是求和变量，v省略时使用系统的默认变量。a和b是求和的开始项和末项。 例 求 clear syms n x s1=symsum(1/n^2,n,1,inf) s1 = 1/6*pi^2 s2=symsum((-1)^(n-1)*x^n/n,n,1,inf) s2 = log(1+x) 泰勒级数 MATLAB提供了taylor函数将函数展开为幂级数，其调用格式为： taylor(f):计算符号表达式f在默认自变量等于0处的5阶Taylor级数展开式。 taylor(f,n,v):计算符号表达式f在自变量v=0处的n-1阶Taylor级数展开式。 taylor(f,n,v,a):计算符号表达式f在自变量v=a处的n-1阶Taylor级数展开式 例 求函数在指定点的泰勒展开式。 x=sym(x); %syms x; f1=(1+x+x^2)/(1-x+x^2); f2=sqrt(1-2*x+x^3)-(1-3*x+x^2)^(1/3); taylor(f1,5,x,3)%求f1在x=3处的4阶Taylor展式 ans = 139/49-16/49*x+38/343*(x-3)^2-78/2401*(x-3)^3+124/16807*(x-3)^4 taylor(f2,3,x) %求f2在x=0处的2阶Taylor展式 ans = 1/6*x^2 根据前面讲的内容，下面做个练习： 定义符号函数F=e2xcos(3x+a)，并对其求x=pi/2时的函数值、当x趋于0时的极限、 一阶导数、不定积分。 syms x a F=exp(2*x)*cos(3*x+a) 或 sym(exp(2*x)*cos(3*x+a)) subs(F,x,pi/2) limit(F,x,0) diff(F,x) int(F,x) 1. 图示化符号函数计算器(funtool) 对于习惯使用计算器或者只是想做一些简单的符号运算与图形处理的，Matlab7.0提供的图示化符号函数计算器是一个较好的选择。功能简单，操作方便，可视性强。 Matlab的可视化数学分析界面 2. 泰勒级数逼近分析器（taylortool） A1=factor(sym(‘243’)) 上次课到此 直接在命令窗口键入funtool即可打开此界面 f=cos(x^3)/(1+x^2) g=1/(5+4*cos(x)) 直接在命令窗口键入taylortool即可打开此界面 MATLAB提供了符号表达式的因式分解与展开的函数，函数的调用格式为： factor(s)：对符号表达式s分解因式。 说明：s可以是正整数、符号整数、符号表达式或其矩阵。当s是正整数时，则分解为质数分解式；当s为符号表达式，则分解为乘积形式；当s为符号整数或符号表达式矩阵，则分解每一个元素为质数分解式或乘积形式。 expand(s)：对符号表达式s进行展开。 collect(s)：对符号表达式s合并同类项。 collect(s,v)：对符号表达式s按变量v合并同类项。 符号表达式的因式分解与展开 A1=factor(sym(243)) A1 = (3)^5 syms x y; A3=factor(x^5+y^5) A3 = (x+y)*(x^4-x^3*y+x^2*y^2-x*y^3+y^4) syms x y a b; A4=factor([a*x-a*b,x^2-b*x;(x-b)^2,x^2-b^2]) A4 = [ a*(x-b), x*(x-b)] [ (x-b)^2, (x-b)*(x+b)] A2=factor(243) A2 =3 3 3 3 3 syms x; f=(x-3)^2*(x-1)-(x+1)^2; expand(f) ans = x^3-8*x^2+13*x-10 ff=(x-y)^3+(x-y)^2+y^4; f1=collect(ff,y) f1 = y^4-y^3+(3*x+1)*y^2+(-3*x^2-2*x)*y+x^3+x^2 MATLAB提供的对符号表达式化简的函数有： simplify(s):对表达式s进行化简； [r,how]=simple(S) :通过对表达式尝试多种不同的算法进行化简，并显示化简过程，寻求符号表达式S的最简形式；r