home *** CD-ROM | disk | FTP | other *** search
/ PC World 2007 January / PCWorld_2007-01_cd.bin / temacd / reaper / reaper147-install.exe / Effects / SStillwell / rbj4eq < prev    next >
Text File  |  2006-10-11  |  7KB  |  272 lines
  1. // Copyright 2006, 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.  
  25. desc:4-band semi-parametric EQ, based on RBJ Filter Cookbook
  26.  
  27. slider1:3<1,5,1{40,80,160,315,500}>Band 1 Freq
  28. slider2:0<-15,15,0.1>Boost/Cut
  29. slider3:3<1,5,1{125,250,500,1k,2k}>Band 2 Freq
  30. slider4:0<-15,15,0.1>Boost/Cut
  31. slider5:3<1,5,1{315,630,1.2k,2.5k,5k}>Band 3 Freq
  32. slider6:0<-15,15,0.1>Boost/Cut
  33. slider7:3<1,5,1{1.6k,3.2k,6.4k,9k,12k}>Band 4 Freq
  34. slider8:0<-15,15,0.1>Boost/Cut
  35.  
  36. @init
  37.   gain1 = 0;
  38.   freq1 = 160;
  39.   a1 = 10^(gain1/40);
  40.   q1 = 0.8;
  41.   w01 = 2 * $pi * freq1/srate;
  42.   cosw01 = cos(w01);
  43.   sinw01 = sin(w01);
  44.   alpha1 = sinw01 / (2 * q1);
  45.  
  46.   b01 = 1 + alpha1 * a1;
  47.   b11 = -2 * cosw01;
  48.   b21 = 1 - alpha1 * a1;
  49.   a01 = 1 + alpha1 / a1;
  50.   a11 = -2 * cosw01;
  51.   a21 = 1 - alpha1 / a1;
  52.   b01 /= a01;
  53.   b11 /= a01;
  54.   b21 /= a01;
  55.   a11 /= a01;
  56.   a21 /= a01;
  57.  
  58.   gain3 = 0;
  59.   freq3 = 500;
  60.   a3 = 10^(gain3/40);
  61.   q3 = 0.8;
  62.   w03 = 2 * $pi * freq3/srate;
  63.   cosw03 = cos(w03);
  64.   sinw03 = sin(w03);
  65.   alpha3 = sinw03 / (2 * q3);
  66.  
  67.   b03 = 1 + alpha3 * a3;
  68.   b13 = -2 * cosw03;
  69.   b23 = 1 - alpha3 * a3;
  70.   a03 = 1 + alpha3 / a3;
  71.   a13 = -2 * cosw03;
  72.   a23 = 1 - alpha3 / a3;
  73.   b03 /= a03;
  74.   b13 /= a03;
  75.   b23 /= a03;
  76.   a13 /= a03;
  77.   a23 /= a03;
  78.  
  79.   gain5 = 0;
  80.   freq5 = 1200;
  81.   a5 = 10^(gain5/40);
  82.   q5 = 0.8;
  83.   w05 = 2 * $pi * freq5/srate;
  84.   cosw05 = cos(w05);
  85.   sinw05 = sin(w05);
  86.   alpha5 = sinw05 / (2 * q5);
  87.  
  88.   b05 = 1 + alpha5 * a5;
  89.   b15 = -2 * cosw05;
  90.   b25 = 1 - alpha5 * a5;
  91.   a05 = 1 + alpha5 / a5;
  92.   a15 = -2 * cosw05;
  93.   a25 = 1 - alpha5 / a5;
  94.   b05 /= a05;
  95.   b15 /= a05;
  96.   b25 /= a05;
  97.   a15 /= a05;
  98.   a25 /= a05;
  99.  
  100.   gain7 = 0;
  101.   freq7 = 6400;
  102.   a7 = 10^(gain7/40);
  103.   q7 = 0.8;
  104.   w07 = 2 * $pi * freq7/srate;
  105.   cosw07 = cos(w07);
  106.   sinw07 = sin(w07);
  107.   alpha7 = sinw07 / (2 * q7);
  108.  
  109.   b07 = 1 + alpha7 * a7;
  110.   b17 = -2 * cosw07;
  111.   b27 = 1 - alpha7 * a7;
  112.   a07 = 1 + alpha7 / a7;
  113.   a17 = -2 * cosw07;
  114.   a27 = 1 - alpha7 / a7;
  115.   b07 /= a07;
  116.   b17 /= a07;
  117.   b27 /= a07;
  118.   a17 /= a07;
  119.   a27 /= a07;
  120.  
  121.  
  122. @slider
  123.   freq1 = (slider1 == 0 ? 40 : (slider1 == 1 ? 80 : (slider1 == 2 ? 160 : (slider1 == 3 ? 315 : 500))));
  124.   gain1 = slider2;
  125.   freq3 = (slider3 == 0 ? 125 : (slider3 == 1 ? 250 : (slider3 == 2 ? 500 : (slider3 == 3 ? 1000 : 2000))));
  126.   gain3 = slider4;
  127.   freq5 = (slider5 == 0 ? 315 : (slider5 == 1 ? 630 : (slider5 == 2 ? 1200 : (slider5 == 3 ? 2500 : 5000))));
  128.   gain5 = slider6;
  129.   freq7 = (slider7 == 0 ? 1600 : (slider7 == 1 ? 3200 : (slider7 == 2 ? 6400 : (slider7 == 3 ? 9000 : 12000))));
  130.   gain7 = slider8;
  131.  
  132.   a1 = 10^(gain1/40);
  133.   q1 = 0.8;
  134.   w01 = 2 * $pi * freq1/srate;
  135.   cosw01 = cos(w01);
  136.   sinw01 = sin(w01);
  137.   alpha1 = sinw01 / (2 * q1);
  138.  
  139.   b01 = 1 + alpha1 * a1;
  140.   b11 = -2 * cosw01;
  141.   b21 = 1 - alpha1 * a1;
  142.   a01 = 1 + alpha1 / a1;
  143.   a11 = -2 * cosw01;
  144.   a21 = 1 - alpha1 / a1;
  145.   b01 /= a01;
  146.   b11 /= a01;
  147.   b21 /= a01;
  148.   a11 /= a01;
  149.   a21 /= a01;
  150.  
  151.   a3 = 10^(gain3/40);
  152.   q3 = 0.8;
  153.   w03 = 2 * $pi * freq3/srate;
  154.   cosw03 = cos(w03);
  155.   sinw03 = sin(w03);
  156.   alpha3 = sinw03 / (2 * q3);
  157.  
  158.   b03 = 1 + alpha3 * a3;
  159.   b13 = -2 * cosw03;
  160.   b23 = 1 - alpha3 * a3;
  161.   a03 = 1 + alpha3 / a3;
  162.   a13 = -2 * cosw03;
  163.   a23 = 1 - alpha3 / a3;
  164.   b03 /= a03;
  165.   b13 /= a03;
  166.   b23 /= a03;
  167.   a13 /= a03;
  168.   a23 /= a03;
  169.  
  170.   a5 = 10^(gain5/40);
  171.   q5 = 0.8;
  172.   w05 = 2 * $pi * freq5/srate;
  173.   cosw05 = cos(w05);
  174.   sinw05 = sin(w05);
  175.   alpha5 = sinw05 / (2 * q5);
  176.  
  177.   b05 = 1 + alpha5 * a5;
  178.   b15 = -2 * cosw05;
  179.   b25 = 1 - alpha5 * a5;
  180.   a05 = 1 + alpha5 / a5;
  181.   a15 = -2 * cosw05;
  182.   a25 = 1 - alpha5 / a5;
  183.   b05 /= a05;
  184.   b15 /= a05;
  185.   b25 /= a05;
  186.   a15 /= a05;
  187.   a25 /= a05;
  188.  
  189.   a7 = 10^(gain7/40);
  190.   q7 = 0.8;
  191.   w07 = 2 * $pi * freq7/srate;
  192.   cosw07 = cos(w07);
  193.   sinw07 = sin(w07);
  194.   alpha7 = sinw07 / (2 * q7);
  195.  
  196.   b07 = 1 + alpha7 * a7;
  197.   b17 = -2 * cosw07;
  198.   b27 = 1 - alpha7 * a7;
  199.   a07 = 1 + alpha7 / a7;
  200.   a17 = -2 * cosw07;
  201.   a27 = 1 - alpha7 / a7;
  202.   b07 /= a07;
  203.   b17 /= a07;
  204.   b27 /= a07;
  205.   a17 /= a07;
  206.   a27 /= a07;
  207.  
  208. @sample
  209.   gain1 != 0 ? (
  210.   ospl0 = spl0;
  211.   spl0 = b01 * spl0 + b11 * xl11 + b21 * xl21 - a11 * yl11 - a21 * yl21;
  212.   xl21 = xl11;
  213.   xl11 = ospl0;
  214.   yl21 = yl11;
  215.   yl11 = spl0;
  216.  
  217.   ospl1 = spl1;
  218.   spl1 = b01 * spl1 + b11 * xr11 + b21 * xr21 - a11 * yr11 - a21 * yr21;
  219.   xr21 = xr11;
  220.   xr11 = ospl1;
  221.   yr21 = yr11;
  222.   yr11 = spl1;
  223.   );
  224.  
  225.   gain3 != 0 ? (
  226.   ospl0 = spl0;
  227.   spl0 = b03 * spl0 + b13 * xl13 + b23 * xl23 - a13 * yl13 - a23 * yl23;
  228.   xl23 = xl13;
  229.   xl13 = ospl0;
  230.   yl23 = yl13;
  231.   yl13 = spl0;
  232.  
  233.   ospl1 = spl1;
  234.   spl1 = b03 * spl1 + b13 * xr13 + b23 * xr23 - a13 * yr13 - a23 * yr23;
  235.   xr23 = xr13;
  236.   xr13 = ospl1;
  237.   yr23 = yr13;
  238.   yr13 = spl1;
  239.   );
  240.  
  241.   gain5 != 0 ? (
  242.   ospl0 = spl0;
  243.   spl0 = b05 * spl0 + b15 * xl15 + b25 * xl25 - a15 * yl15 - a25 * yl25;
  244.   xl25 = xl15;
  245.   xl15 = ospl0;
  246.   yl25 = yl15;
  247.   yl15 = spl0;
  248.  
  249.   ospl1 = spl1;
  250.   spl1 = b05 * spl1 + b15 * xr15 + b25 * xr25 - a15 * yr15 - a25 * yr25;
  251.   xr25 = xr15;
  252.   xr15 = ospl1;
  253.   yr25 = yr15;
  254.   yr15 = spl1;
  255.   );
  256.  
  257.   gain7 != 0 ? (
  258.   ospl0 = spl0;
  259.   spl0 = b07 * spl0 + b17 * xl17 + b27 * xl27 - a17 * yl17 - a27 * yl27;
  260.   xl27 = xl17;
  261.   xl17 = ospl0;
  262.   yl27 = yl17;
  263.   yl17 = spl0;
  264.  
  265.   ospl1 = spl1;
  266.   spl1 = b07 * spl1 + b17 * xr17 + b27 * xr27 - a17 * yr17 - a27 * yr27;
  267.   xr27 = xr17;
  268.   xr17 = ospl1;
  269.   yr27 = yr17;
  270.   yr17 = spl1;
  271.   );
  272.