From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 758/CH8/EX8.11/Ex_8_11.sce | 38 ++++++++++++++++++++++++++++++++++++++
 1 file changed, 38 insertions(+)
 create mode 100755 758/CH8/EX8.11/Ex_8_11.sce

(limited to '758/CH8/EX8.11/Ex_8_11.sce')

diff --git a/758/CH8/EX8.11/Ex_8_11.sce b/758/CH8/EX8.11/Ex_8_11.sce
new file mode 100755
index 000000000..2ece43744
--- /dev/null
+++ b/758/CH8/EX8.11/Ex_8_11.sce
@@ -0,0 +1,38 @@
+//Example 8.11
+clc;clear;close;
+rp=0.707                       //passband ripple
+rs=0.2                        //stopband ripple
+wp=%pi/2;                     //passband frequency
+ws=3*%pi/4;                   //stopband frequency
+T=1;
+fp=wp/T;
+fs=ws/T;
+s=poly(0,'s');
+z=poly(0,'z');
+hs=1;
+//Calculating the order of filter
+num=log((rs^-2 -1)/(rp^-2 -1));
+den=2*log(fs/fp);
+N=ceil(num/den);
+
+//Calculation of cut-off frequency
+fc=fp/(rp^-2 -1)^(0.5/N);
+
+//Calculating filter response
+if modulo(N,2)==1 then
+    b=-2*sin(%pi/(2*N));
+    hs=hs*fc/(s+fc);
+end
+for k=1:N/2
+    b=2*sin((2*k-1)*%pi/(2*N));
+    hs=hs*fc^2/(s^2+b*fc*s+fc^2);
+end
+hs=clean(hs);
+sys=syslin('c',hs);                
+hz=horner(ss2tf(dscr(sys,T)),1/z);        //converting H(s) to H(z)
+
+//Displaying filter response
+[hzm,fr]=frmag(hz,256);
+disp(hz,'Filter Transfer function: ');
+plot(fr,hzm);
+title('Lowpass Butterworth Filter Response');ylabel('Amplitude-->');xlabel('Normalised frequency f/fs-->');
\ No newline at end of file
-- 
cgit