实验四MATLAB符号运算.docx

《实验四MATLAB符号运算.docx》由会员分享，可在线阅读，更多相关《实验四MATLAB符号运算.docx（9页珍藏版）》请在三一办公上搜索。

1、实验四 MATLAB符号运算实验四 MATLAB符号运算 一、 实验目的： 1、 掌握定义符号对象的方法； 2、 掌握符号表达式的运算法则以及符号矩阵运算。 3、 掌握求符号函数极限及导数的方法。 4、 掌握求符号函数定积分和不定积分的方法。 二、 实验原理 1、 符号常量、符号变量、符号表达式的创建 (1) 使用 sym( )创建 输入以下命令，观察 Workspace 中 A、B、f是什么类型的数据，占用多少字节的内存空间。 A=sym(1) %符号常量 B=sym(x) %符号变量 f=sym(2*x2+3y-1) %符号表达式 clear f1=sym(1+2) %有单引号，表示字符串

2、 f2=sym(1+2) %无单引号 f3=sym(2*x+3) f4=sym(2*x+3) %为什么会出错 x=1 f4=sym(2*x+3) 通过看 MATLAB 的帮助可知，sym( )的参数可以是字符串或数值类型，无论是哪种类型都会生成符号类型数据。 (2) 使用 syms 创建 clear syms x y z %注意观察x,y,z都是什么类型的，它们的内容是什么 x,y,z f1=x2+2*x+1 f2=exp(y)+exp(z)2 f3=f1+f2 通过以上实验，知道生成符号表达式的第二种方法：由符号类型的变量经过运算(加减乘除等)得到。又如： f1=sym(x2+y +sin(

3、2) syms x y f2=x2+y+sin(2) x=sym(2) , y=sym(1) f3=x2+y+sin(2) y=sym(w) f4=x2+y+sin(2) 符号矩阵创建 syms a1 a2 a3 a4 A=a1 a2;a3 a4 A(1),A(3) 或者 B=sym( b1 b2 ;b3 b4 ) c1=sym(sin(x) ) c2=sym(x2) c3=sym(3*y+z) c4=sym(3 ) C=c1 c2; c3 c4 2、 符号算术运算 (1) 符号量相乘、相除 符号量相乘运算和数值量相乘一样，分成矩阵乘和数组乘。 a=sym(5);b=sym(7); c1=a*

4、b c2=a/b a=sym(5);B=sym(3 4 5); C1=a*B, C2=aB syms a b A=5 a;b 3; B=2*a b;2*b a; C1=A*B, C2=A.*B C3=AB, C4=A./B (2) 符号数值任意精度控制和运算 任意精度的 VPA 运算可以使用命令 digits(设定默认的精度)和 vpa(对指定对象以新的精度进行计算)来实现。 a=sym(2*sqrt(5)+pi) b=sym(2*sqrt(5)+pi) digits vpa(a) digits(15) vpa(a) c1=vpa(a,56) c2=vpa(b,56) 注意：观察 c1 和 c

5、2 的数据类型，c1 和 c2 是否相等。 3、 符号表达式的操作和转换 符号表达式化简主要包括表达式美化(pretty)、合并同类项(collect)、多项式展开(expand)、因式分解(factor)、化简(simple 或 simplify)等函数。 合并同类项(collect)。分别按x的同幂项和e指数同幂项合并表达式：syms x t; f=(x2+x*exp(-t)+1)*(x+exp(-t); f1=collect(f) f2=collect(f,exp(-t) 对显示格式加以美化(pretty)。针对上例，用格式美化函数可以使显示出的格式更符合数学书写习惯。 pretty(f

6、1) pretty(f2) 注意：与直接输出的 f1 和 f2 对比。 多项式展开(expand)。展开 (x-1)12 成 x 不同幂次的多项式。 clear all syms x; f=(x-1)12; pretty(expand(f) 因式分解(factor)。将表达式 x121作因式分解。 clear all syms x; f=x12-1; pretty(factor(f) 化简(simple 或 simplify)。 将函数 化简。 clear all, syms x; f=(1/x3+6/x2+12/x+8)(1/3); g1=simple(f) g2=simplify(f) 4

7、、 符号极限、符号积分与微分 (1) 求极限函数的调用格式 limit(F,x,a) %返回符号对象F当xa时的极限 limit(F,a) %返回符号对象F当独立变量*a时的极限 limit(F) %返回符号对象F当独立变量0(a=0)时的极限 limit(F,x,a,right) %返回符号对象F当xa时的右极限 limit(F,x,a,left) %返回符号对象F当xa时的左极限 (2) 求积分函数的调用格式 int(F) %求符号对象F关于默认变量的不定积分 int(F,v) %求符号对象F关于指定变量v的不定积分 int(F,a,b) %求符号对象F关于默认变量的从a到b的定积分 in

8、t(F,v,a,b) %求符号对象F关于指定变量v的从a到b的定积分 (3) 求微分函数的调用格式 diff(F) %求符号对象F关于默认变量的微分 diff(F,v) %求符号对象F关于指定变量v的微分 diff(F,n) %求符号对象F关于默认变量的n次微分，n为自然数1、2、3 diff(F, v,n) %求符号对象F关于指定变量v的n次微分 5、 符号方程的求解 (1) 常规方程求解函数的调用格式 g = solve(eq) %求方程(或表达式或字串)eq关于默认变量的解 g = solve(eq,var) %求方程(或表达式或字串)eq关于指定变量var的解 g = solve(eq

9、1,eq2,.,eqn,var1,var2,.,varn) %求方程(或表达式或字串)组 eq1,eq2,.,eqn关于指定变量组var1,var2,.,varn的解 (2) 常微分方程求解 求解常微分方程的函数是 dsolve。应用此函数可以求得常微分方程(组)的通解，以及给定边界条件(或初始条件)后的特解。 常微分方程求解函数的调用格式： r = dsolve(eq1,eq2,.,cond1,cond2,., v) r = dsolve(eq1,eq2,.,cond1,cond2,.,v) 说明： 以上两式均可给出方程 eq1、eq2 .对应初始条件 cond1、cond2 .之下的以 v

10、作为 解变量的各微分方程的解。 常微分方程解的默认变量为 t。 第二式中最多可接受的输入式是 12 个。 微分方程的表达方法。 2 用大写字母Dy表示微分符号dy，在用MATLAB求解常微分方程时，用D2y表示dy，dxdx2依次类推。 边界条件以类似于 y(a) = b 或 Dy(a) = b 的等式给出。其中 y为因变量，a、b 为常数。如果初始条件给得不够，求出的解则为含有 C1、C2 等待定常数的通解。 例如：求微分方程为y=2x 的通解。y=dsolve(Dy=2*x,x) 四、实验要求 1. 已知x=6,y=5,利用符号表达式求z=x+13+x-yZ=9.1631 syms x y

11、 x=6 y=5 z=(x+1)/(sqrt(3+x)-sqrt(y) 截图： 2. 分解因式。 x4-y4； syms x y z=x4-y4 pretty(factor(z) 截图： 125x+75x+15x+1 syms x y=125*x6+75*x4+15*x2+1 pretty(factor(y) 截图： 6423. 化简表达式 sinb1cosb2-cosb1sinb2； syms beta1 beta2 f=sin(beta1)*cos(beta2)-cos(beta1)*sin(beta2) r=simple(f) 截图： 4x2+8x+3 2x+1syms x y=(4*x

12、2+8*x+3)/(2*x+1) r=simple(y) 截图： 4、求syms x f=(x2-1)/(x2-3*x+2) w=limit(f,x,2) 截图： 5、求函数syms x f=cos(2*x)-sin(2*x) int(f,x) 截图： 的积分。 计算定积分syms x f=sin(x)+2 int(f,x,0,pi/6) 截图： 。 6、求函数1-cos(2x)，求y和y。 y=xsyms x y=(1-cos(2*x)/x diff(y,1) 截图： syms x y=(1-cos(2*x)/x diff(y,2) 截图： 7、求下列线性方程组的解syms x y z x y z=solve(x+y+z=10,3*x+2*y+z=14,2*x+3*y-z=1) 截图： 8、求解当y(0)=2，z(0)=7时，微分方程组的解 实验报告提交格式： 1、 实验题目 2、 实验目的 3、 实验内容