function skjdemo4
% SKJDEMO3 demoaa signaalin jaksollisuutta.
%
% Käytetään signaalina kitaran ääntä.

% Jukka Parviainen, parvi@hut.fi, 15.1.2002, 17.1.2001.
% Tiedostot saatavilla 
% http://www.cis.hut.fi/Opinnot/T-61.140/  -> Vinkkejä ja linkkejä

% Alkutoimenpiteet: tarkista että esityskoneessasi on 
% kaiuittimet tai salin kyseessä ollen audioliitäntä.
% Kutista Matlab-työikkuna suhteellisen pieneksi ruudun alareunaan,
% sijoita piirtoikkuna, esimerkiksi ensimmäisen vaiheen jälkeen,
% sopivaan paikkaan, että se näkyy.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

close all;      % sulkee kaikki ikkunat
clc;            % puhdistaa ruudun

fprintf(1,'****************************************\n');
fprintf(1,'****************************************\n');
fprintf(1,'************  SKJDEMO4.m  **************\n');
fprintf(1,'****************************************\n');
fprintf(1,'*******  Jaksolliset signaalit  ********\n');
fprintf(1,'*************  kitara  *****************\n');
fprintf(1,'****************************************\n');
fprintf(1,'\n\n');

[yb, fsb, nbits] = wavread('skjdemo3b.wav');
[yc, fsb, nbits] = wavread('skjdemo3c.wav');
ycF = fft(yc);


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

fprintf(1,'Mikä merkitys harmonisilla komponenteilla on?\n\n');

fprintf(1,'Kuunnellaan seuraavaksi kitaraesitystä\n ja yhden kielen värähdystä.\n\n');

fprintf(1,'** Paina ENTER...\n\n');
pause;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

soundsc(yb, fsb);

fprintf(1,'** Paina ENTER musiikkikappaleen kuuntelemisen jälkeen...\n\n');
pause;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clc;

fprintf(1,'Otetaan YKSI kielen värähdys \n');
fprintf(1,' ja piirretään siitä signaali (kesto n. 0.035 s)\n');
fprintf(1,' ja sen spektri.\n\n');

w = linspace(0, fsb, length(ycF));
t = linspace(0, length(yc)/fsb, length(yc));

clf;
subplot(211),
plot(t, yc);
axis([0 t(end) 1.1*min(yc) 1.1*max(yc)]);
grid on;
xlabel('Aika (sekunnit)');
ylabel('Voimakkuus');
title('Yhden kitarankielen värähdys');

subplot(212);
%plot(w, abs(ycF));
%axis([0 4300 0 1.1*max(abs(ycF))]);
plot(w, 20*log2(abs(ycF)));
axis([0 4300 0 1.1*max(20*log2(abs(ycF)))]);
xlabel('Taajuus (Hz)');
ylabel('logaritminen');
%ylabel('Voimakkuus');
grid on;

fprintf(1,'Värähdys soitetaan seuraavaksi 10 kertaa peräkkäin.\n\n');

fprintf(1,'** Paina ENTER...\n\n');
pause;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for i = 1:10
  soundsc(yc, fsb);
  pause(0.3);
end;

clc;

fprintf(1,'Miltä yksinäinen kitarankielen värähdys kuulosti?\n\n');

fprintf(1,'Mikä merkitys harmonisilla komponenteilla on?\n\n');

fprintf(1,'** LOPPU.\n\n');

