[工学]32符号运算.pptVIP

（1）利用单引号来生成符号表达式 f = sin(x)+5x f —— 符号变量名 sin(x)+5x—— 符号表达式 —— 符号标识 符号表达式一定要用 单引 号括起来matlab才能识别。 （2）用sym函数建立符号表达式 数值矩阵A=[1,2;3,4] A=[a,b;c,d] —— 不识别 ?用matlab函数sym创建矩阵 命令格式：A=sym([ ]) ※ 符号矩阵内容同数值矩阵 ※ 需用sym指令定义 ※ 需用 标识 例如：A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0] 这就完成了一个符号矩阵的创建。 注意：符号矩阵的每一行的两端都有方 括号，这是与 matlab数值矩阵 的一个重要区别。 ?用字符串直接创建矩阵 模仿matlab数值矩阵的创建方法 需保证同一列中各元素字符串有相 同的长度。 ?将数值矩阵转化为符号矩阵 函数调用格式：sym(A) A=[1/3,2.5;1/0.7,2/5] A = 0.3333 2.5000 1.4286 0.4000 sym(A) ans = [ 1/3, 5/2] [10/7, 2/5] 3.2.2 符号运算中的运算符号 和基本函数 1、符号矩阵运算 数值运算中，所有矩阵运算操作指 令都比较直观、简单。例如：a=b+c; a=a*b ；A=2*a^2+3*a-5等。 所有涉及符号运算的操作同数值运算几乎完全相同。差别仅为符号运算中无逻辑运算功能、关系运算符号仅有“==”和“ ~=”。事实为真时，结果为1。 2.2 符号积分变换 常见的积分变换有傅立叶变换、拉普拉斯变换和Z变换。 1．傅立叶(Fourier)变换 在MATLAB中，进行傅立叶变换的函数是： fourier(f,x,t)：求函数f(x)的傅立叶像函数F(t)。 ifourier(F,t,x)：求傅立叶像函数F(t)的原函数f(x)。 例1.已知系统矩阵A= 计算矩阵指数函数eAt。 用拉氏反变换法计算eAt的公式为： eAt = L-1[(SI-A)-1] a=[0 1;-2 -3]; syms s b=(s*eye(2)-a) b = [ s, -1] [ 2, s+3] B=inv(b) B= [ (s+3)/(s^2+3*s+2), 1/(s^2+3*s+2)] [ -2/(s^2+3*s+2), s/(s^2+3*s+2)] E= ilaplace(B) E= [ -exp(-2*t)+2*exp(-t), exp(-t)-exp(-2*t)] [ -2*exp(-t)+2*exp(-2*t), 2*exp(-2*t)-exp(-t)] B= ilaplace(E) B= [ (s+3)/(s^2+3*s+2), 1/(s^2+3*s+2)] [ -2/(s^2+3*s+2), s/(s^2+3*s+2)] 3．Z变换 当函数f(x)呈现为一个离散的时域数列f(n）时，对数列f(n)进行z变换的MATLAB函数是： ztrans(fn)：求fn的Z变换像函数F(z)。 iztrans(Fz)：求Fz的z变换原函数f(n)。 syms x f=x*exp(-x*10); F=ztrans(f) F=ztrans(x*exp(-x*10)) F = z*exp(-10)/(z-exp(-10))^2 fn= iztrans(F) fn= x*exp(-x*10) 3.2.3 符号代数方程求解 例1. f = ax2+bx+c 求解 f=a*x^2+b*x+c; solve(f) —— 对缺省变量x求解 ans = [1/2/a*(-b+(b^2-4*a*c)^(1/2))] [1/2/a*(-b-(b^2-4*a*c)^(1/2))] 例2. 符号方程cos(x)=sin(x) 求解 f1=solv