From 8ac15bc5efafa2afc053c293152605b0e6ae60ff Mon Sep 17 00:00:00 2001
From: Siddharth Agarwal
Date: Tue, 3 Sep 2019 18:27:40 +0530
Subject: Xcos examples from textbooks and for blocks

---
 Working_Examples/154/DEPENDENCIES/ch15_9.sce | 33 ++++++++++++++++++++++++++++
 1 file changed, 33 insertions(+)
 create mode 100755 Working_Examples/154/DEPENDENCIES/ch15_9.sce

(limited to 'Working_Examples/154/DEPENDENCIES/ch15_9.sce')

diff --git a/Working_Examples/154/DEPENDENCIES/ch15_9.sce b/Working_Examples/154/DEPENDENCIES/ch15_9.sce
new file mode 100755
index 0000000..0d9f2fa
--- /dev/null
+++ b/Working_Examples/154/DEPENDENCIES/ch15_9.sce
@@ -0,0 +1,33 @@
+clc
+disp("Example 15.9")
+printf("\n")
+
+s=%s;
+printf("Given")
+disp("L1=0.2H L2=0.1H")
+disp("M=0.1H R=10ohm")
+disp("v1=u(t) a unit step function")
+L1=0.2;L2=0.1
+M=0.1;R=10;
+v1=1;
+w=100;
+//Let Input impedance be Z1 and can be calculated as
+//From the equations in 15.10
+disp("Z1(s)=L1*s-((M*s)^2)/(R+L2*s)")
+Z1=L1*s-(((M*s)^2)/(R+L2*s))
+//Proper rearranging of co-efficients
+Num=Z1('num')/0.01
+Den=Z1('den')*100
+
+disp(Num/Den,"Z1(s)")
+Y1=1/Z1
+disp(Den/Num,"Y1(s)")
+
+//As the input is unit step function the value is 1V for t>0
+//In exponential form the value is represented as exp(s*t) with s=0 as the pole of Y1(s)
+
+//Therefore forced response
+k=1/L1;
+printf("Forced response i1,f=(%d*t) (A)\n",k);
+
+
-- 
cgit