[理学]第二讲 2010年数学建模竞赛培训matlab符号计算.pptVIP

Matlab与科学计算 第二讲 符号推演与计算 引言 符号表达式的定义和数据转换 在MATLAB中为进行符号推演就必须区别数值表达式和符号表达式的区别。 MATLAB提供了一系列按照符号推演规则下的符号命令，如定义符号、定义一个多元函数所有参数中那些是自变量那些是参数等等。当我们在符号推演的过程或结束后，希望计算函数在某些点处的值，我们有一些有关的辅助命令如vpa和subst 定义符号对象的指令 一个符号表达式中，所有的参数都是符号变量或符号常量，因此必须和在计算机语言中的那些表达式有所区别。我们可以用命令sym和syms来定义符号变量和符号表达式。sym 和 syms 常用的格式如下： f = sym(expr) 把表达式expr 转换为符号对象。 syms(‘arg1’, ‘arg2’, ‘arg3’) 将他们定义为符号变量。 syms arg1 arg2 arg3 为上面命令的简写形式 符号表达式中自由变量的确定 在一个符号表达式中可能有多个符号变量，那么哪一个是自变量，哪些是符号常量则对某种运算是非常重要的. MATLAB 提供了一个 findsym 命令，可以对所有符号变量指定一定数量的自变量进行自动认定。 基本语法： findsym（expr） 确定表达式 expr 中所有自由符号变量 findsym（expr，n） 确定expr 靠 x 最近的n个独立变量 符号的数值化和替代 在进行公式推导的过程中，有时我们需要将公式按某些实际参数来进行数值化表示，或在自变量取某些值时求符号公式的值，这时我们可以使用符号数值化和替换等命令。 对符号求值的命令为vpa，即（Variable precision arithmetic）。语法为： R = vpa(A) R = vpa(A,d) 这里 A：对符号表达式A求给定精度的值。 d：输出数值的有效位数 微积分中的符号运算 微积分中可以进行极限、导数、微分、积分、级数展开等解析运算，也可以进行多元函数的微积分运算，结合图形的显示可以更好地帮助我们理解空间微积分的概念和计算。 求极限 求极限即对表达式进行求极限，首先是定义符号表达式，然后对表达式进行极限运算。极限运算的命令语法为： 导数 diff(S) 求函数 S的导数 diff(S,‘v’) 对函数 S中的符号v进行求导运算 diff(S,n) 求函数 S的n阶导数 diff(S,‘v’,n) 对函数 S中的符号v求n阶导数 syms x t diff(sin(x^2)+x*t,x) diff(sin(x^2)+x*t,t) ans = 2*cos(x^2)*x+t ans = x 符号积分 符号积分包括不定积分和定积分，其语法分别为： R = int(S) R = int(S,v) R = int(S,a,b) R = int(S,v,a,b) Taylor展式 泰勒展式命令语法为： r = taylor(f) r = taylor(f,n,v) r = taylor(f,n,v,a) 其中 f：被展开函数，如没有其他参数表示系统内定在0点展开 v：对变量v进行泰勒展开 n：取展式的前n项,系统默认阶数为5 a：在a点展开，若不给出，系统内定为在0展开 符号求和 符号级数求和的命令语法为： r = symsum(s) r = symsum(s,v) r = symsum(s,a,b) r = symsum(s,v,a,b) 其中： s：为含变量k的表达式 v：对表达式s中的变量v求和 b：对表达式中变量k从a到b求和 quiver3函数 向量作图函数，语法为： quiver3(X,Y,Z,U,V,W) quiver3(Z,U,V,W) quiver3(...,scale) quiver3(...,LineSpec) quiver3(...,LineSpec,filled) 这里X，Y，Z为空间的点的坐标，它们为同维数组。 U，V，W为相应点的方向向量， Scale：为方向向量的长度。 矩阵代数中的符号计算 矩阵代数中的符号计算包括对矩阵进行各种四则运算、求矩阵的行列式、矩阵的逆、矩阵的谱分解以及求线性方程组的解等。部分命令和含义如下 符号表达式的操作命令 对符号表达式的操作，就象我们平时在数学课上对多项式的推演一样。这里包括多项式的加、减、乘、除，对多项式的化简等等。 多项式的基本操作 求函数的反函数 求函数的反函数命令为finverse，其语法为： g = finverse(f)，求以x为自变量的函数f的反函数 g = finverse(f,v),求以v为自变量的函数f的反函数 例：finverse(1/tan(x)) A