From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 1430/CH13/EX13.6/exa13_6.jpg | Bin 0 -> 21318 bytes
 1430/CH13/EX13.6/exa13_6.sce |  25 +++++++++++++++++++++++++
 2 files changed, 25 insertions(+)
 create mode 100644 1430/CH13/EX13.6/exa13_6.jpg
 create mode 100644 1430/CH13/EX13.6/exa13_6.sce

(limited to '1430/CH13/EX13.6')

diff --git a/1430/CH13/EX13.6/exa13_6.jpg b/1430/CH13/EX13.6/exa13_6.jpg
new file mode 100644
index 000000000..64df3be12
Binary files /dev/null and b/1430/CH13/EX13.6/exa13_6.jpg differ
diff --git a/1430/CH13/EX13.6/exa13_6.sce b/1430/CH13/EX13.6/exa13_6.sce
new file mode 100644
index 000000000..e0456e7c4
--- /dev/null
+++ b/1430/CH13/EX13.6/exa13_6.sce
@@ -0,0 +1,25 @@
+// Example 13.6
+// Inversion with complex Poles
+s=%s;
+t=0:0.001:10
+num=15*s^2-16*s-7;
+den=(s+2)*(s^2+6*s+25);
+F_s=num/den;
+pfe=pfss(F_s); // partial fraction of the transfer function
+// from pfe(1) we get
+B=10;
+C=-66;
+alpha=3;// from pfe(1)
+beta=sqrt(25-9);//Comparing the denominator of pfe(1) with standard 2nd orderequation
+// Now
+K=B+(%i*(alpha*B-C))/beta;
+// From inverse Laplace Transfrom of pfe(2) we get 
+f1=5*exp(-2*t)
+K_m=abs(K); // Magnitude of K
+phase_K=atan(imag(K),real(K));
+g=K_m*exp(-alpha*t).*cos(beta*t+phase_K);
+f=f1+g;
+plot(t,f)
+xlabel('t')
+ylabel('f(t)')
+title('Function Waveform')
-- 
cgit