2009级MATLAB实验3.docVIP

实验3 MATLAB的符号使用 目的和要求 1.熟练掌握MATLAB符号表达式的创建。 2.掌握符号和数值的替换。 3.熟练掌握符号矩阵的创立及运算。 4.熟练掌握符号微积分。 5.掌握符号积分变换。 6.熟悉符号方程的求解。 内容和步骤 符号运算允许在运算对象和运算过程中出现非数值的符号对象，利用符号对象进行运算。在工程时间和科学研究等各个方面，经常会遇到数值运算无法进行描述的问题，存在非数值问题，引入符号运算就可以解决这方面的问题。 创建符号对象 MATLAB提供了两个建立符号对象的函数：sym和syms，两个函数的用法不同。 sym函数用来建立单个符号量，一般调用格式为： 符号量名=sym(符号字符串) 该函数可以建立一个符号量，符号字符串可以是常量、变量、函数或表达式。应用sym函数还可以定义符号常量，使用符号常量进行代数运算时和数值常量进行的运算不同。 函数sym一次只能定义一个符号变量，使用不方便。MATLAB提供了另一个函数syms，一次可以定义多个符号变量。syms函数的一般调用格式为： syms 符号变量名1 符号变量名2 … 符号变量名n 用这种格式定义符号变量时不要在变量名上加字符串分界符(‘)，变量间用空格而不要用逗号分隔。 2．创建符号表达式和符号表达式的操作 含有符号对象的表达式称为符号表达式。建立符号表达式有以下3种方法： (1)利用单引号来生成符号表达式。 (2)用sym函数建立符号表达式。 (3) 使用已经定义的符号变量组成符号表达式。 对符号表达式f=sinx,g= 进行操作。 (1).创建符号变量。 ①使用sym命令创建符号表达式： f=sym(‘sin(x)’) f= sin(x) g=sym(‘y/exp(-2*t)’) g= y/exp(-2*t) ②syms命令创建符号表达式f,g: syms x y t f=sin(x) f= sin(x) g=y/exp(-2*t) g= y/exp(-2*t) ⑵自由变量的确定。使用findsym确定符号表达式g的自由变量： findsym(g) %得出所有符号变量 Ans= t,y fingsym(g,1) %得出第1个符号变量 ans= y findsym(g,2) %按顺序得出两个符号变量 Ans= y,t ⑶用常数替换符号变量。用行向量替换x,使符号对象f转变为行向量： x=0:10; f=’sin(x)’; y=subs(f,x) y= columns1 through 6 0 0.8415 0.9093 0.1411 -0.7568 -0.9598 Columns 7 through 11 -0.2794 0.6570 0.9894 0.4123 0.4121 -0.5440 x,y都为双精度型数值。 ⑷符号对象与数值的转换.采用double,numeric 和eval将符号对象转换为数值： f=’sin(x)’; f1=subs(f,’5’) f1= Sin((5)) y1=double(f1) y1= -0.9589 y2=eval(f1) y2= -0.9589 ⑸求反函数和复合函数。 ①用finverse函数求f,g的反函数： finverse(f) ans= asin(x) finverse(g) %对默认独立变量y求反函数 ans= y*exp(-2*t) finverse(g,’t’) %对符号变量t求反函数 ans= -1/2*log(y/t) ②用compose 函数求f,g的复合函数： compose(f,g) %计算f(g(x)) ans= sin(y/exp(-2*t)) compose(f,g,’z’) %计算f(g(z)) ans= sin(z/exp(-2*t)) ⑹符号微积分和极限。 ①对f和g用diff 求微分： diff(f) ans= cos(x) diff(g) %对默认自由变量y求微分 ans= 1/exp(t^2) diff(g,’t’) %对符号变量t求微分 Ans= -2*y/exp(t^2)*t ②对f和g用int 求积分：： int(t)