DSPRelated.com
Forums

DSP: basic signal detection in Matlab (help!) :)

Started by theColdest March 2, 2005
Hi, I'm trying to write my first signal detection algorithm in
Matlab... but it's been several years since I've done any DSP... so
it's slow going.  if anyone has a little free-time, can you look at
this and see if you can spot what I'm doing wrong?

the data files can be acquired here:
http://www.podoboq.com/podoboq/basic_detect.zip

the two examples are simple (Class 1 - everything known a priori)
versions of what I am trying to do.  Unfortunately, I only have 299
data points to work with, and a sampling rate of 1Hz (Nyquist ==
0.5Hz).  I am using a Wiener estimation filter for pre-processing, and
for the detection I am using an uncorrelated Wiener-style
high-resolution detection filter.  "run_me1.m" uses a single sine pulse
centered at 100 seconds for the signal I'm trying to detect, and
"run_me2.m" has two sine pulses, one centered at 100 seconds and the
other centered at 250 seconds.

I'm not sure if I'm making a DSP-related error or just an error in the
way I'm manipulating the data in Matlab, but it doesn't seem to work.
after the estimation phase, you can see the pulse coming out at 100
seconds, but then the final output seems like garbage.  Aren't you
supposed to get a narrow impulse at each place where a detection
occurs???

What about SNR?  In this case the ratio is 1:1.  Do I need to
RMS-normalize something before filtering?  I've experimented with all
sorts of things, but maybe I'm just forgetting something dumb!
 
I'd hugely appreciate any help.  thnx!!!

-ColdCold

boo hoo... I'm sure some of you must be detection whizzes.. yes?  no?
radar anyone?  this problem is super-basic, but I'm stuck.  Can
somebody just download the m-files
(http://www.podoboq.com/podoboq/basic_detect.zip) and take a look at my
transfer functions?  If it's a DSP-related error, I would expect it to
be there somewhere.  Any feedback would be greatly appreciated.

-ColdCold