东华大学matlab数学实验课件chapter7.pptVIP

第七章 MATLAB符号计算 7.1 符号对象 7.2 符号矩阵和符号函数 7.3 符号微积分 7.4 符号方程和符号微分方程 7.6 符号计算局限性和Mupad的调用 7.7 习题 7.5 便捷函数作图 7.1 符号对象 1、符号对象的定义 s=sym(arg) 将数值或字符串arg转化为符号对象s,数值为有理表示; s=sym(num,’d’) 将数值表达式转化为符号表达式，数值用十进制表示； syms v1 v2 ? 定义v1,v2,?为符号变量； subs(s,old,new) 将符号表达式s中的符号变量old用new代替 c=sym(pi^2); c = pi^2 syms x y z; %定义x,y,z为符号， d=x^3+2*y^2+c d = x^3+2*y^2+pi^2 A=subs(d,x,c) A = pi^6+2*y^2+pi^2 A=subs(A,y,0.1) A = 971.2788 中间不能加逗号 将表达式d中的x用c替代 将A中的y用0.1替代 符号计算默认精度为32位，是MATLAB数值计算精度的2倍。 digits(n) 将数值计算精度设为n位； x=vpa(s) 求s的数值结果； x=vpa(s,n) 求出s的数值结果，精确到n位数； double(s) 将符号对象转化为双精度数值； char(s) 将符号对象转化为字符串； 2、计算精度和数据类型转换 format long;a=pi^2;format short ; a = 9.86960440108936 vpa(a,5) ans = 9.8696 syms x;f=sym(x^3-1) f = x^3-1 g=char(f) fplot(g,[0 2]) syms x y; sym(…) f= x^3-1 g=‘x^3-1’ 7.2 符号矩阵和符号函数 1、矩阵 clear;A=sym([a,b;c,d]) A = [ a, b] [ c, d] b=inv(A) b = [ d/(a*d-b*c), -b/(a*d-b*c)] [ -c/(a*d-b*c), a/(a*d-b*c)] 符号矩阵运算 factor(expr) 对expr作因式分解 expand(expr) 将expr展开 collect(expr,v) 将expr按变量v合并同类项 simplyfy(expr),simple(expr)将expr化简 g=finverse(f,v) 求函数f(v)的反函数g(v) fg=compose(f,g) 求函数f(v)和g(v)的复合函数f(g(v)) [n,d]=numden(expr)将分式通分，n返回分子，d返回分母 funtool 函数分析图形界面 例：已知f(x)=(x-1)3,g(x)=(x+1)3 求h= f*g 并展开。 解：syms x;f=(x-1)^3;g=(x+1)^3; h=f*g h = (x-1)^3*(x+1)^3 s=expand(h) s = x^6-3*x^4+3*x^2-1 7.3 符号微积分? 1．极限和级数 limit(s,x,a) 返回符号表达式s当x-a时的极限 symsum(s,n,a,b) 返回符号表达式s表示的通项当自变量n由a到b的和。 例： syms x n s1=(1+x/n)^n, s2=(-1)^n*(x^n)/n limit(s1,n,inf), symsum(s2,n,1,inf) diff(s,x) 返回符号表达式s对x的导函数 diff(s,x,n) 返回s对x的n阶导函数。 taylor(s,n,a) 返回s在a点TAYLOR展开到n-1次式. taylortool TAYLOR分析图形界面 jacobian(f,x) 返回向量函数f的Jacobian矩阵 2. 微分 Jacobian矩阵 (3)函数g(x)=lnxsinx在x=1的Taylor展开 例 计算：(1) (2) 的Jacobian矩阵； clear;syms x y; s=x^2*exp(-y);diff(s,x,2) dx=diff(s,x),dxdy=diff(dx,y) t=subs(dxdy,x,1),t=subs(t,y,2) syms x1 x2;f=[x1*exp(x2);cos(x1)]; J=jacobian(f,[x1 x2]) syms x;g=log(x)*sin(x); gt=taylor(g,6,1), gt=vpa(gt,5)