%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% SMS - 11.cvičení: Fyzikální modelování
% R.Čmejla, 11.prosince 2017
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
1. Napište funkci pro základní Karplusův-Strongův algoritmus
2. Demonstrujte závislost algoritmu na fs
3. Napište modifikovaný Karplusův-Strongův algoritmus 
   se zvýšenou přesností generování základní frekvence
4. Modifikovaný KS algoritmus upravte jako funkci fyzikálního modelu
     - délka tónu
     - požadovaná výstupní frekvence
     - typ drnknutí
     - tlumení kobylky
     - tuhost strun
     - hlasitost

   ------------------------------------------------------------------
5. Experimenty s K-S algoritmem a fyzikální modelování  
   KLAVÍR: (KARPLUS STRONG)

    ----
   | x1 |----->-----   obálka
    ----            |     |      ----
        (k*x1 + x2) + ->- x -->-| DP |---> klavír
    ----            |            ----
   | x2 |----->-----
    ---- 

    x1  ... harmonický signál    (vstup fs, f0, T)
    x2  ... výstup K-S algoritmu (vstup fs, f0, T)
    k    ~  0.4
    obálka  t = 0,      50,  150,  500 [ms]
            A = 0.95, 1.00, 0.25, 0.00 [-]
    DP   ... (rezonátor); fm = 100 Hz; B  = 660 Hz

   ------------------------------------------------------------------
   BICÍ: (KARPLUS STRONG)
   https://ccrma.stanford.edu/~sdill/220A-project/drum_KS.m       % matlab code
   https://ccrma.stanford.edu/~sdill/220A-project/drums.html#ks

   https://math.aalto.fi/~fdahl/signals/matlab.html               % matlab code
   https://math.aalto.fi/~fdahl/signals/samples.pdf

   https://www.youtube.com/watch?v=WYsE74odhbs

   ------------------------------------------------------------------
   ELEKTRICKÁ KYTARA
   http://www.music.mcgill.ca/~gary/courses/projects/618_2009/NickDonaldson/
   https://ccrma.stanford.edu/~hugo/220a/project/

   ------------------------------------------------------------------
   HOUSLE (experiment - výtup KS v konvoluci s (vhodnou) obálkou)
   https://people.ece.cornell.edu/land/courses/ece4760/FinalProjects/s2009/jvt6_th389/jvt6_th389/finalproject.html

   ------------------------------------------------------------------
   KLARINET: Fyzikální modelování
   http://www.cs.sfu.ca/~tamaras/cmpt889assignments/clarinet889.m % matlab code
   https://www.cs.sfu.ca/~tamaras/cmpt889projects/evgeny/
   http://www.cs.sfu.ca/~tamaras/cmpt889assignments/PhysModDigitalWaveGuides.pdf