xdel(winsid()) 
  
  // Podstawowe Test 
f1=25; 
A1=50; 
f2=27.5; 
A2=100; 
f3=50; 
A3=50; 
dt=1/1024; 
T=1; 
 //wektor czasu 
 t=0:dt:1-dt; 
 // Utworzenie sygnałów skłądowych oraz sygnału wypadkowego poliharmonicznego 
 s1=A1*sin(2*%pi*f1*t); 
 s2=A2*cos(2*%pi*f2*t); 
 s3=A3*sin(2*%pi*f3*t); 
 sygnal=s1+s2+s3;  
  
  //Wykreslenie sygnalów w dzienicznie czasu  
  scf()
  subplot(311) 
  plot(t,s1) 
  title(' sygnal A1*sin(2*%pi*f1*t)') 
  xlabel('Czestotliwosc [Hz]') 
  ylabel('Amplituda') 
  xgrid 
  subplot(312) 
  plot(t,s2) 
  title(' sygnal A2*cos(2*%pi*f2*t)') 
  xlabel('Czestotliwosc [Hz]') 
  ylabel('Amplituda') 
  xgrid 
   subplot(313) 
  plot(t,sygnal) 
  title('sygnał poliharmoniczny A1*sin(2*%pi*f1*t) + A2*cos(2*%pi*f2*t)') 
  xlabel('Czestotliwosc [Hz]') 
  ylabel('Amplituda') 
  xgrid 
   
   //Utworzenie okna 
   okno=window('hn',length(t)); 
   sw=sygnal.*okno; 
    
    //Szybka transformata Fourera 
    ffsignal=fft(sw,-1); 
     
  // Skalowanie Osi  
 N = T/dt;
incf =1/(N*dt);  
wektor_f=(0:N-1)
wektor_f(((N/2)+1):N)=(wektor_f(((N/2)+1):N))-N 
wektor_f=wektor_f*incf// 
 
 //Wyrysowanie Widma Dwustronnego przeskalowanego
 scf()  
  subplot(211)
title('przeskalowana os amplitud - widmo dwustronne')
plot2d(wektor_f,abs(ffsignal/N));
ylabel('amplituda')  
xlabel('częstotliwośc [Hz]') 
xgrid 
   //jednostronne 
   subplot(212)
title('widmo amplitudowe jednostronne')
plot2d(wektor_f(1:length(ffsignal)/2),abs((1/N)*ffsignal(1,1:length(ffsignal)/2)));
xlabel('częstotliwośc [Hz]')
ylabel('amplituda') 
xgrid  
 
 //Wykreslenie Widma Mocy dla okna Kaisera oraz Hanninga 
  sm1=pspect(100,256,'hn',sygnal); 
  sm2=pspect(100,256,'kr',sygnal,4); 
  smize1=length(sm1); 
  smize2=length(sm2);  
  fr1=((1:smize1/2)/(smize1/length(sygnal)))*incf; 
  fr2=((1:smize2/2)/(smize2/length(sygnal)))*incf; 
  scf()
  subplot(211) 
  plot(fr1,sm1(1:smize1/2))
 xtitle('Widmo Mocy sygnalu poliharmonicnzego dla okna Hanninga') 
 xlabel('Czestotliwosc [Hz]') 
 ylabel('Amplituda') 
 xgrid 
  subplot(212) 
  plot(fr2,sm2(1:smize2/2))
 xtitle('Widmo Mocy sygnalu poliharmonicnzego dla okna Kaisera') 
 xlabel('Czestotliwosc [Hz]') 
 ylabel('Amplituda') 
 xgrid