优·第四讲 MATLAB的符号运算.pptVIP

第三讲 MATLAB的符号运算 —— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号代数方程求解 符号微积分 符号微分方程 一、符号运算的基本操作 什么是符号运算 与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算。 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 3.符号矩阵的创建 数值矩阵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.直接修改 可用?、 ?键找到所要修改的矩阵，直接修改 b.指令修改 用A1=subs(A, new, old)来修改 A1=subs(S, old ,new) 例如：A =[ a, 2*b] [3*a, 0] A(2,2)=4*b A3 = [ a, 2*b] [3*a, 4*b] ? 符号矩阵与数值矩阵的转换 将数值矩阵转化为符号矩阵 函数调用格式：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] 将符号矩阵转化为数值矩阵 函数调用格式： numeric(A) A = [ 1/3, 5/2] [10/7, 2/5] numeric(A) ans = 0.3333 2.5000 1.4286 0.4000 二、符号运算 基本运算符 例1：f= 2*x^2+3*x-5; g= x^2+x-7; syms x f g f=2*x^2+3*x-5; g= x^2+x-7; h=f+g h = 3*x^2+4*x-12 例2：f=cos(x);g= sin(2*x); syms x f=cos(x);g=sin(2*x); f/g+f*g ans = cos(x)/sin(x)+cos(x)*sin(x) 查找符号变量 查找符号表达式中的符号变量 findsym 举例 符号表达式的替换 用给定的数据替换符号表达式中的指定的符号变量 subs 举例 六类常见符号运算 因式分解 函数展开 合并同类项 函数简化 函数简化 函数简化举例 分式通分 horner 多项式 计算极限 计算导数 计算积分 int（被积表达式，积分变量，积分上限， 积分下限）—— 定积分 例1.计算二重不定积分 syms x y F=int(int(x*exp(-x*y),x),y) F = 1/y*exp(-x*y) syms x f=x*exp(-x*10); F=ztrans(f) F=ztrans(x*exp(-x*10); F = z*exp(-10)/(z-exp(-10))^2 例3. 计算指数函数eAt。 用拉氏反变换法计算eAt的公式为： eAt = L-1[(SI-A)-1] 系统矩阵A= a=[0 1;-2 -3]; syms s b=(s*eye(2)-a) b = [ s,