X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?p=CreaPhase.git;a=blobdiff_plain;f=octave_packages%2Ftsa-4.2.4%2Fbisdemo.m;fp=octave_packages%2Ftsa-4.2.4%2Fbisdemo.m;h=d8613e11645cbe1173e70cd47d852a5c5e783f83;hp=0000000000000000000000000000000000000000;hb=c880e8788dfc484bf23ce13fa2787f2c6bca4863;hpb=1705066eceaaea976f010f669ce8e972f3734b05 diff --git a/octave_packages/tsa-4.2.4/bisdemo.m b/octave_packages/tsa-4.2.4/bisdemo.m new file mode 100644 index 0000000..d8613e1 --- /dev/null +++ b/octave_packages/tsa-4.2.4/bisdemo.m @@ -0,0 +1,34 @@ +% BISDEMO (script) Shows BISPECTRUM of eeg8s.mat + +% $Id: bisdemo.m 5090 2008-06-05 08:12:04Z schloegl $ +% Copyright (C) 1997, 1998,2008 by Alois Schloegl +% +% 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 . + +load eeg8s.mat; +[BISPEC,BICOV,ACF]=bispec(eeg8s,30); +[s1,s2]=size(BISPEC); +t1=(1:s1)/max(s1); +t2=(1:s2)/max(s2); +subplot(211); +mesh(t1,t2,abs(BISPEC)); +title('Bispectrum - mesh plot'); + +subplot(212); +if exist('OCTAVE_VERSION')>5 + contour(abs(BISPEC),10,t1,t2); +else + contour(t1,t2,abs(BISPEC),10); +end; +title('Bispectrum - contour plot');