第8章+符号运算.pptVIP

第八章 MATLAB的符号运算 —— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号代数方程求解 符号微积分 符号微分方程 一、符号运算的基本操作 1.什么是符号运算 数值运算中必须先对变量赋值 符号运算无须事先对独立变量赋值 2.符号对象的生成 f = sym(arg) f = sym(argn,flagn) syms(‘argv1’,’argv2’,’argvk’) syms argv1 argv2 argvk flagn: ‘d’ 最接近的十进制浮点精确表示 ‘r’ 最接近有理表示，缺省。 a1=[1/3,pi/7,sqrt(5),pi+sqrt(5)] a2=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) a3=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) a23=a2-a3 a1 =0.3333 0.4488 2.2361 5.3777 a2 =[1/3, pi/7, sqrt(5), 6054707603575008*2^(-50)] a3 =[1/3, pi/7, sqrt(5), pi+sqrt(5)] a23 =[0, 0, 0, 189209612611719/35184372088832-pi-5^(1/2)] 把字符表达式转换为符号变量 方法2 syms x 或者 syms(‘x’) f = sin(x)+5*x 方法1 f =sym(‘[sin(x)+5x]‘)或者 f =sym(‘sin(x)+5*x’) f —— 符号变量名 sin(x)+5x—— 符号表达式 的内容可以是符号表达式，也可以是符号方程。 例： f1=sym(a*x^2+b*x+c)二次三项式 f2=sym(a*x^2+b*x+c=0)方程 ※符号表达式或符号方程可以赋给符号变量，以后调用方便；也可以不赋给符号变量直接参与运算 例1 y=sym(2*sin(x)*cos(x)) y =2*sin(x)*cos(x) y1=simple(y) y1 =sin(2*x) 例2 syms fai1 fai2 y=simple(sin(fai1)*cos(fai2)-cos(fai1)*sin(fai2)) y =sin(fai1-fai2) ?用matlab函数sym创建矩阵 命令格式：A=sym([ ]) ※ 符号矩阵内容同数值矩阵 ※ 需用sym指令定义 ※ 需用标识 3.符号矩阵的创建 数值矩阵A=[1,2;3,4] A=[a,b;c,d] —— 不识别 注意：符号矩阵的每一行的两端都有方括号，这是与 matlab数值矩阵的一个重要区别。 例如： A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0] size(A) ans = 2 2 ② 根据数值数组的方法创建 syms a b; A=[a,2*b;3*a,0] A = [ a, 2*b] [ 3*a, 0] ③ 符号矩阵的修改 b.指令修改 用A1=subs(A,’new’,’old’)来修改 a.直接修改 可用?、 ?键找到所要修改的矩阵，直接修改 例如： A=sym([a,2*b;3*a,0]) A(2,2)=4*b A = [ a, 2*b] [3*a, 4*b] A1=subs(A, c, b) A1 =[ a, 2*c] [3*a, 4*c] 将数值矩阵转化为符号矩阵 函数调用格式：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] 将符号矩阵转化为数值矩阵 char to string do