Code covered by the BSD License

### Highlights from calc_lz_complexity

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

# calc_lz_complexity

### Quang Thai (view profile)

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 1.1

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

24 Sep 2012 1.2

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 1.3

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 1.4

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

15 Mar 2013 1.7

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.