% Experimente zu Spektogram und FFT
global NFFT NSig NBit SSig FSig;
% Anzahl der Werte zur FFT (muss Potenz von 2 sein! also: 128, 256, ... 65536)
NFFT=256;
% Anzahl der Wert des auszuwertenden Signals (min. NFFT)
NSig=1.5*NFFT;
% Aufloesung der Amplitude (ueblich 16, 24, 32)
NBit=16;
% Samplingrate (uebkich 44100, 48000, 88200, 96000 1/sec)
SSig=48000;
% Frequenz des Signals (in Hz)
FSig=440;
%
Pi=asin(1);
% Funktion des Signals
%
tSig=zeros(1,NSig);
ASig=zeros(1,NSig);
FTSig=zeros(1,NFFT);
FFT1=zeros(1,NFFT);
FFT1abs=zeros(1,NFFT);
FFT=zeros(1,NFFT);
FTSpec=zeros((NSig-NFFT),NFFT);

for i=1:NSig
  tSig(i)=(i-1)/SSig;
  ASig(i)=sin(2*Pi*tSig(i)*FSig);
endfor

fig = figure(1);
plot(ASig)
grid

for j=1:(NSig-NFFT)
  for k=1:NFFT
    FTSig(k)=ASig(j+k-1);
  endfor
%  FFT1=fft(ASig,NFFT);
%  FFT1abs=abs(FFT1);
%  FFT=fftshift(FFT1abs/NFFT);
endfor

fig = figure(fig+1);
plot(FTSig)
grid