Code covered by the BSD License

### Highlights from calc_lz_complexity

Be the first to rate this file! 23 Downloads (last 30 days) File Size: 6.28 KB File ID: #38211

# calc_lz_complexity

17 Sep 2012 (Updated )

Calculates the Lempel-Ziv complexity of binary sequence - a measure of its "randomness"

File Information
Description

This function calculates the complexity of a finite binary sequence, according to the work presented by Abraham Lempel and Jacob Ziv in the paper "On the Complexity of Finite Sequences", published in "IEEE Transactions on Information Theory", Vol. IT-22, no. 1, January 1976.

From that perspective, the algorithm could be referred to as "LZ76".

The function supports two methods of evaluating sequence complexity:
1. Decomposition into an exhaustive production process
2. Decomposition into a primitive production process
Exhaustive complexity can be considered a lower limit of the complexity measurement approach proposed in LZ76, and primitive complexity an upper limit.

Currently, only sequences with binary alphabets (0, 1) are supported.

Feel free to email me if you find this function useful, find bugs with it, or have any suggestions for improvements.

Required Products MATLAB
MATLAB release MATLAB 7.13 (R2011b)
18 Sep 2012

Forgot to include a helper function that's used by the main calc_lz_complexity() function.

24 Sep 2012

Changed the way normalized complexity was calculated. It is now normalized against (n/log2(n)), rather than n, where n is the length of the sequence.

Upon reviewing Lempel-Ziv's paper, this normalization seems to make more sense.

12 Nov 2012

Major speed improvement. The old version was actually very inefficient and slow. A rethink of the eigenfunction calculation was warranted.

On my PC, this new version runs about 140x faster than the previous version!!!

13 Nov 2012

Another optimisation tweak that should result in this implementation running faster again for most sequences than the previous version.

15 Mar 2013

The eigenfunction that corresponds with the input sequence, and calculated internally, can now be obtained by the caller of calc_lz_complexity() as a return value.