home *** CD-ROM | disk | FTP | other *** search
/ PC World 2007 March / PCWorld_2007-03_cd.bin / audio-video / reaper / reaper174-install.exe / Effects / SStillwell / rbj7eq < prev    next >
Text File  |  2006-12-11  |  11KB  |  510 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:7-band + HPF EQ, based on RBJ Filter Cookbook
  26.  
  27. slider1:10<10,120.0.1>HPF Freq
  28. slider2:0<-15,15,0.1>100 Hz
  29. slider3:0<-15,15,0.1>200 Hz
  30. slider4:0<-15,15,0.1>400 Hz
  31. slider5:0<-15,15,0.1>800 Hz
  32. slider6:0<-15,15,0.1>2.5 kHz
  33. slider7:0<-15,15,0.1>6 kHz
  34. slider8:0<-15,15,0.1>12 kHz
  35.  
  36. @init
  37.   gain0 = 0;
  38.   freq0 = 20;
  39.   q0 = 1;
  40.   a0 = 1;
  41.   w00 = 2 * $pi * freq0/srate;
  42.   cosw00 = cos(w00);
  43.   sinw00 = sin(w00);
  44.   alpha0 = sinw00 / (2 * q0);
  45.  
  46.   b00 = (1 + cosw00) / 2;
  47.   b10 = -(1 + cosw00);
  48.   b20 = (1 + cosw00) / 2;
  49.   a00 = 1 + alpha0;
  50.   a10 = -2 * cosw00;
  51.   a20 = 1 - alpha0;
  52.   b00 /= a00;
  53.   b10 /= a00;
  54.   b20 /= a00;
  55.   a10 /= a00;
  56.   a20 /= a00;
  57.  
  58.   gain1 = 0;
  59.   freq1 = 100;
  60.   q1 = 0.8;
  61.   a1 = 10^(gain1/40);
  62.   w01 = 2 * $pi * freq1/srate;
  63.   cosw01 = cos(w01);
  64.   sinw01 = sin(w01);
  65.   alpha1 = sinw01 / (2 * q1);
  66.  
  67.   b01 = 1 + alpha1 * a1;
  68.   b11 = -2 * cosw01;
  69.   b21 = 1 - alpha1 * a1;
  70.   a01 = 1 + alpha1 / a1;
  71.   a11 = -2 * cosw01;
  72.   a21 = 1 - alpha1 / a1;
  73.   b01 /= a01;
  74.   b11 /= a01;
  75.   b21 /= a01;
  76.   a11 /= a01;
  77.   a21 /= a01;
  78.  
  79.   gain2 = 0;
  80.   freq2 = 200;
  81.   q2 = 0.8;
  82.   a2 = 10^(gain2/40);
  83.   w02 = 2 * $pi * freq2/srate;
  84.   cosw02 = cos(w02);
  85.   sinw02 = sin(w02);
  86.   alpha2 = sinw02 / (2 * q2);
  87.  
  88.   b02 = 1 + alpha2 * a2;
  89.   b12 = -2 * cosw02;
  90.   b22 = 1 - alpha2 * a2;
  91.   a02 = 1 + alpha2 / a2;
  92.   a12 = -2 * cosw02;
  93.   a22 = 1 - alpha2 / a2;
  94.   b02 /= a02;
  95.   b12 /= a02;
  96.   b22 /= a02;
  97.   a12 /= a02;
  98.   a22 /= a02;
  99.  
  100.   gain3 = 0;
  101.   freq3 = 400;
  102.   q3 = 0.8;
  103.   a3 = 10^(gain3/40);
  104.   w03 = 2 * $pi * freq3/srate;
  105.   cosw03 = cos(w03);
  106.   sinw03 = sin(w03);
  107.   alpha3 = sinw03 / (2 * q3);
  108.  
  109.   b03 = 1 + alpha3 * a3;
  110.   b13 = -2 * cosw03;
  111.   b23 = 1 - alpha3 * a3;
  112.   a03 = 1 + alpha3 / a3;
  113.   a13 = -2 * cosw03;
  114.   a23 = 1 - alpha3 / a3;
  115.   b03 /= a03;
  116.   b13 /= a03;
  117.   b23 /= a03;
  118.   a13 /= a03;
  119.   a23 /= a03;
  120.  
  121.   gain4 = 0;
  122.   freq4 = 800;
  123.   q4 = 0.8;
  124.   a4 = 10^(gain4/40);
  125.   w04 = 2 * $pi * freq4/srate;
  126.   cosw04 = cos(w04);
  127.   sinw04 = sin(w04);
  128.   alpha4 = sinw04 / (2 * q4);
  129.  
  130.   b04 = 1 + alpha4 * a4;
  131.   b14 = -2 * cosw04;
  132.   b24 = 1 - alpha4 * a4;
  133.   a04 = 1 + alpha4 / a4;
  134.   a14 = -2 * cosw04;
  135.   a24 = 1 - alpha4 / a4;
  136.   b04 /= a04;
  137.   b14 /= a04;
  138.   b24 /= a04;
  139.   a14 /= a04;
  140.   a24 /= a04;
  141.  
  142.   gain5 = 0;
  143.   freq5 = 2500;
  144.   q5 = 0.8;
  145.   a5 = 10^(gain5/40);
  146.   w05 = 2 * $pi * freq5/srate;
  147.   cosw05 = cos(w05);
  148.   sinw05 = sin(w05);
  149.   alpha5 = sinw05 / (2 * q5);
  150.  
  151.   b05 = 1 + alpha5 * a5;
  152.   b15 = -2 * cosw05;
  153.   b25 = 1 - alpha5 * a5;
  154.   a05 = 1 + alpha5 / a5;
  155.   a15 = -2 * cosw05;
  156.   a25 = 1 - alpha5 / a5;
  157.   b05 /= a05;
  158.   b15 /= a05;
  159.   b25 /= a05;
  160.   a15 /= a05;
  161.   a25 /= a05;
  162.  
  163.   gain6 = 0;
  164.   freq6 = 6000;
  165.   q6 = 0.8;
  166.   a6 = 10^(gain6/40);
  167.   w06 = 2 * $pi * freq6/srate;
  168.   cosw06 = cos(w06);
  169.   sinw06 = sin(w06);
  170.   alpha6 = sinw06 / (2 * q6);
  171.  
  172.   b06 = 1 + alpha6 * a6;
  173.   b16 = -2 * cosw06;
  174.   b26 = 1 - alpha6 * a6;
  175.   a06 = 1 + alpha6 / a6;
  176.   a16 = -2 * cosw06;
  177.   a26 = 1 - alpha6 / a6;
  178.   b06 /= a06;
  179.   b16 /= a06;
  180.   b26 /= a06;
  181.   a16 /= a06;
  182.   a26 /= a06;
  183.  
  184.   gain7 = 0;
  185.   freq7 = 12000;
  186.   q7 = 0.8;
  187.   a7 = 10^(gain7/40);
  188.   w07 = 2 * $pi * freq7/srate;
  189.   cosw07 = cos(w07);
  190.   sinw07 = sin(w07);
  191.   alpha7 = sinw07 / (2 * q7);
  192.  
  193.   b07 = 1 + alpha7 * a7;
  194.   b17 = -2 * cosw07;
  195.   b27 = 1 - alpha7 * a7;
  196.   a07 = 1 + alpha7 / a7;
  197.   a17 = -2 * cosw07;
  198.   a27 = 1 - alpha7 / a7;
  199.   b07 /= a07;
  200.   b17 /= a07;
  201.   b27 /= a07;
  202.   a17 /= a07;
  203.   a27 /= a07;
  204.  
  205. cDcAdd = 10^-30;
  206. cDenorm = 10^-30;
  207.  
  208. @slider
  209.   freq0 = slider1;
  210.   gain1 = slider2;
  211.   gain2 = slider3;
  212.   gain3 = slider4;
  213.   gain4 = slider5;
  214.   gain5 = slider6;
  215.   gain6 = slider7;
  216.   gain7 = slider8;
  217.  
  218.   gain0 = 0;
  219.   q0 = 1;
  220.   a0 = 1;
  221.   w00 = 2 * $pi * freq0/srate;
  222.   cosw00 = cos(w00);
  223.   sinw00 = sin(w00);
  224.   alpha0 = sinw00 / (2 * q0);
  225.  
  226.   b00 = (1 + cosw00) / 2;
  227.   b10 = -(1 + cosw00);
  228.   b20 = (1 + cosw00) / 2;
  229.   a00 = 1 + alpha0;
  230.   a10 = -2 * cosw00;
  231.   a20 = 1 - alpha0;
  232.   b00 /= a00;
  233.   b10 /= a00;
  234.   b20 /= a00;
  235.   a10 /= a00;
  236.   a20 /= a00;
  237.  
  238.   freq1 = 100;
  239.   q1 = 0.8;
  240.   a1 = 10^(gain1/40);
  241.   w01 = 2 * $pi * freq1/srate;
  242.   cosw01 = cos(w01);
  243.   sinw01 = sin(w01);
  244.   alpha1 = sinw01 / (2 * q1);
  245.  
  246.   b01 = 1 + alpha1 * a1;
  247.   b11 = -2 * cosw01;
  248.   b21 = 1 - alpha1 * a1;
  249.   a01 = 1 + alpha1 / a1;
  250.   a11 = -2 * cosw01;
  251.   a21 = 1 - alpha1 / a1;
  252.   b01 /= a01;
  253.   b11 /= a01;
  254.   b21 /= a01;
  255.   a11 /= a01;
  256.   a21 /= a01;
  257.  
  258.   freq2 = 200;
  259.   q2 = 0.8;
  260.   a2 = 10^(gain2/40);
  261.   w02 = 2 * $pi * freq2/srate;
  262.   cosw02 = cos(w02);
  263.   sinw02 = sin(w02);
  264.   alpha2 = sinw02 / (2 * q2);
  265.  
  266.   b02 = 1 + alpha2 * a2;
  267.   b12 = -2 * cosw02;
  268.   b22 = 1 - alpha2 * a2;
  269.   a02 = 1 + alpha2 / a2;
  270.   a12 = -2 * cosw02;
  271.   a22 = 1 - alpha2 / a2;
  272.   b02 /= a02;
  273.   b12 /= a02;
  274.   b22 /= a02;
  275.   a12 /= a02;
  276.   a22 /= a02;
  277.  
  278.   freq3 = 400;
  279.   q3 = 0.8;
  280.   a3 = 10^(gain3/40);
  281.   w03 = 2 * $pi * freq3/srate;
  282.   cosw03 = cos(w03);
  283.   sinw03 = sin(w03);
  284.   alpha3 = sinw03 / (2 * q3);
  285.  
  286.   b03 = 1 + alpha3 * a3;
  287.   b13 = -2 * cosw03;
  288.   b23 = 1 - alpha3 * a3;
  289.   a03 = 1 + alpha3 / a3;
  290.   a13 = -2 * cosw03;
  291.   a23 = 1 - alpha3 / a3;
  292.   b03 /= a03;
  293.   b13 /= a03;
  294.   b23 /= a03;
  295.   a13 /= a03;
  296.   a23 /= a03;
  297.  
  298.   freq4 = 800;
  299.   q4 = 0.8;
  300.   a4 = 10^(gain4/40);
  301.   w04 = 2 * $pi * freq4/srate;
  302.   cosw04 = cos(w04);
  303.   sinw04 = sin(w04);
  304.   alpha4 = sinw04 / (2 * q4);
  305.  
  306.   b04 = 1 + alpha4 * a4;
  307.   b14 = -2 * cosw04;
  308.   b24 = 1 - alpha4 * a4;
  309.   a04 = 1 + alpha4 / a4;
  310.   a14 = -2 * cosw04;
  311.   a24 = 1 - alpha4 / a4;
  312.   b04 /= a04;
  313.   b14 /= a04;
  314.   b24 /= a04;
  315.   a14 /= a04;
  316.   a24 /= a04;
  317.  
  318.   freq5 = 2500;
  319.   q5 = 0.8;
  320.   a5 = 10^(gain5/40);
  321.   w05 = 2 * $pi * freq5/srate;
  322.   cosw05 = cos(w05);
  323.   sinw05 = sin(w05);
  324.   alpha5 = sinw05 / (2 * q5);
  325.  
  326.   b05 = 1 + alpha5 * a5;
  327.   b15 = -2 * cosw05;
  328.   b25 = 1 - alpha5 * a5;
  329.   a05 = 1 + alpha5 / a5;
  330.   a15 = -2 * cosw05;
  331.   a25 = 1 - alpha5 / a5;
  332.   b05 /= a05;
  333.   b15 /= a05;
  334.   b25 /= a05;
  335.   a15 /= a05;
  336.   a25 /= a05;
  337.  
  338.   freq6 = 6000;
  339.   q6 = 0.8;
  340.   a6 = 10^(gain6/40);
  341.   w06 = 2 * $pi * freq6/srate;
  342.   cosw06 = cos(w06);
  343.   sinw06 = sin(w06);
  344.   alpha6 = sinw06 / (2 * q6);
  345.  
  346.   b06 = 1 + alpha6 * a6;
  347.   b16 = -2 * cosw06;
  348.   b26 = 1 - alpha6 * a6;
  349.   a06 = 1 + alpha6 / a6;
  350.   a16 = -2 * cosw06;
  351.   a26 = 1 - alpha6 / a6;
  352.   b06 /= a06;
  353.   b16 /= a06;
  354.   b26 /= a06;
  355.   a16 /= a06;
  356.   a26 /= a06;
  357.  
  358.   freq7 = 12000;
  359.   q7 = 0.8;
  360.   a7 = 10^(gain7/40);
  361.   w07 = 2 * $pi * freq7/srate;
  362.   cosw07 = cos(w07);
  363.   sinw07 = sin(w07);
  364.   alpha7 = sinw07 / (2 * q7);
  365.  
  366.   b07 = 1 + alpha7 * a7;
  367.   b17 = -2 * cosw07;
  368.   b27 = 1 - alpha7 * a7;
  369.   a07 = 1 + alpha7 / a7;
  370.   a17 = -2 * cosw07;
  371.   a27 = 1 - alpha7 / a7;
  372.   b07 /= a07;
  373.   b17 /= a07;
  374.   b27 /= a07;
  375.   a17 /= a07;
  376.   a27 /= a07;
  377.  
  378. @sample
  379.  
  380.   freq0 > 10 ? (
  381.   ospl0 = spl0;
  382.   spl0 = b00 * spl0 + b10 * xl10 + b20 * xl20 - a10 * yl10 - a20 * yl20;
  383.   xl20 = xl10;
  384.   xl10 = ospl0;
  385.   yl20 = yl10;
  386.   yl10 = abs(spl0) < cDenorm ? 0 : spl0 ;
  387.  
  388.   ospl1 = spl1;
  389.   spl1 = b00 * spl1 + b10 * xr10 + b20 * xr20 - a10 * yr10 - a20 * yr20;
  390.   xr20 = xr10;
  391.   xr10 = ospl1;
  392.   yr20 = yr10;
  393.   yr10 = abs(spl1) < cDenorm ? 0 : spl1 ;
  394.   );
  395.  
  396. spl0 += cDcAdd;
  397. spl1 += cDcAdd;
  398.  
  399.   gain1 != 0 ? (
  400.   ospl0 = spl0;
  401.   spl0 = b01 * spl0 + b11 * xl11 + b21 * xl21 - a11 * yl11 - a21 * yl21;
  402.   xl21 = xl11;
  403.   xl11 = ospl0;
  404.   yl21 = yl11;
  405.   yl11 = spl0;
  406.  
  407.   ospl1 = spl1;
  408.   spl1 = b01 * spl1 + b11 * xr11 + b21 * xr21 - a11 * yr11 - a21 * yr21;
  409.   xr21 = xr11;
  410.   xr11 = ospl1;
  411.   yr21 = yr11;
  412.   yr11 = spl1;
  413.   );
  414.  
  415.   gain2 != 0 ? (
  416.   ospl0 = spl0;
  417.   spl0 = b02 * spl0 + b12 * xl12 + b22 * xl22 - a12 * yl12 - a22 * yl22;
  418.   xl22 = xl12;
  419.   xl12 = ospl0;
  420.   yl22 = yl12;
  421.   yl12 = spl0;
  422.  
  423.   ospl1 = spl1;
  424.   spl1 = b02 * spl1 + b12 * xr12 + b22 * xr22 - a12 * yr12 - a22 * yr22;
  425.   xr22 = xr12;
  426.   xr12 = ospl1;
  427.   yr22 = yr12;
  428.   yr12 = spl1;
  429.   );
  430.  
  431.   gain3 != 0 ? (
  432.   ospl0 = spl0;
  433.   spl0 = b03 * spl0 + b13 * xl13 + b23 * xl23 - a13 * yl13 - a23 * yl23;
  434.   xl23 = xl13;
  435.   xl13 = ospl0;
  436.   yl23 = yl13;
  437.   yl13 = spl0;
  438.  
  439.   ospl1 = spl1;
  440.   spl1 = b03 * spl1 + b13 * xr13 + b23 * xr23 - a13 * yr13 - a23 * yr23;
  441.   xr23 = xr13;
  442.   xr13 = ospl1;
  443.   yr23 = yr13;
  444.   yr13 = spl1;
  445.   );
  446.  
  447.   gain4 != 0 ? (
  448.   ospl0 = spl0;
  449.   spl0 = b04 * spl0 + b14 * xl14 + b24 * xl24 - a14 * yl14 - a24 * yl24;
  450.   xl24 = xl14;
  451.   xl14 = ospl0;
  452.   yl24 = yl14;
  453.   yl14 = spl0;
  454.  
  455.   ospl1 = spl1;
  456.   spl1 = b04 * spl1 + b14 * xr14 + b24 * xr24 - a14 * yr14 - a24 * yr24;
  457.   xr24 = xr14;
  458.   xr14 = ospl1;
  459.   yr24 = yr14;
  460.   yr14 = spl1;
  461.   );
  462.  
  463.   gain5 != 0 ? (
  464.   ospl0 = spl0;
  465.   spl0 = b05 * spl0 + b15 * xl15 + b25 * xl25 - a15 * yl15 - a25 * yl25;
  466.   xl25 = xl15;
  467.   xl15 = ospl0;
  468.   yl25 = yl15;
  469.   yl15 = spl0;
  470.  
  471.   ospl1 = spl1;
  472.   spl1 = b05 * spl1 + b15 * xr15 + b25 * xr25 - a15 * yr15 - a25 * yr25;
  473.   xr25 = xr15;
  474.   xr15 = ospl1;
  475.   yr25 = yr15;
  476.   yr15 = spl1;
  477.   );
  478.  
  479.   gain6 != 0 ? (
  480.   ospl0 = spl0;
  481.   spl0 = b06 * spl0 + b16 * xl16 + b26 * xl26 - a16 * yl16 - a26 * yl26;
  482.   xl26 = xl16;
  483.   xl16 = ospl0;
  484.   yl26 = yl16;
  485.   yl16 = spl0;
  486.  
  487.   ospl1 = spl1;
  488.   spl1 = b06 * spl1 + b16 * xr16 + b26 * xr26 - a16 * yr16 - a26 * yr26;
  489.   xr26 = xr16;
  490.   xr16 = ospl1;
  491.   yr26 = yr16;
  492.   yr16 = spl1;
  493.   );
  494.  
  495.   gain7 != 0 ? (
  496.   ospl0 = spl0;
  497.   spl0 = b07 * spl0 + b17 * xl17 + b27 * xl27 - a17 * yl17 - a27 * yl27;
  498.   xl27 = xl17;
  499.   xl17 = ospl0;
  500.   yl27 = yl17;
  501.   yl17 = spl0;
  502.  
  503.   ospl1 = spl1;
  504.   spl1 = b07 * spl1 + b17 * xr17 + b27 * xr27 - a17 * yr17 - a27 * yr27;
  505.   xr27 = xr17;
  506.   xr17 = ospl1;
  507.   yr27 = yr17;
  508.   yr17 = spl1;
  509.   );
  510.