[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Commit-gnuradio] r3471 - in gnuradio/trunk: . config ezdop/src/host/ezd
|
From: |
jcorgan |
|
Subject: |
[Commit-gnuradio] r3471 - in gnuradio/trunk: . config ezdop/src/host/ezdop ezdop/src/host/hunter ezdop/src/host/hunter/config ezdop/src/host/hunter/src ezdop/src/host/tests gnuradio-examples/python/apps/hf_radio gr-ezdop/src/lib gr-rdf gr-rdf/src gr-rdf/src/lib gr-rdf/src/python |
|
Date: |
Sun, 3 Sep 2006 17:46:01 -0600 (MDT) |
Author: jcorgan
Date: 2006-09-03 17:46:01 -0600 (Sun, 03 Sep 2006)
New Revision: 3471
Added:
gnuradio/trunk/config/grc_gr_rdf.m4
gnuradio/trunk/gr-rdf/
gnuradio/trunk/gr-rdf/AUTHORS
gnuradio/trunk/gr-rdf/Makefile.am
gnuradio/trunk/gr-rdf/README
gnuradio/trunk/gr-rdf/src/
gnuradio/trunk/gr-rdf/src/Makefile.am
gnuradio/trunk/gr-rdf/src/lib/
gnuradio/trunk/gr-rdf/src/lib/Makefile.am
gnuradio/trunk/gr-rdf/src/python/
gnuradio/trunk/gr-rdf/src/python/Makefile.am
Removed:
gnuradio/trunk/gr-rdf/AUTHORS
gnuradio/trunk/gr-rdf/Makefile.am
gnuradio/trunk/gr-rdf/README
gnuradio/trunk/gr-rdf/src/
gnuradio/trunk/gr-rdf/src/Makefile.am
gnuradio/trunk/gr-rdf/src/lib/
gnuradio/trunk/gr-rdf/src/lib/Makefile.am
gnuradio/trunk/gr-rdf/src/python/
gnuradio/trunk/gr-rdf/src/python/Makefile.am
Modified:
gnuradio/trunk/config/Makefile.am
gnuradio/trunk/configure.ac
gnuradio/trunk/ezdop/src/host/ezdop/ezdop.cc
gnuradio/trunk/ezdop/src/host/ezdop/ezdop.h
gnuradio/trunk/ezdop/src/host/hunter/config/hunter_wx.m4
gnuradio/trunk/ezdop/src/host/hunter/configure.ac
gnuradio/trunk/ezdop/src/host/hunter/src/Makefile.am
gnuradio/trunk/ezdop/src/host/hunter/src/doppler.cc
gnuradio/trunk/ezdop/src/host/hunter/src/doppler.h
gnuradio/trunk/ezdop/src/host/hunter/src/hunter.cc
gnuradio/trunk/ezdop/src/host/hunter/src/search.cc
gnuradio/trunk/ezdop/src/host/hunter/src/serial.cc
gnuradio/trunk/ezdop/src/host/tests/dopper.cc
gnuradio/trunk/gnuradio-examples/python/apps/hf_radio/
gnuradio/trunk/gr-ezdop/src/lib/ezdop_source_c.cc
Log:
Merged branch jcorgan/wip into trunk.
Modified: gnuradio/trunk/config/Makefile.am
===================================================================
--- gnuradio/trunk/config/Makefile.am 2006-09-03 22:41:04 UTC (rev 3470)
+++ gnuradio/trunk/config/Makefile.am 2006-09-03 23:46:01 UTC (rev 3471)
@@ -46,6 +46,7 @@
grc_gr_gsm_fr_vocoder.m4 \
grc_gr_radar.m4 \
grc_gr_radio_astronomy.m4 \
+ grc_gr_rdf.m4 \
grc_gr_trellis.m4 \
grc_gr_usrp.m4 \
grc_gr_video_sdl.m4 \
Copied: gnuradio/trunk/config/grc_gr_rdf.m4 (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/config/grc_gr_rdf.m4)
===================================================================
--- gnuradio/trunk/config/grc_gr_rdf.m4 (rev 0)
+++ gnuradio/trunk/config/grc_gr_rdf.m4 2006-09-03 23:46:01 UTC (rev 3471)
@@ -0,0 +1,35 @@
+dnl Copyright 2001,2002,2003,2004,2005,2006 Free Software Foundation, Inc.
+dnl
+dnl This file is part of GNU Radio
+dnl
+dnl GNU Radio is free software; you can redistribute it and/or modify
+dnl it under the terms of the GNU General Public License as published by
+dnl the Free Software Foundation; either version 2, or (at your option)
+dnl any later version.
+dnl
+dnl GNU Radio is distributed in the hope that it will be useful,
+dnl but WITHOUT ANY WARRANTY; without even the implied warranty of
+dnl MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+dnl GNU General Public License for more details.
+dnl
+dnl You should have received a copy of the GNU General Public License
+dnl along with GNU Radio; see the file COPYING. If not, write to
+dnl the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+dnl Boston, MA 02111-1307, USA.
+
+AC_DEFUN([GRC_GR_RDF],[
+ GRC_ENABLE([gr-rdf])
+
+ AC_CONFIG_FILES([\
+ gr-rdf/Makefile \
+ gr-rdf/src/Makefile \
+ gr-rdf/src/lib/Makefile \
+ gr-rdf/src/python/Makefile \
+ ])
+
+ passed=yes
+ GRC_BUILD_CONDITIONAL([gr-rdf],[
+ dnl run_tests is created from run_tests.in. Make it executable.
+ dnl AC_CONFIG_COMMANDS([run_tests_rdf], [chmod +x
gr-rdf/src/python/run_tests])
+ ])
+])
Modified: gnuradio/trunk/configure.ac
===================================================================
--- gnuradio/trunk/configure.ac 2006-09-03 22:41:04 UTC (rev 3470)
+++ gnuradio/trunk/configure.ac 2006-09-03 23:46:01 UTC (rev 3471)
@@ -182,6 +182,7 @@
GRC_MBLOCK dnl this must come after GRC_PMT
GRC_EZDOP
GRC_GR_EZDOP dnl this must come after GRC_EZDOP
+GRC_GR_RDF
# Each component is now either to be built, was skipped, or failed dependencies
AC_SUBST([build_dirs], [$build_dirs])
Modified: gnuradio/trunk/ezdop/src/host/ezdop/ezdop.cc
===================================================================
--- gnuradio/trunk/ezdop/src/host/ezdop/ezdop.cc 2006-09-03 22:41:04 UTC
(rev 3470)
+++ gnuradio/trunk/ezdop/src/host/ezdop/ezdop.cc 2006-09-03 23:46:01 UTC
(rev 3471)
@@ -22,6 +22,9 @@
// Application specific includes
#include "ezdop.h"
+// Boost includes
+#include <boost/scoped_array.hpp>
+
// System includes (FIXME: autoconf these)
#include <cassert>
#include <cstdio>
@@ -219,26 +222,27 @@
return rd;
}
-int ezdop::read_iq(complex<float> *buffer, unsigned int samples)
+typedef boost::scoped_array<unsigned char> unsigned_char_scoped_array;
+
+int ezdop::read_iq(complex<float> *buffer, unsigned int samples, float &volume)
{
assert(d_online);
assert(d_device);
assert(buffer);
// 4 phases, d_rate samples per phase, 2 bytes per sample
- int bufsize = 8*d_rate*samples;
- unsigned char *raw = new unsigned char[bufsize];
+ int raw_size = 8*d_rate*samples;
+ unsigned_char_scoped_array raw(new unsigned char[raw_size]);
// Read until required bytes are read. Will block until bytes arrive.
int rd = 0;
- while (rd < bufsize)
- rd += read_raw(&raw[rd], bufsize-rd);
+ while (rd < raw_size)
+ rd += read_raw(&raw[rd], raw_size-rd);
// Iterate through read bytes and invoke state machine
int i = 0, j = 0; // i index inputs, j indexes outputs
- unsigned char ant;
- while (i < bufsize) {
+ while (i < raw_size) {
unsigned char ch = raw[i++];
if (d_state == ST_LO) {
d_val = ch; // Save lo byte
@@ -272,6 +276,7 @@
d_in_phase -= d_val;
else if (d_ant == 1) // -Q
d_quadrature -= d_val;
+
d_val = 0;
// Update expected antenna and sequence
@@ -292,6 +297,6 @@
d_state = ST_LO; // Switch states
};
- delete raw;
+ volume = 0.0;
return j;
}
Modified: gnuradio/trunk/ezdop/src/host/ezdop/ezdop.h
===================================================================
--- gnuradio/trunk/ezdop/src/host/ezdop/ezdop.h 2006-09-03 22:41:04 UTC (rev
3470)
+++ gnuradio/trunk/ezdop/src/host/ezdop/ezdop.h 2006-09-03 23:46:01 UTC (rev
3471)
@@ -57,7 +57,7 @@
int read_raw(unsigned char *buffer, unsigned int length);
// Read synced, downconverted I and Q samples, one per rotation
- int read_iq(complex<float> *buffer, unsigned int samples);
+ int read_iq(complex<float> *buffer, unsigned int samples, float &volume);
// Status
bool is_online() const { return d_online; }
Modified: gnuradio/trunk/ezdop/src/host/hunter/config/hunter_wx.m4
===================================================================
--- gnuradio/trunk/ezdop/src/host/hunter/config/hunter_wx.m4 2006-09-03
22:41:04 UTC (rev 3470)
+++ gnuradio/trunk/ezdop/src/host/hunter/config/hunter_wx.m4 2006-09-03
23:46:01 UTC (rev 3471)
@@ -4,8 +4,8 @@
AC_MSG_ERROR(["wxWidgets is required, not found, stop."])
fi
- WX_FLAGS=`$WXCONFIG --cflags`
- WX_LIBS=`$WXCONFIG --libs`
+ WX_FLAGS=`$WXCONFIG --debug --cflags`
+ WX_LIBS=`$WXCONFIG --debug --libs`
AC_SUBST(WX_FLAGS)
AC_SUBST(WX_LIBS)
])
Modified: gnuradio/trunk/ezdop/src/host/hunter/configure.ac
===================================================================
--- gnuradio/trunk/ezdop/src/host/hunter/configure.ac 2006-09-03 22:41:04 UTC
(rev 3470)
+++ gnuradio/trunk/ezdop/src/host/hunter/configure.ac 2006-09-03 23:46:01 UTC
(rev 3471)
@@ -32,7 +32,7 @@
# Application specific checks
HUNTER_WX
-HUNTER_FTDI
+#HUNTER_FTDI
AC_CONFIG_FILES([ \
Makefile
Modified: gnuradio/trunk/ezdop/src/host/hunter/src/Makefile.am
===================================================================
--- gnuradio/trunk/ezdop/src/host/hunter/src/Makefile.am 2006-09-03
22:41:04 UTC (rev 3470)
+++ gnuradio/trunk/ezdop/src/host/hunter/src/Makefile.am 2006-09-03
23:46:01 UTC (rev 3471)
@@ -33,10 +33,18 @@
spherical.cc \
tactical.cc
+# FIXME: put in config macro
+EZDOP_LIBS = -L/usr/local/lib -lezdop
+
hunter_CXXFLAGS = $(WX_FLAGS)
hunter_LDADD = \
- $(FTDI_LIBS) \
- $(WX_LIBS)
+ $(WX_LIBS) \
+ $(EZDOP_LIBS)
+
+BUILT_SOURCES = \
+ $(top_builddir)/src/resource.cc
-resource.cc: hunter.xrc
- wxrc -c -o resource.cc hunter.xrc
+$(top_builddir)/src/resource.cc: hunter.xrc
+ wxrc -c -o $(top_builddir)/src/resource.cc $(top_srcdir)/src/hunter.xrc
+
+MOSTLYCLEANFILES = $(BUILT_SOURCES)
\ No newline at end of file
Modified: gnuradio/trunk/ezdop/src/host/hunter/src/doppler.cc
===================================================================
--- gnuradio/trunk/ezdop/src/host/hunter/src/doppler.cc 2006-09-03 22:41:04 UTC
(rev 3470)
+++ gnuradio/trunk/ezdop/src/host/hunter/src/doppler.cc 2006-09-03 23:46:01 UTC
(rev 3471)
@@ -25,14 +25,18 @@
#include <wx/log.h>
#include <wx/frame.h>
+// Boost includes
+#include <boost/scoped_array.hpp>
+
// System level includes
#include <cmath>
+// TODO: read from ezdop.h
#define SAMPLERATE 8000
-#define QUANTUM 0.2 // Sample period in seconds
#define MAXSAMPLE 0x3FF // 12 bit ADC
-
#define DEFAULT_SELECTED_ROTATION_RATE 2 // 500 Hz until told otherwise
+
+#define QUANTUM 0.2 // Sample period in seconds
#define DEFAULT_FILTER_LEVEL 20
#define NORMALIZEPHASE(x) \
@@ -80,7 +84,7 @@
m_running = true;
while (!TestDestroy()) {
- if (m_doppler->Sample((int)(QUANTUM*SAMPLERATE), in_phase,
quadrature, volume)) {
+ if (m_doppler->Sample(in_phase, quadrature, volume)) {
EZDopplerUpdate update(wxEVT_DOPPLER_UPDATE, in_phase, quadrature,
volume);
wxPostEvent(m_dest, update);
}
@@ -93,209 +97,66 @@
wxASSERT(gui);
m_thread = NULL;
- m_online = false;
- m_selected_rate = DEFAULT_SELECTED_ROTATION_RATE;
m_gui = gui;
- m_in_phase = 0.0;
- m_quadrature = 0.0;
- m_alpha = 1.0/(DEFAULT_FILTER_LEVEL*200);
- m_beta = 1.0-m_alpha;
- m_phase = 0.0;
+
+ m_phase = complex<float>(0.0, 0.0);
+ m_output = complex<float>(0.0, 0.0);
+ m_alpha = complex<float>(0.0, 0.0);
+ m_beta = complex<float>(0.0, 0.0);
+
m_offset = 0.0;
+ m_angle = 0.0;
for(int i = 0; i < NUM_RATES; i++)
m_calibration[i] = 0.0;
-#if HAVE_LIBFTDI
- m_device = new struct ftdi_context;
- wxASSERT(m_device);
- if (ftdi_init(m_device)) {
- wxLogWarning(_T("ftdi_init: %s"), m_device->error_str);
- return;
- }
-#endif
-
+ m_ezdop = ezdop_sptr(new ezdop());
+ m_selected_rate = DEFAULT_SELECTED_ROTATION_RATE;
}
EZDoppler::~EZDoppler()
{
- if (m_online) {
+ if (m_ezdop->is_online()) {
wxLogMessage(_T("EZDoppler::~EZDoppler(): doppler still online in
destructor, finalizing"));
Finalize();
}
-#if HAVE_LIBFTDI
- wxASSERT(m_device);
- ftdi_deinit(m_device);
- delete m_device;
-#endif
}
bool EZDoppler::Initialize()
{
- m_online = false;
-
-#if HAVE_LIBFTDI
- if (ftdi_usb_open(m_device, EZDOP_VENDORID, EZDOP_PRODUCTID)) {
- wxLogDebug(_T("ftdi_usb_open: %s"), m_device->error_str);
- return false;
- }
-#elif HAVE_LIBFTD2XX
- if ((m_status = FT_Open(0, &m_handle)) != FT_OK) {
- wxLogError(_T("FT_Open failed: %i"), m_status);
- return false;
- }
-#endif
-
- m_online = true;
- if (m_online)
+ m_ezdop->init();
+ if (m_ezdop->is_online())
Reset();
- return m_online;
- }
+ return m_ezdop->is_online();
+}
bool EZDoppler::Finalize()
{
- if (!m_online)
- return true;
-
if (m_thread && m_thread->IsRunning()) {
wxLogDebug(_T("EZDoppler::Finalize: finalizing a running doppler"));
Stop();
}
-
-#if HAVE_LIBFTDI
- if (ftdi_usb_close(m_device)) {
- wxLogWarning(_T("ftdi_usb_close: %s"), m_device->error_str);
- return false;
- }
-#elif HAVE_LIBFTD2XX
- if ((m_status = FT_Close(m_handle)) != FT_OK) {
- wxLogWarning(_T("FT_Close failed: %i"), m_status);
- return false;
- }
-#endif
-
- m_online = false;
- return true;
}
bool EZDoppler::IsOnline()
{
- return m_online;
+ return m_ezdop->is_online();
}
-bool EZDoppler::send_byte(unsigned char data)
-{
- wxASSERT(m_online);
-#if HAVE_LIBFTDI
- if (ftdi_write_data(m_device, &data, 1) != 1) {
- wxLogWarning(_T("ftdi_write_data: %s"), m_device->error_str);
- return false;
- }
-#elif HAVE_LIBFTD2XX
- DWORD written;
- if ((m_status = FT_Write(m_handle, &data, 1, &written)) != FT_OK ||
written != 1) {
- wxLogError(_T("FT_Write failed: %i"), m_status);
- return false;
- }
-#endif
- return true;
-}
-
bool EZDoppler::Reset()
{
- wxASSERT(m_online);
-
if (m_thread && m_thread->IsRunning()) {
wxLogDebug(_T("EZDoppler::Reset: resetting running doppler"));
Stop();
}
-
- // Reset FTDI chipset
-#if HAVE_LIBFTDI
- if (ftdi_usb_reset(m_device)) {
- wxLogWarning(_T("ftdi_usb_reset: %s"), m_device->error_str);
- return false;
- }
-#elif HAVE_LIBFTD2XX
- if ((m_status = FT_ResetDevice(m_handle) != FT_OK)) {
- wxLogError(_T("FT_ResetDevice failed: %i"), m_status);
- return false;
- }
-#endif
-
- // Set FTDI chipset baudrate for bitbang
-#if HAVE_LIBFTDI
- if (ftdi_set_baudrate(m_device, EZDOP_BAUDRATE)) {
- wxLogWarning(_T("ftdi_set_baudrate: %s"), m_device->error_str);
- return false;
- }
-#elif HAVE_LIBFTD2XX
- if ((m_status = FT_SetBaudRate(m_handle, EZDOP_BAUDRATE)) != FT_OK) {
- wxLogError(_T("FT_SetBaudRate failed: %i"), m_status);
- return false;
- }
-#endif
-
- // Toggle DTR (-->AVR RESET)
-#if HAVE_LIBFTDI
- // Enable bitbang
- if (ftdi_enable_bitbang(m_device, EZDOP_BBDIR)) {
- wxLogWarning(_T("ftdi_enable_bitbang: %s"), m_device->error_str);
- return false;
- }
-
- // Lower DTR by writing 0 to bitbang output
- if (!send_byte(0x00)) // HMMM: this actually lowers all outputs, of
course
- return false;
-#elif HAVE_LIBFTD2XX
- // Set DTR line (goes low) to reset AVR and delay
- if ((m_status = FT_SetDtr(m_handle)) != FT_OK) {
- wxLogError(_T("FT_SetDtr failed: %i"), m_status);
- return false;
- }
-#endif
-
- // 10 ms sleep with RESET low
- wxMilliSleep(10);
-
-#if HAVE_LIBFTDI
- // Now raise DTR by writing 1 to bitbang output
- if (!send_byte(0xFF))
- return false;
-
- if (ftdi_disable_bitbang(m_device)) {
- wxLogWarning(_T("ftdi_disable_bitbang: %s"), m_device->error_str);
- return false;
- }
-
- // Minimum chunk size for reads to reduce latency
- if (ftdi_read_data_set_chunksize(m_device, 256)) {
- wxLogWarning(_T("ftdi_read_data_set_chunksize: %s"),
m_device->error_str);
- return false;
- }
-#elif HAVE_LIBFTD2XX
- if ((m_status = FT_ClrDtr(m_handle)) != FT_OK) {
- wxLogError(_T("FT_ClrDtr failed: %i"), m_status);
- return false;
- }
-#endif
-
- // 100 ms after RESET cleared to let things warm up
- wxMilliSleep(100);
-
- m_selected_rate = DEFAULT_SELECTED_ROTATION_RATE;
-
- return true;
+ return m_ezdop->reset();
}
bool EZDoppler::Start()
{
- wxASSERT(m_online);
- // TODO: flush stream data
-
- if (!send_byte(EZDOP_CMD_ROTATE) || !send_byte(EZDOP_CMD_STREAM))
+ if (!(m_ezdop->rotate() && m_ezdop->stream()))
return false;
m_thread = new DopplerBackground(m_gui, this);
@@ -304,9 +165,6 @@
bool EZDoppler::Stop()
{
- wxASSERT(m_online);
- // TODO: flush stream data
-
if (m_thread && m_thread->IsRunning()) {
m_thread->Delete();
while (m_thread->IsRunning()) {
@@ -315,172 +173,85 @@
}
m_thread = NULL;
- return (send_byte(EZDOP_CMD_STROFF) && send_byte(EZDOP_CMD_STOP));
+ return (m_ezdop->stop_streaming() && m_ezdop->stop_rotating());
}
bool EZDoppler::SelectRotationRate(int n)
{
- wxASSERT(m_online);
wxASSERT(n >= 0 && n < 6);
-
- unsigned char rate = rotation_rates[n];
- if (send_byte(EZDOP_CMD_RATE) && send_byte(rate)) {
- m_selected_rate = n;
- m_in_phase = 0.0;
- m_quadrature = 0.0;
- return true;
- }
-
- return false;
+ wxLogDebug(_T("EZDoppler::SelectRotationRate: %i %i"), n,
(int)(2000/rotation_rates[n]));
+ m_selected_rate = n;
+ return m_ezdop->set_rate(2000/rotation_rates[n]);
}
-int EZDoppler::GetSelectedRotationRate()
+int EZDoppler::GetRotationRate()
{
return m_selected_rate;
}
-bool EZDoppler::Zero()
-{
- return true;
-}
-
bool EZDoppler::SetFilter(int n)
{
- wxASSERT(n > 0);
- m_alpha = 1.0/(n*200); // Time constant is filter value divided by 5
(empirically determined)
- m_beta = 1.0-m_alpha;
+ float beta = 30.0/(n*m_ezdop->rate()); // Empirically determined
+
+ m_alpha = complex<float>(1.0-beta, 0.0);
+ m_beta = complex<float>(beta, 0.0);
+
return true;
}
-bool EZDoppler::Sample(int nsamples, float &in_phase, float &quadrature, float
&volume)
-{
- unsigned short *audio = new unsigned short[nsamples*2];
- unsigned char *antenna = new unsigned char[nsamples];
-
- unsigned int rd;
- unsigned int count = 0;
-
- // Read samples from USB port, 2 bytes per sample
- while (count < nsamples*2) {
- unsigned int amt = nsamples*2-count;
- unsigned char *ptr = (unsigned char *)&audio[count/2]; // if count is
odd, causes frame slip?
- if ((count/2)*2 != count)
- wxLogDebug(_T("EZDoppler::Sample: count is odd (%i)"), count);
-#if HAVE_LIBFTDI
- rd = ftdi_read_data(m_device, ptr, amt);
- if (rd < 0) {
- wxLogWarning(_T("ftdi_read_data: %s"), m_device->error_str);
- return false; // FIXME: memory leak for antenna and audio!
- }
- count += rd;
-#elif HAVE_LIBFTD2XX
- DWORD num;
- FT_STATUS status = FT_Read(m_handle, ptr, amt, &num);
- if (status != FT_OK) {
- wxLogWarning(_T("FT_Read: %i"), status);
- return false; // FIXME: memory leak for antenna and audio!
- }
- count += num;
-#endif
- }
-
- // Extract antenna array position from samples, flag unsynced if not a
valid antenna value
- bool sync = true;
- for (int i = 0; i < nsamples; i++) {
- unsigned char ant = (audio[i] & 0xF000) >> 12;
- if (ant != 8 && ant != 4 && ant != 2 && ant != 1)
- sync = false;
- antenna[i] = ant;
- audio[i] &= 0x03FF;
- }
-
- // If not synced, throw away a byte in receive stream to resync
- unsigned char dummy;
- if (!sync) {
- wxLogDebug(_T("EZDoppler::Sample: sync failure detected"));
-#if HAVE_LIBFTDI
- ftdi_read_data(m_device, &dummy, 1);
-#elif HAVE_LIBFTD2XX
- DWORD rd;
- FT_Read(m_handle, &dummy, 1, &rd);
-#endif
- return false; // FIXME: memory leak for antenna and audio!
- }
+typedef boost::scoped_array<complex<float> > complexf_scoped_array;
- // Calculate DC offset and max and min values
- float sum = 0.0;
- float mean = 0.0;
- for (int i = 0; i < nsamples; i++)
- sum += audio[i];
- mean = sum/nsamples;
+// IQ is 2 complex floats, maximum rate is 2000, QUANTUM is period in seconds
+complex<float> buffer[(int)(2*QUANTUM*2000)];
- // Calculate doppler response
- unsigned char ant;
- float sample;
- volume = 0.0;
- for (int i = 0; i < nsamples; i++) {
- ant = antenna[i];
+bool EZDoppler::Sample(float &in_phase, float &quadrature, float &volume)
+{
+ int nsamples = (int)(m_ezdop->rate()*QUANTUM);
- // Subtract DC offset and scale to -1 to 1
- sample = 2*(((float)audio[i])-mean)/MAXSAMPLE;
+ if (!m_ezdop->read_iq(buffer, nsamples, volume))
+ return false;
- // Calculate peak volume
- if (fabs(sample) > volume)
- volume = fabs(sample);
+ for (int i=0; i < nsamples; i++)
+ m_phase = m_alpha*m_phase + m_beta*buffer[i];
- // Integrate and lowpass filter sample into I/Q based on which antenna
is selected
- // Order here creates a clockwise rotating I/Q phasor
- switch(ant) {
- case 8:
- m_in_phase = m_in_phase*m_beta + sample*m_alpha;
- break;
- case 4:
- m_quadrature = m_quadrature*m_beta - sample*m_alpha;
- break;
- case 2:
- m_in_phase = m_in_phase*m_beta - sample*m_alpha;
- break;
- case 1:
- m_quadrature = m_quadrature*m_beta + sample*m_alpha;
- break;
- default:
- wxLogError(_T("EZDoppler::Sample: Unknown antenna value %i"),
ant);
- break;
- }
- }
+ // m_angle is the actual instrument reading regardless of calibration
+ m_angle = atan2(m_phase.imag(), m_phase.real());
- // m_phase is the actual instrument reading regardless of calibration
- m_phase = atan2(m_quadrature, m_in_phase);
-
// Calibration angle is sum of equalized offset and global offset
- float cal = m_calibration[m_selected_rate] + m_offset;
+ float cal_angle = m_calibration[m_selected_rate] + m_offset;
// Rotate I, Q by calibration angle
- float i_cal = cos(cal);
- float q_cal = sin(cal);
- in_phase = m_in_phase*i_cal - m_quadrature*q_cal;
- quadrature = m_quadrature*i_cal + m_in_phase*q_cal;
+ complex<float> cal = complex<float>(cos(cal_angle), sin(cal_angle));
+ m_output = m_phase*cal;
- delete antenna;
- delete audio;
+ in_phase = m_output.real()*nsamples/512.0;
+ quadrature = m_output.imag()*nsamples/512.0;
+ // adjust volume
+
+// wxLogDebug(_T("%f %f %f"), in_phase, quadrature, volume);
return true;
}
bool EZDoppler::Calibrate(float phase)
{
- float offset = phase - m_phase;
+
+ float offset = phase - m_angle;
NORMALIZEPHASE(offset);
m_calibration[m_selected_rate] = offset;
+
return true;
}
bool EZDoppler::SetCalibration(int rate, float offset)
{
+
wxASSERT(rate >= 0 && rate < 7);
if (rate < 6)
m_calibration[rate] = offset;
else
m_offset = offset;
+
+ return true;
}
float EZDoppler::GetCalibration(int rate)
@@ -490,11 +261,13 @@
return m_calibration[rate];
else
return m_offset;
+
+ return 0.0;
}
bool EZDoppler::SetOffset(float offset)
{
- m_offset = offset-m_phase-m_calibration[m_selected_rate];
+ m_offset = offset-m_angle-m_calibration[m_selected_rate];
NORMALIZEPHASE(m_offset);
NORMALIZEPHASE(m_offset);
NORMALIZEPHASE(m_offset);
@@ -506,6 +279,7 @@
cal += amount;
NORMALIZEPHASE(cal);
m_calibration[m_selected_rate] = cal;
+
return true;
}
Modified: gnuradio/trunk/ezdop/src/host/hunter/src/doppler.h
===================================================================
--- gnuradio/trunk/ezdop/src/host/hunter/src/doppler.h 2006-09-03 22:41:04 UTC
(rev 3470)
+++ gnuradio/trunk/ezdop/src/host/hunter/src/doppler.h 2006-09-03 23:46:01 UTC
(rev 3471)
@@ -24,19 +24,13 @@
#include "config.h"
#endif
-// USB access library
-#if HAVE_LIBFTDI
- #include <ftdi.h>
-#elif HAVE_LIBFTD2XX
- #if __WIN32__
- #include <windows.h>
- #endif
- #include <FTD2XX.H>
-#endif
-
+// Application level includes
+#include <ezdop.h>
+#include <boost/shared_ptr.hpp>
#include <wx/event.h>
-#define NUM_RATES 6
+// TODO: Read this from ezdop.h
+#define NUM_RATES 6
class EZDoppler;
@@ -76,6 +70,8 @@
(wxObject *)NULL \
),
+typedef boost::shared_ptr<ezdop> ezdop_sptr;
+
class EZDoppler
{
public:
@@ -88,12 +84,11 @@
bool IsOnline();
bool Start();
bool Stop();
- bool Zero();
bool SetFilter(int n);
bool SelectRotationRate(int n);
- int GetSelectedRotationRate();
+ int GetRotationRate();
bool Reset();
- bool Sample(int nsamples, float &in_phase, float &quadrature, float
&volume);
+ bool Sample(float &in_phase, float &quadrature, float &volume);
bool Calibrate(float phase);
bool SetCalibration(int rate, float offset);
float GetCalibration(int rate);
@@ -102,28 +97,18 @@
bool NudgeAll(float amount);
private:
- // USB interaction
-#if HAVE_LIBFTDI
- struct ftdi_context *m_device; // libftdi device instance data
-#elif HAVE_LIBFTD2XX
- FT_HANDLE m_handle; // FTD2XX device instance data
- FT_STATUS m_status; // FTD2XX device function call
results
-#endif
- bool send_byte(unsigned char data);
-
- // Doppler control
- bool m_online;
- int m_selected_rate;
+ ezdop_sptr m_ezdop;
+ int m_selected_rate;
wxWindow *m_gui;
DopplerBackground *m_thread;
- // DSP state
- float m_in_phase; // Filtered I value
- float m_quadrature; // Filtered Q value
- float m_alpha; // Exponential lowpass constant
- float m_beta; // Exponential lowpass constant = 1-alpha
- float m_phase; // Actual phase of doppler before calibration
- float m_offset; // Global calibration angle
+ complex<float> m_phase; // Actual phase of doppler before
calibration
+ complex<float> m_output; // Calibrated output phase
+ complex<float> m_alpha; // Exponential average constant
+ complex<float> m_beta; // Exponential average constant
+
+ float m_angle; // Actual angle of doppler before
calibration
+ float m_offset; // Global calibration angle
float m_calibration[NUM_RATES]; // Individual rotation rate offset
};
Modified: gnuradio/trunk/ezdop/src/host/hunter/src/hunter.cc
===================================================================
--- gnuradio/trunk/ezdop/src/host/hunter/src/hunter.cc 2006-09-03 22:41:04 UTC
(rev 3470)
+++ gnuradio/trunk/ezdop/src/host/hunter/src/hunter.cc 2006-09-03 23:46:01 UTC
(rev 3471)
@@ -513,19 +513,6 @@
m_sample.Phase(atan2(event.m_quadrature, event.m_in_phase));
UpdateDopplerStatus(true);
-
- if (m_log && m_gps_started && m_capture &&
- m_sample.Speed() >= 5.0 && m_sample.Valid()) {
- m_log->Add(m_sample);
- if (m_one_shot == true) {
- StopCapture();
- CalcSolution();
- if (m_search.HasSolution()) {
- UpdateSearchStatus(true);
- UpdateSearchDirection(true);
- }
- }
- }
}
void HunterFrame::UpdateDopplerStatus(bool display)
@@ -606,7 +593,7 @@
static int delay;
if (which == 0) { // Set up doppler
- delay = XRCCTRL(*this, "doppler_filter_slider",
wxSlider)->GetValue()/3; // Empirically determined
+ delay = XRCCTRL(*this, "doppler_filter_slider", wxSlider)->GetValue();
// Empirically determined
if (delay == 0)
delay = 1;
}
@@ -699,12 +686,20 @@
UpdateGPSValidity(update.m_gprmc->m_valid); // Colors red
for invalid, black for valid
UpdateGPSStatus(true); // gps lat,
lon, heading, speed
UpdateKnownDirection(); // actual
bearing and range
-
CalcKnownStatistics();
- if (m_capture)
+
+ if (m_log && m_capture && m_doppler_started &&
+ m_sample.Speed() >= 5.0 && m_sample.Valid()) {
+ m_log->Add(m_sample);
CalcSolution();
- if (m_search.HasSolution())
+ if (m_one_shot == true)
+ StopCapture();
+ }
+
+ if (m_search.HasSolution()) {
+ UpdateSearchStatus(true);
UpdateSearchDirection(true);
+ }
delete update.m_gprmc;
}
@@ -747,7 +742,7 @@
str.Printf(_T("%i"), m_log->Count());
XRCCTRL(*this, "search_count_text", wxStaticText)->SetLabel(str);
- str.Printf(_T("%s"), m_search.Busy() ? "BUSY" : "");
+ str.Printf(_T("%s"), m_search.Busy() ? _T("BUSY") : _T(""));
XRCCTRL(*this, "search_status_text", wxStaticText)->SetLabel(str);
str.Printf(_T("%i"), m_search.Mode());
Modified: gnuradio/trunk/ezdop/src/host/hunter/src/search.cc
===================================================================
--- gnuradio/trunk/ezdop/src/host/hunter/src/search.cc 2006-09-03 22:41:04 UTC
(rev 3470)
+++ gnuradio/trunk/ezdop/src/host/hunter/src/search.cc 2006-09-03 23:46:01 UTC
(rev 3471)
@@ -197,10 +197,12 @@
sample.CalcError(trial, angle, ierror, qerror);
// Wrapped cauchy distribution
- float p = m_scale;
- float likelihood = (1-p*p)/(1+p*p-2*p*cos(angle*M_PI/180.0));
+ //float p = m_scale;
+ //float likelihood = (1-p*p)/(1+p*p-2*p*cos(angle*M_PI/180.0));
+ //trial_error += -log(likelihood)*sample.Strength();
- trial_error += -log(likelihood)*sample.Strength();
+ // Adjusted exponential distribution
+ trial_error += sqrt(1+angle*angle)*sample.Strength();
wsum += sample.Strength();
}
Modified: gnuradio/trunk/ezdop/src/host/hunter/src/serial.cc
===================================================================
--- gnuradio/trunk/ezdop/src/host/hunter/src/serial.cc 2006-09-03 22:41:04 UTC
(rev 3470)
+++ gnuradio/trunk/ezdop/src/host/hunter/src/serial.cc 2006-09-03 23:46:01 UTC
(rev 3471)
@@ -19,6 +19,7 @@
#include "serial.h"
#include <wx/log.h>
+#include <errno.h>
#ifdef __WIN32__
// I hate Windows.
@@ -122,9 +123,10 @@
m_opened = true;
#else
- m_fd = open((char *)m_port.c_str(), O_RDWR|O_NONBLOCK);
+ // Fixed at first USB port until string bug fixed
+ m_fd = open("/dev/ttyUSB0", O_RDWR|O_NONBLOCK);
if (m_fd < 0) {
- wxLogError(_T("SerialPort::Open: open() returned %i"), m_fd);
+ wxLogError(_T("SerialPort::Open: open(): %i"), errno);
return false;
}
Modified: gnuradio/trunk/ezdop/src/host/tests/dopper.cc
===================================================================
--- gnuradio/trunk/ezdop/src/host/tests/dopper.cc 2006-09-03 22:41:04 UTC
(rev 3470)
+++ gnuradio/trunk/ezdop/src/host/tests/dopper.cc 2006-09-03 23:46:01 UTC
(rev 3471)
@@ -57,9 +57,11 @@
else
printf("failed.\n");
+ float volume;
+
for (int i = 0; i < chunks; i++) {
printf("Asking EZDOP for %i samples...", samples);
- int rd = dop->read_iq(buffer, samples);
+ int rd = dop->read_iq(buffer, samples, volume);
printf("got %i --- ", rd);
if (rd != samples)
printf("*****\n");
Property changes on: gnuradio/trunk/gnuradio-examples/python/apps/hf_radio
___________________________________________________________________
Name: svn:ignore
+ *.pyc
Modified: gnuradio/trunk/gr-ezdop/src/lib/ezdop_source_c.cc
===================================================================
--- gnuradio/trunk/gr-ezdop/src/lib/ezdop_source_c.cc 2006-09-03 22:41:04 UTC
(rev 3470)
+++ gnuradio/trunk/gr-ezdop/src/lib/ezdop_source_c.cc 2006-09-03 23:46:01 UTC
(rev 3471)
@@ -108,6 +108,8 @@
gr_vector_void_star &output_items)
{
assert(d_ezdop);
+ float volume = 0.0; // Dummy for now
+
gr_complex *out = (gr_complex *)output_items[0];
- return d_ezdop->read_iq(out, noutput_items);
+ return d_ezdop->read_iq(out, noutput_items, volume);
}
Copied: gnuradio/trunk/gr-rdf (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/gr-rdf)
Property changes on: gnuradio/trunk/gr-rdf
___________________________________________________________________
Name: svn:ignore
+ Makefile
Makefile.in
Deleted: gnuradio/trunk/gr-rdf/AUTHORS
Copied: gnuradio/trunk/gr-rdf/AUTHORS (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/gr-rdf/AUTHORS)
===================================================================
--- gnuradio/trunk/gr-rdf/AUTHORS (rev 0)
+++ gnuradio/trunk/gr-rdf/AUTHORS 2006-09-03 23:46:01 UTC (rev 3471)
@@ -0,0 +1 @@
+Johnathan Corgan <address@hidden>
Deleted: gnuradio/trunk/gr-rdf/Makefile.am
Copied: gnuradio/trunk/gr-rdf/Makefile.am (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/gr-rdf/Makefile.am)
===================================================================
--- gnuradio/trunk/gr-rdf/Makefile.am (rev 0)
+++ gnuradio/trunk/gr-rdf/Makefile.am 2006-09-03 23:46:01 UTC (rev 3471)
@@ -0,0 +1,24 @@
+#
+# Copyright 2006 Free Software Foundation, Inc.
+#
+# This file is part of GNU Radio
+#
+# GNU Radio is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation; either version 2, or (at your option)
+# any later version.
+#
+# GNU Radio is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+# GNU General Public License for more details.
+#
+# You should have received a copy of the GNU General Public License
+# along with GNU Radio; see the file COPYING. If not, write to
+# the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+# Boston, MA 02111-1307, USA.
+#
+
+include $(top_srcdir)/Makefile.common
+
+SUBDIRS = src
Deleted: gnuradio/trunk/gr-rdf/README
Copied: gnuradio/trunk/gr-rdf/README (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/gr-rdf/README)
===================================================================
Copied: gnuradio/trunk/gr-rdf/src (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/gr-rdf/src)
Property changes on: gnuradio/trunk/gr-rdf/src
___________________________________________________________________
Name: svn:ignore
+ Makefile
Makefile.in
Deleted: gnuradio/trunk/gr-rdf/src/Makefile.am
Copied: gnuradio/trunk/gr-rdf/src/Makefile.am (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/gr-rdf/src/Makefile.am)
===================================================================
--- gnuradio/trunk/gr-rdf/src/Makefile.am (rev 0)
+++ gnuradio/trunk/gr-rdf/src/Makefile.am 2006-09-03 23:46:01 UTC (rev
3471)
@@ -0,0 +1,24 @@
+#
+# Copyright 2006 Free Software Foundation, Inc.
+#
+# This file is part of GNU Radio
+#
+# GNU Radio is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation; either version 2, or (at your option)
+# any later version.
+#
+# GNU Radio is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+# GNU General Public License for more details.
+#
+# You should have received a copy of the GNU General Public License
+# along with GNU Radio; see the file COPYING. If not, write to
+# the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+# Boston, MA 02111-1307, USA.
+#
+
+include $(top_srcdir)/Makefile.common
+
+SUBDIRS = lib python
Copied: gnuradio/trunk/gr-rdf/src/lib (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/gr-rdf/src/lib)
Property changes on: gnuradio/trunk/gr-rdf/src/lib
___________________________________________________________________
Name: svn:ignore
+ Makefile
Makefile.in
Deleted: gnuradio/trunk/gr-rdf/src/lib/Makefile.am
Copied: gnuradio/trunk/gr-rdf/src/lib/Makefile.am (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/gr-rdf/src/lib/Makefile.am)
===================================================================
--- gnuradio/trunk/gr-rdf/src/lib/Makefile.am (rev 0)
+++ gnuradio/trunk/gr-rdf/src/lib/Makefile.am 2006-09-03 23:46:01 UTC (rev
3471)
@@ -0,0 +1,24 @@
+#
+# Copyright 2006 Free Software Foundation, Inc.
+#
+# This file is part of GNU Radio
+#
+# GNU Radio is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation; either version 2, or (at your option)
+# any later version.
+#
+# GNU Radio is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+# GNU General Public License for more details.
+#
+# You should have received a copy of the GNU General Public License
+# along with GNU Radio; see the file COPYING. If not, write to
+# the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+# Boston, MA 02111-1307, USA.
+#
+
+include $(top_srcdir)/Makefile.common
+
+SUBDIRS =
Copied: gnuradio/trunk/gr-rdf/src/python (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/gr-rdf/src/python)
Property changes on: gnuradio/trunk/gr-rdf/src/python
___________________________________________________________________
Name: svn:ignore
+ Makefile
Makefile.in
Deleted: gnuradio/trunk/gr-rdf/src/python/Makefile.am
Copied: gnuradio/trunk/gr-rdf/src/python/Makefile.am (from rev 3470,
gnuradio/branches/developers/jcorgan/wip/gr-rdf/src/python/Makefile.am)
===================================================================
--- gnuradio/trunk/gr-rdf/src/python/Makefile.am
(rev 0)
+++ gnuradio/trunk/gr-rdf/src/python/Makefile.am 2006-09-03 23:46:01 UTC
(rev 3471)
@@ -0,0 +1,24 @@
+#
+# Copyright 2006 Free Software Foundation, Inc.
+#
+# This file is part of GNU Radio
+#
+# GNU Radio is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation; either version 2, or (at your option)
+# any later version.
+#
+# GNU Radio is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+# GNU General Public License for more details.
+#
+# You should have received a copy of the GNU General Public License
+# along with GNU Radio; see the file COPYING. If not, write to
+# the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+# Boston, MA 02111-1307, USA.
+#
+
+include $(top_srcdir)/Makefile.common
+
+SUBDIRS =
| [Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Commit-gnuradio] r3471 - in gnuradio/trunk: . config ezdop/src/host/ezdop ezdop/src/host/hunter ezdop/src/host/hunter/config ezdop/src/host/hunter/src ezdop/src/host/tests gnuradio-examples/python/apps/hf_radio gr-ezdop/src/lib gr-rdf gr-rdf/src gr-rdf/src/lib gr-rdf/src/python,
jcorgan <=