From c600ebcb67961fe6007ba41fd5ad987da3af7f6e Mon Sep 17 00:00:00 2001 From: Brijeshcr Date: Thu, 6 Jul 2017 15:48:47 +0530 Subject: Revert "LinearAlgebra Function Added" --- .../src/c/signalProcessing/filt_sinc/dfilt_sincs.c | 49 ++++++++++++++++++++++ 1 file changed, 49 insertions(+) create mode 100644 2.3-1/src/c/signalProcessing/filt_sinc/dfilt_sincs.c (limited to '2.3-1/src/c/signalProcessing/filt_sinc/dfilt_sincs.c') diff --git a/2.3-1/src/c/signalProcessing/filt_sinc/dfilt_sincs.c b/2.3-1/src/c/signalProcessing/filt_sinc/dfilt_sincs.c new file mode 100644 index 00000000..1b7d1b18 --- /dev/null +++ b/2.3-1/src/c/signalProcessing/filt_sinc/dfilt_sincs.c @@ -0,0 +1,49 @@ +/* Copyright (C) 2017 - IIT Bombay - FOSSEE + + This file must be used under the terms of the CeCILL. + This source file is licensed as described in the file COPYING, which + you should have received as part of this distribution. The terms + are also available at + http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt + Author: Ankit Raj + Organization: FOSSEE, IIT Bombay + Email: toolbox@scilab.in + */ +#include +#include +#include "filt_sinc.h" +void dfilt_sincs(double N,double fc,double* oup) +{ + double no2=(N-1)/2; + int ino2=(int)no2; + double wl=fc*2*M_PI; + int sz=N; + double xn[sz]; + double i; + int l,m; + for(i=-no2;i<=no2;i++) + { + l=i+no2; + xn[l]=sin(wl*i); + } + double xd[sz]; + double j; + printf("\n"); + for(j=-no2;j<=no2;j++) + { + m=j+no2; + xd[m]=M_PI*j; + } + if(ino2==no2) + { + xn[(int)no2]=2*fc; + xd[(int)no2]=1; + } + double x[sz]; + int k; + for(k=0;k