MATLAB第7讲符号运算ppt课件.ppt

第七讲 MATLAB的符号运算—— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包SymbolicMath Toolbox 符号运算的功能ß符号表达式、符号矩阵的创立ß符号线性代数ß因式分解、展开和简化ß符号代数方程求解ß符号微积分ß符号微分方程一、符号运算的根本操作1.什么是符号运算ß与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达ß特点： 运算对象可以是没赋值的符号变量 可以获得任意精度的解ßSymbolic Math Toolbox——符号运算工具包通过调用Maple软件实现符号计算的ßmaple软件——主要功能是符号运算，占据符号软件的主导地位　　 2. 符号变量与符号表达式f = 'sin〔x〕+5x'f —— 符号变量名sin〔x〕+5x—— 符号表达式' '—— 符号标识v符号表达式一定要用' ' 单引 号括起来matlab才能识别。

' ' 的内容可以是符号表达式，也可以是符号方程例： f1='ax^2+bx+c' —— 二次三项式 f2= 'ax^2+bx+c=0' —— 方程 f3='Dy+y^2=1' ——微分方程※符号表达式或符号方程可以赋给符号变量，以后调用方便；也可以不赋给符号变量直接参与运算3.符号矩阵的创立 数值矩阵A=[1,2;3,4] A=[a,b;c,d] —— 不识别用matlab函数sym创立矩阵〔symbolic 的缩写〕命令格式：A=sym〔'[ ]'〕 ※ 符号矩阵内容同数值矩阵 ※ 需用sym指令定义 ※ 需用' '标识例如：A = sym〔'[a , 2*b ; 3*a , 0]'〕 A = [ a, 2*b] [3*a, 0] 这就完成了一个符号矩阵的创立注意：符号矩阵的每一行的两端都有方 括号，这是与 matlab数值矩阵的 一个重要区别。

用字符串直接创立矩阵v 模拟matlab数值矩阵的创立方法v 需保证同一列中各元素字符串有相 同的长度例：A =['[ a,2*b]'; '[3*a, 0]'] A = [ a, 2*b] [3*a, 0] 符号矩阵的修改 a.直接修改 可用、 键找到所要修改的矩阵，直接修改 b.指令修改v用A1=subs〔A, {old}, [new]〕来修改subs〔A, 'a','3*a'〕 A=sym〔'[a b;c d]'〕二、符号运算1.符号矩阵运算 数值运算中，所有矩阵运算操作指令都比较直观、简单例如：a=b+c; a=a*b ；A=2*a^2+3*a-5等 所有涉及符号运算的操作都需要先定义符号然后进展运算.例1：f= 2*x^2+3*x-5; g= x^2+x-7;>> syms x>> f=2*x^2+3*x-5; g= x^2+x-7;>> h=f+gh = 3*x^2+4*x-12例2：f=cos〔x〕;g= sin〔x〕;>> syms x>> f=cos〔x〕;g=sin〔x〕;>> f/g+f*gans =cos〔x〕/sin〔x〕+cos〔x〕*sin〔x〕 2. 任意精度的数学运算￿在symbolic中有三种不同的算术运算：1.数值类型 matlab的浮点算术运算2.有理数类型 maple的准确符号运算3.vpa类型 maple的任意精度算术 运算 ß浮点算术运算1/2+1/3 －－〔定义输出格式format long〕ans =0.83333333333333ß符号运算sym〔1/2〕+〔1/3〕ans = 5/6 －－准确解 ß任意精度算术运算digits〔n〕 —— 设置可变精度，缺省16位vpa〔x,n〕 —— 显示可变精度计算digits〔25〕vpa〔1/2+1/3〕ans =.8333333333333333333333333vpa〔5/6,40〕 ans =.8333333333333333333333333333333333333333 a=sym〔〔'[1/4,exp〔〔1〕〕;log〔〔3〕〕,3/7]'〕〕a =[ 1/4,exp〔〔1〕〕][log〔〔3〕〕, 3/7]vpa〔〔a,10〕〕ans =[.2500000000, 2.718281828][1.098612289, .4285714286]3. 符号微积分与积分变换ßdiff〔f〕 — 对缺省变量求微分ßdiff〔f,v〕 — 对指定变量v求微分ßdiff〔f,v,n〕 —对指定变量v求n阶微分ßint〔f〕 — 对f表达式的缺省变量求积分ßint〔f,v〕 — 对f表达式的v变量求积分ßint〔f,v,a,b〕 — 对f表达式的v变量在〔a,b〕 区间求定积分int〔'被积表达式'，'积分变量'，'积分上限'，￿￿￿￿￿￿￿￿￿'积分下限'〕—— 定积分——缺省时为不定积分taylor〔f,n〕 —— 泰勒级数展开ztrans〔f〕 —— Z变换iztrans〔f〕 —— 反Z变换Laplace〔f〕 —— 拉氏变换ilaplace〔f〕 —— 反拉氏变换fourier〔f〕 —— 付氏变换ifourier〔f〕 —— 反付氏变换例1.计算二重不定积分 F=int〔〔int〔〔'x*exp〔〔-x*y〕〕','x'〕〕,'y'〕〕F= 1/y*exp〔〔-x*y〕〕例2.计算 f='x*exp〔-x*10〕'的Z变换 F=ztrans〔〔f〕〕 F= z*exp〔〔-10〕〕/〔〔z-exp〔〔-10〕〕〕〕^2>> syms x y>> F=int〔int〔x*exp〔-x*y〕,x〕,y〕F =1/y*exp〔-x*y〕>> syms x>> f=x*exp〔-x*10〕;>> F=ztrans〔f〕>> F=ztrans〔x*exp〔-x*10〕;F =z*exp〔-10〕/〔z-exp〔-10〕〕^2 例3. 计算指数函数eAt。

用拉氏反变换法计算eAt的公式为： eAt = L-1[〔SI-A〕-1]系统矩阵A= eAt =结果：>> a=[0 1;-2 -3];>> syms s>> b=〔s*eye〔2〕-a〕b =[ s, -1][ 2, s+3]>> B=inv〔b〕 [ 〔s+3〕/〔s^2+3*s+2〕, 1/〔s^2+3*s+2〕] [ -2/〔s^2+3*s+2〕, s/〔s^2+3*s+2〕] >> b11=ilaplace〔sym〔b,1,1〕〕;b〔1,1〕=b11;>> b12=ilaplace〔sym〔b,1,2〕〕;b〔1,2〕=b12;>> b21=ilaplace〔sym〔b,2,1〕〕;b〔2,1〕=b21;>> b22=ilaplace〔sym〔b,2,2〕〕;b〔2,2〕=b22;>> bb = [ -exp〔-2*t〕+2*exp〔-t〕, exp〔-t〕-exp〔-2*t〕] [ -2*exp〔-t〕+2*exp〔-2*t〕, 2*exp〔-2*t〕-exp〔-t〕] 4.符号代数方程求解 matlab符号运算可以解一般的线性方程、非线性方程及一般的代数方程、代数方程组。

当方程组不存在符号解时，又无其他自由参数，那么给出数值解命令格式：solve〔f〕 —— 求一个方程的解Solve〔f1,f2, …fn〕 —— 求n个方程的解 例1. f = ax2+bx+c 求解f='a*x^2+b*x+c';ß solve〔f〕 —— 对缺省变量x求解ans =[1/2/a*〔-b+〔b^2-4*a*c〕^〔1/2〕〕][1/2/a*〔-b-〔b^2-4*a*c〕^〔1/2〕〕]计算机格式一般格式例2. 符号方程cos〔x〕=sin〔x〕 tan〔2*x〕=sin〔x〕求解f1=solve〔'cos〔x〕=sin〔x〕'〕,f1 =1/4*pi• solve〔f , 'b' 〕 —— 对指定变量b求解ans =-〔a*x^2+c〕/xf2=solve〔'tan〔2*x〕=sin〔x〕'〕f2 = atan〔1/2*〔-2*3^〔1/2〕〕^〔1/2〕,1/2+1/2*3^〔1/2〕〕 atan〔-1/2*〔-2*3^〔1/2〕〕^〔1/2〕,1/2+1/2*3^〔1/2〕〕 atan〔1/2*2^〔1/2〕*3^〔1/4〕/〔1/2-1/2*3^〔1/2〕〕〕+pi -atan〔1/2*2^〔1/2〕*3^〔1/4〕/〔1/2-1/2*3^〔1/2〕〕〕-pi 0 pi例3. 解方程组 x+y+z=1 x-y+z=2 2x-y-z=1g1='x+y+z=1',g2='x-y+z=2',g3='2*x-y-z=1'f=solve〔g1,g2,g3〕 [x,y,z]=solve〔'x+y+z=1','x-y+z=2','2*x-y-z=1'〕f=solve〔'x+y+z=1','x-y+z=2','2*x-y-z=1'〕f = x: [1x1 sym] f.x ans =2/3 y: [1x1 sym] f.y ans =-1/2 z: [1x1 sym] f.z ans =5/6 [x,y,z]=solve〔'x+y+z=1','x-y+z=2','2*x-y-z=1'〕 x = 2/3 y =-1/2 z =5/65. 符号微分方程求解 —— 用一个函数可以方便地得到微 分方程的符号解符号微分方程求解指令：dsolve命令格式：dsolve〔f,g〕ßf —— 微分方程，可多至12个微分方程的求 解；g为初始条件ß默认自变量为 'x',可任意指定自变量't', 'u'等ß微分方程的各阶导数项以大写字母D表示 或或或y的一阶导数—— Dyy的二阶导数—— D2yy的 n 阶导数—— Dny[y1,y2…]=dsolve〔x1,x2,…xn〕 —— 返回 微分方程的解一阶微分方程dsolve〔'Dx=y','Dy=x','x〔0〕=0','y〔0〕=1'〕ans =x〔t〕 = sin〔t〕, y〔t〕 = cos〔t〕二阶微分方程dsolve〔'D2y=-a^2*y','y〔0〕=1','Dy〔pi/a〕=0'〕ans =cos〔a*x〕例3.y=dsolve〔'D2y+2*Dy+2*y=0','y〔0〕=1','Dy〔0〕=0'〕ans =exp〔-x〕*cos〔x〕+exp〔-x〕*sin〔x〕ezplot〔y〕 —— 方程解y〔t〕的时间曲线图求该方程的解三、maple函数——符号运算的扩展maple——是专门进展数学运算的软件工具， 具有超强的符号运算才能，提供了 几乎包括所有数学领域的专用函数matlab——依赖于maple的内核与函数库，扩 展了自己的符号运算功能。

matlab还设计了对maple库函数的调用功能使得已有的maple数学功能，可以扩大matlab中,作为自身符号运算才能的扩展1. maple内核访问函数可以访问maple内核的matlab函数: maple ——— 访问maple内核函数 mapleinit —— maple函数初始化 mpa ———— maple函数定义 mhelp ——— maple函数帮助命令 procread —— maple函数程序安装. maple 的调用格式maple〔'表达式'〕 —— 将表达式送至maple内核， 返回符号表达式结果maple 〔'函数'，变量1，变量2〕 ——调用maple函数，传递给定 变量 例1. 展开5阶 bernoulli 多项式，计算 x=3 时bernoulli 数a=maple〔'bernoulli〔5,x〕'〕a =-1/6*x+5/3*x^3+x^5-5/2*x^4a=maple〔'bernoulli〔5,3〕'〕a =85例2. 化简三角函数式sin2x+cos2xa=maple〔'simplify〔sin〔x〕^2+cos〔x〕^2〕;'〕a =1例4. 求f〔t〕=e-3tsint的拉式变换f=maple〔'laplace〔exp〔-3*t〕*sin〔t〕,t,s〕;'〕f =1/〔〔s+3〕^2+1〕例4. 寻找二次多项式的完全平方 f 〔x〕 = x2+2x+2a=maple〔'completesquare〔x^2+2*x+2〕'〕a =completesquare〔x^2+2*x+2〕 将工具包装入内存maple〔'with〔student〕;'〕a=maple〔'completesquare〔x^2+2*x+2〕'〕a =〔x+1〕^2+1ßmaple软件中的所有函数，在初始化时并没有完全装入内存，可用readlib指令把库函数读入内存，或用with指令将应用工具包装入内存。

ß调用格式maple〔'readlib〔函数名〕;'〕maple〔'with〔工具包名〕;'〕例5.求sin〔x2+y2〕在x=0,y=0处泰勒级数展开式，8阶截断maple〔'mtaylor〔sin〔x^2+y^2〕,[x=0,y=0],8〕'〕ans =mtaylor〔sin〔x^2+y^2〕,[x = 0, y = 0],8〕maple〔'readlib〔mtaylor〕;'〕maple〔'mtaylor〔sin〔x^2+y^2〕,[x=0,y=0],8〕'〕ans =x^2+y^2-1/6*x^6-1/2*y^2*x^4-1/2*y^4*x^2-1/6*y^62. mpa —— maple变量定义ß任何一个matlab定义的函数f，可使用mpa语句直接调用，还可把 f 定义成maple变量vßmaple的工作空间与matlab工作空间是互相独立的， 所以f 与v是属于不同工作空间中的变量ßmpa的调用格式： mpa〔'v',f〕 mpa v ff为matlab工作空间中已存在的变量例. 电磁力计算公式为试I=0.5，x=0.1邻域展开泰勒级数，3阶截断，令常数 ，1.直接调用maple〔'readlib〔mtaylor〕;'〕maple〔'mtaylor〔k*I^2/x^2,[I=0.5,x=0.1],3〕;'〕2.定义符号函数f〔matlab6.1无map函数〕f='k*I^2/x^2';maple〔'mtaylor〔f,[I=0.5,x=0.1],3〕;'〕ans =mtaylor〔f,[I = .5, x = .1],3〕mpa〔'u',f〕maple〔'mtaylor〔u,[I=0.5,x=0.1],3〕;'〕ans =25.*k-.50e3*k*〔x-.1〕+.10e3*k*〔I-.5〕+7500.000000000000*k*〔x-.1〕^2+.1e3*k*〔I-.5〕^2-.20e4*k*〔I-.5〕*〔x-.1〕注意：matlab符号运算时，可以识别matlab定义的符号变量，但在调用 maple 函数时，需将matlab变量定义为maple变量后，所调用的函数方可识别和执行3.mhelp —— maple函数帮助命令 mhelp 是协助检索maple库函数的专用命令调用格式：mhelp 相关词条例如： mhelp intro — maple介绍 mhelp maple — maple命令格式 mhelp tutorial —maple入门 mhelp index —maple检索工具词条函数词条mhelp index 用于工具包检索library ——maple标准库函数packages —— 应用工具包libmisc —— 其它库函数statements —— maple语句描绘expressions —— maple表达式datatypes —— maple数据格式tables —— maple表格和阵列procedures —— maple程序misc —— maple其它应用一般帮助文本主要包括以下部分ßFUNCTION—— 函数功能说明ßCALLING SEQUENCE—— 调用格式ßPARAMETERS —— 调用参数说明ßSYNOPSIS —— 语法说明ßEXAMPLES —— 应用举例ßSEE ALSO —— 相关词条4.maple库函数maple库函数共分四类maple内部函数：驻留函数任何条件下都可调用 mhelp index[internal]  maple的外部函数—读库定义部分： 调用时先执行读库命令，因此与内部函 数一样可直接调用 mhelp index[external]  maple的外部函数—读库装入部分 maple其余外部函数需要在使用前执行maple〔'readlib〔函数名〕;'〕命令将其装入内存 mhelp index[libmisc]  maple的惰性函数—不能直接调用，还需一些函数如mod,evala,evalf 等才能调用 mhelp index[intert] 小 结 本节介绍了matlab语言的符号运算功能，通过学习应该掌握：ß掌握如何创立、修改符号矩阵ß掌握符号运算功能ßmaple函数调用ßmhelp检索。