summaryrefslogtreecommitdiff
path: root/2.3-1/src/c/signalProcessing/zpch2
diff options
context:
space:
mode:
Diffstat (limited to '2.3-1/src/c/signalProcessing/zpch2')
-rw-r--r--2.3-1/src/c/signalProcessing/zpch2/dzpch2s.c63
-rw-r--r--2.3-1/src/c/signalProcessing/zpch2/int_zpch2.h18
-rw-r--r--2.3-1/src/c/signalProcessing/zpch2/zpch2.h28
3 files changed, 0 insertions, 109 deletions
diff --git a/2.3-1/src/c/signalProcessing/zpch2/dzpch2s.c b/2.3-1/src/c/signalProcessing/zpch2/dzpch2s.c
deleted file mode 100644
index 22a46a64..00000000
--- a/2.3-1/src/c/signalProcessing/zpch2/dzpch2s.c
+++ /dev/null
@@ -1,63 +0,0 @@
-/* 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<stdio.h>
-#include<math.h>
-#include "zpch2.h"
-#include "multiplication.h"
-#include "abs.h"
-#include "division.h"
-#define PI 3.14159265
-double dzpch2s(double N, double A, double omegar,doubleComplex* zeros,doubleComplex* poles)
-{
- int k,j;
- double e,xk,realVal,gain;
- doubleComplex accumulate1,accumulate2,tp1,tp2;
- accumulate1=DoubleComplex(1,0);
- accumulate2=DoubleComplex(1,0);
- int ct=0;
- for(k=1;k<=N;k++)
- {
- if(k!=((N+1)/2))
- {
- ct++;
- double num=PI*(2*k-1);
- double den=2*N;
- xk=cos(num/den);
- double sk=omegar/xk;
- zeros[k-1]=DoubleComplex(0,sk);
- tp1=DoubleComplex(0,sk);
- accumulate1=zmuls(accumulate1,tp1);
- }
- }
- for(j=1;j<=N;j++)
- {
- double num=PI*(2*j-1);
- double den=2*N;
- double xk1=num/den;
- double Gamma=pow((A+sqrt(A*A-1)),(1/N));
- double alpha=-((Gamma-1/Gamma)/2)*sin(xk1);
- double Beta=((Gamma+1/Gamma)/2)*cos(xk1);
- double normal=alpha*alpha+Beta*Beta;
- poles[j-1]=DoubleComplex((omegar*alpha)/normal,(omegar*Beta)/normal);
- tp2=DoubleComplex((omegar*alpha)/normal,(omegar*Beta)/normal);
- accumulate2=zmuls(accumulate2,tp2);
- }
- double qt;
- double ra2=zreals(accumulate2);
- double ra1=zreals(accumulate1);
- qt=dldivs(ra1,ra2);
- gain=dabss(qt);
-
- return gain;
-
-}
-
diff --git a/2.3-1/src/c/signalProcessing/zpch2/int_zpch2.h b/2.3-1/src/c/signalProcessing/zpch2/int_zpch2.h
deleted file mode 100644
index 683a7614..00000000
--- a/2.3-1/src/c/signalProcessing/zpch2/int_zpch2.h
+++ /dev/null
@@ -1,18 +0,0 @@
-/* Copyright (C) 2016 - 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
- */
-
-#ifndef __INT_ZPCH2_H__
-#define __INT_ZPCH2_H__
-
-#define d0d0d0zpch2z2z2(N,A,omegar,zeros,poles) dzpch2s(N,A,omegar,zeros,poles)
-
-#endif /* !INT_ZPCH2_H__! */
diff --git a/2.3-1/src/c/signalProcessing/zpch2/zpch2.h b/2.3-1/src/c/signalProcessing/zpch2/zpch2.h
deleted file mode 100644
index 0d191d5c..00000000
--- a/2.3-1/src/c/signalProcessing/zpch2/zpch2.h
+++ /dev/null
@@ -1,28 +0,0 @@
-/* Copyright (C) 2016 - 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
- */
-
-#ifndef __ZPCH2_H__
-#define __ZPCH2_H__
-#include "types.h"
-#include "doubleComplex.h"
-
-#ifdef __cplusplus
-extern "C" {
-#endif
-
-void dzpch2s(double N, double A, double omegar,doubleComplex* zeros,doubleComplex* poles);
-
-#ifdef __cplusplus
-} /* extern "C" */
-#endif
-
-#endif /* __ZPCH2_H__ */