第三章_Malab符号计算.pptVIP

第3章 MATLAB符号计算 主讲人：李良群 ; 所谓符号计算是指在运算时，无须事先对变量赋值，而将所得到结果以标准的符号形式来表示。;符号计算是精确计算，可以获得任意精度的解； Symbolic Math Toolbox——符号运算工具包通过调用Maple软件实现符号计算的。 Maple软件——主要功能是符号运算，它占据符号软件的主导地位。　　 ;一、符号对象的创建;如: f1=a?x^2+b?x+c —— 二次三项式 f2= a?x^2+b?x+c=0 —— 方程 f3=Dy+y^2=1 ——微分方程; 符号表达式由符号变量、函数、算术运算符等组成。符号表达式的书写格式与数值表达式相同。; syms函数;（二） 符号函数的生成; 符号表达式中变量的确定;例3 查询符号函数 f=x^n g=sin(at+b) 中的系统缺省变量。;（三） 符号矩阵的创建;用sym函数创建符号矩阵;将数值矩阵转化为符号矩阵;A =sym([1/3,5/2;10/7,2/5]) eval(A);（四）任意精度的数学运算 ;digits(25) vpa(1/2+1/3);二、符号表达式运算;例9 符号表达式的四则运算示例。;f=cos(x);g= sin(2*x); symop(f,/,g,+,f,*,g) ans = cos(x)/sin(2*x)+cos(x)*sin(2*x);syms x f=2*x^2+3*x-5; g= x^2+x-7; h=f+g;该函数将符号表达式s转换为分子和分母都是整系数的最佳多项式，并分别存放在n与d中。;factor(S)——对S分解因式，S是符号表达式或符号矩阵; pretty(S)——将符号表达式S按照类似书写习惯的方式显示; expand(S)——对S进行展开，S是符号表达式或符号矩阵; collect(S)——对S合并同类项，S是符号表达式或符号矩阵; collect(S,v)——对S按变量v合并同类项，S是符号表达式或符号矩阵。;syms x f=(x+1)^5; expand(f);syms x t; f=x*(x*(x-6)+12)*t; collect(f) collect(f,t);（四）符号表达式的化简;（五）符号表达式嵌套形式重写;三、符号微积分;例18 求极限 ;例19 求极限。 syms a m x; f=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/(x+a); limit(f,x,a) ans = (1/2*a*exp(sin(a))+1/2*a-exp(tan(a))+1)/a syms x t; limit((1+2*t/x)^(3*x),x,inf) ans = exp(6*t) syms x; f=x*(sqrt(x^2+1)-x); limit(f,x,inf,left) ans = 1/2 syms x; f=(sqrt(x)-sqrt(2)-sqrt(x-2))/sqrt(x*x-4); limit(f,x,2,right) ans = -1/2;（二）符号导数;例21 求函数的导数。;f1=a*cos(t);f2=b*sin(t); diff(f2)/diff(f1) %按参数方程求导公式求y对x的导数 (diff(f1)*diff(f2,2)-diff(f1,2)*diff(f2))/(diff(f1))^3 %求y对x的二阶导数;（三）符号积分;例22 求不定积分。;x=sym(x);t=sym(t); int(abs(1-x),1,2) %求定积分 f=1/(1+x^2); int(f,-inf,inf) %求定积分 int(4*t*x,x,2,sin(t)) %求定积分 f=x^3/(x-1)^100; I=int(f,2,3