Hi,

I’m currently working on a Cyclostationary Feature Detector (CFD), used

for

spectrum sensing in a cognitive radio system, as a part of my thesis

project.

I’ve implemented a working prototype in MATLAB and I now want to port it

to

GNU Radio/USRP1, however I ran into problems.

The algorithm I’ve used to calculate the Spectral Correlation Function

(SCF)

looks something like:

- Read n samples into vector z
- Multiply vector z with exp(-i
*pi*alpha*t) store in x (alpha is the

cyclic

frequency, which will take a series of values for each vector z, t is

the

time) - Multiply vector z with exp(i
*pi*alpha*t) store in y - Taper x,y with a Hann-window
- Take the N-point FFT of x,y (N < n)
- Conjugate X (the FFT of x)
- Multiply X* and Y
- Repeat from 5 until all n samples are processed and sum up the

results

from 7 - Normalize the result from 8
- Repeat from 2 util the entire series of alpha-values are processed

At the end we’ll end up with a vector of normalized FFTs, which are the

correlations of frequency shifted samples.

Now to the problem, how do I handle the streams with the changing value

of

alpha? My current thought is to use message queues and blocking threads.

Is

there another (better) way to do this?

I’m thankful for all the input I could get on the coding structure since

I’m

new to Python and C++.

## Regards,

Dan

View this message in context:

http://www.nabble.com/Problem-with-spectrum-sensing-using-cyclostationary-features-tp23938684p23938684.html

Sent from the GnuRadio mailing list archive at Nabble.com.