home *** CD-ROM | disk | FTP | other *** search
/ PC World 2007 March / PCWorld_2007-03_cd.bin / audio-video / reaper / reaper174-install.exe / Effects / SStillwell / exciter < prev    next >
Text File  |  2007-01-14  |  5KB  |  174 lines
  1. // Copyright 2007, Thomas Scott Stillwell
  2. // All rights reserved.
  3. //
  4. //Redistribution and use in source and binary forms, with or without modification, are permitted 
  5. //provided that the following conditions are met:
  6. //
  7. //Redistributions of source code must retain the above copyright notice, this list of conditions 
  8. //and the following disclaimer. 
  9. //
  10. //Redistributions in binary form must reproduce the above copyright notice, this list of conditions 
  11. //and the following disclaimer in the documentation and/or other materials provided with the distribution. 
  12. //
  13. //The name of Thomas Scott Stillwell may not be used to endorse or 
  14. //promote products derived from this software without specific prior written permission. 
  15. //
  16. //THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR 
  17. //IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND 
  18. //FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS 
  19. //BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES 
  20. //(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 
  21. //PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, 
  22. //STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF 
  23. //THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  24. desc:exciter - treble enhancer
  25.  
  26. slider1:0<0,100,0.1>mix%
  27. slider2:0<0,100,0.1>drive%
  28. slider3:5000<2000,10000,100>frequency
  29.  
  30. @init
  31.  
  32.   gain1 = 0;
  33.   freq1 = freq;
  34.   a1 = 1;
  35.   s1 = 1;
  36.   q1 = 1 / (sqrt((a1 + 1/a1)*(1/s1 - 1) + 2));
  37.   w01 = 2 * $pi * freq1/srate;
  38.   cosw01 = cos(w01);
  39.   sinw01 = sin(w01);
  40.   alpha1 = sinw01 / (2 * q1);
  41.  
  42.   b01 = (1 + cosw01)/2;
  43.   b11 = -(1 + cosw01);
  44.   b21 = (1 + cosw01)/2;
  45.   a01 = 1 + alpha1;
  46.   a11 = -2 * cosw01;
  47.   a21 = 1 - alpha1;
  48.   b01 /= a01;
  49.   b11 /= a01;
  50.   b21 /= a01;
  51.   a11 /= a01;
  52.   a21 /= a01;
  53.  
  54.   gain3 = 0;
  55.   freq3 = freq * 2;
  56.   a3 = 10^(gain3/40);
  57.   s3 = 1;
  58.   q3 = 1 / (sqrt((a3 + 1/a3)*(1/s3 - 1) + 2));
  59.   w03 = 2 * $pi * freq3/srate;
  60.   cosw03 = cos(w03);
  61.   sinw03 = sin(w03);
  62.   alpha3 = sinw03 / (2 * q3);
  63.  
  64.   b03 = (1 + cosw03)/2;
  65.   b13 = -(1 + cosw03);
  66.   b23 = (1 + cosw03)/2;
  67.   a03 = 1 + alpha3;
  68.   a13 = -2 * cosw03;
  69.   a23 = 1 - alpha3;
  70.   b03 /= a03;
  71.   b13 /= a03;
  72.   b23 /= a03;
  73.   a13 /= a03;
  74.   a23 /= a03;
  75.  
  76. @slider
  77.  
  78.   mix=(10 * slider1) / 100;
  79.   drive=slider2 / 100;
  80.   mix1 = 1 - mix;
  81.   drive1 = 1 / (1 - (drive / 2));
  82.   drive2 = drive / 2;
  83.   freq = slider3;
  84.  
  85.   gain1 = 0;
  86.   freq1 = freq;
  87.   a1 = 1;
  88.   s1 = 1;
  89.   q1 = 1 / (sqrt((a1 + 1/a1)*(1/s1 - 1) + 2));
  90.   w01 = 2 * $pi * freq1/srate;
  91.   cosw01 = cos(w01);
  92.   sinw01 = sin(w01);
  93.   alpha1 = sinw01 / (2 * q1);
  94.  
  95.   b01 = (1 - cosw01)/2;
  96.   b11 = (1 - cosw01);
  97.   b21 = (1 - cosw01)/2;
  98.   a01 = 1 + alpha1;
  99.   a11 = -2 * cosw01;
  100.   a21 = 1 - alpha1;
  101.   b01 /= a01;
  102.   b11 /= a01;
  103.   b21 /= a01;
  104.   a11 /= a01;
  105.   a21 /= a01;
  106.  
  107.   gain3 = 0;
  108.   freq3 = freq * 2;
  109.   a3 = 10^(gain3/40);
  110.   s3 = 1;
  111.   q3 = 1 / (sqrt((a3 + 1/a3)*(1/s3 - 1) + 2));
  112.   w03 = 2 * $pi * freq3/srate;
  113.   cosw03 = cos(w03);
  114.   sinw03 = sin(w03);
  115.   alpha3 = sinw03 / (2 * q3);
  116.  
  117.   b03 = (1 + cosw03)/2;
  118.   b13 = -(1 + cosw03);
  119.   b23 = (1 + cosw03)/2;
  120.   a03 = 1 + alpha3;
  121.   a13 = -2 * cosw03;
  122.   a23 = 1 - alpha3;
  123.   b03 /= a03;
  124.   b13 /= a03;
  125.   b23 /= a03;
  126.   a13 /= a03;
  127.   a23 /= a03;
  128.  
  129.  
  130. @sample
  131.  
  132.   dry0 = spl0;
  133.   dry1 = spl1;
  134.  
  135.   //hpf to extract HF signal
  136.   ospl0 = dry0;
  137.   dry0 = b01 * dry0 + b11 * xl11 + b21 * xl21 - a11 * yl11 - a21 * yl21;
  138.   xl21 = xl11;
  139.   xl11 = ospl0;
  140.   yl21 = yl11;
  141.   yl11 = dry0;
  142.  
  143.   ospl1 = dry1;
  144.   dry1 = b01 * dry1 + b11 * xr11 + b21 * xr21 - a11 * yr11 - a21 * yr21;
  145.   xr21 = xr11;
  146.   xr11 = ospl1;
  147.   yr21 = yr11;
  148.   yr11 = dry1;
  149.  
  150.   //enhance HF signal
  151.   wet0 = drive1 * (1.5 * dry0) * (1 - abs((1.5 * dry0) * drive2));
  152.   wet1 = drive1 * (1.5 * dry1) * (1 - abs((1.5 * dry1) * drive2)); 
  153.   wet0 /= 1.5;
  154.   wet1 /= 1.5;
  155.  
  156.   //hpf to isolate enhance signal
  157.   ospl0 = wet0;
  158.   wet0 = b03 * wet0 + b13 * xl13 + b23 * xl23 - a13 * yl13 - a23 * yl23;
  159.   xl23 = xl13;
  160.   xl13 = ospl0;
  161.   yl23 = yl13;
  162.   yl13 = wet0;
  163.  
  164.   ospl1 = wet1;
  165.   wet1 = b03 * wet1 + b13 * xr13 + b23 * xr23 - a13 * yr13 - a23 * yr23;
  166.   xr23 = xr13;
  167.   xr13 = ospl1;
  168.   yr23 = yr13;
  169.   yr13 = wet1;
  170.  
  171.  
  172.   spl0 = spl0 + (mix) * wet0;
  173.   spl1 = spl1 + (mix) * wet1;
  174.