MATLAB原理应用实验报告第三章标记运算.doc

《MATLAB原理及应用》实验报告 第三章 MATLAB的符号运算 一．实验目的 1、掌握符号对象的命名方法 2、掌握符号表达式的基本运算 3、掌握符号级数的求法 二．实验设备 计算机、MATLAB软件 三．实验内容 1.确定符号表达式的变量 为了简化符号对象的操作和计算，MATLAB为用户提过了findsym命令。 r=findsym(S)确定符号表达式或者矩阵S中自由符号变量 r=findsym(S，n)确定符号表达式或者矩阵S中靠近x最近的n个独立符号变量。 【实验3-1】使用MATLAB的命令确定符号表达式的变量。 在MATLAB的命令窗口中输入下例内容： syms a x y z t 确定下面简单符号表达式中的符号变量信息： findsym(sin(pi*t)) ans = t 确定下面简单符号表达式中的符号变量信息： findsym(x+i*y-j*z) ans = x, y, z 确定下面简单符号表达式中的符号变量信息： findsym(a+y,1) ans = y 符号表达式元算 符号表达式的四则运算 表达式的四则运算与数字运算一样，用+、-、/、运算符实现，其运算结果依然是一个符号表达式。 【实验3-2】 在MATLAB的命令窗口中输入下例内容： f=sym(2*x^2+3*x-5);%定义符号表达式 g=sym(x^2-x+7); f+g ans = 3*x^2+2*x+2 ans = 3*x^2+2*x+2 f^g ans = (2*x^2+3*x-5)^(x^2-x+7) 3.符号表达式的提取分子和分母运算 如果符号表达式是一个有理分式或可以展开为有理分式，可以可利用numden函数来提取符号表达式的分子或分母。期一般调用格式为 [n,d]=numden函数来提取符号表达式 该函数提取的符号表达式s的分子和分母，分别将它们存放在n和d中。 【实验3-3】 在MATLAB的命令窗口中输入下例内容： f= sym(a*x/(b+x)); [n,d]=numden(f) n = a*x d = b+x numden函数在提取各部分之前，将符号表达式有利化后返回所得分子和分母 g=sym((x^2+3)/(2*x-1)+3*x/(x+1)); [n,d]=numden(g) n = x^3+7*x^2+3 d = (2*x-1)*(x+1) 如果符号表达式是一个符号矩阵，numden返回两个新矩阵n和d，其中n是分子矩阵，d是分母矩阵。 h=sym([3/2,(2*x+1)/3;a/x+a/y,x+4]) h = [ 3/2, (2*x+1)/3] [ a/x+a/y, x+4] [n,d]=numden(h) n = [ 3, 2*x+1] [ a*(y+x), x+4] d = [ 2, 3] [ x*y, 1] 4.符号表达式的因式分解与展开 MATLAB提供了符号表达式分解与展开的函数，函数的调用格式为 ①factor(s) :对符号表达式s分解因式。 ②expand(s):对符号表达式s进行展开。 ③collect(s):对符号表达式s进行合并同类型。 ④collect(S，v) 将表达式S中相同次幂的合并，v的默认值是x 【实验3-4】 syms x y 下面简单符号表达式s1因式分解 s1=x^3-y^3; factor(s1) ans = (x-y)*(x^2+x*y+y^2) 下面简单符号表达式s2进行展开 s2=(-7*x^2-8*y^2)*(-x^2+3*y^2); expand(s2) ans = 7*x^4-13*x^2*y^2-24*y^4 下面简单符号表达式s3按变量y合并同类型 s3=(x+y)*(x^2+y^2+1); collect(s3,y) ans = y^3+x*y^2+(x^2+1)*y+x*(x^2+1) 下面简单符号整数分解因式 factor(sym(630)) ans = (2)*(3)^2*(