第三章 matlab语言的符号运算 (2).pptVIP

例题：计算拉普拉斯反变换 例题：计算Z变换 计算反Z变换。 * * 第三章 MATLAB语言的符号运算 —— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号代数方程求解 符号微积分 符号微分方程 3.1 符号运算的基本操作 什么是符号运算 与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算。 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 特点： ? 运算对象可以是没赋值的符号变量 可以获得任意精度的解 　　 2. 符号变量与符号表达式 f = sin(x)+5x f —— 符号变量名 sin(x)+5x—— 符号表达式 —— 符号标识 符号表达式一定要用 单引 号括起来matlab才能识别。 的内容可以是符号表达式，也可以是符号方程。 例： f1=a?x^2+b?x+c —— 二次三项式 f2= a?x^2+b?x+c=0 —— 方程 f3=Dy+y^2=1 ——微分方程 ※符号表达式或符号方程可以赋给符号变量，以后调用方便；也可以不赋给符号变量直接参与运算 3.符号变量定义 R = SYM(‘(-1 + sqrt(5))/2’) R = (-1 + sqrt(5))/2 vpa(R,4) ans = .6180 syms w t F = sin(w*t); syms s G=(s+1)/((s+2)(s+3)) G = (s+1)/(s+2)/(s+3) E = sym(‘Dy + 5 * y = t’); A = sym(‘[a11 a12;a21 a22]’); 例题：f = sin(w*t + pi/3), 求该函数对时间t的导数，并在x=0邻域展开泰勒级数的前4项。 f = sym(‘sin(w*t + pi/3)’); f1 = diff(f,’t’) f2 = taylor(f,’t’,4) 例题：已知函数y = sin(wt-pi/6),求其对时间t的导数，并展开泰勒级数的前4项。 4.符号矩阵的创建 数值矩阵A=[1,2;3,4] A=[a,b;c,d] —— 不识别 ?用matlab函数sym创建矩阵（symbolic 的缩写) 命令格式：A=sym([ ]) ※ 符号矩阵内容同数值矩阵 ※ 需用sym指令定义 ※ 需用 标识 例如：A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0] 这就完成了一个符号矩阵的创建。 注意：符号矩阵的每一行的两端都有方 括号，这是与 matlab数值矩阵的 一个重要区别。 ?用字符串直接创建矩阵 模仿matlab数值矩阵的创建方法 需保证同一列中各元素字符串有相 同的长度。 例：A =[[ a,2*b]; [3*a, 0]] A = [ a, 2*b] [3*a, 0] ? 符号矩阵的修改 a.直接修改 可用?、 ?键找到所要修改的矩阵，直接修改 b.指令修改 用A1=subs(A, ‘new’, ‘old’)来修改 例如：A =[ a, 2*b] [3*a, 0] A1=subs(A, c, b) A1 =[ a, 2*c] [3*a, 0] 函数调用格式：sym(A) A=