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

---
 40/CH9/EX9.10a/Exa_9_10a.jpg | Bin 0 -> 16092 bytes
 40/CH9/EX9.10a/Exa_9_10a.sce |  23 +++++++++++++++++++++++
 2 files changed, 23 insertions(+)
 create mode 100755 40/CH9/EX9.10a/Exa_9_10a.jpg
 create mode 100755 40/CH9/EX9.10a/Exa_9_10a.sce

(limited to '40/CH9/EX9.10a')

diff --git a/40/CH9/EX9.10a/Exa_9_10a.jpg b/40/CH9/EX9.10a/Exa_9_10a.jpg
new file mode 100755
index 000000000..5c242b942
Binary files /dev/null and b/40/CH9/EX9.10a/Exa_9_10a.jpg differ
diff --git a/40/CH9/EX9.10a/Exa_9_10a.sce b/40/CH9/EX9.10a/Exa_9_10a.sce
new file mode 100755
index 000000000..3921304f2
--- /dev/null
+++ b/40/CH9/EX9.10a/Exa_9_10a.sce
@@ -0,0 +1,23 @@
+//Bilinear design of second order filters
+s=%s;z=%z;
+fo=6;Wo=2*%pi*fo/25;
+delf=5;S=25;
+B=cos(2*%pi*fo/25)
+C=tan(%pi*delf/25)
+HS=1/(s+1);
+HZ=horner(HS,(z^2-(2*B*z)+1)/(C*(z^2)-C))
+f=0:0.5:12.5;
+HZ1=horner(HZ,exp(%i*2*%pi*f'/25));
+HZ1=abs(HZ1);
+W2=(%pi*delf/25)+acos(cos(Wo)*cos(%pi*delf/25))
+W1=W2-(2*%pi*delf/25)
+f1=S*W1/(2*%pi),f2=S*W2/(2*%pi)
+f3=[f1;fo;f2];
+HZf=abs(horner(HZ,exp(%i*2*%pi*f3'/25)));
+a=gca();
+a.x_location="origin";
+plot2d(f,HZ1,rect=[0 0 13 1]);
+plot2d3('gnn',f3,HZf);
+xlabel('Analog Frequency f[kHZ]');
+ylabel('Magnitude');
+xtitle('Band pass filter fo=6kHZ,delf=5kHZ');
-- 
cgit