Rank: 210046 based on 0 downloads (last 30 days) and 0 file submitted
photo

Harry Commin

E-mail
Company/University
Imperial College London
Lat/Long
51.49839019775391, -0.176586002111435

Personal Profile:
Professional Interests:
MIMO Radar

 

Watch this Author's files

 

Comments and Ratings by Harry View all
Updated File Comments Rating
14 Dec 2012 Bell Polynomials of the Second Kind Recursive algorithm for computing Bell polynomials of the second kind Author: Moysey Brio

I needed to alter the code to allow for vector inputs (producing a 3D array output). It may be useful for others if this functionality were added.

14 Dec 2012 Bell Polynomials of the Second Kind Recursive algorithm for computing Bell polynomials of the second kind Author: Moysey Brio

I have found this very useful for implementing Faa di Bruno's formula.

I thought the function inputs were a bit confusing for this, since I don't think it's useful to restrict Nin or Kin in this case. Therefore, for the nth derivative, we just use:

B = IncompleteBellPoly(n,n,g_derivatives);

and then use B(n,k) as required.

13 Dec 2012 Restricted Integer Composition Generates restricted and unrestricted integer compositions Author: Theophanes Raptis

This works correctly on every set of inputs I've given it so far (some thousands).

One slight issue is that the function prints compositions to the console, rather than returning them in an array.

A quick and ugly fix is to pass some variable "comps" into and out of the function N2N(). Then, replace the two disp() calls with "comps = [comps;cell]" and "comps = [comps;n]", respectively. I'm sure there's a much neater way of doing this, though.

13 Nov 2012 Khatri-Rao product An efficient implementation of the Khatri-Rao product. Author: Laurent S

The function help states that the result is an I*J-by-K matrix. This is not always true. Example:

kr(zeros(10000,0),zeros(2,0))

ans =

Empty matrix: 2-by-0

21 Mar 2010 avilib: a pedestrian AVI creator wrapper creates avi files using syntax akin to fopen / fwrite / fclose Author: us

This is absolutely excellent! I'll never use Windows Movie Maker again!

Just one thing to note is that a few more usage examples might be useful. For example, framerate is set using: "aid=aopen('-f','foo','-d',32);" ... not "aid=aopen('-f','foo','-d','32');" or "aid=aopen('-f','foo','-d 32');" etc.

This may be obvious, but it could save some people a little time.

Contact us