（新）4 matlab的符号运算.pptVIP

如果Z为复矩阵，则Z为它的复数共轭转置，非共轭转置使用Z. 或conj(Z)求得． * 奇异值分解是线性代数中一种重要的矩阵分解，在信号处理、统计学等领域有重要应用。 1、奇异值分解非常有用，对于矩阵A(m*n)，存在U(m*m)，V(n*n)，S(m*n)，满足A = U*S*V’。U和V中分别是A的奇异向量，而S是A的奇异值。AA的正交单位特征向量组成U，特征值组成SS，AA的正交单位特征向量组成V，特征值（与AA相同）组成SS。因此，奇异值分解和特征值问题紧密联系。 2、奇异值分解提供了一些关于A的信息，例如非零奇异值的数目（S的阶数）和A的秩相同，一旦秩r确定，那么U的前r列构成了A的列向量空间的正交基。 矩阵的奇异值分解为: A=USV, 其中U,V是正交阵(所谓B为正交阵是指B=B-1, 即BB=I), S为对角阵. * * syms x; f1=cos(x*x); diff(f1) %未指定求导变量和阶数，按默认规则处理 ans = (-2)*x*sin(x^2) diff(f1,x,2) %求f1对x的二阶导数 ans = - 2*sin(x^2) - 4*x^2*cos(x^2) diff(f1,x,3) %求f1对x的三阶导数 ans = 8*x^3*sin(x^2) - 12*x*cos(x^2) 求(1) syms a b t; f21=a*(t-sin(t)); f22=b*(1-cos(t)); diff(f22,t)/diff(f21,t) %求y对x的一阶导数 ans = -(b*sin(t))/(a*(cos(t) - 1)) 求(2) syms x y; f3= x^6-3*y^4+2*x^2*y^2; diff(f3,x) %求z对x的偏导数 ans = 6*x^5 + 4*x*y^2 diff(f3,y) %求z对y的偏导数 ans = 4*x^2*y - 12*y^3 diff(diff(f3,x),y) ans = 8*x*y * 求(3) 符号积分由函数int来实现，一般调用格式如下。 int(s)：没有指定积分变量和积分阶数时，系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分。 int(s,v)：以v为自变量，对被积函数或符号表达式s求不定积分。 int(s,a,b):求符号表达式s对默认符号变量从a到b得定积分。a和b为双精度或符号数量。 * 符号积分 §3.2 符号微积分 int(s,v,a,b)：求定积分运算。a、b分别表示定积分的下限和上限。该函数求被积函数在区间[a，b]上的定积分。a和b可以是两个具体的数，也可以是一个符号表达式，还可以是无穷（inf）。 当函数f关于变量x在闭区间[a，b]上可积时，函数返回一个定积分结果。 当a、b中有一个是inf时，函数返回一个广义积分。 当a、b中有一个符号表达式时，函数返回一个符号函数。 * 【例】分别求下列积分 （1） （2） （3） （4）求 的不定积分。 * f=sym(1/(1+x^2)); f1=int(f) %求不定积分 f1 = atan(x) f2=int(f,a,b) %求定积分 f2 = atan(b) - atan(a) f3=int(f,1,2) %求定积分 f3 = atan(2)-1/4*pi eval(f3) %计算积分值 ans = 0.3218 （1） （2） （3） syms x t; %求符号矩阵的不定积分 A=[cos(x*t),sin(x*t);-sin(x*t),cos(x*t)]; int(A,t) ans = [ 1/x*sin(x*t), -cos(x*t)/x] [ cos(x*t)/x, 1/x*sin(x*t)] * 求（4） 的不定积分。 以t为自变量 级数 掌握级数求和的方法和将函数展开为泰勒级数的方法。 级数 级数符号求和 函数的泰勒级数 * 求无穷级数的和需要符号表达式求和函数symsum，其调用格式为： symsum(s):