9 files changed, 245 insertions, 0 deletions
diff --git a/2414/CH13/EX13.1/Ex13_1.sce b/2414/CH13/EX13.1/Ex13_1.sce
new file mode 100755
index 000000000..0d223b246
--- /dev/null
+++ b/2414/CH13/EX13.1/Ex13_1.sce
@@ -0,0 +1,31 @@
+clc;
+close();
+clear();
+//page no 442
+//prob no. 13.1
+Gb=1;
+mprintf('(a) SSB:  Gb=%i  \n',Gb);
+mprintf('   GbdB=%i  dB\n',10*log10(Gb));
+mprintf('(b) DSB:  Gb=%i  \n',Gb);
+mprintf('   GbdB=%i  dB\n',10*log10(Gb));
+m=0.5;
+Gb=m^2/(2+m^2);
+mprintf('(c) AM(m=.5):  Gb=%.3f  \n',Gb);
+mprintf('   GbdB=%.3f  dB\n',10*log10(Gb));
+m=1;
+Gb=m^2/(2+m^2);
+mprintf('(d) AM(m=1):  Gb=%.3f  \n',Gb);
+mprintf('   GbdB=%.3f  dB\n',10*log10(Gb));
+delta_phi=5;
+Gb=delta_phi^2/2;
+mprintf('(e) FM(delta phi=5)):  Gb=%.1f  \n',Gb);
+mprintf('   GbdB=%.3f  dB\n',10*log10(Gb));
+D=5;
+Gb=3*D^2/2;
+mprintf('(f) FM(D=5):  Gb=%.1f  \n',Gb);
+mprintf('   GbdB=%.3f  dB\n',10*log10(Gb));
+Wf1=7.07;
+Gb=3/2*D^2*%pi/6*Wf1;
+mprintf('(g) FM(D=5, W/f1=7.07):  Gb=%.1f  \n',Gb);
+mprintf('   GbdB=%.2f  dB\n',10*log10(Gb));
+
diff --git a/2414/CH13/EX13.2/Ex13_2.sce b/2414/CH13/EX13.2/Ex13_2.sce
new file mode 100755
index 000000000..c353d8401
--- /dev/null
+++ b/2414/CH13/EX13.2/Ex13_2.sce
@@ -0,0 +1,32 @@
+clc;
+close();
+clear();
+//page no 443
+//prob no. 13.2
+GR=1;
+mprintf('(a) SSB:  GR=%i  \n',GR);
+mprintf('   GRdB=%i  dB\n',10*log10(GR));
+GR=2;
+mprintf('(b) DSB:  GR=%i  \n',GR);
+mprintf('   GRdB=%.2f  dB\n',10*log10(GR));
+m=0.5;
+GR=2*m^2/(2+m^2);
+mprintf('(c) AM(m=.5):  GR=%.4f  \n',GR);
+mprintf('   GRdB=%.3f  dB\n',10*log10(GR));
+m=1;
+GR=2*m^2/(2+m^2);
+mprintf('(d) AM(m=1):  GR=%.3f  \n',GR);
+mprintf('   GRdB=%.2f  dB\n',10*log10(GR));
+delta_phi=5;
+GR=(1+delta_phi)*delta_phi^2;
+mprintf('(e) FM(delta phi=5)):  GR=%.1f  \n',GR);
+mprintf('   GRdB=%.3f  dB\n',10*log10(GR));
+D=5;
+GR=3*D^2*(1+D);
+mprintf('(f) FM(D=5):  GR=%.1f  \n',GR);
+mprintf('   GRdB=%.3f  dB\n',10*log10(GR));
+Wf1=7.07;
+GR=3*(1+D)*D^2*%pi/6*Wf1;
+mprintf('(g) FM(D=5, W/f1=7.07):  GR=%.1f  \n',GR);
+mprintf('   GRdB=%.2f  dB\n',10*log10(GR));
+
diff --git a/2414/CH13/EX13.3/Ex13_3.sce b/2414/CH13/EX13.3/Ex13_3.sce
new file mode 100755
index 000000000..337ad9ad9
--- /dev/null
+++ b/2414/CH13/EX13.3/Ex13_3.sce
@@ -0,0 +1,31 @@
+clc;
+close();
+clear();
+//page no 447
+//prob no. 13.3
+k=1.38*10^-23;   //Boltzmann's const
+//Temperatures in K
+Ti=150;
+Te=325;
+Tsys=Ti+Te;
+mprintf(' Tsys=%i K \n',Tsys);
+D=5;
+W=15;   //kHz
+B=2*(1+D)*W;
+mprintf('  B=%i  kHz\n',B);
+Nsys=k*Tsys*B*10^3;   //W
+mprintf('  Nsys=%.3f fW\n',Nsys*10^15);
+PR=50*10^-12;   //W
+SNsys=PR/Nsys; 
+mprintf('  (S/N)sys=%i  \n',SNsys);
+GR=3*(1+D)*D^2
+mprintf('  GR=%.0f  \n',GR);
+SNoutput=GR*SNsys;
+mprintf('  (S/N)output=%.0f  \n',SNoutput);
+mprintf('  (S/N)out,dB=%.2f  dB\n',10*log10(SNoutput));
+mprintf('  (S/N)sys,dB=%.2f  dB\n',10*log10(SNsys));
+GRdb=10*log10(GR);
+mprintf('  GR,dB=%.2f dB \n',GRdb);
+mprintf('  (S/N)output,dB=%.2f  dB\n',10*log10(SNoutput));
+
+
diff --git a/2414/CH13/EX13.4/Ex13_4.sce b/2414/CH13/EX13.4/Ex13_4.sce
new file mode 100755
index 000000000..584b404de
--- /dev/null
+++ b/2414/CH13/EX13.4/Ex13_4.sce
@@ -0,0 +1,8 @@
+clc;
+close();
+clear();
+//page no 450
+//prob no. 13.4
+N=16;  //bit
+SNoutdB=1.76+6.02*N;
+mprintf(' (S/N)output,dB=%.2f dB \n',SNoutdB);
diff --git a/2414/CH13/EX13.5/Ex13_5.sce b/2414/CH13/EX13.5/Ex13_5.sce
new file mode 100755
index 000000000..d5f958e0b
--- /dev/null
+++ b/2414/CH13/EX13.5/Ex13_5.sce
@@ -0,0 +1,11 @@
+clc;
+close();
+clear();
+//page no 450
+//prob no. 13.5
+SNoutdB=53;
+N=(SNoutdB-1.76)/6.02;
+mprintf(' N=%.2f bits \n',N);
+N=9;   //roundup
+mprintf(' N=%i bits \n',N);
+
diff --git a/2414/CH13/EX13.6/Ex13_6.sce b/2414/CH13/EX13.6/Ex13_6.sce
new file mode 100755
index 000000000..9bb572c1d
--- /dev/null
+++ b/2414/CH13/EX13.6/Ex13_6.sce
@@ -0,0 +1,39 @@
+clc;
+close();
+clear();
+//page no 453
+//prob no. 13.6
+N=6;   //bits per word
+M=2^N;
+mprintf(' M=%i  \n',M);
+Pr=200*10^-15;   //W
+R=2*10^6;   //bits/s
+Eb=Pr/R;
+mprintf(' Bit energy ,Eb=%.0f*10^-21  \n',Eb*10^21);
+k=1.38*10^-23;   //Boltzmann cons
+Ti=300;  //K
+Te=425;  //K
+Tsys=Ti+Te;
+nsys=k*Tsys;                 
+mprintf('  Noise power spectral density ,nsys=%.0f*10^-20 W/Hz \n',nsys*10^20);
+rho=Eb/nsys;
+mprintf(' Bit energy , rho=%.0f  \n',rho);
+rhodB=10*log10(rho);
+mprintf(' Bit energy in db, rho,dB=%.0f dB \n',rhodB);
+//part a
+Pe=4*10^-6;
+SNout=1.5*M^2/(1+4*M^2*Pe);
+mprintf(' \n(a) (S/N)output=%.0f (or %0.2f  dB) \n',SNout,10*log10(SNout));
+//part b
+Pe=2.3*10^-5;
+SNout=1.5*M^2/(1+4*M^2*Pe);
+mprintf(' \n(b) (S/N)output=%.0f (or %0.2f  dB) \n',SNout,10*log10(SNout));
+//part c
+Pe=8*10^-4;
+SNout=1.5*M^2/(1+4*M^2*Pe);
+mprintf(' \n(c) (S/N)output=%.1f (or %0.2f  dB) \n',SNout,10*log10(SNout));
+//part d
+Pe=3.5*10^-3;
+SNout=1.5*M^2/(1+4*M^2*Pe);
+mprintf(' \n(d) (S/N)output=%.1f (or %0.2f  dB) \n',SNout,10*log10(SNout));
+
diff --git a/2414/CH13/EX13.7/Ex13_7.sce b/2414/CH13/EX13.7/Ex13_7.sce
new file mode 100755
index 000000000..4c6882d17
--- /dev/null
+++ b/2414/CH13/EX13.7/Ex13_7.sce
@@ -0,0 +1,25 @@
+clc;
+close();
+clear();
+//page no 455
+//prob no. 13.7
+//data from ex 13.6
+M=2^6;
+Pr=200*10^-15;   //W
+R=8*10^6;   //bits/s    (changed)
+Eb=Pr/R;
+
+k=1.38*10^-23;   //Boltzmann cons
+Ti=300;  //K
+Te=425;  //K
+Tsys=Ti+Te;
+nsys=k*Tsys;                 
+//mprintf('  Noise power spectral density ,nsys=%.0f*10^-20 W/Hz \n',nsys*10^20);
+rho=Eb/nsys;
+mprintf(' Bit energy , rho=%.1f  \n',rho);
+rhodB=10*log10(rho);
+mprintf(' Bit energy in db, rho,dB=%.2f dB \n',rhodB);
+
+Pe=1.3*10^-2;
+SNout=1.5*M^2/(1+4*M^2*Pe);
+mprintf(' \n (S/N)output=%.2f (or %0.2f  dB) \n',SNout,10*log10(SNout));
diff --git a/2414/CH13/EX13.8/Ex13_8.sce b/2414/CH13/EX13.8/Ex13_8.sce
new file mode 100755
index 000000000..6c7b6cac5
--- /dev/null
+++ b/2414/CH13/EX13.8/Ex13_8.sce
@@ -0,0 +1,41 @@
+clc;
+close();
+clear();
+//page no 455
+//prob no. 13.8
+Pe=10^-5;
+R=1*10^6;  //bits/s
+k=1.38*10^-23;   //Boltzmann cons
+Ti=475;  //K
+Te=250;  //K
+Tsys=Ti+Te;
+nsys=k*Tsys;  //W/Hz  
+function Eb=E(rhodb)    //function for Eb
+    rho=10^(rhodb/10);
+    Eb=nsys*rho;
+endfunction
+function Pr=P(E)      //function for Pr
+    Pr=R*Eb;
+endfunction
+function display(rhodb,pt)
+    Eb=E(rhodb);
+    Pr=P(E);
+mprintf('\n(%c)Bit energy , Eb=%.2f*10^-21 J \n',pt,Eb*10^21);
+mprintf(' Required reciver carrier power , Pr=%.2f fW \n',Pr*10^15);
+
+endfunction
+//Part a
+rhodb=9.6;
+display(rhodb,'a');
+
+//Part b
+rhodb=10.3;
+display(rhodb,'b');
+
+//Part c
+rhodb=12.6;
+display(rhodb,'c');
+
+//Part d
+rhodb=13.4;
+display(rhodb,'d');
diff --git a/2414/CH13/EX13.9/Ex13_9.sce b/2414/CH13/EX13.9/Ex13_9.sce
new file mode 100755
index 000000000..5e8eb4a5c
--- /dev/null
+++ b/2414/CH13/EX13.9/Ex13_9.sce
@@ -0,0 +1,27 @@
+clc;
+close();
+clear();
+//page no 456
+//prob no. 13.9
+
+//Data form ex13.8
+Pe=10^-5;
+R=2*10^6;  //bits/s  (changed)
+k=1.38*10^-23;   //Boltzmann cons
+Ti=475;  //K
+Te=250;  //K
+Tsys=Ti+Te;
+nsys=k*Tsys;  //W/Hz  
+function Eb=E(rhodb)    //function for Eb
+    rho=10^(rhodb/10);
+    Eb=nsys*rho;
+endfunction
+function Pr=P(E)      //function for Pr
+    Pr=R*Eb;
+endfunction
+
+rhodb=9.6;
+ Eb=E(rhodb);
+ Pr=P(E);
+mprintf('\nBit energy , Eb=%.2f*10^-21 J \n',Eb*10^21);
+mprintf(' Required reciver carrier power , Pr=%.2f fW \n',Pr*10^15);