MATLAB与数学建模2014_第6章符号计算.ppt

第6章 MATLAB符号计算;6.1 符号计算基础;6.1 符号计算基础;例6.1 考察符号变量和数值变量的差别。 在 MATLAB命令窗口，输入命令： a=sym(a); b=sym(b); c=sym(c); d=sym(d); %定义4个符号变量 w=10; x=5; y=-8; z=11; %定义4个数值变量 A=[a, b; c, d] %建立符号矩阵A B=[w, x; y, z] %建立数值矩阵B det(A) %计算符号矩阵A的行列式 det(B) %计算数值矩阵B的行列式 A = [ a, b] [ c, d] B = 10 5 -8 11;例6.2 比较符号常数与数值在代数运算时的差别。 pi1=sym(pi); k1=sym(8); k2=sym(2); k3=sym(3); % 定义符号变量 pi2=pi; r1=8; r2=2; r3=3; % 定义数值变量 sin(pi1/3) % 计算符号表达式值 sin(pi2/3) % 计算数值表达式值 sqrt(k1) % 计算符号表达式值 sqrt(r1) % 计算数值表达式值 sqrt(k3+sqrt(k2)) % 计算符号表达式值 sqrt(r3+sqrt(r2)) % 计算数值表达式值;6.1 符号计算基础;clear; syms a b c d; A1=[a,b;c,d]; A2=det(A1); A3=subs(A2,a,1); A3=subs(A3,b,2); A3=subs(A3,c,3); A3=subs(A3,d,4) ;6.1 符号计算基础;例6.5 计算3阶范得蒙矩阵行列式的值。设A是一个由符号变量a,b,c确定的范得蒙矩阵。 syms a b c U=[a,b,c] A=[[1,1,1];U;U.^2] %建立范得蒙符号矩阵 det(A) %计算A的行列式值 ;6.1 符号计算基础;6.1.2 基本的符号运算 1. 符号表达式运算 (1) 符号表达式的四则运算 例6.7 符号表达式的四则运算示例。 在 MATLAB命令窗口，输入命令： syms x y z; f=2*x+x^2*x-5*x+x^3 %符号表达式的结果为最简形式 f=2*x/(5*x) %符号表达式的结果为最简形式 f=(x+y)*(x-y) %符号表达式的结果不是x^2-y^2，而是 %(x+y)*(x-y);6.1 符号计算基础;例6.8 对符号矩阵A的每个元素分解因式。 命令如下： syms a b x y; A=[2*a^2*b^3*x^2-4*a*b^4*x^3+10*a*b^6*x^4, 3*x*y-5*x^2; 4, a^3-b^3]; factor(A) %对A的每个元素分解因式;例6.9 计算表达式S的值。 syms x y; s=(-7*x^2-8*y^2)*(-x^2+3*y^2); expand(s) %对s展开 collect(s,x) %对s按变量x合并同类项(无同类项) factor(ans) % 对ans分解因式 expand(cos(x+y));6.1 符号计算基础;6.1 符号计算基础;6.1.3 符号表达式中变量的确定 MATLAB中的符号可以表示符号变量和符号常数。findsym可以帮助用户查找一个符号表达式中的符号变量。该函数的调用格式为：findsym(S,n) 函数返回符号表达式S中的n个符号变量，若没有指定