matlab2补充.pptVIP

MATLAB的符号运算 —— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号代数方程求解 符号微积分 符号微分方程 一、符号运算的基本操作 什么是符号运算 与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算。 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 特点： ? 运算对象可以是没赋值的符号变量 ? 可以获得任意精度的解　　 2. 符号变量与符号表达式 f =sym( sin(x)+5x‘) f —— 符号变量名 sin(x)+5x—— 符号表达式 —— 符号标识 符号表达式一定要用 单引 号括起来matlab才能识别。 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数值矩阵的 一个重要区别。 ? 符号矩阵的修改 a.直接修改 可用?、 ?键找到所要修改的矩阵，直接修改 b.指令修改 用A1=subs(A, old, new)来修改 A1=subs(S) ? 符号矩阵与数值矩阵的转换 将数值矩阵转化为符号矩阵 函数调用格式：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] 求极限 Limit（f，x，a） Limit（f，x，a，‘left’）左趋近a Limit（f，x，a,’right’）右趋近a 3. 符号微积分与积分变换 diff(f) — 对缺省变量求微分 diff(f,v) — 对指定变量v求微分 diff(f,n) — 对缺省变量求n阶微分 diff(f,v,n) —对指定变量v求n阶微分 int(f) — 对f表达式的缺省变量求积分 int(f,v) — 对f表达式的v变量求积分 int(f,v,a,b) — 对f表达式的v变量在（a,b) 区间求定积分，a，b也可以是符号 符号序列的求和 S＝symsum（f，v，a，b）求通式在指定的v变量取遍【a，b】中所有整数的和。 符号变换 例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 4.符号代数方程求解 例1. f = ax2+bx+c 求解 f=sym(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)