《MATLAB 8.X程序设计及典型应用》课件第5章.pptVIP

5.4符号方程的求解

5.4.1符号代数方程的求解

MATLAB求解代数方程或代数方程组的指令为solve()，调用格式为

S=solve(eq，v)

计算单个方程eq=0关于变量v的解S。eq可以是含等号的符号方程。v缺省时MATLAB将利用指令findsym()找寻第一自由变量作为求解变量。

[y1，y2，…，yn]=solve(‘eq1’，‘eq2’，…，‘eqn’，‘v1’，‘v2’，…，‘vn’)

计算方程组eq1=0，eq2=0，…，eqn=0关于变量v1，v2，…，vn的解y1，y2，…，yn。如果只有一个输出变量，则该输出变量类型为构架数组，MATLAB把各个解分别存放在此构架数组的不同域中，每个域值为对应域名的变量的解。v缺省时MATLAB将利用指令findsym()找寻n个自由变量作为求解变量。5.4.2符号微分方程的求解

MATLAB求解符号微分方程和符号微分方程组的指令为dsolve()，调用格式为：

dsolve(‘eq1’，‘eq2’，…，‘eqn’，‘condl’，‘cond2’，…，‘condn’，‘v’)

计算带有初始条件cond1，cond2，…，condn的微分方程组eq1，eq2，…，eqn的解。v为微分变量，缺省时MATLAB将以t作为微分变量。【例5-35】求二阶微分方程

在初始条件下的特解。5．符号表达式的化简

符号表达式可以用许多等价的形式来表示。为了让用户得到最简洁的表达形式，MATLAB提供了两个指令进行符号化简：simplify()和simple()。其调用格式相同：

simplify(S)或simple(S)

对符号表达式S化简，输出最简化形式。二者的区别在于对于较为复杂的表达式化简，不带输出变量的simple()往往给出符号表达式的化简过程，而simplify()只给出化简后的最简形式。【例5-17】化简指令simplify()和simple()使用实例：化简表达式。【例5-19】已知矩阵A=，计算A的平方，以及A中各元素的平方，并化简。5.2.4符号表达式的置换

MATLAB符号数学工具箱中提供了用于符号表达式的置换函数，用于简化表达式的输出，方便用户阅读。

subs(ES，OLD，NEW)：用变量NEW替换符号表达式ES中的变量OLD。

[Y，X]=subexpr(ES，X)：自动查找表达式中最长的子表达式，并用符号变量X来置换，输出置换后的符号表达式Y和子表达式X。【例5-21】置换函数subexpr()使用实例：计算符号矩阵的逆矩阵。5.3符?号?微?积?分

MATLAB符号计算工具箱中提供了很多关于符号微积分的运算函数，用户只要正确调用这些函数就可以准确地进行符号微积分的运算，减轻了手工运算的劳苦，更避免了因疏忽引起的差错。5.3.1符号极限和符号微分

1.符号极限

MATLAB计算符号极限的函数为limit()。其格式如下：由执行结果可知：2.符号微分

MATLAB中可用来计算符号表达式微分的指令形式多样，主要有diff()、jacobian()和taylor()。diff()的调用格式为：

dfdtn=diff(f，t，n)

计算符号表达式f对符号变量t的n阶微分(即计算

)。n缺省时默认为计算一阶微分1；t缺省时按照findsym()找寻的第一自变量求微分。jacobian()的调用格式：

fj=jacobian(f(t))

计算多元符号表达式f(t)的一阶微分雅克比(Jacobian)矩阵。雅可比矩阵定义为向量对向量的一阶微分矩阵：taylor()的调用格式如下：

ftlr=taylor(f，n，t，a)

计算符号表达式f(t)在t=a处的(n-1)阶泰勒级数展开

(即计算)。【例5-23】已知函数矩阵5.3.2符号级数/序列求和与符号积分

1．级数/序列的符号求和

级数或序列求和是高等数学中常见的运算，MATLAB中的实现函数为symsum()，格式如下：

Sum=symsum(f，t，a，b)

计算序列/级数在指定变量t取区间[a，b]中所有整数时的和(即计算)。如果t缺省，MATLAB按照findsym()找寻第一自变量。t可以取有限值，也可以是无穷大。用户如果不指定t的区间，MATLAB将按照默认计算区间求和。【例5-26】序列/级数求和指令sy