2013 MATLAB实验五.docVIP

实验五 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 findsym(g,1) %得出第1个符号变量 ans= y findsym(g,2) %按顺序得出两个符号变量 Ans= y,t ⑶用常数替换符号变量。用行向量替换x,使符号对象f转变为行向量： x=0:10; f=sin(x); y=subs(f,x) ⑷符号对象与数值的转换.采用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(f) %求不定积分 Ans= -cos(x) int(g) Ans= 1/2*y^2/exp(-2*t) int(g, t) %对t求不定积分 Ans= 1/2*y/exp(-2*t) int(g, t,0,10) %对t 求定积分 Ans= 1/2*y*exp(20)-1/2*y