1 function sbispec(BISPEC)
2 % SBISPEC show BISPECTRUM
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>
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.
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.
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/>.
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,:);
29 mesh(t1,t2,abs(BISPEC));
30 title('Bispectrum - mesh plot');
33 contour(t1,t2,abs(BISPEC));
34 %contour(t1,t2,log(abs(BISPEC)));
35 title('Bispectrum - contour plot');