第5章matlab符号运算.pptVIP

(1) 第一步：在命令窗口输入如下程序。 sym x y=(x^4+1)/(x^2+1) f=int(y) 第二步：按“enter”键，得到的结果如下。 f=2*atan(x) - x + x^3/3 (2) 第一步：在命令窗口输入如下程序。 sym x y=1/(1+x+x^2) f=int(y) 第二步：按“enter”键，得到的结果如下。 f=(2*3^(1/2)*atan((2*3^(1/2)*x)/3 + 3^(1/2)/3))/3 (3) 第一步：在命令窗口输入如下程序。 syms a x y=a*x^3 f=int(y,0,2) 第二步：按“enter”键，得到的结果如下。 f=4*a (4) 第一步：在命令窗口输入如下程序。 syms x y=sqrt(x+1)/(sqrt(x+1)+sqrt(3-x)) f=int(y,x,0,2) 第二步：按“enter”键，得到的结果如下。 y =(x + 1)^(1/2)/((x + 1)^(1/2) + (3 - x)^(1/2)) 5.3.4 级数求和 symsum函数用于级数的求和。该函数的调用格式为： r=symsum(s) 自变量为findsym函数所确定的符号变量，设其为k，则该表达式计算 s从0到k－1的和； r=symsum(s,v) 计算表达式s从0到v－1的和； r=symsum(s,a,b) 计算自变量从a到b之间s的和； r=symsum(s,v,a,b) 计算v从a到b之间的s的和。其中s表示一个级数的通项，是一个符号表达式。v是求和变量，v省略时使用系统的默认变量。a和b是求和的开始项和末项。 【例5-17】求级数前6项的部分和s。 解： 第一步：在命令窗口输入如下程序。 syms x s=symsum(1/(1+x^2),1,6) 第二步：按“enter”键，得到的结果如下。 s =75581/81770 MATLAB提供了taylor函数将函数展开为幂级数，所以函数 taylor 可以用于实现 Taylor级数的计算。该函数的调用格式为： r=taylor(f) 计算表达式f的Taylor级数，自变量由findsym函数确定，计算f的在0的6阶Taylor级数； r=taylor(f,n,v) 指定自变量v和阶数n； r=taylor(f,n,v,a) 指定自变量v、结束n，计算f在a的级数。 该函数将函数f按变量v展开为泰勒级数，展开到第n项(即变量v的n-1次幂)为止，n的缺省值为6。v的缺省值与diff函数相同。参数a指定将函数f在自变量v＝a处展开，a的缺省值是0。 【例5-18】求函数的泰勒级数展开式。 (1) 求y＝sin(sinx)的三阶泰勒级数展开式。 (2) 求y＝tan(x)的2阶泰勒级数展开式。 (1) 第一步：在命令窗口输入如下程序。 syms x y=sin(sin(x)) f=taylor(y,x,4) 第二步：按“enter”键，得到的结果如下。 f=x - x^3/3 (2) 第一步：在命令窗口输入如下程序。 syms x y=tan(x) f=taylor(y,x,3) 第二步：按“enter”键，得到的结果如下。 f=x 5.4　符号线性代数 5.4.1　基本符号运算 在Matlab数学运算中，基本符号运算符包括＋、-、*、.*、\、.\、/、^、.^。 (1) 加减运算。 A＋B，A-B。如果A、B为同类型的矩阵时，A＋B，A-B分别对对应分量进行加减。 (2)??A*B符号矩阵乘法。 A*B为线性代数中定义的矩阵乘法。按乘法定义要求必须有矩阵A的列数等于矩阵B的行数，否则该函数返回错误信息。 (3)??A.*B符号数组的乘法。 A.*B为按参量A与B对应的分量进行相乘。A与B必须为同型阵列，或至少有一个为标量。 (4)??A\B矩阵的左除法。 X=A\B为符号线性方程组A*X=B的解。我们指出的是，A\B近似地等于inv(A)*B。若X不存在或者不唯一，则产生一警告信息。矩阵A可以是矩形矩阵(即非正方形矩阵)，但此时要求方程组必须是相容的。 (5)??A.\B数组的左除法。 A.\B为按对应的分量进行相除。若A与B为同型阵列时。若A与B中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行操作。 (6)??B/A矩阵的右除法。 X=B/A为符号线性方程组X*A=B的解。需要指出的是，B/A粗略地等于B*inv(A)。若X不存在或者不唯一，则产生一警告信息。矩阵A可以是矩形矩阵(即非正方形矩阵)，但此时要求方程组必须是相容的。 (7)??A./B数组的右除法。 A./B为按对应