Trying to have high-pass simple filter

2 files changed, 22 insertions, 0 deletions

diff --git a/scilab2c/tests/unit_tests/test004_Filter/scilabcode/high_pass.sci b/scilab2c/tests/unit_tests/test004_Filter/scilabcode/high_pass.sci
new file mode 100644
index 00000000..a3af8322
--- /dev/null
+++ b/scilab2c/tests/unit_tests/test004_Filter/scilabcode/high_pass.sci
@@ -0,0 +1,11 @@
+// Return RC high-pass filter output_signal samples, 

+// given input_signal samples,

+// time interval dt,

+// R and C

+function output_signal = high_pass(input_signal, dt, R, C)

+  alpha = R * C / (R * C + dt)

+  output_signal(1) = input_signal(1)

+  for i = 2 : size(input_signal, "*")

+    output_signal(i) = alpha * output_signal(i-1) + alpha * (input_signal(i) - input_signal(i-1))

+  end

+endfunction
\ No newline at end of file
diff --git a/scilab2c/tests/unit_tests/test004_Filter/scilabcode/main.sci b/scilab2c/tests/unit_tests/test004_Filter/scilabcode/main.sci
new file mode 100644
index 00000000..9b32b2bd
--- /dev/null
+++ b/scilab2c/tests/unit_tests/test004_Filter/scilabcode/main.sci
@@ -0,0 +1,11 @@
+function main()
+  sup = 5000;
+  M = 50 * rand(sup,1);
+  //scf(1);
+  //plot(1:sup, M);
+  
+  
+  N = high_pass(M, 0.1, 1, 1);
+  //scf(2);
+  //plot(1:sup, N);
+endfunction
\ No newline at end of file