From efd4b2645e3fd3c3cdd7b65501a0e7ae26d03c59 Mon Sep 17 00:00:00 2001
From: Ankit Raj
Date: Wed, 21 Jun 2017 11:00:07 +0530
Subject: Functions added - string related and signal processing

---
 src/c/signalProcessing/zpch1/dzpch1s.c | 53 ++++++++++++++++++++++++++++++++++
 1 file changed, 53 insertions(+)
 create mode 100644 src/c/signalProcessing/zpch1/dzpch1s.c

(limited to 'src/c/signalProcessing/zpch1/dzpch1s.c')

diff --git a/src/c/signalProcessing/zpch1/dzpch1s.c b/src/c/signalProcessing/zpch1/dzpch1s.c
new file mode 100644
index 0000000..ab7b633
--- /dev/null
+++ b/src/c/signalProcessing/zpch1/dzpch1s.c
@@ -0,0 +1,53 @@
+/* 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 "zpch1.h"
+#include "multiplication.h"
+#include "conj.h"
+#include "abs.h"
+#define PI 3.14159265
+double dzpch1s(double N,double e,double wc,doubleComplex* out)
+{
+	doubleComplex accumulate,tp;
+      double B,r,R,gain,realVal;
+      double temp=sqrt(1+e*e);
+    
+    B=pow((temp+1)/e,1/N);
+    
+    r=wc*((B*B-1)/(2*B));
+    R=wc*((B*B+1)/(2*B));
+    
+    int k;
+	double t1=1;
+	double t2=0;
+	accumulate=DoubleComplex(t1,t2);
+    for(k=0;k<N;k++)
+    {
+      double theta;
+    theta=(PI/2)+((2*k+1)*PI)/(2*N);
+      double xk,yk;
+    xk=r*cos(theta);
+    yk=R*sin(theta);
+    out[k]=DoubleComplex(xk,yk);
+	tp=DoubleComplex(xk,yk);
+    accumulate=zmuls(accumulate,tp);
+    }
+	
+    realVal=zreals(accumulate);
+	gain=dabss(realVal);
+    if (N==2*(int)(N/2)) 
+            gain=gain/sqrt(1+e*e);
+return gain;
+		
+}
+
-- 
cgit