KroneckerDelta in matlab

前端 未结 4 1893
心在旅途
心在旅途 2021-01-21 02:01

This link shows that there is a kronecker delta function in matlab. However:

>> help kroneckerDelta

    kroneckerDelta not found

I am u

相关标签:
4条回答
  • 2021-01-21 02:11

    The Kronecker delta returns 1 if j==k...

    So you could simplify the expression with:

    function d=kronDel(j,k)
    
    d=j==k
    
    end
    

    Luckily, MATLAB represents booleans as (0,1)

    0 讨论(0)
  • 2021-01-21 02:12

    Your link to is to the MuPAD function kroneckerDelta -note the URL and the funky typography of the examples. You won't see it in any version of Matlab because it's only available through MuPAD (type mupad in your command window and try in the window that launches). I have no idea when it was added to MuPAD, I know that it's at least in R2012b. You may have it even if the help command returns nothing.

    If you have kroneckerDelta in R2011b, you won't be able to run it from the regular command window or Editor in the normal fashion.

    evalin(symengine,'kroneckerDelta(1,1)')
    

    or the more flexible

    feval(symengine,'kroneckerDelta',1,1)
    

    See more here. However, if you're not working with symbolic math, there's really no reason to use this function that I can see -it's not even vectorized! I'd go with a solution that fully emulates kroneckerDelta's behavior in double precision:

    function d=kronDel(m,n)
    if nargin == 1
        d = double(m==0);
    else
        d = double(m==n);
    end
    
    0 讨论(0)
  • 2021-01-21 02:25

    I don't see it in my R2012b so perhaps not. Unless you need the symbolic math, you could always write your own. Something as simple as

    function d = kronDel(j,k)
    
    if j == k
        d = 1;
    else
        d = 0;
    end
    
    0 讨论(0)
  • 2021-01-21 02:33

    You can also do this inline, like

    ( a == b )
    

    However, using an anonymous function is a nice way to convert a one liner like this into something that is more readable

    kronDel = @(j, k) j==k ;
    
    kronDel( 2, 1 )
    kronDel( 2, 2 )
    
    0 讨论(0)
提交回复
热议问题