No BSD License  

4.49333

4.5 | 76 ratings Rate this file 177 Downloads (last 30 days) File Size: 2.92 MB File ID: #3013

HOSA - Higher Order Spectral Analysis Toolbox

by

 

12 Feb 2003 (Updated )

Spectral and polyspectral analysis, and time-frequency distributions.

| Watch this File

File Information
Description

There is much more information in a stochastic non-Gaussian or deterministic signal than is conveyed by its autocorrelation and power
spectrum. Higher-order spectra which are defined in terms of the higher-order moments or cumulants of a signal, contain this additional information. The Higher-Order Spectral Analysis (HOSA) Toolbox provides comprehensive higher-order spectral analysis capabilities for signal processing applications. The toolbox is an excellent resource for the advanced researcher and the practicing engineer, as well as the novice
student who wants to learn about concepts and algorithms in statistical signal processing.

The HOSA Toolbox is a collection of M-files that implement a variety of advanced signal processing algorithms for the estimation of cross- and auto-cumulants (including correlations), spectra and olyspectra,bispectrum, and bicoherence, and omputation of time-frequency
distributions. Based on these, algorithms for parametric and non-parametric blind system identification, time-delay estimation, harmonic retrieval, phase-coupling, direction of arrival estimation, parameter estimation of Volterra (non-linear) models, and adaptive linear prediction are implemented. Also included are algorithms for testing of Gaussianity and Linearity of a time series. A full tutorial and demo set are included in the toolbox.

Acknowledgements

This file inspired Fast Third Order Cross Cumulant Calculation.

Required Products Signal Processing Toolbox
MATLAB release MATLAB 6.5 (R13)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (116)
22 Sep 2014 yue

rate

14 Jul 2014 Khujasta

unable to install HOSA on MATLAB2012b..
can anybody help..i have followed all the steps in the read me text..but to no avail

19 Jun 2014 Tiffany

Anybody else having issues getting the 'Hosa Demo' to work? It tells me
'Undefined function 'choices' for input arguments of type 'char'. Granted, I can go into the folder and open any demo I need, but I'm not really sure what they all are.

Second question (and maybe I just need to read on this), but how does the window work? From the demo, it seemed he chose an arbitrary size? I'm looking to calculate the bi-spectrum. I'm used to making my own hanning window when I do a regular FFT; this appears to be different.

28 Aug 2013 yangyang  
21 Jul 2013 zhongguo

very good

08 Dec 2012 Damon Bradley

Also, to remove the annoyance of inexact case matches when you use the any of the HOSA functions, I'd recommend putting this line of code at the top of your script:

warning('off','MATLAB:dispatcher:InexactCaseMatch').

For some strange reason, all of the HOSA function file names are written in all caps.

13 Nov 2012 Lyuboslav Petrov

Convert from upper to lower all files in both folder using this script:

files = dir;
for i = 3:length(files)
system(['REN ',lower(files(i).name),' ',lower(files(i).name)]);
end

13 Nov 2012 Lyuboslav Petrov  
25 Sep 2012 Wolfgang

Hi,

i have tried as described in the readme.txt. Unfortunately all file names and extentions are in uppercase after extracting !

Because of that no chance to use anything (like hosaver or hosademo).

Please advice what to do.

03 Jul 2012 Russell

I have the HOSA pack working, but I run out of memory with the WIG2C function when using a array of 2500 rows and 1 column. Can the function not handle a array like this?

07 Mar 2012 MIAO DU

Hi,all, I tried to use the function bispecd.m on my data, but it always shows
??? Subscript indices must either be real positive integers or logicals.

Error in ==> BISPECD at 123
xseg = y(locseg);

How can I fix this problem?
Thanks.

18 Jan 2012 Keith Thornicroft

Hi All,

I can't seem to get the toolbox installed having followed the text file to the T. I am using MATLAB Student Version 7.12.0.635 (R2011a) with an additional wavelet toolbox. Any ideas? Does it even work for a student licence?

Cheers!

19 Sep 2011 Jered Wells

For those who get an error due to the info.xml file: to fix the error in hosa\info.xlm, the line containing <area>matlab</area> should be changed to <type>matlab</type>. Now, to try HOSA...

27 Apr 2011 XIAYANG xiayang

It is good and helpful.

22 Apr 2011 MD. TANJIMUDDIN

The installation seems to have worked incompletely. I get the following error:
>> HOSADEMO
??? Undefined function or method 'choices' for input arguments of type 'char'.
Error in ==> HOSADEMO at 68
choices('HoSatDemo', 'HOSA Demo', labels, callbacks, interruptible);
My Matlab cannot find the command/function "choices".
I have MATLAB Version 7.8.0.347 (R2009a)
Signal Processing Toolbox Version 6.11 (R2009a)
Can you please help me solve this?
Thanks
please tell me how to use these toolboxes.please reply on my email tanjii023@yahoo.com

21 Apr 2011 Evan

For those in need of choices.m

http://www.maths.lth.se/matstat/wafo/documentation/wafodoc/wafo/misc/choices.html

has the source code.

Create a new .m, paste in code, delete all code line numbers :)

06 Mar 2011 henry gorx

sorry, I was wrong. The comment was for other function (MATS).

05 Mar 2011 henry gorx

Hi, I tested the function Bicorrelation(xV,tauV,flag), how I must to interpret the results? for example, bicV=-0.18738 and cumbicV=0.70683. Many thanks for your help.

27 Jan 2011 leo

it is grateful

13 Jan 2011 wu

File "choices.m" is lost. Is the problem of version of matlab, and mine is 7.1?

12 Jan 2011 wu

great tool of researchers whose major is signal processing

07 Jan 2011 andrea ruju

Why the window specification for frequency-domain smoothing is differetn in bispectrum and bicoherence estimation?

Hi, i'm a phd student in coastal oceanography. My e-mail is matumaloa@gmail.com

04 Dec 2010 asandu asandu

Help: why the bispectrum of the gaussian noise does not equal zero in HOSA
Hi :
I am a student from Wuhan, China. I have download the toolbox of HOSA. Thank you very much of your job. BUt I have met a problem in compute the bispectrum of the gaussian noise, because the result of it does not equal to zero. I can not find the solution of it. can you help me? and thanks.
The signal of gaussian noise is given by the commad 'rand'.
my e_mail: xhdu1981@yahoo.com.cn

08 Nov 2010 Sayan

having errors please help
>> hosademo
??? Undefined function or method 'choices' for input arguments of type 'char'.

Error in ==> HOSADEMO at 68
choices('HoSatDemo', 'HOSA Demo', labels, callbacks, interruptible);

02 Nov 2010 wang ?

thank you

12 Oct 2010 Rafael J.

Great toolbox, but there are clearly some typos and errors. Capitalized filenames causes warnings (or the fact that they don't match the functions calls with respective to caps).

Also, why does bicoherx and bispecdx have three inputs when most people only have two inputs for a cross-bicoherence or cross-bispectrum analysis? Anyone figure this issue out? Which two inputs should we use? (eg x-x-y, x-y-x, etc)

19 May 2010 Salai Selvam V

Excellent!!!

18 May 2010 Ali Ali

Thanks alot. fantastic tool

06 Nov 2009 zamrud fathima

The installation seems to have worked incompletely. I get the following error:
>> HOSADEMO
??? Undefined function or method 'choices' for input arguments of type 'char'.
Error in ==> HOSADEMO
choices('HoSatDemo', 'HOSA Demo', labels, callbacks, interruptible);
My Matlab cannot find the command/function "choices".
Can you please help me solve this?
Thanks.
please send reply to my e-mail

03 Nov 2009 Rae

Is there a simple way to edit bicoher to normalize it such that it's giving the bounded (0-1) version rather than the current skewness version?

I can't access the Kim and Powers paper and if others have sucessfully done this, I'd love to know how!

11 Aug 2009 Paul

The installation seems to have worked incompletely. I get the following error:
>> HOSADEMO
??? Undefined function or method 'choices' for input arguments of type 'char'.
Error in ==> HOSADEMO at 68
choices('HoSatDemo', 'HOSA Demo', labels, callbacks, interruptible);
My Matlab cannot find the command/function "choices".
I have MATLAB Version 7.8.0.347 (R2009a)
Signal Processing Toolbox Version 6.11 (R2009a)
Can you please help me solve this?
Thanks,
Paul Axon, Sydney, Australia
paul.axon@thalesgroup.com.au

11 Jun 2009 David Gagnon

I meant to submit this post first...
(Both posts refer to BICOHER.m)

Beware that there are multiple ways to normalize the bispectrum, and they are all referred to as the "bicoherence." The formula that this m-file uses is NOT the one that bounds the bicoherence between 0 and 1. If you are looking for that one, check out Kim and Powers' 1979 paper "Digital Bispectral Analysis and Its Applications to Nonlinear Wave Interactions." They provide an explanation of how to code that version, and it has worked reasonably well for me.

11 Jun 2009 David Gagnon

By the way, "segsamp" (actually implemented as "nsamp") is the number of samples per segment, that is, the number of points in each chunk of data that you're throwing into the bicoherence. You can use datasets of any length, and the program will typically calculate fft's with a length that is a power of two and is greater than or equal to "segsamp." It does this by adding zeros to the data record. However, you can specify "nfft" to tell the program how many points you want the fft to use. If you set "segsamp" and "nfft" to be equal, then you can avoid zero-padding. Powers of two are calculated faster, but MATLAB can handle fft's of any length.

As for the frequency axis, "waxis" is based on the number of points that are used in the fft. As a result, "0.5" on the axis corresponds to the highest frequency that can be resolved with that number of points, also known as the Nyquist frequency. Exactly how you want to scale the axis will depend on the nature of your data.

On an unrelated note, the demo gives me errors, too.

04 Jun 2009 Rakesh

I have tried WIG2.m function on a times series on 28800 samples (288 sec signal sampled at 100Hz) but function doesnt seem to handle this big signal. However, I have also tried it on chunks of this signal but it takes forever to complete. Is it my computer or the program itself? any help?

08 Apr 2009 Gabriel Akira Schreiber

But, however, lowercase filenames would be a great improvement.

08 Apr 2009 Gabriel Akira Schreiber

This toolbox is greatly appreciated. Is this still developed/maintained?

27 Oct 2008 Edwin Wu

Does anybody have the experience to use TDE.M? I encounter some problem when I use it as a way to estimate the time delay of 2 channels of signal. I would be so appreciated for your helps.My msn is wjr0529@gmail.com,it is also my E-mail.I would like to make friends with anyone of you.join me , share our knowledge and be brilliant together.

21 Oct 2008 Bessie Bao

EEG

19 Oct 2008 elahe alavi  
27 Aug 2008 Yue ZHOU

To identify whether one time serial is norm-distribution

04 Aug 2008 Femi Osunnuyi

I've had some problems using the 'bicoher' function. With most industrial data sets, i obtain bicoherence values greater than 1, which in theory is not meant to be possible.
Any ideas?
Also is 'segasamp' meant to be square root of no. of data set (which must be even powers of 2 [4,16,64,256,1024], right?)?

28 Jul 2008 Alex Yakubovich

Does anybody what the frequency units are in BICOHER and BISPECD? Some of the examples in the tutorial suggest it is in Hz, but that doesn't make sense given the plot I get with my signal.

17 Jun 2008 Guy Incognito  
26 May 2008 jason yun  
19 May 2008 Kacimi Abderrahmane  
01 May 2008 Vinayak swarnkar

I downloaded it..but facing problem in running it. Even 'hosademo.m' dosent work...it gives error saying 'Undefined function or method 'choices' for input arguments of type 'char'.'. Can some please help in this.
Thanks you

26 Jan 2008 as as  
04 Dec 2007 bahman abc

thx!

07 Nov 2007 Alhade Mohamed

hi

11 Oct 2007 guang yang

it's good

30 Sep 2007 Morteza Dardel

I want some more description about “nlpow” command (estimation with volterra serises). If this command can be used for more than one input and output cases?

11 Sep 2007 navya showry

I have downloaded this file, but .exe file is missing from the extracted one.. plz advice how to proceed further or email me the latest updated one.

07 Sep 2007 Ali Ozgur Argunsah

hi,
?s there any linux matlab version of this toolbox?

23 Aug 2007 ahmed badawy

pls I hope to make good relation in working
in HOSA and Signal Separation Using HOSA

13 Jun 2007 yhaofly yhaofly

your work is execellent,thank you very much,
I will support you for ever.

01 Jun 2007 song jianjiang

good article

17 May 2007 Chang David

I want to know if the software is only compliant with Matlab 6.0(R12 version)

27 Jan 2007 Florencio Calle  
03 Jan 2007 Vishnuvenkatesh Dhage  
04 Dec 2006 liu xiangle

the help system is not covinent.

17 Sep 2006 Frederico Melo  
04 Jun 2006 Scott Miller

How do you open this thing? I have tried to extract it and nothing happens. It seems corrupted.

to reply take out chinese food

02 Jun 2006 M Nage

not see yet

27 May 2006 ali dashti

needed box

20 Apr 2006 Mohsen Afrasiabi  
19 Apr 2006 Sivaganesh Mullapudi  
23 Mar 2006 Jon Hobbs

Great Box

23 Dec 2005 meng yusong  
22 Dec 2005 Ghassan Maalouli

Excellent toolbox. Very good documentation.
Brilliant author.

05 Dec 2005 Jon-Fredrik Hopperstad

Excellent toolbox with very good documentation.

01 Dec 2005 essam anwer

i can not run this toolbox how can download and who can help me

27 Nov 2005 xu guisen

how to download? who teach me!

15 Nov 2005 Neil L Scott

Demo fails.
Crashes with GUI callback issues.
Uppercase file names???

09 Oct 2005 Bruce Lena

great tools for me

03 Aug 2005 xunkai wei

great

07 May 2005 xjtuljp Lin  
19 Apr 2005 chen xushuang

very good

10 Apr 2005 ZHONGCHAO HUANG

It is really good and very helpful to me

11 Mar 2005 Marcin J...

I would be thankfull if you can send me more informations about high order moments. Thank for help.

08 Mar 2005 farzaneh shayegh  
15 Jan 2005 hawk wei

Anyway but can you share it ?
It's really great

10 Jan 2005 Chuanshuang Hu

It looks great. I need HOSA. Also, could you tell me whether it is possible to compute Auto-and cross-spectra or not using HOSA?

11 Dec 2004 Wei Zhang

why this toolbox is not concluded in Matlab 7 ?

22 Nov 2004 Nej Nej

Kindly send me a HOSA tool box for study purpose

09 Nov 2004 Kim Gunnium

I need HOSA for learning.

05 Nov 2004 Kianoush Nazarpour  
03 Nov 2004 Jim Green

I need much higher order spectral analysis toolbox

14 Aug 2004 sina niki

Dear Sir,
Problem:
We want simulate by matlab "Passive Radar localization by time difference of arrival esimator".
For this purpose we have two site on the ground and we want by two widwband receivers estimate two same received pulse at separate station.
could you help me to send m.files for this purpose.
thanks your helps.
Regards.

09 Aug 2004 Feng Changjian

I have like this toolbox because it can help me to complete my school work.

04 Aug 2004 abdelali belabassi

i'm etudiant ,i need it to complet my thesis , thanks

27 Jun 2004 Jianwen Luo  
18 Jun 2004 Marcus Kruger  
15 Jun 2004 xnkj xnkj

hao dong dong

11 Jun 2004 zheng qing  
30 May 2004 wang shenghui

ok

29 Apr 2004 Dr. Abdul-Bary Sulaiman

A very useful toolbox, but i\I surprised that it is only found in V5.3 Matlab and has been cancelled in the uper versions !!!

28 Feb 2004 Camilo Londoño

Nice work!!!!

22 Dec 2003 Brian Powell

Excellent toolbox. The only negative is that it is distributed as an EXE. If you distribute as a tar, it should work everywhere. After unpacking on a borrowed windows box, I use it on Solaris and MacOSX.

18 Dec 2003 Norbert Marwan

Why this toolbox is not provided as a zip-package or Matlab-installation file (eg. created with Makeinstall). Moreover, all filenames are upperletter. The exe-package as well as the upperletter filenames make it hard for a Unix user to use this toolbox.

11 Nov 2003 Jaime Aviles  
10 Sep 2003 qiang xueming

excellent

26 Aug 2003 Dale Clifford

This IS a matlab toolbox and installs the same even as far as the icon in the view window of Matlab. Previously not available as a student edition, but now free.This toolbox provides much information on signals, especially if the researcher is looking for non-linear dynamics. Works well when used along the side of wavelets.

12 Aug 2003 w w

very good

05 Aug 2003 Moh Al-jebori  
06 Jul 2003 Mokrane Malek  
16 Jun 2003 dylan mckay  
01 Jun 2003 extra ll  
13 May 2003 Ma Yongfeng

It is good and helpful.

10 Apr 2003 x yg  
20 Mar 2003 guojun yang  
19 Mar 2003 Usama Zaher

Excellent application

If there are more examples it would be nice.Especially if it is concerned with some application of sensor validation

05 Mar 2003 D. Gansawat  
28 Feb 2003 Dr. Cornelis Los

This the best non-wavelet, higher-order time-frequency analysis of time series and the only one freely available. Use it closely together with the excellent book for the best results and for easy and accurate interpretation.

21 Feb 2003 George Barakos

Pity it is not part of the main MatLab release!

17 Feb 2003 - -

Is this the same as TMW's HOSA Toolbox?

16 Feb 2003 John Brown  

Contact us