MATLAB第4章 符号运算要点.pptVIP

第4章 MATLAB的符号运算 4.1 符号运算的基本操作 §4.1 符号运算的基本操作 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阶 bernoulli 多项式，计算 x=3 时bernoulli 数。 a=maple(bernoulli(5,x)) a = -1/6*x+5/3*x^3+x^5-5/2*x^4 a=maple(bernoulli(5,3)) a = 85 例：寻找二次多项式的完全平方 f (x) = x2+2x+2 a=maple(completesquare(x^2+2*x+2)) a = completesquare(x^2+2*x+2) maple(with(student);) a=maple(completesquare(x^2+2*x+2)) a = (x+1)^2+1 例：求sin(x2+y2)在x=0,y=0处泰勒级数展开式，8阶截断。 maple(mtaylor(sin(x^2+y^2),[x=0,y=0],8)) ans = mtaylor(sin(x^2+y^2),[x = 0, y = 0],8) maple(readlib(mtaylor);) maple(mtaylor(sin(x^2+y^2),[x=0,y=0],8)) ans = x^2+y^2-1/6*x^6-1/2*y^2*x^4-1/2*y^4*x^2-1/6*y^6 * * * 4.2 符号运算 4.3 maple函数—符号运算的扩展 MATLAB不仅具有数值运算功能，还开发了在MATLAB环境下实现符号计算的工具包： Symbolic Math Toolbox 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号代数方程求解 符号微积分 符号微分方程 一、什么是符号运算 数值运算中必须先对变量赋值，然后才能参与运算。符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 1、概念 2、特点 ? 运算对象可以是没赋值的符号变量； ? 可以获得任意精度的解。 ? Symbolic Math Toolbox——符号运算工具包通过调用Maple软件实现符号计算的。 ? maple软件——主要功能是符号运算， 它占据符号软件的主导地位。 3、途径 二、 符号变量与符号表达式 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 ——微分方程 符号表达式或符号方程可以赋给符号变量，以后调用方便；也可以不赋给符号变量直接参与运算 三、符号矩阵的创建 数值矩阵 A=[1,2;3,4] A=[a,b;c,d] — 不识别 1、用MATLAB函数sym创建矩阵（symbolic的缩写) 命令格式：A=sym([ ]) ※ 符号矩阵内容同数值矩阵 ※ 需用sym指令定义 ※ 需用 标识 例如：A = sym ([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0] 这就完成了一个符号矩阵的创建。 注意：符号矩阵的每一行的两端都有方括号，这是与 MATLAB数值矩阵的一个重要区别。 2、用字符串直接创建矩阵 模仿matlab数值矩阵的创建方法 需保证同一列中各元素字符串有相同的长度。 例：A =[[a,2*b]; [3*a, 0]] A = [ a, 2*b] [3*a, 0] 3、符号矩阵的修改 a.直接修改 可用?、 ?键找到所要修改的矩阵，直接修改 b.指令修改 用A1=subs(S, old ,