第五章MATLAB符号运算.pptVIP

第5章 MATLAB符号运算;5.1 符号对象及其表达方式 ;定义符号常量 ;例5.1 将一组数值常量定义成符号常量。 log(2) %数值常量 f1=sym(log(2))　　　　　　 　 %符号常量，注意f1在工作空间中的类型 ;定义符号变量 ;;5.1.2 符号表达式 ;; f4=r^2+sin(x)+cos(y)+log(s)+exp(t); g=sin(a)+cos(b) ; % f4和g均为字符串变量。 ;5.1.3 符号矩阵 ;;例5.8 用size()函数求符号矩阵的大小 s=size(M) [s_r,s_c]=size(M) s_r=size(M,1) s_c=size(M,2);例5.9 抽取符号矩阵中的元素 B=sym([a,b,c;d,e,f;g,h,k]) B(2,3),B(6) ;5.2 符号算术运算 ;;A*B、A/B可分别用来求A和B两个符号矩阵的乘法与除法。A.*B则用来实现两个符号数组的乘法。其中矩阵除法也可用来求解符号线性方程组的解。 syms a b c d e f g h; A = [a b; c d]; B = [e f; g h]; C1 = A.*B D1=A*B;;5.3.1 表达式中的独立变量 ;findsym( )查找MATLAB的符号表达式或矩阵中哪一个变量是独立变量 r = findsym(S)：以字母表的顺序返回表达式S中的所有符号变量（注：符号变量由除了i与j的字母与数字构成的、字母打头的字符串）。若S中没有任何的符号变量，则findsym返回一空字符串。 r = findsym(S,n)：返回表达式S中依接近x的顺序而排列的n个符号变量 ;syms a x y z t alpha beta S1 = findsym(sin(pi*t*alpha+beta),1) S2 = findsym(x+i*y-j*z+eps-nan) S2 = NaN, x, y, z S3 = findsym(a+y,2) S3 = y,a ;5.3.2 表达式化简 ;1. 合并同类项（collect） ;;2. 多项式展开（expand） ;3. 因式分解（factor） ;;4. 一般化简（simplify） ;5. 不定化简（simple） ;6. 通分（numden） ;7. 书写格式美化（pretty） ;5.4.1 符号极限 ;例5-22 求极限示例。 syms x a t h n; L1 = limit((cos(x)-1)/x) L2 = limit(1/x^2,x,0,right) 例5-23 求： 、 的极限。 syms x y f=sin(x)/x; %表达式赋值 g=sin(x+2*y); %表达式赋值 fx=limit(f) %求f(x)的极限 gx=limit(g,y,0) %求f(x)的极限 ;5.4.2 符号微分 ;例5-24 求一次符号微分示例。 syms x n y=sin(x)^n*cos(n*x); Xd=diff(y) ，Nd=diff(y, n) 例5-25 求二次符号微分示例。 syms t f=exp(-t)*sin(t); diff(f,t,2) ;5.4.3 符号积分 ;例5-27 求积分示例。 syms x z t alpha INT1 = int(-2*x/(1+x^3)^2) INT2 = int(x/(1+z^2),z) 积分比微分复杂得多。积分不一定是以封闭形式存在，或许存在但软件也许找不到，或者软件可明显地求解，但超过内存或时间限制。当MATLAB不能找到积分表达时，它将返回未经计算的函数形式。 ;5.4.4 符号Taylor级数展开 ;解析函数f(x)在点x=a的Taylor级数定义为： ;;