第三章MATLAB符号运算功能36new.pptVIP

第三章 符号运算功能 —— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号代数方程求解 符号微积分 符号微分方程 一、符号运算的基本操作 什么是符号运算 与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算。 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 2. 符号表达式与方程 （1）单引号法 f = sin(x)+5*x f —— 符号变量名 sin(x)+5x—— 符号表达式 —— 符号标识 符号表达式一定要用 单引 号括起来matlab才能识别。 3.符号矩阵的创建 （1）函数sym创建矩阵（symbolic 的缩写) 命令格式：A=sym([ ]) ※ 符号矩阵内容同数值矩阵 （2）用字符串直接创建矩阵 模仿matlab数值矩阵的创建方法 需保证同一列中各元素字符串有相 同的长度。 ? 符号矩阵的修改 a.直接修改 可用?、 ?键找到所要修改的矩阵，直接修改 b.指令修改 用A1=subs(S, ‘old’ ,‘new’)来修改 注：多个参量需替代时，应放入{}中 例如：A=sym([a,2*b;3*a,0]) A =[ a, 2*b] [ 3*a, 0] A1=subs(A,0,’4*b’) （4） 符号矩阵与数值矩阵的转换 将数值矩阵转化为符号矩阵 函数调用格式：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 二、符号运算 符号矩阵运算 与数值运算中的运算符和基本函数几乎完全相同。 例：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) 2. 符号微积分与积分变换 symsum(f,x,a,b) — 求表达式f在自变量x取[a,b]范围的级数和 diff(f) — 对缺省变量求微分 diff(f,v) — 对指定变量v求微分 diff(f,v,n) —对指定变量v求n阶微分 int(f) — 对f表达式的缺省变量求积分 int(f,v) — 对f表达式的v变量求积分 int(f,v,a,b) — 对f表达式的v变量在（a,b) 区间求定积分 5.符号代数方程求解 例： 解方程组 x+y+z=1 x-y+z=2 2x-y-z=1 f=solve(x+y+z=1,x-y+z=2,2*x-y-z=1) f = x: [1x1 sym] f.x ans =2/3 y: [1x1 sym] f.y ans =-1/2 z: [1x1 sym] f.z ans =5/6 5. 符号微分方程求解 一阶微分方程 [x,y]=dsolve(Dx=y,Dy=x,x(0)=0,y(0)=1, t) 结果：x= sin(t) y= cos(t) * 的内容可以是符号表达式，也可以是符号方程。 例： f1=a?x^2+b?x+c —— 二次三项式 f2= a?x^2+b?x+c=0 —— 方程 f3=Dy+y^2=1 ——微分方程 ※符号表达式或符号方程可以赋给符号变量，以后调用方便；也可以不赋给符号变量直