MATLAB语言及应用---第7章符号计算.pptVIP

第七章 符号计算 主要内容 7.1 符号对象和符号表达式 1、符号表达式的生成 7.1 符号对象和符号表达式 例题开讲 例题开讲 7.1 符号对象和符号表达式 2、符号计算中的算符和基本函数 7.1 符号对象和符号表达式 3、识别对象类别的指令 7.1 符号对象和符号表达式 4、符号表达式中自由变量的确定 例题开讲 例题开讲 7.2 符号微积分 1、符号极限 7.2 符号微积分 2、符号积分 7.2 符号微积分 3、符号微分 测验 1）syms x; f=sin(x^2) , b=diff(f); 则 b= ________ 2）A=sym(‘[cos(x),sin(x);x^2+x+1 tan(x)]’) B=det(A);C=diff(A(3)); B=______________ C=______________ * *杨惠--《matlab语言及应用》--第五章 1、符号对象和符号表达式 2、符号微积分 规定： （1）在定义符号运算时，首先要定义基本的符号对象，然后 用这些基本符号对象去构成新的表达式，进而从事所需 的符号运算。 （2）在运算中，凡是包含符号对象的表达式所生成的衍生对 象也都是符号对象。 生成方式： 指令 sym、syms f=sym(arg) f=sym(argn,flagn) argv=sym(‘argv’,flagv) syms(‘argv1’,’argv2’,’argvk’) syms argv1 argv2 argvk 把数字、字符串或表达式arg转换成符号对象f 把字符argv1 argv2 argvk定义为基本符号对象 把数值或数值表达式argn转换成flagn格式的符号对象 按flagv指定的要求把字符串’argv’定义为符号对象 例一： 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)]) 例二： syms x; y=sym(2*sin(x)*cos(x)) y=simple(y) 按规则把已有的y符号表达式化为最简形式 例：求以下矩阵的行列式值、逆和特征根。 syms a11 a12 a21 a22; A=[a11,a12;a21,a22] DA=det(A) IA=inv(A) EA=eig(A) 是空格！ 不能用逗号 说明： 符号计算表达式的算符和基本函数，无论在形状、名称还是 使用方法上，都与数值计算中的算符和基本函数几乎完全相同。 异同： 1）在关系运算符中，只有是否等于的概念。“==”和“~=” 2）函数atan2仅能用于数值计算 3）符号计算中没有log2,log10计算，只有自然对数的运算。 4）符号计算中，没有提供求相角的指令。 5）svd（奇异值分解）指令在符号运算中与数值计算不同。 %（1）生成三种不同类型的矩阵，给出不同的显示形式 a=1;b=2;c=3;d=4; Mn=[a,b;c,d] Mc=[a,b;c,d] Ms=sym(Mc) %（2）三种矩阵的大小不同 SizeMn=size(Mn) SizeMc=size(Mc) SizeMs=size(Ms) %（3）用class获得每种矩阵的类别 CMn=class(Mn) CMc=class(Mc) CMs=class(Ms) %（4）用isa判断每种矩阵的类别 isa(Mn,double) isa(Mc,char) isa(Ms,sym) %（5）利用whos观察内存变量的类别和其它属性 whos Mn Mc Ms findsym(expr) findsym(expr,N) 确认表达式expr中所有的自由符号变量 从表达式expr中确认出靠x最近的N个独立自变量 说明：如果表达式是矩阵形式， 指令是对整个矩阵进行的 %（1）生成符号变量 syms a b x X Y k=sym(3); z=sym(c*sqrt(delta)+y*sin(theta)); EXPR=a*z*X+(b*x^2+k)*Y; %（2）找出EXPR中的全部自由符号变量 findsym(EXPR) %（3）在EXPR中确定一个自由符号变量 findsym(EXPR,1) %（4）在EXP