COMPGROUPS.NET | Search | Post Question | Groups | Stream | About | Register

### Cepstrum Help

• Email
• Follow

```Hi Guys,

I'm trying to calculate the cepstrum of a signal to see if I can find
some similarities between it and another signal.  I am doing the
following, but I'm not sure the data coming out looks quite right.  I
have an 8000hz signal and 512 sample window.  I am passing this to
FFTW and then on each result, doing...

float sample = (float)Math.Sqrt(Math.Pow(complex1,
2) + Math.Pow(complex2, 2));

to get the magnitude, and then

sample = (float)Math.Log(Math.Pow(sample, 2));

to get the log (as per what I've been reading on cepstrums in
wikipedia).  I then run the resulting array again through FFTW this
time with of course a 256 sample window, and then on each of the
results I again do...

float sample = (float)Math.Sqrt(Math.Pow(complex1,
2) + Math.Pow(complex2, 2));

to get the magnitude.

Am I missing something here, or do I have all the steps?

Thanks
Ray

```
 0

See related articles to this posting

```On Nov 23, 11:01=A0pm, Raeldor <rael...@gmail.com> wrote:
> Hi Guys,
>
> I'm trying to calculate the cepstrum of a signal to see if I can find
> some similarities between it and another signal. =A0

> Am I missing something here, or do I have all the steps?

As I recall, computing the cepstrum is a 3-stage procedure:

1) Compute the DFT of the data
2) Compute some logs of the DFT
3) Compute the IDFT of the LogDFT.

You will of course need to verify this with a textbook.

Apart from that, be aware that the cepstrum is a mess.
A simplified log, like you seem to be doing, where the
phase of the DFT is ignored, will cause you to mis-represent
the signal and so the cepstra you compare might not
correspond quite the way you think or suppose they do.

Computing the 'exact' log of the complex DFT coefficients
is an even greater mess, as you will be dealing with all
kinds of uniqueness issues caused by the complex log not
being well-behaved.

There are very good reasons why the cepstrum is not used
more than it actually is. You need to be aware of these
reasons, and deal with every single one of them, to have
any chance of succeeding.

Rune
```
 0

```On Nov 23, 9:55=A0pm, Rune Allnor <all...@tele.ntnu.no> wrote:
> On Nov 23, 11:01=A0pm, Raeldor <rael...@gmail.com> wrote:
>
> > Hi Guys,
>
> > I'm trying to calculate the cepstrum of a signal to see if I can find
> > some similarities between it and another signal. =A0
> > Am I missing something here, or do I have all the steps?
>
> As I recall, computing the cepstrum is a 3-stage procedure:
>
> 1) Compute the DFT of the data
> 2) Compute some logs of the DFT
> 3) Compute the IDFT of the LogDFT.
>
> You will of course need to verify this with a textbook.
>
> Apart from that, be aware that the cepstrum is a mess.
> A simplified log, like you seem to be doing, where the
> phase of the DFT is ignored, will cause you to mis-represent
> the signal and so the cepstra you compare might not
> correspond quite the way you think or suppose they do.
>
> Computing the 'exact' log of the complex DFT coefficients
> is an even greater mess, as you will be dealing with all
> kinds of uniqueness issues caused by the complex log not
> being well-behaved.
>
> There are very good reasons why the cepstrum is not used
> more than it actually is. You need to be aware of these
> reasons, and deal with every single one of them, to have
> any chance of succeeding.
>
> Rune

Thank you for your informative (and honest) reply.  I'll try and look
for other ways to find features first, but I may come back to this
later if I can't see any other avenues.  Thank you!
```
 0

```On Nov 24, 4:33=A0pm, Raeldor <rael...@gmail.com> wrote:
> On Nov 23, 9:55=A0pm, Rune Allnor <all...@tele.ntnu.no> wrote:

> > There are very good reasons why the cepstrum is not used
> > more than it actually is. You need to be aware of these
> > reasons, and deal with every single one of them, to have
> > any chance of succeeding.
>
> > Rune
>

It's not about honesty, but reality. An idea is not good just
because it is new or cool. Or arcane. Once you read up on
a subject, you will find that most textbooks cover more or
less the same subjects in more or less the same order - that's
certainly the case with DSP.

As you keep reading, you find all kinds of seemingly nifty
tricks and ideas. Like the cepstrum. If you track the publishing
history, you very often find that the ideas are old. I know the
xcepstrum was extensively studied in the late '60s / early '70s,
with some papers and textbooks published. But the material is
nowhere to be seen in the textbooks, except as asides.

Once you discover these kinds of things, start asking 'why':

- Why isn't this idea covered in more depth in more places?
- Why isn't this idea used in more applications?

and so on.

That way you very soon find reasons like I outlined about the
cepstrum, and you know why to be cautious about the technique,
and when it might be used.

Rune
```
 0

```On Nov 24, 5:55=A0am, Rune Allnor <all...@tele.ntnu.no> wrote:
> On Nov 23, 11:01=A0pm, Raeldor <rael...@gmail.com> wrote:
>
> > Hi Guys,
>
> > I'm trying to calculate the cepstrum of a signal to see if I can find
> > some similarities between it and another signal. =A0
> > Am I missing something here, or do I have all the steps?
>
> As I recall, computing the cepstrum is a 3-stage procedure:
>
> 1) Compute the DFT of the data
> 2) Compute some logs of the DFT
> 3) Compute the IDFT of the LogDFT.

No, stage (3) is just a DFT.

Cepstrum is the DFT of the log of the DFT.

For implementation (real, complex, power etc) then yes you must look
for the many variants.

In case anyone is interested, I invented a similarly weird use of the
DFT in ultrasound prostate (rectal) scanning that is called the
'sprectum'). :-)

Chris
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
Chris Bore
BORES Signal Processing
www.bores.com

>
> You will of course need to verify this with a textbook.
>
> Apart from that, be aware that the cepstrum is a mess.
> A simplified log, like you seem to be doing, where the
> phase of the DFT is ignored, will cause you to mis-represent
> the signal and so the cepstra you compare might not
> correspond quite the way you think or suppose they do.
>
> Computing the 'exact' log of the complex DFT coefficients
> is an even greater mess, as you will be dealing with all
> kinds of uniqueness issues caused by the complex log not
> being well-behaved.
>
> There are very good reasons why the cepstrum is not used
> more than it actually is. You need to be aware of these
> reasons, and deal with every single one of them, to have
> any chance of succeeding.
>
> Rune

```
 0

4 Replies
359 Views

Similar Articles

12/11/2013 1:34:15 PM
[PageSpeed]

Similar Artilces:

Help!! Help!! Help!!
I wish to vectorize the following code: s = tf('s'); W = logspace(-1,2,50); COMBINATION = combn(W, 2); K = cell(length(COMBINATION), 1); K0 = 1 + (1/(2*s)) + 0.5*s/(0.2*s + 1); % PID controller to be approximated as PI gapValue = zeros(length(COMBINATION), 1); for i = 1:length(COMBINATION) K{i} = (COMBINATION(i,1)*(1 + (1/(COMBINATION(i,2)*s)))); % PI controller end for i = 1:length(COMBINATION) gapValue(i,1) = gapmetric(K0, K{i}); end where K0 is the Nominal transfer function and K is a cell array including transfer functions whose gap metric from K0 are...

Help... Help
3500.0&nbsp; store&nbsp; in the&nbsp; chip cup is&nbsp; 0x45,0x5A,0xC0,0x00 that&nbsp; chip cup send&nbsp; 0x45,0x5A,0xC0,0x00 to computer with&nbsp; NI-CAN how can i change the&nbsp; 0x45,0x5A,0xC0,0x00 to 3500.0 in the&nbsp; labview? how can i change the&nbsp; 3500.0 to 0x45,0x5A,0xC0,0x00 in the&nbsp; labview? ...

Re: Help ! Help!! Help!!! revisited
In the spirit of the new year, I decided this morning to review some of my earliest posts in this newsgroup In the process, I came across the following: % Newsgroups: comp.lang.prolog % From: t...@ecrc.de (Thom Fruehwirth) % Subject: Re: Help ! Help!! Help!!! % Message-ID: <Cot5rE.Cwu@ecrc.de> % Sender: n...@ecrc.de % Reply-To: t...@ecrc.de % Organization: European Computer-Industry Research Centre GmbH. % References: <2o49fk\$ab9@email.tuwien.ac.at> % Date: Mon, 25 Apr 1994 09:16:25 GMT % Lines: 22 % % Somebody [namely me] wrote: % % : Given *only* the facts % % : a(1). % : ...

HELP!! HELP!! HELP!! failure memory address
Por favor, si alguien me puede ayudar, tengo un notebook dell latitude c640 y desde ayer me arroja un error al encenderce Memory write/read failure at 1FFF0000, read FF107F10 expecting FF10FF10 Memory address line failure at 1FFF0000, read FF10F10 expecting 00000000 Strike the F1 key to continue, F2 to run the setup utility Y no puedo entrar al pc... intent� reparar el sistema operativo, formatear y no me deja :( tengo windows xp professional instalado Por favor, si alguien me puede ayudar, guiar, orientar... se lo agradeceria mucho!!! ���por favor!!! de verdad estoy muy urgida con...

Missing Help topics in Help
I was trying to look up Access Help for the GetSetting, SaveSetting commands to implement in my DB code when i press F1, up pops the Microsoft Access Help, if i click Index and type in Get, all that comes up in the box below is : GetChunk GetHiddenAttribute GetOption The Help Menu in Access looks like this Microsoft Access Help F1 Show the Office Assistant What's this? Office on the Web Detect and Repair About Microsoft Access I had a look in the Access VB Help too , but that was the same too? Am i missing something, i checked that Access Help was fully installed using Add/Remove...

Help
hi Gurus, Here my data structure 1234 => { 'open' => '1234-1', 'analyzed' => '1234-2', 'feedback' => '1234-3', 'closed' => '1234-4' } 3455 => { 'open' => '3455-1', 'analyzed' => '34552', 'feedback' => '3455-3', 'closed' => '3455-4' } My task is to display numbers [like 1234 or 3455] if they have ('open&...

Help
How do you insert an image (logo) into an existing PDF? sjb schrieb: > How do you insert an image (logo) into an existing PDF? This is a FAQ. Please use Google Groups to find answers. Use "pdfpages" as a keyword. Short answer: commercial: Adobe Acrobat, Touch-Up tool Foxit PDF Editor vector graphics programs with PDF import and export free: pdflatex and pdfpages and eso-pic packages details: http://www-user.tu-chemnitz.de/~ralk/latex/pdfoverlay.php Ralf -- Ralf Koenig Wissenschaftlicher Mitarbeiter an der Professur Rechnernetze und verteilte Systeme TU Chemnitz, Zi. 1/...

Help
I would like some help about how to programme a neural network in fortran 90 or 95.Any piece of advice will be welcome (books,web pages,examples) Thanks a lot Vega Vega wrote: > I would like some help about how to programme a neural network in > fortran 90 or 95.Any piece of advice will be welcome (books,web > pages,examples) There is Fortran 90 code for neural networks at http://mensch.org/neural/nnet.tar.gz and http://www.mscs.mu.edu/~ffritz/User_guide/Examples/Least_Squares/NN_train_3/train.f90 ,and Fortran 77 codes at http://dphs10.saclay.cea.fr/Spp/Experiences/OPAL/opalcern/nn...

help
hi, i have a problem since last week. i can not log in me hotmail email at home, but at office every thing is OK. (I use the same lablop). i use verizon DSL at home and already called them. they told me every thing is fine after they checked the signal. so i donot know how to do. please help. thanks This is a MIME GnuPG-signed message. If you see this text, it means that your E-mail or Usenet software does not support MIME signed messages. --=_mimegpg-commodore.email-scan.com-26410-1119042322-0001 Content-Type: text/plain; format=flowed; charset="US-ASCII" Content-Disposition: in...

help
Hi! I am trying to design a dc motor controller in a Spartan 3. Does somebody know how i can generate the speed reference for this issue? I have designed with a vhdl module reading values of an array but it doesn't work, any way? thanks in advance What kind of dc motor? Step motor? Hi! It is a Maxom Motor of direct current but it is not a step motor. It is the Maxon Re40 Graphite Brushes of 150 Watt. Do you know it? ...

HELP
Hello everyone, Being on the brink of purchasing a new laptop, I have a query concerning hp's offerings as mentioned in the title. the machines seem identical but for a few items, namely : dv5166ea dv5140eu -------- -------- Proc : 1.66Ghz Intel Core AMD TURION 64 Mobile Processor ML-32 Duo Processor T2300 with PowerNow Technology (comment : no idea what differences this makes) Cache : Level 2 cache 2MB Level 2 cache 512KB (comment : much more high speed memory for 5166 right ? faster machine ?) Memory :...

HELP
I'm Italian and I'm sorry for my English!!!I'm new of this ns and so I don't know if my question is suitable. I have a problem: using the C++ language,the QT library and the OpenGL library realize a CAD that can: a) load a glut three-dimensional model b) load a three-dimensional model in dxf format Thank you very much "Mj23" <snowteo13@libero.it> wrote... > I'm Italian and I'm sorry for my English!!!I'm new of this ns and so I don't > know if my question is suitable. The problem I see is that there _is_ no question. > I have a probl...

help me.......
Hello, excused my English, I am writing from Italy I have a problem, I have realized a program in vb6 for changed fields in rows pdf. in some cases the my file product hey give an error to me 20 in the opening file, but i view files only contained for fields but no grafica. help me....... ...

Help!!!!!!!!!!!!!!!
Hey i am a newbie in d feild of VHDL. Can any1 tell me how to interface FPGA'S and CPLD'S with d system(computer) and how a VHDL code is fed into d FPGA or CPLD..Thx for d kind help in advance.. In article <1110265377.226619.190750@g14g2000cwa.googlegroups.com>, celebration.of.being@gmail.com says... > Hey i am a newbie in d feild of VHDL. Can any1 tell me how to interface > FPGA'S and CPLD'S with d system(computer) and how a VHDL code is fed > into d FPGA or CPLD..Thx for d kind help in advance.. > Since this is a newsgroup read also by people not so fam...

HELP ME
Hello All As part of our Company Expansion we are insatalling a LAN in our network also we are looking for a solution to connect our remote office to remote site. the current configuration in Remote office and the HO Is as follows Remote Site ISP<dynamic IP> LINKSYS WIRELSEE ROUTER<Ethernet connection to>Desktop PC Head Office ISP<Dynamic IP> Lucent ADSL modem< Ethernet cable connection to> 3com office connect firewall<from firewall to Dlinkswitch connected by ethernet cable>Dlink Switch All other PC's are connected to the Switch Now We are lookin...

help ?
(sorry on my bad english) today I installed solaris 10 for a first time a everything went ok,I have installed them after 6 hours trying... A than I've run into problems. after I run solaris I have been promted for login name and password and that is a problem, I dont know them. For a time of installation I was promted to write domain name , I did, and for login passw I left empty space... so please tell me what do I do or at least where to find answer. p.s I know how to use windows very good and programming under them so i am try to learn something about programming under linux so plea...

help me
I got a difficult problem about javascript. Can someone please help me? If I use the <form> tag, and I use the <input type ="text" name="input"> tag, I can get the value by using document.formname.input.value , but if I am using <input type="file" name="file"> to browse for files, how can I get its value? Anyone who are interested to answer my question please email to me at cyx.thunderfoot@gmai.com , thanks a lot. cyx.thunderfoot@gmail.com wrote: > I got a difficult problem about javascript. Can someone please help me? > If I ...

Help!
I have a four page document in Adobe used for travel expenses and we have 300 to do! The information need for the forms has been kept in an MSAccess Database. At the beginning of this process, I exported information into an .xml file. When I import the .xml file in to the claim only one record appears. When I scroll through the document only the first record is completed in the form. How do I get the form to display all records? I am a rookie with Adobe so please have patience with me!~ carla.burton@gmail.com wrote: >I have a four page document in Adobe used for travel expenses and ...

Help me
Who can give me(A fresh man) some advice about the FPGA learning,and how to become an expert in it.Please recommend some materials to me. Thanks! zqhpnp@gmail.com wrote: > Who can give me(A fresh man) some advice about the FPGA learning,and > how to become an expert in it.Please recommend some materials to me. > Thanks! > First, even though vendors and some professors may think you can do complex things easily with the modern tools without a whole lot of knowledge about the underlying hardware (a.k.a System Generator), Learning basic logic is ******VERY******* important. I ...

help me!
i want get a gmail accont, E-mail to:aheadache@163.com thank you aheadache@163.com said the following on 10/23/2005 2:01 AM: > i want get a gmail accont, E-mail to:aheadache@163.com > thank you Apply for it. Sheesh. But, what in the name of Hades does this have to do with Javascript? -- Randy comp.lang.javascript FAQ - http://jibbering.com/faq & newsgroup weekly ...