bisprectum_test2.asv

高阶谱工具箱 Version 2.0.3 (R12 compliant) 27 Dec 2

%EX_SUNS: Sunspot data analysis
%
echo off

% A. Swami April 15, 1995 
% Copyright (c) 1991-2001 by United Signals & Systems, Inc. 
%       $Revision: 1.6 $

%     RESTRICTED RIGHTS LEGEND
% Use, duplication, or disclosure by the Government is subject to
% restrictions as set forth in subparagraph (c) (1) (ii) of the 
% Rights in Technical Data and Computer Software clause of DFARS
% 252.227-7013. 
% Manufacturer: United Signals & Systems, Inc., P.O. Box 2374, 
% Culver City, California 90231. 
%
%  This material may be reproduced by or for the U.S. Government pursuant 
%  to the copyright license under the clause at DFARS 252.227-7013. 

clear, 
echo on
echo on
nt=1400;
f0=100;
Fs=4*f0;
nf=128;
t=1:nt;
x=sin(2*pi*f0*t)+rand(1,nt);
%[R,f1,f2]=bispectr(x,nf,nt,fs)

 %load laff.dat 		       % LAUGHTER data
 sp =x;
 %Fs = 8192;
 %load sunspot.dat             %  SUNSPOT data 
 %sp=sunspot(:,2); 
 figno(1) = gcf; 
 
 clc
 subplot(211),  plot( 1:length(sp),sp), grid on 
 subplot(212), hist(sp), 
 set(gcf,'Name','sunspot data') 

%Hit any key to continue
pause  

% These are the annual sunspot counts for years 1700-1987
%   sunspot data are all positive valued;
%   differencing the data is useful in such cases 

 sp = diff(sp);   
% sp_txt = 'sunspot data (differenced)'; 
  sp_txt =sp;
  ex_eda

echo on
%  End of sunspot demo
%  Hit any key to clear all plots, and return to previous menu
pause

echo off
for k=1:length(figno)
    close(figno(k))      % figno has figure numbers returned by ex_eda
end   
clear figno

return