二讲Matlab符号计算12.pptVIP

第二讲 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 ——微分方程;例2 定义一个符号函数 fxy=(a*x^2+b*y^2)/c^2 ，分别求该函数对x、y的导数和对x的积分。;例3 比较符号常数与数值在代数运算时的差别。; 符号表达式由符号变量、函数、算术运算符等组成。符号表达式的书写格式与数值表达式相同。;例4 已知一复数表达式 z=x+i*y，试求其共轭复数，并求该表达式与其共轭复数乘积的多项式。 为了使乘积表达式 x^2+y^2非负，这里，把变量x和y定义为实数。 x=sym(x,real);y=sym(x,real);; syms函数;（二） 符号函数的生成; 符号表达式中变量的确定;例6 查询符号函数 f=x^n g=sin(at+b) 中的系统缺省变量。;（三） 符号矩阵的创建;用sym函数创建符号矩阵;将数值矩阵转化为符号矩阵;A =sym([1/3,5/2;10/7,2/5]) eval(A); 符号矩阵的修改;syms a x; f=a*sin(x)+5 f = a*sin(x)+5 f1=subs(f,sin(x),sym(y)) f1 = a*y+5 f2=subs(f,{a,x},{2,pi/3}) f2 = 3^(1/2)+5 f3=subs(f,{a,x},{2,pi/3}) f3 = 6.7321 f4=subs(subs(f,a,2),x,0:pi/6:pi) f4 = 5.0000 6.0000 6.7321 7.0000 6.7321 6.0000 5.0000 f5=subs(f,{a,x},{0:6,0:pi/6:pi}) f5 = 5.0000 5.5000 6.7321 8.0000 8.4641 7.5000 5.0000; 实际上，Matlab的符号数学工具箱提供了两个符号表达式的替换函数subexpr和subs，可以通过符号替换使表达式的输出形式简化，以得到一个简单的表达式。;例12 求解并化简 3 次方程 x3+ax+1=0 的符号解。;[r,s]=subexpr(t,s);数值运算中，所有矩阵运算操作指令都比较直观、简单。而符号运算就不同了，所有涉及符号矩阵运算的操作都有专用符号矩阵函数来进行，这些函数作用于单个的数据无意义。;例13 计算3阶范得蒙矩阵行列式的值。设A是一个由符号变量a,b,c确定的范得蒙矩阵。 ;二、符号表达式运算;例14 符号表达式的四则运算示例。;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);（四）符号表达式的化简;（五）符号表达式嵌套形式重写;（六）任意精度的数学运算 ;digits(25) vpa