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/3432/CH3/EX3.18/Ex3_18.sce    |  71 +++++++++++++++++++++++++
 Working_Examples/3432/CH3/EX3.18/Ex3_18_f1.pdf | Bin 0 -> 28707 bytes
 Working_Examples/3432/CH3/EX3.18/Ex3_18_f3.pdf | Bin 0 -> 26282 bytes
 3 files changed, 71 insertions(+)
 create mode 100755 Working_Examples/3432/CH3/EX3.18/Ex3_18.sce
 create mode 100755 Working_Examples/3432/CH3/EX3.18/Ex3_18_f1.pdf
 create mode 100755 Working_Examples/3432/CH3/EX3.18/Ex3_18_f3.pdf

(limited to 'Working_Examples/3432/CH3/EX3.18')

diff --git a/Working_Examples/3432/CH3/EX3.18/Ex3_18.sce b/Working_Examples/3432/CH3/EX3.18/Ex3_18.sce
new file mode 100755
index 0000000..311eb6f
--- /dev/null
+++ b/Working_Examples/3432/CH3/EX3.18/Ex3_18.sce
@@ -0,0 +1,71 @@
+//Example 3.18 Satellite Transfer Function
+
+xdel(winsid())//close all graphics Windows
+clear;
+clc;
+//------------------------------------------------------------------
+//(a)
+//Given
+d=1 //meters
+I=5000 //Kg-meter^2
+
+//Coefficients of numerator and denominator of the transfer function
+// of satellite
+numG=[d/I 0];
+denG=[0 0 1];
+
+// Transfer function
+Ns=poly(numG,'s','coeff');
+Ds=poly(denG,'s','coeff');
+sysG=syslin('c',Ns/Ds);
+t=0:0.01:10;
+[i j]=size(t);
+
+//------------------------------------------------------------------
+//(b)
+// Thrust input after 5 sec.
+u=zeros(1,j);
+w=find(t>=5 & t<=5+0.1);
+u(w)=25;
+plot(t,u);
+exec .\fig_settings.sci; //custom script for setting figure properties
+title("Transient response of the satellite...
+ (a) Thrust input",'fontsize',3);
+xlabel('Time t (sec.)','fontsize',2)
+ylabel('Fc','fontsize',2)
+
+//Transient response of the satellite to the thrust input as a pulse 
+sysd=dscr(sysG,0.01); //sample data system model
+y=flts(u,sysd);        //impulse response 
+figure, plot(t,y*180/%pi);
+exec .\fig_settings.sci; //custom script for setting figure properties
+title("Transient response of the satellite(double-pulse)...
+ (b) satellite attitude",'fontsize',3);
+xlabel('Time t (sec.)','fontsize',2)
+ylabel('$\theta(deg)$','fontsize',2)
+//------------------------------------------------------------------
+// Thrust input double-pulse.
+u=zeros(1,j);
+w1=find(t>=5 & t<=5+0.1);
+u(w1)=25;
+w2=find(t>=6.1 & t<=6.1+0.1);
+u(w2)=-25;
+figure,
+plot(t,u);
+exec .\fig_settings.sci; //custom script for setting figure properties
+title("Transient response of the satellite (double-pulse)...
+ (a) Thrust input",'fontsize',3);
+xlabel('Time t (sec.)','fontsize',2)
+ylabel('Fc','fontsize',2)
+
+//Transient response of the satellite to the thrust input as a pulse 
+sysd=dscr(sysG,0.01); //sample data system model
+y=flts(u,sysd);        //impulse response 
+figure, plot(t,y*180/%pi);
+exec .\fig_settings.sci; //custom script for setting figure properties
+title("Transient response of the satellite(double-pulse)...
+ (b) satellite attitude",'fontsize',3);
+xlabel('Time t (sec.)','fontsize',2)
+ylabel('$\theta(deg)$','fontsize',2)
+
+//------------------------------------------------------------------
diff --git a/Working_Examples/3432/CH3/EX3.18/Ex3_18_f1.pdf b/Working_Examples/3432/CH3/EX3.18/Ex3_18_f1.pdf
new file mode 100755
index 0000000..b79b5d1
Binary files /dev/null and b/Working_Examples/3432/CH3/EX3.18/Ex3_18_f1.pdf differ
diff --git a/Working_Examples/3432/CH3/EX3.18/Ex3_18_f3.pdf b/Working_Examples/3432/CH3/EX3.18/Ex3_18_f3.pdf
new file mode 100755
index 0000000..ef8fa00
Binary files /dev/null and b/Working_Examples/3432/CH3/EX3.18/Ex3_18_f3.pdf differ
-- 
cgit