MATLAB符号运算运用幻灯片.ppt

第3章 MATLAB的符号运算 3.1 符号表达式的建立 3.1.2 创建符号变量和表达式 3.1.3 符号矩阵 3.2 符号表达式的代数运算 3.2.1 符号表达式的代数运算 【例】求矩阵的行列式值、非共轭转置和特征值。 【例】符号表达式f=2x2+3x+4与g=5x+6的代数运算。 3.2.2 符号表达式的操作和转换 1. 自由变量的确定原则 小写字母i和j不能作为自由变量。 符号表达式中如果有多个字符变量，则按照以下顺序选择自由变量：首先选择x作为自由变量；如果没有x，则选择在字母顺序中最接近x的字符变量；如果与x相同距离，则在x后面的优先。 大写字母比所有的小写字母都靠后。 2. findsym函数 如果不确定符号表达式中的自由符号变量，可以用findsym函数来自动确定。 语法： findsym(f,n) %确定自由符号变量 说明：f可以是符号表达式或符号矩阵；n为按顺序得出符号变量的个数，当n省略时，则不按顺序得出f中所有的符号变量。 3、符号表达式的替换 MATLAB 中，可以通过符号替换使表达式的形式简化。符号工具箱中提供了两个函数用于表达式的替换：1．subexpr 该函数自动将表达式中重复出现的比较长的子表达式或字符串用变量替换，该函数的调用格式为： subexpr(s,s1)，指定用符号变量 s1 来代替符号表达式s（可以是矩阵）中重复出现的字符串。替换后的结果由 ans 返回，被替换的字符串由 s1返回； [Y,s1] = subexpr(X,‘s1’)，该命令与上面的命令不同之处在于第二个参数为字符串，该命令用来替换表达式中重复出现的字符串。 2. subs 函数 subs 可以用指定符号替换表达式中的某一特定符号。 subs(s) subs(s, new) subs(s,old,new) 5、 符号表达式与多项式的转换 3.3 符号极限、微积分和级数求和 3.3.2 符号微分 3.3.3 符号积分 3.3.4 符号级数 【例】求级数  1+x+x2+…+xk+…的和。 syms x k s1=symsum(1/k^2,1,10) %计算级数的前10项和 s1 = 1968329/1270080 s2=symsum(1/k^2,1,inf) %计算级数和 s2 = 1/6*pi^2 s3=symsum(x^k,k,0,inf) %计算对k为自变量的级数和 s3 = -1/(x-1) 【例】求ex的泰勒展开式 syms x s1=taylor(exp(x),8) %展开前8项 s1 = 1+x+1/2*x^2+1/6*x^3+1/24*x^4+1/120*x^5+1/720*x^6+1/5040*x^7 s2=taylor(exp(x)) %默认展开前5项 s2 = 1+x+1/2*x^2+1/6*x^3+1/24*x^4+1/120*x^5 3.4 符号方程的求解 3.4.2 符号常微分方程 3. 5 符号积分变换 【例】计算f(t)= 的fourier变换F以及F的fourier反变换。 syms t w F=fourier(1/t,t,w) %fourier变换 F = i*pi*(Heaviside(-w)-Heaviside(w)) f=ifourier(F,t) %fourier反变换 f = 1/t f=ifourier(F) %fourier反变换默认x为自变量 f = 1/x 程序分析：其中Heaviside(t)是单位阶跃函数 ，函数名为数学家Heaviside的名字。 3.5.2 拉普拉斯(Laplace)变换 及其反变换 【例】求sin(at)和阶跃函数的Laplace变换。 syms a t s F1=laplace(sin(a*t),t,s) %求sinat的Laplace变换 F1 = a/(s^2+a^2) F2=laplace(sym(Heaviside(t))) %求阶跃函数的Laplace变换 F2 = 1/s 3.6 符号数值任意精度控制和运算 符号计算是对未赋值的符号对象(可以是常数、变量、表达式)进行运算和处理。MATLAB具有符号数学工具箱(Symbolic Math Toolbox)，将符号运算结合到MATLA