zoukankan      html  css  js  c++  java
  • MATLAB符号计算

    1.  建立符号对象

    建立符号变量和符号常量
    MATLAB提供了两个建立符号对象的函数:sym和syms,两个函数的用法不同。
      (1) sym函数
      sym函数用来建立单个符号量,一般调用格式为:
    符号量名=sym('符号字符串')
    该函数可以建立一个符号量,符号字符串可以是常量、变量、函数或表达式。
    应用sym函数还可以定义符号常量,使用符号常量进行代数运算时和数值常量进行的运算不同。下面的命令用于比较符号常量与数值常量在代数运算时的差别。

      (2) syms函数
      函数sym一次只能定义一个符号变量,使用不方便。MATLAB提供了另一个函数syms,一次可以定义多个符号变量。syms函数的一般调用格式为:
    syms  符号变量名1 符号变量名2 … 符号变量名n
    用这种格式定义符号变量时不要在变量名上加字符串分界符(‘),变量间用空格而不要用逗号分隔。

    2. 符号表达式运算

      1 符号表达式的四则运算
      符号表达式的加、减、乘、除运算可分别由函数symadd、symsub、symmul和symdiv来实现,幂运算可以由sympow来实现。
      2 符号表达式的提取分子和分母运算
      如果符号表达式是一个有理分式或可以展开为有理分式,可利用numden函数来提取符号表达式中的分子或分母。其一般调用格式为:
      [n,d]=numden(s)
      该函数提取符号表达式s的分子和分母,分别将它们存放在n与d中。

      3 符号表达式的因式分解与展开
      MATLAB提供了符号表达式的因式分解与展开的函数,函数的调用格式为:
      factor(s):对符号表达式s分解因式。
      expand(s):对符号表达式s进行展开。
      collect(s):对符号表达式s合并同类项。
      collect(s,v):对符号表达式s按变量v合并同类项。

      4 符号表达式的化简
      MATLAB提供的对符号表达式化简的函数有:
      simplify(s):应用函数规则对s进行化简。
      simple(s):调用MATLAB的其他函数对表达式进行综合化简,并显示化简过程。

      5 符号表达式与数值表达式之间的转换
      利用函数sym可以将数值表达式变换成它的符号表达式。
      函数numeric或eval可以将符号表达式变换成数值表达式。

    3 . 符号表达式中变量的确定

      MATLAB中的符号可以表示符号变量和符号常量。findsym可以帮助用户查找一个符号表达式中的的符号变量。该函数的调用格式为:
      findsym(s,n),函数返回符号表达式s中的n个符号变量,若没有指定n,则返回s中的全部符号变量。

    4. 符号矩阵

      符号矩阵也是一种符号表达式,所以前面介绍的符号表达式运算都可以在矩阵意义下进行。但应注意这些函数作用于符号矩阵时,是分别作用于矩阵的每一个元素。
    由于符号矩阵是一个矩阵,所以符号矩阵还能进行有关矩阵的运算。MATLAB还有一些专用于符号矩阵的函数,这些函数作用于单个的数据无意义。例如
    transpose(s):返回s矩阵的转置矩阵。
    determ(s):返回s矩阵的行列式值。
    其实,曾介绍过的许多应用于数值矩阵的函数,如diag、triu、tril、inv、det、rank、eig等,也可直接应用于符号矩阵。

    5 . 符号微积分

       1 符号极限

    limit函数的调用格式为:
    (1) limit(f,x,a):求符号函数f(x)的极限值。即计算当变量x趋近于常数a时,f(x)函数的极限值。
    (2) limit(f,a):求符号函数f(x)的极限值。由于没有指定符号函数f(x)的自变量,则使用该格式时,符号函数f(x)的变量为函数findsym(f)确定的默认自变量,即变量x趋近于a。

    (3) limit(f):求符号函数f(x)的极限值。符号函数f(x)的变量为函数findsym(f)确定的默认变量;没有指定变量的目标值时,系统默认变量趋近于0,即a=0的情况。
    (4) limit(f,x,a,'right'):求符号函数f的极限值。'right'表示变量x从右边趋近于a。
    (5) limit(f,x,a,‘left’):求符号函数f的极限值。‘left’表示变量x从左边趋近于a。

     1 %限1:
     2 syms a m x;
     3 f=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/(x+a);
     4 limit(f,x,a)
     5 ans =
     6 (1/2*a*exp(sin(a))+1/2*a-exp(tan(a))+1)/a
     7 %限2:
     8 syms x t;
     9 limit((1+2*t/x)^(3*x),x,inf)
    10 ans =
    11 exp(6*t)

      2  符号导数

    diff函数用于对符号表达式求导数。该函数的一般调用格式为:
    diff(s):没有指定变量和导数阶数,则系统按findsym函数指示的默认变量对符号表达式s求一阶导数。
    diff(s,'v'):以v为自变量,对符号表达式s求一阶导数。
    diff(s,n):按findsym函数指示的默认变量对符号表达式s求n阶导数,n为正整数。
    diff(s,'v',n):以v为自变量,对符号表达式s求n阶导数。

      3  符号积分

    符号积分由函数int来实现。该函数的一般调用格式为:
    int(s):没有指定积分变量和积分阶数时,系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分。
    int(s,v):以v为自变量,对被积函数或符号表达式s求不定积分。
    int(s,v,a,b):求定积分运算。a,b分别表示定积分的下限和上限。该函数求被积函数在区间[a,b]上的定积分。a和b可以是两个具体的数,也可以是一个符号表达式,还可以是无穷(inf)。当函数f关于变量x在闭区间[a,b]上可积时,函数返回一个定积分结果。当a,b中有一个是inf时,函数返回一个广义积分。当a,b中有一个符号表达式时,函数返回一个符号函数。

      4  积分变换

      常见的积分变换有傅立叶变换、拉普拉斯变换和Z变换。

        1.傅立叶(Fourier)变换
        在MATLAB中,进行傅立叶变换的函数是:
        fourier(f,x,t):求函数f(x)的傅立叶像函数F(t)。
        ifourier(F,t,x):求傅立叶像函数F(t)的原函数f(x)。

        2.拉普拉斯(Laplace)变换
        在MATLAB中,进行拉普拉斯变换的函数是:
        laplace(fx,x,t):求函数f(x)的拉普拉斯像函数F(t)。
        ilaplace(Fw,t,x):求拉普拉斯像函数F(t)的原函数f(x)。

        3.Z变换
        当函数f(x)呈现为一个离散的数列f(n)时,对数列f(n)进行z变换的MATLAB函数是:
        ztrans(fn,n,z):求fn的Z变换像函数F(z)。
        iztrans(Fz,z,n):求Fz的z变换原函数f(n)。

    6.  级  数

      1   级数符号求和

      求无穷级数的和需要符号表达式求和函数symsum,其调用格式为:
      symsum(s,v,n,m)
      其中s表示一个级数的通项,是一个符号表达式。v是求和变量,v省略时使用系统的默认变量。n和m是求和的开始项和末项。

      2  函数的泰勒级数
      MATLAB提供了taylor函数将函数展开为幂级数,其调用格式为:
      taylor(f,v,n,a)
      该函数将函数f按变量v展开为泰勒级数,展开到第n项(即变量v的n-1次幂)为止,n的缺省值为6。v的缺省值与diff函数相同。参数a指定将函数f在自变量v=a处展开,a的缺省值是0。

      

  • 相关阅读:
    121.买卖股票 求最大收益1 Best Time to Buy and Sell Stock
    409.求最长回文串的长度 LongestPalindrome
    202.快乐数 Happy Number
    459.(KMP)求字符串是否由模式重复构成 Repeated Substring Pattern
    326.是否为3的平方根 IsPowerOfThree
    231.是否为2的平方根 IsPowerOfTwo
    461.求两个数字转成二进制后的“汉明距离” Hamming Distance
    206.反转单链表 Reverse Linked List
    448. 数组中缺少的元素 Find All Numbers Disappeared in an Array
    常见表单元素处理
  • 原文地址:https://www.cnblogs.com/hxsyl/p/2787577.html
Copyright © 2011-2022 走看看