updated the code

2 files changed, 53 insertions, 68 deletions

diff --git a/605/CH11/EX11.3/11_3.sce b/605/CH11/EX11.3/11_3.sce
index a721d2b24..518c81de5 100755
--- a/605/CH11/EX11.3/11_3.sce
+++ b/605/CH11/EX11.3/11_3.sce
@@ -1,30 +1,25 @@
-

-S11=0.614*exp(-%i*167.4*%pi/180)

-S21=2.187*exp(%i*32.4*%pi/180)

-S12=0.046*exp(%i*65*%pi/180)

-S22=0.716*exp(-%i*83*%pi/180)

-del=(S11*S22-S12*S21)

-B1=1+mod(S11)^2-mod(S22)^2-mod(del)^2

-B2=1+mod(S22)^2-mod(S11)^2-mod(del)^2

-C1=S11-Ss22*del

-C2=S22-Ss11*del

-function x=mod(n)

-    r=real(n)

-    i=imag(n)

-    x=sqrt(r^2+i^2)

-endfunction

-

-disp(mod(del),"del=")

-k=(1-mod(S11)^2-mod(S22)^2+mod(del)^2)/(2*mod(S12*S21))

-disp(k,"k=")

-FMS=(B1-sqrt(B1^2-4*mod(C1)^2))/(2*C1)

-disp(FMS,"FMS=")

-FML= (B2-sqrt(B2^2-4*mod(C2)^2))/(2*C2)

-disp(FML,"FML=")

-GTmax=mod(S21)/mod(S12)*(k-sqrt(k^2-1))

-printf("\nGTmax=%.3f",GTmax)

-printf("\nGTmax in dB=%.2f dB",10*log10(GTmax))

-

-

-

-

+S11=0.614*exp(-%i*167.4*%pi/180)
+S21=2.187*exp(%i*32.4*%pi/180)
+S12=0.046*exp(%i*65*%pi/180)
+S22=0.716*exp(-%i*83*%pi/180)
+del=(S11*S22-S12*S21)
+B1=1+abs(S11)^2-abs(S22)^2-abs(del)^2
+B2=1+abs(S22)^2-abs(S11)^2-abs(del)^2
+C1=S11-S22*del
+C2=S22-S11*del
+function x=mod(n)
+    r=real(n)
+    i=imag(n)
+    x=sqrt(r^2+i^2)
+endfunction
+
+disp(mod(del),"del=")
+k=(1-mod(S11)^2-mod(S22)^2+mod(del)^2)/(2*mod(S12*S21))
+disp(k,"k=")
+FMS=(B1-sqrt(B1^2-4*mod(C1)^2))/(2*C1)
+disp(FMS,"FMS=")
+FML= (B2-sqrt(B2^2-4*mod(C2)^2))/(2*C2)
+disp(FML,"FML=")
+GTmax=mod(S21)/mod(S12)*(k-sqrt(k^2-1))
+printf("\nGTmax=%.3f",GTmax)
+printf("\nGTmax in dB=%.2f dB",10*log10(GTmax))
\ No newline at end of file
diff --git a/605/CH9/EX9.7/9_7.sce b/605/CH9/EX9.7/9_7.sce
index fc55c51bc..7b455df66 100755
--- a/605/CH9/EX9.7/9_7.sce
+++ b/605/CH9/EX9.7/9_7.sce
@@ -1,38 +1,28 @@
-

-Gr=0.01

-m=acosh(sqrt(10^(0.1*L)-1)/sqrt(10^(0.1*Gr-1)))/acosh(W/Wc)

-printf("\nm=%.0f\n",m)

-m=3

-E=log(coth(Gr/17.37))

-X=sinh(E/2/m)

-n=3

-gp=1

-

-for p=1:1:n

-    ap=sin((2*p-1)*%pi/2/m)

-    bp=X^2+sin(p*%pi/m)^2

-    printf("\nap=%.4f\nbp=%.4f\n",ap,bp)

-end

-gp=0.62425

-printf("\ng0=g4=1")

-printf("\np=1\tgp=0.62425")

-for p=2:1:n

-    gp=4*sin((2*(p-1)-1)*%pi/2/m)*sin((2*p-1)*%pi/2/m)/(X^2+sin((p-1)*%pi/m)^2)/gp

-     printf("\np=%.0f\tgp=%.5f",p,gp)

-end

-

-printf("\nL1=L3=%.4e H\nC1=%.4e F",75*0.62425/(2*%pi*10^8),0.9662/(75*2*%pi*10^8))

-

-

-

-

-

-

-

-

-

-

-

-

-

-

+clc
+clear
+L = 5;
+W  =400;
+ Wc = 100;
+Gr=0.01
+m=acosh(sqrt(10^(0.1*L)-1)/sqrt(10^(0.1*Gr-1)))/acosh(W/Wc)
+printf("\nm=%.0f\n",m)
+m=3
+E=log(coth(Gr/17.37))
+X=sinh(E/2/m)
+n=3
+gp=1
+
+for p=1:1:n
+    ap=sin((2*p-1)*%pi/2/m)
+    bp=X^2+sin(p*%pi/m)^2
+    printf("\nap=%.4f\nbp=%.4f\n",ap,bp)
+end
+gp=0.62425
+printf("\ng0=g4=1")
+printf("\np=1\tgp=0.62425")
+for p=2:1:n
+    gp=4*sin((2*(p-1)-1)*%pi/2/m)*sin((2*p-1)*%pi/2/m)/(X^2+sin((p-1)*%pi/m)^2)/gp
+     printf("\np=%.0f\tgp=%.5f",p,gp)
+end
+
+printf("\nL1=L3=%.4e H\nC1=%.4e F",75*0.62425/(2*%pi*10^8),0.9662/(75*2*%pi*10^8))
\ No newline at end of file