octave_packages/tsa-4.2.4/sbispec.m

   1 function sbispec(BISPEC)
   2 % SBISPEC show BISPECTRUM
   3
   4 %       $Id: sbispec.m 5090 2008-06-05 08:12:04Z schloegl $
   5 %       Copyright (C) 1997, 1998, 2008 by Alois Schloegl <a.schloegl@ieee.org>
   6 %
   7 %    This program is free software: you can redistribute it and/or modify
   8 %    it under the terms of the GNU General Public License as published by
   9 %    the Free Software Foundation, either version 3 of the License, or
  10 %    (at your option) any later version.
  11 %
  12 %    This program is distributed in the hope that it will be useful,
  13 %    but WITHOUT ANY WARRANTY; without even the implied warranty of
  14 %    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15 %    GNU General Public License for more details.
  16 %
  17 %    You should have received a copy of the GNU General Public License
  18 %    along with this program.  If not, see <http://www.gnu.org/licenses/>.
  19
  20
  21 [s1,s2]=size(BISPEC);
  22 t1=(0:s1-1)/max(s1);
  23 t2=(0:s2/2-1)/max(s2);
  24 tmp=tril(NaN*ones(s1,s2),-1);
  25 BISPEC=BISPEC+tmp+rot90(tmp);
  26 BISPEC=BISPEC(1:s1/2,:);
  27
  28 subplot(211);
  29 mesh(t1,t2,abs(BISPEC));
  30 title('Bispectrum - mesh plot');
  31
  32 subplot(212);
  33 contour(t1,t2,abs(BISPEC));
  34 %contour(t1,t2,log(abs(BISPEC)));
  35 title('Bispectrum - contour plot');