7符号计算matlab.pptVIP

七：MATLAB符号计算 MATLAB 课堂 content 前言 与数值运算的区别 7.1 符号计算基础 7.2符号函数的极限 7.3符号微积分 7.4 级数 7.5 代数方程的符号求解 7.6 常微分方程的符号求解 7.7符号分析可视化 7.8符号计算局限性和Maple调用 符号运算的特点： ? 运算对象可以是没赋值的符号变量 ? 可以获得任意精度的解 Symbolic Math Toolbox——符号运算工具包通过调用Maple软件实现符号计算的。 maple软件——主要功能是符号运算， 它占据符号软件的主导地位。　　 7.1 符号计算基础 例: a=sym(‘a’)建立符号量a， k2=sym(2); 建立单个符号量k2 用户可以在表达式中使用变量a和k2进行各种运算。 f = sin(x)+5x f —— 符号变量名 sin(x)+5x—— 符号表达式 —— 符号标识 数值矩阵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.指令修改 用R=subs(S, old, new)来修改 例如： ? 符号矩阵与数值矩阵的转换 将数值矩阵转化为符号矩阵 函数调用格式：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 在symbolic中有三种不同的算术运算： 数值类型 matlab的浮点算术运算 有理数类型 maple的精确符号运算 vpa类型 maple的任意精度算术 运算 浮点算术运算 1/2+1/3 －－(定义输出格式format long) ans = 0.83333333333333 符号运算 sym(1/2)+(1/3) ans = 5/6 －－精确解 任意精度算术运算 digits(n) —— 设置可变精度，缺省16位 vpa(x,n) —— 显示可变精度计算 digits(25) vpa(1/2+1/3) %在digits指定精度下，给出结果 ans = .8333333333333333333333333 vpa(5/6,40) %在40位精度下给出计算结果 ans = .8333333333333333333333333333333333333333 a=sym([1/4,exp(1);log(3),3/7]) a = [ 1/4,exp(1)] [log(3), 3/7] vpa(a,10) ans = [.2500000000, 2.718281828] [1.098612289, .4285714286] 7.2 符号函数的极限 limit函数的调用格式为： 7.3 符号微积分 7.3.1不定积分 在MATLAB中，求不定积分