第五讲 MATLAB的符号计算功能.pptVIP

MATLAB第五讲 MATLAB的符号计算功能 一、符号表达式和符号矩阵的创建 1、 符号变量的赋值 2、符号常量 3、符号矩阵的创建 4、符号方程的创建 5、将数字矩阵转化为符号矩阵 6、符号表达式的升幂 7、符号表达式的合并 8、变量代换 二、符号函数的微积分 1、求导数 2、求积分 3、 求泰勒级数taylortool(f) 三、求符号方程的精确解和近似解 1、求解单个符号方程 2 、对代数方程组求解 四、符号矩阵的基本运算 1、加、减、乘 2、求逆运算 3、求符号矩阵的除法 4、求矩阵的特征值和特征向量 5、求符号矩阵的行列式 6、求符号矩阵的约当矩阵 7、求符号矩阵的奇异值 五、符号函数画图 六、求符号函数的零点 七、求微分方程的解 一、入门 Symbolic工具包 MATLAB有一个符号计算工具包叫作 Symbolic Math Toolbox 其中有60多个专用函数。包括微积分、线性代数、方程求解、多项式的简约与展开、特殊数学函数等。 sym的指令在C:\MATLAB6p1\toolbox\symbolic，如果在搜索路境中没有设定则要添加。 符号表达式和符号矩阵的创建 符号表达式是数字、函数、变量的MATLAB字符串，或字符串数组。 符号运算是指使用已知的规则和给定的符号恆等式求解符号方程。 生成符号表达式用引号或sym函数 例 1、符号变量的赋值 f1=sin(x)^2’；f2=exp(-x^2/2)’；f3=1/(1+x^2)’； f1 = sin(x)^2 f2 = exp(-x^2/2) f3 = 1/(1+x^2) 符号常量 没有变量的符号表达式叫作符号常量 5、isstr( )用来检测变量是否符号变量 例 6、符号表达式的升幂 7、符号表达式的合并 8、变量代换 例1计算符号函数的导数 f=‘sin(x)^2’ %定义函数的符号表达式 f = sin(x)^2 diff(f) ans = 2*sin(x)*cos(x) diff(f,2) ans = 2*cos(x)^2-2*sin(x)^2 f=‘sin(x)^2’ ； int(f,x) ans = -1/2*sin(x)*cos(x)+1/2*x int(1/(1+x^2)) ans = atan(x) 在框中可以交互作用，给出所需的阶数，立即返回表达式和图形。 几个常用命令 Solve(‘方程’）%求精确解 vpa(S,n) %求n位有效数字的近似解 numeric(S) %将不含自由自变量的近似解转 化为数值解 digits(n) % 设定近似解的有效位数 subs(S,Dsym,Fsym) %将数值Fsym带入自变 量Dsym R1=solve(x^2-x-1=0) R1 = [ 1/2*5^(1/2)+1/2] [ 1/2-1/2*5^(1/2)] RV=vpa(R1) RV = [ 1.6180339887498948482045868343657] [ -.61803398874989484820458683436570] RV4=vpa(R1,4) RV4 = [ 1.618] [ -.6180] RV30=vpa(R1,16) RV30 = [ 1.618033988749895] [ -.6180339887498950] numeric(R1) ans = 1.6180 -0.6180 1、 加、减、乘 symadd(A,B) %符号加 symsub(A,B） %符号减 symmul(A,B) %符号乘 例 求二阶符号矩阵的逆 M=sym([a,b;c,d]) M = [ a, b] [ c, d] inv(M) ans = [ d/(a*d-b*c), -b/(a*d-b*c)] [ -c/(a*d-b*c), a/(a*d-b*c)] fzero(fun,[a,b],[ ],options) 在求零点之前先要确定一个根的隔离区间，也就是说，[a,b]区间内有且仅有函数的一个零点，这里采用的是对分法求零点。 fzero(f