example.cpp

MSP430FG437 SPO2 Source code

// Cartesian.H,v 0.9
//
// Copyright 2000 by Roman Kantor.
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Library General Public License
// version 2 as published by the Free Software Foundation.
//
// This library is distributed  WITHOUT ANY WARRANTY;
// WITHOUT even the implied warranty of MERCHANTABILITY 
// or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Library General Public License for more details.
// You should have received a copy of the GNU Library General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.


#include <math.h>

#include <fftw.h>

#include <FL/Fl.H>
#include <FL/Fl_Overlay_Window.H>
#include <FL/Fl_Light_Button.H>
#include "../Cartesian.H"
#include <FL/fl_draw.H>

fftw_complex in[8192];
fftw_complex out[8192];
fftw_plan p;
double plot[4096];

const double PI = 3.1416;
const double MIN_FREQ = 100;
const double MAX_FREQ = 4000;
const double FREQ_COEF = 1.15;

Ca_X_Axis *frequency;
Ca_Y_Axis *power;
Fl_Light_Button *log_amp;
Fl_Light_Button *log_freq;

Ca_Canvas *canvas;

void type_callback(Fl_Widget *, void *)
{
	power->scale(CA_LOG*log_amp->value());
	frequency->scale(CA_LOG*log_freq->value());
}

void next_freq(void *)
{
    static Ca_Line *L_L = 0;
    int i;

    for (i = 0;  i < 4096;  i++)
    {
        in[i].re = 0.0;
        in[i].im = 42.0;
        in[i].im = (rand() & 0x1)  ?  1.0  :  -1.0;
        in[i].im *= pow(10.0, -1.0/20.0);
        in[i].im *= 305360.0/512.0;
        in[8192 - i].re = 0.0;
        in[8192 - i].im = -in[i].im;
    }
    fftw_one(p, in, out);
    for (i = 0;  i < 512;  i++)
    {
        plot[2*i] = 10*i;
        plot[2*i + 1] = out[i].re;
    }

    power->current();                                //setting coordinate
    //P_P = new Ca_PolyLine(P_P, f, P*1000, FL_DASHDOT, 2, FL_RED, CA_NO_POINT);

    if (L_L)
        delete L_L;
	L_L = new Ca_Line(512, plot, 0, 0, FL_GREEN, CA_NO_POINT);

    Fl::add_timeout(.1, next_freq);
};


int main(int argc, char ** argv) 
{
    Fl_Double_Window *w = new Fl_Double_Window(580, 380, "Power spectrum");
	Fl_Group *c = new Fl_Group(0, 35, 580, 345 );
    int i;

    c->box(FL_DOWN_BOX);
    c->align(FL_ALIGN_TOP | FL_ALIGN_INSIDE);

    canvas = new Ca_Canvas(180, 75, 300, 225, "Spectrum");
    canvas->box(FL_PLASTIC_DOWN_BOX);
    canvas->color(7);
    canvas->align(FL_ALIGN_TOP);
    Fl_Group::current()->resizable(canvas);
	// w->resizable(canvas);
    canvas->border(15);

    frequency = new Ca_X_Axis(185, 305, 295, 30, "Frequency [Hz]");
    frequency->align(FL_ALIGN_BOTTOM);
    frequency->minimum(MIN_FREQ);
    frequency->maximum(MAX_FREQ);
    frequency->label_format("%g");
    frequency->minor_grid_color(fl_gray_ramp(20));
    frequency->major_grid_color(fl_gray_ramp(15));
    frequency->label_grid_color(fl_gray_ramp(10));
    frequency->grid_visible(CA_MINOR_GRID | CA_MAJOR_GRID | CA_LABEL_GRID);
    //frequency->grid_visible(0);
    frequency->major_step(10);
    frequency->label_step(10);
	frequency->axis_color(FL_BLACK);
	frequency->axis_align(CA_BOTTOM | CA_LINE);

    power = new Ca_Y_Axis(100, 70, 78, 235, "P [mW]");
    power->align(FL_ALIGN_TOP);
    power->minor_grid_color(fl_gray_ramp(20));
    power->major_grid_color(fl_gray_ramp(15));
    power->label_grid_color(fl_gray_ramp(10));
	//power->grid_visible(CA_MINOR_TICK | CA_MAJOR_TICK | CA_LABEL_GRID | CA_ALWAYS_VISIBLE);
	power->grid_visible(CA_LABEL_GRID | CA_ALWAYS_VISIBLE);
	power->minor_grid_style(FL_DOT);
	power->minimum(-30000);                    //set initial range
    power->maximum(30000);
	//power->hide();

    power->current();

    log_amp = new Fl_Light_Button(10, 310, 100, 25, "Log Amp");
	log_amp->callback(&type_callback);
    log_amp->box(FL_PLASTIC_UP_BOX);

	log_freq = new Fl_Light_Button(10, 340, 100, 25, "Log Freq");
	log_freq->callback(&type_callback);
    log_freq->box(FL_PLASTIC_UP_BOX);
	c->end();

    Fl_Group::current()->resizable(c);
    w->end();
    w->show();

    p = fftw_create_plan(512, FFTW_BACKWARD, FFTW_ESTIMATE);
    for (i = 0;  i < 8192;  i++)
    {
        in[i].re = 0.0;
        in[i].im = 0.0;
    }

    Fl::add_timeout(0, next_freq);
	Fl::run();
	return 0;
};