updated the code

3 files changed, 38 insertions, 35 deletions

diff --git a/773/CH17/EX17.03/17_03.sci b/773/CH17/EX17.03/17_03.sci
index 3681df84b..66b3bbdd7 100755
--- a/773/CH17/EX17.03/17_03.sci
+++ b/773/CH17/EX17.03/17_03.sci
@@ -1,7 +1,7 @@
-//function//

-s=%s;

-//Creating cont-time transfer function

-TFcont=syslin('c',3/(s^4+(2*s^3)+(3*s)+2)) 

-SScont=tf2ss(TFcont)

-//CCF form

-[Ac,Bc,U,ind]=canon(SScont(2),SScon(3))

+//function//
+s=%s;
+//Creating cont-time transfer function
+TFcont=syslin('c',3/(s^4+(2*s^3)+(3*s)+2)) 
+SScont=tf2ss(TFcont)
+//CCF form
+[Ac,Bc,U,ind]=canon(SScont(2),SScont(3))
\ No newline at end of file
diff --git a/773/CH18/EX18.08.01/18_08_01.sci b/773/CH18/EX18.08.01/18_08_01.sci
index 3a33aec94..2bdf13a3d 100755
--- a/773/CH18/EX18.08.01/18_08_01.sci
+++ b/773/CH18/EX18.08.01/18_08_01.sci
@@ -1,5 +1,7 @@
-//symsum//

-exec ztransfer.sce;

-sequence=[0 2 0 0 -3 0 0 8]

-y=ztransfer(sequence);

-disp(y,"ans=")

+//symsum//
+ function [Ztransfer]= ztransfer_new (sequence,n)
+  z = poly (0, 'z' , 'r' );  Ztransfer = sequence *(1/z )^n'
+endfunction
+sequence=[0 2 0 0 -3 0 0 8]
+y=ztransfer(sequence);
+disp(y,"ans=")
\ No newline at end of file
diff --git a/773/CH8/EX8.17/8_17.sci b/773/CH8/EX8.17/8_17.sci
index 43f56c3f2..4a7f03293 100755
--- a/773/CH8/EX8.17/8_17.sci
+++ b/773/CH8/EX8.17/8_17.sci
@@ -1,23 +1,24 @@
-//coefficient//

-s=%s;

-F=syslin('c',(25)/(s^2+7*s)); //Creates transfer function in forward path

-B=syslin('c',(1+0*s)/(1+0*s)); //Creates transfer function in backward path

-k=20/25;   //k=gain factor

-CL=k*(F/.B) //Calculates closed-loop transfer function

-// compare CL with Wn^2/(s^2+2*zeta*Wn+Wn^2)

-y=denom(CL) //extracting the denominator of CL

-z=coeff(y)  //extracting the coefficients of the denominator polynomial

-//Wn^2=z(1,1) ,comparing the coefficients

-Wn=sqrt(z(1,1))            // Wn=natural frequency

-//2*zeta*Wn=z(1,2)

-zeta=z(1,2)/(2*Wn)             // zeta=damping factor

-Wd=Wn*sqrt(1-zeta^2)

-Tp=%pi/Wd     //Tp=peak time

-Mp=100*exp((-%pi*zeta)/sqrt(1-zeta^2))  //peak overshoot

-Td=(1+0.7*zeta)/Wn      //Td=delay time

-a=atan(sqrt(1-zeta^2)/zeta)

-Tr=(%pi-a)/Wd           //Tr=rise time

-Ts=4/(zeta*Wn)          //Ts=settling time

-//y(t)=expression for output

-y=(20/25)*(1-(exp(-1*zeta*Wn*t)/sqrt(1-zeta^2))*sin(Wd*t+atan(zeta/sqrt(1-zeta^2))));

-disp(y,"Y(t)")

+//coefficient//
+s=%s;
+t = %t;
+F=syslin('c',(25)/(s^2+7*s)); //Creates transfer function in forward path
+B=syslin('c',(1+0*s)/(1+0*s)); //Creates transfer function in backward path
+k=20/25;   //k=gain factor
+CL=k*(F/.B) //Calculates closed-loop transfer function
+// compare CL with Wn^2/(s^2+2*zeta*Wn+Wn^2)
+y=denom(CL) //extracting the denominator of CL
+z=coeff(y)  //extracting the coefficients of the denominator polynomial
+//Wn^2=z(1,1) ,comparing the coefficients
+Wn=sqrt(z(1,1))            // Wn=natural frequency
+//2*zeta*Wn=z(1,2)
+zeta=z(1,2)/(2*Wn)             // zeta=damping factor
+Wd=Wn*sqrt(1-zeta^2)
+Tp=%pi/Wd     //Tp=peak time
+Mp=100*exp((-%pi*zeta)/sqrt(1-zeta^2))  //peak overshoot
+Td=(1+0.7*zeta)/Wn      //Td=delay time
+a=atan(sqrt(1-zeta^2)/zeta)
+Tr=(%pi-a)/Wd           //Tr=rise time
+Ts=4/(zeta*Wn)          //Ts=settling time
+//y(t)=expression for output
+y=(20/25)*(1-(exp(-1*zeta*Wn*t)/sqrt(1-zeta^2))*sin(Wd*t+atan(zeta/sqrt(1-zeta^2))));
+disp(y,"Y(t)")
\ No newline at end of file