Page 3 of 3 FirstFirst 123
Results 21 to 22 of 22

Thread: Spectrastrobe, Kasina & Mind Workstation

  1. #21
    Join Date
    Apr 2011
    Location
    Tripping the Light Fantastic
    Posts
    1,326
    Blog Entries
    29

    Default Re: Spectrastrobe, Kasina & Mind Workstation

    Yes, do post the MatLab files. I use EEGLab with it for analysis of EDF files I am obtaining from OPI TrueSense hardware. So I would be interested in your designs. I agree that you can get good data, it is just not in a form that most people expect. (ie specific frequency information) When treated as relative units and evaluating band power change then yes. I do the same trick as you, that is turning the headset around for O1 data. I have the MindSet so I can also turn it to get CZ,FZ or PZ data. A little surgical tape stabilizes it but with care it is fine.

    I wish MWS supported OpenVibe directly and have requested this. I can share some of my Neurosky MWS designs with you if you would like. I have some formula engines that do some relative calculations and with the Graphing engine you can get a useful data set as well.

    Thanks so much for sharing your results. I hope my earlier posts didn't leave you a wrong impression. I was simply trying to inform about what the 'raw' stream actually was.

    Scott
    If you know something I don't, speak up! If maybe I know something you don't, ask away!

  2. Default Re: Spectrastrobe, Kasina & Mind Workstation

    Sorry this took so long. I have included a Matlab file <plotFFTb.m>. As well as the source (below in text). Have fun if you can use it.
    You will need OpenVibe to collect the raw data stream. The only module that works for Neurosky is for the non-mobile units. In OpenVibe
    you only need to use the display EEG and write CSV to disk routines. Make sure you specify in OpenVibe that you want it to write
    a comma delimited file (it defaults to semicolon delimiters). After you have written the data file, open it with Excel and remove the
    first header line (you must delete it). Then save the file with a .CSV extension. Run Matlab and run the plotFFTb function. When it
    asks for a file name type the saved data file without the .CSV extension. The program will then plot the power/frequency data as well
    as the proportion of the total power in each relevant frequency band.

    John

    ------------------------------------------------------------------------------------

    % plotFFTb.m
    %
    % Christopher DiMattina & John Rossi III
    %
    %
    % Florida Gulf Coast University
    %
    % Brief Description: This program loads an EEG data file and plots the FFT
    % and PSD. It also calculates the signal power in the
    % alpha, theta, delta beta and gamma frequency bands
    % as a proportion of overall power
    %
    % Version: 2.0 Beta
    %
    % Created: 9-24-2014
    %
    % Copyright (C) 2014 Christopher DiMattina & John Rossi III
    %
    % This program 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 3 of the License, or
    % (at your option) any later version.
    %
    % This program 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 this program. If not, see <http://www.gnu.org/licenses/>.
    %

    Fs = 512; % Sampling frequency = 512 Hz
    fmax = Fs/2; % Nyquist frequency

    % Get Filename

    filename = input('Filename to load? ', 's');
    filename = strcat(filename, '.csv');

    M = csvread(filename);

    % M = csvread('test.csv');

    t = M(:,1);
    y = M(:,2);
    y = y-mean(y); % Remove DC component

    N = 4096;
    FFTy = fftshift(abs(fft(y,N)));
    faxis = linspace(-fmax,fmax,N-1);

    plotInd = (N/2):length(faxis);
    FFTy_temp = FFTy(2:end);

    faxis = faxis(plotInd);
    FFTy = FFTy_temp(plotInd);

    FFTpow = (FFTy.^2);

    plot(faxis,FFTy);

    deltaInd = find( (faxis > 0.1) & (faxis <=3));
    alphaInd = find( (faxis > 8 ) & (faxis <= 13) );
    betaInd = find( (faxis > 13) & (faxis <= 30) );
    gammaInd = find( (faxis > 30) & (faxis <=100) );
    thetaInd = find( (faxis > 3 ) & (faxis <= 8) );
    junkloInd = find( (faxis < 0.1) );
    junkhiInd = find( (faxis > 100));

    deltaFrac = sum(FFTpow(deltaInd))/sum(FFTpow);
    alphaFrac = sum(FFTpow(alphaInd))/sum(FFTpow);
    betaFrac = sum(FFTpow(betaInd))/sum(FFTpow);
    gammaFrac = sum(FFTpow(gammaInd))/sum(FFTpow);
    thetaFrac = sum(FFTpow(thetaInd))/sum(FFTpow);
    junkloFrac = sum(FFTpow(junkloInd))/sum(FFTpow);
    junkhiFrac = sum(FFTpow(junkhiInd))/sum(FFTpow);

    disp(sprintf('alpha = %.3f',alphaFrac));
    disp(sprintf('theta = %.3f',thetaFrac));
    disp(sprintf('delta = %.3f',deltaFrac));
    disp(sprintf('beta = %.3f',betaFrac));
    disp(sprintf('gamma = %.3f',gammaFrac));
    disp(sprintf('lowout = %.3f',junkloFrac));
    disp(sprintf('hiout = %.3f',junkhiFrac));
    Attached Files Attached Files

Thread Information

Users Browsing this Thread

There are currently 1 users browsing this thread. (0 members and 1 guests)

Similar Threads

  1. Mind Workstation Update with Enhanced Kasina Supoort
    By neuroasis in forum Kasina & Limina
    Replies: 16
    Last Post: 01-09-2014, 04:09 AM
  2. Mind Workstation Effects
    By Thalassa in forum Kasina & Limina
    Replies: 13
    Last Post: 12-29-2013, 10:35 AM
  3. Article: MindPlace Kasina and Mind Workstation
    By neuroasis in forum General Forum
    Replies: 0
    Last Post: 12-21-2013, 12:44 PM
  4. mind workstation
    By filemon in forum Proteus
    Replies: 2
    Last Post: 02-12-2012, 07:36 AM
  5. The Mind Workstation
    By hseaver in forum General Forum
    Replies: 3
    Last Post: 05-10-2008, 02:25 PM

Tags for this Thread

Posting Permissions

  • You may not post new threads
  • You may not post replies
  • You may not post attachments
  • You may not edit your posts
  •