diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /52 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '52')
202 files changed, 2607 insertions, 0 deletions
diff --git a/52/CH1/EX1.1/Example1_1.jpg b/52/CH1/EX1.1/Example1_1.jpg Binary files differnew file mode 100755 index 000000000..978563b89 --- /dev/null +++ b/52/CH1/EX1.1/Example1_1.jpg diff --git a/52/CH1/EX1.1/Example1_1.sce b/52/CH1/EX1.1/Example1_1.sce new file mode 100755 index 000000000..ee8718380 --- /dev/null +++ b/52/CH1/EX1.1/Example1_1.sce @@ -0,0 +1,19 @@ +//Example 1.1
+//Sketch the continuous time signal x(t)=2*exp(-2t) and also its discrete time equivalent signal with a sampling period T = 0.2sec
+clear;
+clc ;
+close ;
+t=0:0.01:2;
+x1=2*exp(-2*t);
+subplot(1,2,1);
+plot(t,x1);
+xlabel('t');
+ylabel('x(t)');
+title('CONTINUOUS TIME PLOT');
+n=0:0.2:2;
+x2=2*exp(-2*n);
+subplot(1,2,2);
+plot2d3(n,x2);
+xlabel('n');
+ylabel('x(n)');
+title('DISCRETE TIME PLOT');
\ No newline at end of file diff --git a/52/CH1/EX1.11/Example1_11.sce b/52/CH1/EX1.11/Example1_11.sce new file mode 100755 index 000000000..8ddb5dad7 --- /dev/null +++ b/52/CH1/EX1.11/Example1_11.sce @@ -0,0 +1,12 @@ +//Example 1.11
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Testing Stability of Given System
+clear;
+clc ;
+close ;
+syms n;
+x =(1/2)^n
+X= symsum (x,n ,0, %inf );
+//Display the result in command window
+disp (X,"Summation is :");
+disp('Hence Summation < infinity. Given System is Stable');
\ No newline at end of file diff --git a/52/CH1/EX1.12/Example1_12.sce b/52/CH1/EX1.12/Example1_12.sce new file mode 100755 index 000000000..8bc62a789 --- /dev/null +++ b/52/CH1/EX1.12/Example1_12.sce @@ -0,0 +1,10 @@ +//Example 1.12
+//Program to Compute convolution of given sequences
+//x(n)=[3 2 1 2], h(n)=[1 2 1 2];
+clear;
+clc ;
+close ;
+x=[3 2 1 2];
+h=[1 2 1 2];
+y=convol(x,h);
+disp(y);
diff --git a/52/CH1/EX1.13/Example1_13.sce b/52/CH1/EX1.13/Example1_13.sce new file mode 100755 index 000000000..cc281d736 --- /dev/null +++ b/52/CH1/EX1.13/Example1_13.sce @@ -0,0 +1,10 @@ +//Example 1.13
+//Program to Compute convolution of given sequences
+//x(n)=[1 2 1 1], h(n)=[1 -1 1 -1];
+clear;
+clc ;
+close ;
+x=[1 2 1 1];
+h=[1 -1 1 -1];
+y=convol(x,h);
+disp(round(y));
diff --git a/52/CH1/EX1.18/Example1_18.sce b/52/CH1/EX1.18/Example1_18.sce new file mode 100755 index 000000000..c2c2b07d7 --- /dev/null +++ b/52/CH1/EX1.18/Example1_18.sce @@ -0,0 +1,11 @@ +//Example 1.18
+//Program to Compute Cross-correlation of given sequences
+//x(n)=[1 2 1 1], h(n)=[1 1 2 1];
+clear;
+clc ;
+close ;
+x=[1 2 1 1];
+h=[1 1 2 1];
+h1=[1 2 1 1];
+y=convol(x,h1);
+disp(round(y));
diff --git a/52/CH1/EX1.19/Example1_19.sce b/52/CH1/EX1.19/Example1_19.sce new file mode 100755 index 000000000..eeafe48b0 --- /dev/null +++ b/52/CH1/EX1.19/Example1_19.sce @@ -0,0 +1,11 @@ +//Example 1.19
+//To find input x(n)
+//h(n)=[1 2 1], y(n)=[1 5 10 11 8 4 1]
+clear;
+clc ;
+close ;
+z=%z;
+a=z^6+5*(z^(5))+10*(z^(4))+11*(z^(3))+8*(z^(2))+4*(z^(1))+1;
+b=z^6+2*z^(5)+1*z^(4);
+x =ldiv(a,b,5);
+disp (x,"x(n)=");
\ No newline at end of file diff --git a/52/CH1/EX1.2/Example1_2.jpg b/52/CH1/EX1.2/Example1_2.jpg Binary files differnew file mode 100755 index 000000000..efe3593d5 --- /dev/null +++ b/52/CH1/EX1.2/Example1_2.jpg diff --git a/52/CH1/EX1.2/Example1_2.sce b/52/CH1/EX1.2/Example1_2.sce new file mode 100755 index 000000000..73f5a4492 --- /dev/null +++ b/52/CH1/EX1.2/Example1_2.sce @@ -0,0 +1,19 @@ +//Example 1.2
+//Sketch the continuous time signal x=sin(7*t)+sin(10*t) and also its discrete time equivalent signal with a sampling period T = 0.2sec
+clear;
+clc ;
+close ;
+t=0:0.01:2;
+x1=sin(7*t)+sin(10*t);
+subplot(1,2,1);
+plot(t,x1);
+xlabel('t');
+ylabel('x(t)');
+title('CONTINUOUS TIME PLOT');
+n=0:0.2:2;
+x2=sin(7*n)+sin(10*n);
+subplot(1,2,2);
+plot2d3(n,x2);
+xlabel('n');
+ylabel('x(n)');
+title('DISCRETE TIME PLOT');
\ No newline at end of file diff --git a/52/CH1/EX1.3.a/Example1_3_a.sce b/52/CH1/EX1.3.a/Example1_3_a.sce new file mode 100755 index 000000000..a032bd1f7 --- /dev/null +++ b/52/CH1/EX1.3.a/Example1_3_a.sce @@ -0,0 +1,10 @@ +//Example 1.3 (a)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Calculate Following Summations
+clear;
+clc ;
+close ;
+syms n;
+X= symsum (sin(2*n),n ,2, 2);
+//Display the result in command window
+disp (X,"The Value of summation comes out to be:");
\ No newline at end of file diff --git a/52/CH1/EX1.3.b/Example1_3_b.sce b/52/CH1/EX1.3.b/Example1_3_b.sce new file mode 100755 index 000000000..b36910b8c --- /dev/null +++ b/52/CH1/EX1.3.b/Example1_3_b.sce @@ -0,0 +1,10 @@ +//Example 1.3 (b)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Calculate Following Summations
+clear;
+clc ;
+close ;
+syms n;
+X= symsum (%e^(2*n),n ,0, 0);
+//Display the result in command window
+disp (X,"The Value of summation comes out to be:");
\ No newline at end of file diff --git a/52/CH1/EX1.32.a/Example1_32_a.jpg b/52/CH1/EX1.32.a/Example1_32_a.jpg Binary files differnew file mode 100755 index 000000000..77bbfee64 --- /dev/null +++ b/52/CH1/EX1.32.a/Example1_32_a.jpg diff --git a/52/CH1/EX1.32.a/Example1_32_a.sce b/52/CH1/EX1.32.a/Example1_32_a.sce new file mode 100755 index 000000000..bb85f818b --- /dev/null +++ b/52/CH1/EX1.32.a/Example1_32_a.sce @@ -0,0 +1,25 @@ +//Example 1.32
+//Program to Plot Magnitude and Phase Responce
+clear;
+clc ;
+close ;
+w=-%pi:0.01:%pi;
+H=1+2*cos(w)+2*cos(2*w);
+//caluculation of Phase and Magnitude of H
+[phase_H,m]=phasemag(H);
+Hm=abs(H);
+a=gca();
+subplot(2,1,1);
+a.y_location="origin";
+plot2d(w/%pi,Hm);
+xlabel('Frequency in Radians')
+ylabel('abs(Hm)');
+title('MAGNITUDE RESPONSE');
+subplot(2,1,2);
+a=gca();
+a.x_location="origin";
+a.y_location="origin";
+plot2d(w/(2*%pi),phase_H);
+xlabel('Frequency in Radians');
+ylabel('<(H)');
+title('PHASE RESPONSE');
\ No newline at end of file diff --git a/52/CH1/EX1.37/Example1_37.jpg b/52/CH1/EX1.37/Example1_37.jpg Binary files differnew file mode 100755 index 000000000..73ced07fe --- /dev/null +++ b/52/CH1/EX1.37/Example1_37.jpg diff --git a/52/CH1/EX1.37/Example1_37.sce b/52/CH1/EX1.37/Example1_37.sce new file mode 100755 index 000000000..6f140f082 --- /dev/null +++ b/52/CH1/EX1.37/Example1_37.sce @@ -0,0 +1,26 @@ +//Example 1.37
+//Program to Plot Magnitude and Phase Responce
+//y(n)=1/2[x(n)+x(n-2)]
+clear;
+clc ;
+close ;
+w=0:0.01:%pi;
+H=(1+cos(2*w)-%i*sin(2*w))/2;
+//caluculation of Phase and Magnitude of H
+[phase_H,m]=phasemag(H);
+Hm=abs(H);
+a=gca();
+subplot(2,1,1);
+a.y_location="origin";
+plot2d(w/%pi,Hm);
+xlabel('Frequency in Radians')
+ylabel('abs(Hm)');
+title('MAGNITUDE RESPONSE');
+subplot(2,1,2);
+a=gca();
+a.x_location="origin";
+a.y_location="origin";
+plot2d(w/(2*%pi),phase_H);
+xlabel('Frequency in Radians');
+ylabel('<(H)');
+title('PHASE RESPONSE');
\ No newline at end of file diff --git a/52/CH1/EX1.38/Example1_38.jpg b/52/CH1/EX1.38/Example1_38.jpg Binary files differnew file mode 100755 index 000000000..976ecb950 --- /dev/null +++ b/52/CH1/EX1.38/Example1_38.jpg diff --git a/52/CH1/EX1.38/Example1_38.sce b/52/CH1/EX1.38/Example1_38.sce new file mode 100755 index 000000000..8d0106f2e --- /dev/null +++ b/52/CH1/EX1.38/Example1_38.sce @@ -0,0 +1,26 @@ +//Example 1.38
+//Program to Plot Magnitude and Phase Responce
+//0.5delta(n)+delta(n-1)+0.5delta(n-2)
+clear;
+clc ;
+close ;
+w=-%pi:0.01:%pi;
+H=0.5+exp(-%i*w)+0.5*exp(-%i*w);
+//caluculation of Phase and Magnitude of H
+[phase_H,m]=phasemag(H);
+Hm=abs(H);
+a=gca();
+subplot(2,1,1);
+a.y_location="origin";
+plot2d(w/%pi,Hm);
+xlabel('Frequency in Radians')
+ylabel('abs(Hm)');
+title('MAGNITUDE RESPONSE');
+subplot(2,1,2);
+a=gca();
+a.x_location="origin";
+a.y_location="origin";
+plot2d(w/(2*%pi),phase_H);
+xlabel('Frequency in Radians');
+ylabel('<(H)');
+title('PHASE RESPONSE');
\ No newline at end of file diff --git a/52/CH1/EX1.4.a/Example1_4_a.sce b/52/CH1/EX1.4.a/Example1_4_a.sce new file mode 100755 index 000000000..05002e5d6 --- /dev/null +++ b/52/CH1/EX1.4.a/Example1_4_a.sce @@ -0,0 +1,15 @@ +//Example 1.4 (a)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Find Energy and Power of Given Signals
+clear;
+clc ;
+close ;
+syms n N;
+x=(1/3)^n;
+E= symsum (x^2,n ,0, %inf);
+//Display the result in command window
+disp (E,"Energy:");
+p=(1/(2*N+1))*symsum (x^2,n ,0, N);
+P=limit(p,N,%inf);
+disp (P,"Power:");
+//The Energy is Finite and Power is 0. Therefore the given signal is an Energy Signal
\ No newline at end of file diff --git a/52/CH1/EX1.4.d/Example1_4_d.sce b/52/CH1/EX1.4.d/Example1_4_d.sce new file mode 100755 index 000000000..70de2441d --- /dev/null +++ b/52/CH1/EX1.4.d/Example1_4_d.sce @@ -0,0 +1,15 @@ +//Example 1.4 (d)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Find Energy and Power of Given Signals
+clear;
+clc ;
+close ;
+syms n N;
+x=%e^(2*n);
+E= symsum (x^2,n ,0, %inf);
+//Display the result in command window
+disp (E,"Energy:");
+p=(1/(2*N+1))*symsum (x^2,n ,0, N);
+P=limit(p,N,%inf);
+disp (P,"Power:");
+//The Energy andPower is infinite. Therefore the given signal is an neither Energy Signal nor Power Signal
\ No newline at end of file diff --git a/52/CH1/EX1.45/Example1_45.jpg b/52/CH1/EX1.45/Example1_45.jpg Binary files differnew file mode 100755 index 000000000..dee4c49e3 --- /dev/null +++ b/52/CH1/EX1.45/Example1_45.jpg diff --git a/52/CH1/EX1.45/Example1_45.sce b/52/CH1/EX1.45/Example1_45.sce new file mode 100755 index 000000000..eb624a1b8 --- /dev/null +++ b/52/CH1/EX1.45/Example1_45.sce @@ -0,0 +1,22 @@ +//Example 1.45
+//
+clear;
+clc ;
+close ;
+t=0:0.01:10;
+x=2*cos(5*t)+cos(300*t);
+x1=2*cos(5*t);
+b=[0.05 0.05];
+a=[1 -0.9];
+y=filter(b,a,x);
+subplot(2,1,1);
+plot(t,x);
+xlabel('Time in Sec');
+ylabel('Amplitude');
+subplot(2,1,2);
+plot(t,y);
+subplot(2,1,2);
+plot(t,x1,':');
+title('x:SIGNAL WITHOUT NOISE y:SIGNAL WITH NOISE');
+xlabel('Time in Sec');
+ylabel('Amplitude');
\ No newline at end of file diff --git a/52/CH1/EX1.5.a/Example1_5_a.jpg b/52/CH1/EX1.5.a/Example1_5_a.jpg Binary files differnew file mode 100755 index 000000000..34bf2bf73 --- /dev/null +++ b/52/CH1/EX1.5.a/Example1_5_a.jpg diff --git a/52/CH1/EX1.5.a/Example1_5_a.sce b/52/CH1/EX1.5.a/Example1_5_a.sce new file mode 100755 index 000000000..050bff03d --- /dev/null +++ b/52/CH1/EX1.5.a/Example1_5_a.sce @@ -0,0 +1,20 @@ +//Example 1.5 (a)
+//To Determine Whether Given Signal is Periodic or not
+clear;
+clc ;
+close ;
+t=0:0.01:2;
+x1=exp(%i*6*%pi*t);
+subplot(1,2,1);
+plot(t,x1);
+xlabel('t');
+ylabel('x(t)');
+title('CONTINUOUS TIME PLOT');
+n=0:0.2:2;
+x2=exp(%i*6*%pi*n);
+subplot(1,2,2);
+plot2d3(n,x2);
+xlabel('n');
+ylabel('x(n)');
+title('DISCRETE TIME PLOT');
+//Hence Given Signal is Periodic with N=1
\ No newline at end of file diff --git a/52/CH1/EX1.5.c/Example1_5_c.jpg b/52/CH1/EX1.5.c/Example1_5_c.jpg Binary files differnew file mode 100755 index 000000000..476b3bbd0 --- /dev/null +++ b/52/CH1/EX1.5.c/Example1_5_c.jpg diff --git a/52/CH1/EX1.5.c/Example1_5_c.sce b/52/CH1/EX1.5.c/Example1_5_c.sce new file mode 100755 index 000000000..08cd6aa2f --- /dev/null +++ b/52/CH1/EX1.5.c/Example1_5_c.sce @@ -0,0 +1,20 @@ +//Example 1.5 (c)
+//To Determine Whether Given Signal is Periodic or not
+clear;
+clc ;
+close ;
+t=0:0.01:10;
+x1=cos(2*%pi*t/3);
+subplot(1,2,1);
+plot(t,x1);
+xlabel('t');
+ylabel('x(t)');
+title('CONTINUOUS TIME PLOT');
+n=0:0.2:10;
+x2=cos(2*%pi*n/3);
+subplot(1,2,2);
+plot2d3(n,x2);
+xlabel('n');
+ylabel('x(n)');
+title('DISCRETE TIME PLOT');
+//Hence Given Signal is Periodic with N=3
\ No newline at end of file diff --git a/52/CH1/EX1.5.d/Example1_5_d.jpg b/52/CH1/EX1.5.d/Example1_5_d.jpg Binary files differnew file mode 100755 index 000000000..f198e232c --- /dev/null +++ b/52/CH1/EX1.5.d/Example1_5_d.jpg diff --git a/52/CH1/EX1.5.d/Example1_5_d.sce b/52/CH1/EX1.5.d/Example1_5_d.sce new file mode 100755 index 000000000..10dd7274a --- /dev/null +++ b/52/CH1/EX1.5.d/Example1_5_d.sce @@ -0,0 +1,20 @@ +//Example 1.5 (d)
+//To Determine Whether Given Signal is Periodic or not
+clear;
+clc ;
+close ;
+t=0:0.01:50;
+x1=cos(%pi*t/3)+cos(3*%pi*t/4);
+subplot(1,2,1);
+plot(t,x1);
+xlabel('t');
+ylabel('x(t)');
+title('CONTINUOUS TIME PLOT');
+n=0:0.2:50;
+x2=cos(%pi*n/3)+cos(3*%pi*n/4);
+subplot(1,2,2);
+plot2d3(n,x2);
+xlabel('n');
+ylabel('x(n)');
+title('DISCRETE TIME PLOT');
+//Hence Given Signal is Periodic with N=24
\ No newline at end of file diff --git a/52/CH1/EX1.57.a/Example1_57_a.sce b/52/CH1/EX1.57.a/Example1_57_a.sce new file mode 100755 index 000000000..8c48c13e5 --- /dev/null +++ b/52/CH1/EX1.57.a/Example1_57_a.sce @@ -0,0 +1,10 @@ +//Example 1.57 (a)
+//Program to Compute discrete convolution of given sequences
+//x(n)=[1 2 -1 1], h(n)=[1 0 1 1];
+clear;
+clc ;
+close ;
+x=[1 2 -1 1];
+h=[1 0 1 1];
+y=convol(x,h);
+disp(round(y));
diff --git a/52/CH1/EX1.61/Example1_61.sce b/52/CH1/EX1.61/Example1_61.sce new file mode 100755 index 000000000..124b6fb97 --- /dev/null +++ b/52/CH1/EX1.61/Example1_61.sce @@ -0,0 +1,11 @@ +//Example 1.61
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Fourier transform of (3)^n u(n)
+clear;
+clc ;
+close ;
+syms n;
+x =(3) ^n;
+X= symsum (x,n ,0, %inf )
+//Display the result in command window
+disp (X,"The Fourier Transform does not exit as x(n) is not absolutely summable and approaches infinity i.e.");
\ No newline at end of file diff --git a/52/CH1/EX1.62/Example1_62.sce b/52/CH1/EX1.62/Example1_62.sce new file mode 100755 index 000000000..fefc888d9 --- /dev/null +++ b/52/CH1/EX1.62/Example1_62.sce @@ -0,0 +1,10 @@ +//Example 1.62
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Fourier transform of (0.8)^|n| u(n)
+clear;
+clc ;
+close ;
+syms w n;
+X= symsum ((0.8)^n*%e^(%i*w*n),n ,1, %inf )+symsum ((0.8)^n*%e^(-%i*w*n),n ,0, %inf )
+//Display the result in command window
+disp (X,"The Fourier Transform comes out to be:");
\ No newline at end of file diff --git a/52/CH1/EX1.64.a/Example1_64_a.jpg b/52/CH1/EX1.64.a/Example1_64_a.jpg Binary files differnew file mode 100755 index 000000000..4677df5c5 --- /dev/null +++ b/52/CH1/EX1.64.a/Example1_64_a.jpg diff --git a/52/CH1/EX1.64.a/Example1_64_a.sce b/52/CH1/EX1.64.a/Example1_64_a.sce new file mode 100755 index 000000000..47bdcdf6e --- /dev/null +++ b/52/CH1/EX1.64.a/Example1_64_a.sce @@ -0,0 +1,25 @@ +//Example 1.64 (a)
+//Program to Calculate Plot Magnitude and Phase Responce
+clear;
+clc ;
+close ;
+w=0:0.01:%pi;
+H=1/(1-0.5*%e^(-%i*w));
+//caluculation of Phase and Magnitude of H
+[phase_H,m]=phasemag(H);
+Hm=abs(H);
+a=gca();
+subplot(2,1,1);
+a.y_location="origin";
+plot2d(w/%pi,Hm);
+xlabel('Frequency in Radians')
+ylabel('abs(Hm)');
+title('MAGNITUDE RESPONSE');
+subplot(2,1,2);
+a=gca();
+a.x_location="origin";
+a.y_location="origin";
+plot2d(w/(2*%pi),phase_H);
+xlabel('Frequency in Radians');
+ylabel('<(H)');
+title('PHASE RESPONSE');
\ No newline at end of file diff --git a/52/CH1/EX1.64.c/Example1_64_c.jpg b/52/CH1/EX1.64.c/Example1_64_c.jpg Binary files differnew file mode 100755 index 000000000..19bedf0c6 --- /dev/null +++ b/52/CH1/EX1.64.c/Example1_64_c.jpg diff --git a/52/CH1/EX1.64.c/Example1_64_c.sce b/52/CH1/EX1.64.c/Example1_64_c.sce new file mode 100755 index 000000000..5b8f5c9d9 --- /dev/null +++ b/52/CH1/EX1.64.c/Example1_64_c.sce @@ -0,0 +1,25 @@ +//Example 1.64 (c)
+//Program to Calculate Plot Magnitude and Phase Responce
+clear;
+clc ;
+close ;
+w=0:0.01:%pi;
+H=1/(1-0.9*%i*%e^(-%i*w));
+//caluculation of Phase and Magnitude of H
+[phase_H,m]=phasemag(H);
+Hm=abs(H);
+a=gca();
+subplot(2,1,1);
+a.y_location="origin";
+plot2d(w/%pi,Hm);
+xlabel('Frequency in Radians')
+ylabel('abs(Hm)');
+title('MAGNITUDE RESPONSE');
+subplot(2,1,2);
+a=gca();
+a.x_location="origin";
+a.y_location="origin";
+plot2d(w/(2*%pi),phase_H);
+xlabel('Frequency in Radians');
+ylabel('<(H)');
+title('PHASE RESPONSE');
\ No newline at end of file diff --git a/52/CH11/EX11.3/Program11_3.sce b/52/CH11/EX11.3/Program11_3.sce new file mode 100755 index 000000000..4f5338f47 --- /dev/null +++ b/52/CH11/EX11.3/Program11_3.sce @@ -0,0 +1,13 @@ +//Program 11.3
+//Program To Calculate the value of the function
+//Y=A*B
+clear;
+clc;
+close;
+//Input Data
+A=input('Enter Integer Number A =');
+B=input('Enter Integer Number B =');
+//Multiplication Computation
+Y=A*B;
+//Display the result in command window
+disp(Y,"Y = A*B = ");
\ No newline at end of file diff --git a/52/CH11/EX11.5/Program11_5.sce b/52/CH11/EX11.5/Program11_5.sce new file mode 100755 index 000000000..1edd39d51 --- /dev/null +++ b/52/CH11/EX11.5/Program11_5.sce @@ -0,0 +1,17 @@ +//Program 11.5
+//Program To Calculate the value of the function
+//Y=A*X1+B*X2+C*X3
+clear;
+clc;
+close;
+//Data
+A=1;
+B=2;
+C=3;
+X1=4;
+X2=5;
+X3=6;
+//Compute the function
+Y=A*X1+B*X2+C*X3;
+//Display the result in command window
+disp(Y,"Y = A*X1+B*X2+C*X3 = ");
\ No newline at end of file diff --git a/52/CH2/EX2.1/Example2_1.sce b/52/CH2/EX2.1/Example2_1.sce new file mode 100755 index 000000000..bd90f996e --- /dev/null +++ b/52/CH2/EX2.1/Example2_1.sce @@ -0,0 +1,15 @@ +//Example 2.1
+//Z- transform of [1 0 3 -1 2]
+clear;
+clc ;
+close ;
+function[za]=ztransfer(sequence,n)
+z=poly(0,'z','r')
+za=sequence*(1/z)^n'
+endfunction
+x1=[1 0 3 -1 2];
+n=0:length(x1)-1;
+zz=ztransfer(x1,n);
+//Display the result in command window
+disp (zz,"Z-transform of sequence is:");
+disp('ROC is the entire plane except z = 0');
\ No newline at end of file diff --git a/52/CH2/EX2.10/Example2_10.sce b/52/CH2/EX2.10/Example2_10.sce new file mode 100755 index 000000000..9a853b88d --- /dev/null +++ b/52/CH2/EX2.10/Example2_10.sce @@ -0,0 +1,11 @@ +//Example 2.10
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z transform of r^n.cos(Wo*n)
+clc;
+syms r Wo n z;
+x1=(r^n)*exp(sqrt(-1)*Wo*n);
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=(r^n)*exp(-sqrt(-1)*Wo*n);
+X2=symsum(x2*(z^-n),n,0,%inf);
+X=(X1+X2)/2;
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/52/CH2/EX2.11/Example2_11.sce b/52/CH2/EX2.11/Example2_11.sce new file mode 100755 index 000000000..e8bdfd8c3 --- /dev/null +++ b/52/CH2/EX2.11/Example2_11.sce @@ -0,0 +1,12 @@ +//Example 2.11
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of n.a^n u(n)
+clear;
+clc ;
+close ;
+syms a n z;
+x =(a) ^n;
+X= symsum (x*(z^(-n)),n ,0, %inf )
+Y = diff (X,z);
+//Display the result in command window
+disp (Y,"Z-transform of n.a^n u(n) is:");
\ No newline at end of file diff --git a/52/CH2/EX2.13.a/Example2_13_a.sce b/52/CH2/EX2.13.a/Example2_13_a.sce new file mode 100755 index 000000000..42728da95 --- /dev/null +++ b/52/CH2/EX2.13.a/Example2_13_a.sce @@ -0,0 +1,15 @@ +//Example 2.13 (a)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of (-1/5)^n u(n)+5(1/2)^(-n)u(-n-1)
+clear;
+clc ;
+close ;
+syms n z;
+x1 =(-1/5)^n ;
+X1= symsum (x1 *(z^(-n)),n ,0, %inf );
+x2 =(1/2)^(-n);
+X2= symsum (5* x2 *(z^(-n)),n ,0, %inf );
+X = (X1 -X2);
+//Display the result in command window
+disp (X,"Z-transform of [3(3^n)-4(2)^n] u(n) is:");
+disp('ROC is the Region 1/5 < mod(z) < 2');
\ No newline at end of file diff --git a/52/CH2/EX2.13.b/Example2_13_b.sce b/52/CH2/EX2.13.b/Example2_13_b.sce new file mode 100755 index 000000000..30a683da5 --- /dev/null +++ b/52/CH2/EX2.13.b/Example2_13_b.sce @@ -0,0 +1,13 @@ +//Example 2.13 (b)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z transform
+clc;
+syms n z k;
+x1=1;
+X1=symsum(x1*z^(-n),n,0,0);
+x2=1;
+X2=symsum(x2*z^(-n),n,1,1);
+x3=1;
+X3=symsum(x3*z^(-n),n,2,2);
+X=0.5*X1+X2-1/3*X3;
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/52/CH2/EX2.13.c/Example2_13_c.sce b/52/CH2/EX2.13.c/Example2_13_c.sce new file mode 100755 index 000000000..6256d04bf --- /dev/null +++ b/52/CH2/EX2.13.c/Example2_13_c.sce @@ -0,0 +1,11 @@ +//Example 2.13 (c)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of u(n-2)
+clear;
+clc ;
+close ;
+syms n z;
+x =1;
+X= (1/(z^2))*symsum (x*(z^(-n)),n ,0, %inf );
+//Display the result in command window
+disp (X,"Z-transform of u(n-2) is:");
\ No newline at end of file diff --git a/52/CH2/EX2.13.d/Example2_13_d.sce b/52/CH2/EX2.13.d/Example2_13_d.sce new file mode 100755 index 000000000..d9a228551 --- /dev/null +++ b/52/CH2/EX2.13.d/Example2_13_d.sce @@ -0,0 +1,15 @@ +//Example 2.13 (d)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of (n+0.5)((1/3)^n)u(n)
+clear;
+clc ;
+close ;
+syms n z;
+x1 =(1/3)^n;
+X11= symsum (x1*(z^(-n)),n ,0, %inf )
+X1 = diff (X11,z);
+x2 =(1/3) ^(n);
+X2= symsum (0.5* x2 *(z^(-n)),n ,0, %inf );
+X = (X1+X2);
+//Display the result in command window
+disp (X,"Z-transform of (n+0.5)((1/3)^n)u(n) is:");
\ No newline at end of file diff --git a/52/CH2/EX2.16/Example2_16.sce b/52/CH2/EX2.16/Example2_16.sce new file mode 100755 index 000000000..8dd5a5a10 --- /dev/null +++ b/52/CH2/EX2.16/Example2_16.sce @@ -0,0 +1,11 @@ +//Example 2.16
+//To find input h(n)
+//a=[1 2 -4 1], b=[1]
+clear;
+clc ;
+close ;
+z=%z;
+a=z^3+2*(z^(2))-4*(z)+1;
+b=z^3;
+h =ldiv(a,b,4);
+disp (h,"h(n)=");
\ No newline at end of file diff --git a/52/CH2/EX2.17/Example2_17.jpg b/52/CH2/EX2.17/Example2_17.jpg Binary files differnew file mode 100755 index 000000000..c177d3489 --- /dev/null +++ b/52/CH2/EX2.17/Example2_17.jpg diff --git a/52/CH2/EX2.17/Example2_17.sce b/52/CH2/EX2.17/Example2_17.sce new file mode 100755 index 000000000..dbc5d815d --- /dev/null +++ b/52/CH2/EX2.17/Example2_17.sce @@ -0,0 +1,9 @@ +//Example 2.17
+//To draw the pole-zero plot
+clear;
+clc ;
+close ;
+z=%z
+H1Z=((z)*(z-1))/((z-0.25)*(z-0.5));
+xset('window',1);
+plzr(H1Z);
\ No newline at end of file diff --git a/52/CH2/EX2.19/Example2_19.jpg b/52/CH2/EX2.19/Example2_19.jpg Binary files differnew file mode 100755 index 000000000..0b0ee5ab2 --- /dev/null +++ b/52/CH2/EX2.19/Example2_19.jpg diff --git a/52/CH2/EX2.19/Example2_19.sce b/52/CH2/EX2.19/Example2_19.sce new file mode 100755 index 000000000..a3bb6f748 --- /dev/null +++ b/52/CH2/EX2.19/Example2_19.sce @@ -0,0 +1,25 @@ +//Example 2.19
+//Program to Plot Magnitude and Phase Responce
+clear;
+clc ;
+close ;
+w=-%pi:0.01:%pi;
+H=1/(1-0.5*(cos(w)-%i*sin(w)));
+//caluculation of Phase and Magnitude of H
+[phase_H,m]=phasemag(H);
+Hm=abs(H);
+a=gca();
+subplot(2,1,1);
+a.y_location="origin";
+plot2d(w/%pi,Hm);
+xlabel('Frequency in Radians');
+ylabel('abs(Hm)');
+title('MAGNITUDE RESPONSE');
+subplot(2,1,2);
+a=gca();
+a.x_location="origin";
+a.y_location="origin";
+plot2d(w/(2*%pi),phase_H);
+xlabel('Frequency in Radians');
+ylabel('<(H)');
+title('PHASE RESPONSE');
\ No newline at end of file diff --git a/52/CH2/EX2.2/Example2_2.sce b/52/CH2/EX2.2/Example2_2.sce new file mode 100755 index 000000000..1cf7f1f5a --- /dev/null +++ b/52/CH2/EX2.2/Example2_2.sce @@ -0,0 +1,15 @@ +//Example 2.2
+//Z- transform of [-3 -2 -1 0 1]
+clear;
+clc ;
+close ;
+function[za]=ztransfer(sequence,n)
+z=poly(0,'z','r')
+za=sequence*(1/z)^n'
+endfunction
+x1=[-3 -2 -1 0 1];
+n=-(length(x1)-1):0;
+zz=ztransfer(x1,n);
+//Display the result in command window
+disp (zz,"Z-transform of sequence is:");
+disp('ROC is the entire plane except z = %inf');
\ No newline at end of file diff --git a/52/CH2/EX2.20.a/Example2_20_a.sce b/52/CH2/EX2.20.a/Example2_20_a.sce new file mode 100755 index 000000000..23a069f2a --- /dev/null +++ b/52/CH2/EX2.20.a/Example2_20_a.sce @@ -0,0 +1,11 @@ +//Example 2.10 (a)
+//To find input h(n)
+//X(z)=(z+0.2)/((z+0.5)(z-1);
+clear;
+clc ;
+close ;
+z=%z;
+a=(z+0.5)*(z-1);
+b=z+0.2;
+h =ldiv(b,a,4);
+disp (h,"h(n)=");
\ No newline at end of file diff --git a/52/CH2/EX2.22/Example2_22.sce b/52/CH2/EX2.22/Example2_22.sce new file mode 100755 index 000000000..cd24b210e --- /dev/null +++ b/52/CH2/EX2.22/Example2_22.sce @@ -0,0 +1,11 @@ +//Example 2.22
+//To find input x(n)
+//X(z)=1/(2*z^(-2)+2*z^(-1)+1);
+clear;
+clc ;
+close ;
+z=%z;
+a=(2+2*z+z^2);
+b=z^2;
+h =ldiv(b,a,6);
+disp (h,"First six values of h(n)=");
\ No newline at end of file diff --git a/52/CH2/EX2.23/Example2_23.sce b/52/CH2/EX2.23/Example2_23.sce new file mode 100755 index 000000000..776b7f7e9 --- /dev/null +++ b/52/CH2/EX2.23/Example2_23.sce @@ -0,0 +1,11 @@ +//Example 2.23
+//To find input x(n)
+//X(z)=1/(1-2z^(-1))(1-z^(-1))^2;
+clear;
+clc ;
+close ;
+z=%z;
+a=(z-2)*(z-1)^2;
+b=z^3;
+h =ldiv(b,a,6);
+disp (h,"First six values of h(n)=");
\ No newline at end of file diff --git a/52/CH2/EX2.3/Example2_3.sce b/52/CH2/EX2.3/Example2_3.sce new file mode 100755 index 000000000..559f6dd8d --- /dev/null +++ b/52/CH2/EX2.3/Example2_3.sce @@ -0,0 +1,15 @@ +//Example 2.3
+//Z- transform of [2 -1 3 2 1 0 2 3 -1]
+clear;
+clc ;
+close ;
+function[za]=ztransfer(sequence,n)
+z=poly(0,'z','r')
+za=sequence*(1/z)^n'
+endfunction
+x1=[2 -1 3 2 1 0 2 3 -1];
+n=-4:4;
+zz=ztransfer(x1,n);
+//Display the result in command window
+disp (zz,"Z-transform of sequence is:");
+disp('ROC is the entire plane except z = 0 and z = %inf');
\ No newline at end of file diff --git a/52/CH2/EX2.34/Example2_34.jpg b/52/CH2/EX2.34/Example2_34.jpg Binary files differnew file mode 100755 index 000000000..31338edf8 --- /dev/null +++ b/52/CH2/EX2.34/Example2_34.jpg diff --git a/52/CH2/EX2.34/Example2_34.sce b/52/CH2/EX2.34/Example2_34.sce new file mode 100755 index 000000000..818855caa --- /dev/null +++ b/52/CH2/EX2.34/Example2_34.sce @@ -0,0 +1,27 @@ +//Example 2.34
+//To plot the impulse responce of the system analyically and using scilab
+clear;
+clc ;
+close ;
+n=0:1:50;
+x=[1,zeros(1,50)];
+b=[1 2];
+a=[1 -3 -4];
+yanaly=6/5*4.^n-1/5*(-1).^n;//Analytical Solution
+ymat=filter(b,a,x);
+subplot(3,1,1);
+plot2d3(n,x);
+xlabel('n');
+ylabel('x(n)');
+title('INPUT SEQUENCE (IMPULSE FUNCTION)');
+subplot(3,1,2);
+plot2d3(n,yanaly);
+xlabel('n');
+ylabel('y(n)');
+title('OUTPUT SEQUENCE yanaly');
+subplot(3,1,3);
+plot2d3(n,ymat);
+xlabel('n');
+ylabel('y(n)');
+title('OUTPUT SEQUENCE ymat');
+//As the Analtical Plot matches the Scilab Plot hence it is the Responce of the system
\ No newline at end of file diff --git a/52/CH2/EX2.35.a/Example2_35_a.jpg b/52/CH2/EX2.35.a/Example2_35_a.jpg Binary files differnew file mode 100755 index 000000000..981d4af0c --- /dev/null +++ b/52/CH2/EX2.35.a/Example2_35_a.jpg diff --git a/52/CH2/EX2.35.a/Example2_35_a.sce b/52/CH2/EX2.35.a/Example2_35_a.sce new file mode 100755 index 000000000..3c8484922 --- /dev/null +++ b/52/CH2/EX2.35.a/Example2_35_a.sce @@ -0,0 +1,9 @@ +//Example 2.35 (a)
+//To draw the pole-zero plot
+clear;
+clc ;
+close ;
+z=%z
+H1Z=(z)/(z^2-z-1);
+xset('window',1);
+plzr(H1Z);
\ No newline at end of file diff --git a/52/CH2/EX2.35.b/Example2_35_b.jpg b/52/CH2/EX2.35.b/Example2_35_b.jpg Binary files differnew file mode 100755 index 000000000..df5226c81 --- /dev/null +++ b/52/CH2/EX2.35.b/Example2_35_b.jpg diff --git a/52/CH2/EX2.35.b/Example2_35_b.sce b/52/CH2/EX2.35.b/Example2_35_b.sce new file mode 100755 index 000000000..73cc4f1b7 --- /dev/null +++ b/52/CH2/EX2.35.b/Example2_35_b.sce @@ -0,0 +1,27 @@ +//Example 2.35 (b)
+//To plot the responce of the system analyically and using scilab
+clear;
+clc ;
+close ;
+n=0:1:20;
+x=ones(1,length(n));
+b=[0 1];
+a=[1 -1 -1];
+yanaly=0.447*(1.618).^n-0.447*(-0.618).^n;//Analytical Solution
+[ymat,zf]=filter(b, a, x);
+subplot(3,1,1);
+plot2d3(n,x);
+xlabel('n');
+ylabel('x(n)');
+title('INPUT SEQUENCE (STEP FUNCTION)');
+subplot(3,1,2);
+plot2d3(n,yanaly);
+xlabel('n');
+ylabel('y(n)');
+title('OUTPUT SEQUENCE yanaly');
+subplot(3,1,3);
+plot2d3(n,ymat,zf);
+xlabel('n');
+ylabel('y(n)');
+title('OUTPUT SEQUENCE ymat');
+//As the Analtical Plot matches the Scilab Plot hence it is the Responce of the system
\ No newline at end of file diff --git a/52/CH2/EX2.38/Example2_38.jpg b/52/CH2/EX2.38/Example2_38.jpg Binary files differnew file mode 100755 index 000000000..f3bddeada --- /dev/null +++ b/52/CH2/EX2.38/Example2_38.jpg diff --git a/52/CH2/EX2.38/Example2_38.sce b/52/CH2/EX2.38/Example2_38.sce new file mode 100755 index 000000000..075bee4f8 --- /dev/null +++ b/52/CH2/EX2.38/Example2_38.sce @@ -0,0 +1,27 @@ +//Example 2.38
+//To plot the responce of the system analyically and using scilab
+clear;
+clc ;
+close ;
+n=0:1:20;
+x=n;
+b=[0 1 1];
+a=[1 -0.7 0.12];
+yanaly=38.89*(0.4).^n-26.53*(0.3).^n-12.36+4.76*n;//Analytical Solution
+ymat=filter(b,a,x);
+subplot(3,1,1);
+plot2d3(n,x);
+xlabel('n');
+ylabel('x(n)');
+title('INPUT SEQUENCE (RAMP FUNCTION)');
+subplot(3,1,2);
+plot2d3(n,yanaly);
+xlabel('n');
+ylabel('y(n)');
+title('OUTPUT SEQUENCE yanaly');
+subplot(3,1,3);
+plot2d3(n,ymat);
+xlabel('n');
+ylabel('y(n)');
+title('OUTPUT SEQUENCE ymat');
+//As the Analtical Plot matches the Scilab Plot hence it is the Responce of the system
\ No newline at end of file diff --git a/52/CH2/EX2.4/Example2_4.sce b/52/CH2/EX2.4/Example2_4.sce new file mode 100755 index 000000000..bc0802459 --- /dev/null +++ b/52/CH2/EX2.4/Example2_4.sce @@ -0,0 +1,12 @@ +//Example 2.4
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of a^n u(n)
+clear;
+clc ;
+close ;
+syms a n z;
+x =a^n
+X= symsum (x*(z^(-n)),n ,0, %inf );
+//Display the result in command window
+disp (X,"Z-transform of a^n u(n) with is:");
+disp('ROC is the Region mod(z) > a')
\ No newline at end of file diff --git a/52/CH2/EX2.40/Example2_40.sce b/52/CH2/EX2.40/Example2_40.sce new file mode 100755 index 000000000..d92ecb1c0 --- /dev/null +++ b/52/CH2/EX2.40/Example2_40.sce @@ -0,0 +1,11 @@ +//Example 2.40
+//To find input x(n)
+//h(n)=1 2 3 2, y(n)=[1 3 7 10 10 7 2]
+clear;
+clc ;
+close ;
+z=%z;
+a=z^6+3*(z^(5))+7*(z^(4))+10*(z^(3))+10*(z^(2))+7*(z^(1))+2;
+b=z^6+2*z^(5)+3*z^(4)+2*z^(3);
+x =ldiv(a,b,4);
+disp (x,"x(n)=");
\ No newline at end of file diff --git a/52/CH2/EX2.41.a/Example2_41_a.sce b/52/CH2/EX2.41.a/Example2_41_a.sce new file mode 100755 index 000000000..68a491f3d --- /dev/null +++ b/52/CH2/EX2.41.a/Example2_41_a.sce @@ -0,0 +1,12 @@ +//Example 2.41 (a)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of n.(-1)^n u(n)
+clear;
+clc ;
+close ;
+syms a n z;
+x =(-1) ^n;
+X= symsum (x*(z^(-n)),n ,0, %inf )
+Y = diff (X,z);
+//Display the result in command window
+disp (Y,"Z-transform of n.(-1)^n u(n) is:");
\ No newline at end of file diff --git a/52/CH2/EX2.41.b/Example2_41_b.sce b/52/CH2/EX2.41.b/Example2_41_b.sce new file mode 100755 index 000000000..b5640e261 --- /dev/null +++ b/52/CH2/EX2.41.b/Example2_41_b.sce @@ -0,0 +1,12 @@ +//Example 2.41 (b)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of n^2 u(n)
+clear;
+clc ;
+close ;
+syms n z;
+x =1;
+X= symsum (x*(z^(-n)),n ,0, %inf )
+Y = diff(diff (X,z),z);
+//Display the result in command window
+disp (Y,"Z-transform of n^2 u(n) is:");
\ No newline at end of file diff --git a/52/CH2/EX2.41.c/Example2_41_c.sce b/52/CH2/EX2.41.c/Example2_41_c.sce new file mode 100755 index 000000000..3f675ec2e --- /dev/null +++ b/52/CH2/EX2.41.c/Example2_41_c.sce @@ -0,0 +1,12 @@ +//Example 2.41 (c)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z transform of (-1)^n.cos(%pi/3*n)
+clc;
+syms n z;
+Wo=%pi/3;
+x1=exp(sqrt(-1)*Wo*n);
+X1=(-1)^n*symsum(x1*(z^-n),n,0,%inf);
+x2=exp(-sqrt(-1)*Wo*n);
+X2=(-1)^n*symsum(x2*(z^-n),n,0,%inf);
+X=(X1+X2)/2;
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/52/CH2/EX2.45/Example2_45.jpg b/52/CH2/EX2.45/Example2_45.jpg Binary files differnew file mode 100755 index 000000000..1a24182b0 --- /dev/null +++ b/52/CH2/EX2.45/Example2_45.jpg diff --git a/52/CH2/EX2.45/Example2_45.sce b/52/CH2/EX2.45/Example2_45.sce new file mode 100755 index 000000000..ae750c8e1 --- /dev/null +++ b/52/CH2/EX2.45/Example2_45.sce @@ -0,0 +1,9 @@ +//Example 2.45
+//To draw the pole-zero plot
+clear;
+clc ;
+close ;
+z=%z
+H1Z=((z)*(z+1))/(z^2-z+0.5);
+xset('window',1);
+plzr(H1Z);
\ No newline at end of file diff --git a/52/CH2/EX2.5/Example2_5.sce b/52/CH2/EX2.5/Example2_5.sce new file mode 100755 index 000000000..e33226607 --- /dev/null +++ b/52/CH2/EX2.5/Example2_5.sce @@ -0,0 +1,12 @@ +//Example 2.5
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of -b^n u(-n-1)
+clear;
+clc ;
+close ;
+syms b n z;
+x =b^n
+X= symsum (x*(z^(-n)),n ,0, %inf );
+//Display the result in command window
+disp (X,"Z-transform of b^n u(n) with is:");
+disp('ROC is the Region mod(z) < b')
\ No newline at end of file diff --git a/52/CH2/EX2.53.a/Example2_53_a.sce b/52/CH2/EX2.53.a/Example2_53_a.sce new file mode 100755 index 000000000..7bc4991f5 --- /dev/null +++ b/52/CH2/EX2.53.a/Example2_53_a.sce @@ -0,0 +1,14 @@ +//Example 2.53 (a)
+//Z- transform of [3 1 2 5 7 0 1]
+clear;
+clc ;
+close ;
+function[za]=ztransfer(sequence,n)
+z=poly(0,'z','r')
+za=sequence*(1/z)^n'
+endfunction
+x1=[3 1 2 5 7 0 1];
+n=-3:3;
+zz=ztransfer(x1,n);
+//Display the result in command window
+disp (zz,"Z-transform of sequence is:");
\ No newline at end of file diff --git a/52/CH2/EX2.53.b/Example2_53_b.sce b/52/CH2/EX2.53.b/Example2_53_b.sce new file mode 100755 index 000000000..0fede7b16 --- /dev/null +++ b/52/CH2/EX2.53.b/Example2_53_b.sce @@ -0,0 +1,8 @@ +//Example 2.53 (b)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z transform of delta(n)
+clc;
+syms n z;
+x=1;
+X=symsum(x*z^(-n),n,0,0);
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/52/CH2/EX2.53.c/Example2_53_c.sce b/52/CH2/EX2.53.c/Example2_53_c.sce new file mode 100755 index 000000000..cf0d921ae --- /dev/null +++ b/52/CH2/EX2.53.c/Example2_53_c.sce @@ -0,0 +1,8 @@ +//Example 2.53 (c)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z transform of delta(n)
+clc;
+syms n z k;
+x=1;
+X=symsum(x*z^(-n),n,k,k);
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/52/CH2/EX2.53.d/Example2_53_d.sce b/52/CH2/EX2.53.d/Example2_53_d.sce new file mode 100755 index 000000000..92e090818 --- /dev/null +++ b/52/CH2/EX2.53.d/Example2_53_d.sce @@ -0,0 +1,8 @@ +//Example 2.53 (d)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z transform of delta(n)
+clc;
+syms n z kc;
+x=1;
+X=symsum(x*z^(-n),n,-k,-k);
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/52/CH2/EX2.54/Example2_54.sce b/52/CH2/EX2.54/Example2_54.sce new file mode 100755 index 000000000..daeb4b2da --- /dev/null +++ b/52/CH2/EX2.54/Example2_54.sce @@ -0,0 +1,11 @@ +//Example 2.54
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z transform of cos(Wo*n)
+clc;
+syms Wo n z;
+x1=exp(sqrt(-1)*Wo*n);
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=exp(-sqrt(-1)*Wo*n);
+X2=symsum(x2*(z^-n),n,0,%inf);
+X=(X1+X2)/2;
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/52/CH2/EX2.58/Example2_58.jpg b/52/CH2/EX2.58/Example2_58.jpg Binary files differnew file mode 100755 index 000000000..553104730 --- /dev/null +++ b/52/CH2/EX2.58/Example2_58.jpg diff --git a/52/CH2/EX2.58/Example2_58.sce b/52/CH2/EX2.58/Example2_58.sce new file mode 100755 index 000000000..741170cc8 --- /dev/null +++ b/52/CH2/EX2.58/Example2_58.sce @@ -0,0 +1,27 @@ +//Example 2.58
+//To plot the responce of the system analyically and using scilab
+clear;
+clc ;
+close ;
+n=0:1:20;
+x=[1 zeros(1,20)];
+b=[1 -0.5];
+a=[1 -1 3/16];
+yanaly=0.5*(0.75).^n+0.5*(0.25).^n;//Analytical Solution
+ymat=filter(b,a,x);
+subplot(3,1,1);
+plot2d3(n,x);
+xlabel('n');
+ylabel('x(n)');
+title('INPUT SEQUENCE (IMPULSE FUNCTION)');
+subplot(3,1,2);
+plot2d3(n,yanaly);
+xlabel('n');
+ylabel('y(n)');
+title('OUTPUT SEQUENCE yanaly');
+subplot(3,1,3);
+plot2d3(n,ymat);
+xlabel('n');
+ylabel('y(n)');
+title('OUTPUT SEQUENCE ymat');
+//As the Analtical Plot matches the Scilab Plot hence it is the Responce of the system
\ No newline at end of file diff --git a/52/CH2/EX2.6/Example2_6.sce b/52/CH2/EX2.6/Example2_6.sce new file mode 100755 index 000000000..4f1275b00 --- /dev/null +++ b/52/CH2/EX2.6/Example2_6.sce @@ -0,0 +1,12 @@ +//Example 2.6
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of 2^n u(n)
+clear;
+clc ;
+close ;
+syms n z;
+x =(2) ^n
+X= symsum (x*(z^(-n)),n ,0, %inf );
+//Display the result in command window
+disp (X,"Z-transform of 2^n u(n) is:");
+disp('ROC is the Region mod(z) > 2');
\ No newline at end of file diff --git a/52/CH2/EX2.7/Example2_7.sce b/52/CH2/EX2.7/Example2_7.sce new file mode 100755 index 000000000..55142ef4a --- /dev/null +++ b/52/CH2/EX2.7/Example2_7.sce @@ -0,0 +1,14 @@ +//Example 2.7
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of [3(3^n)-4(2)^n] u(n)
+clear;
+clc ;
+close ;
+syms n z;
+x1 =(3) ^(n);
+X1= symsum (3* x1 *(z^(-n)),n ,0, %inf );
+x2 =(4) ^(n);
+X2= symsum (4* x2 *(z^(-n)),n ,0, %inf );
+X = (X1 -X2);
+//Display the result in command window
+disp (X,"Z-transform of [3(3^n)-4(2)^n] u(n) is:");
\ No newline at end of file diff --git a/52/CH2/EX2.8.a/Example2_8_a.sce b/52/CH2/EX2.8.a/Example2_8_a.sce new file mode 100755 index 000000000..a7653ac47 --- /dev/null +++ b/52/CH2/EX2.8.a/Example2_8_a.sce @@ -0,0 +1,11 @@ +//Example 2.8 (a)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z transform of cos(Wo*n)
+clc;
+syms Wo n z;
+x1=exp(sqrt(-1)*Wo*n);
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=exp(-sqrt(-1)*Wo*n);
+X2=symsum(x2*(z^-n),n,0,%inf);
+X=(X1+X2)/2;
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/52/CH2/EX2.9/Example2_9.sce b/52/CH2/EX2.9/Example2_9.sce new file mode 100755 index 000000000..984425f71 --- /dev/null +++ b/52/CH2/EX2.9/Example2_9.sce @@ -0,0 +1,11 @@ +//Example 2.9
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Z- transform of (1/3)^n u(n-1)
+clear;
+clc ;
+close ;
+syms n z;
+x =(1/3)^n;
+X= (1/z)*symsum (x*(z^(-n)),n ,0, %inf );
+//Display the result in command window
+disp (X,"Z-transform of (1/3)^n u(n-1) is:");
\ No newline at end of file diff --git a/52/CH3/EX3.1/Example3_1.sce b/52/CH3/EX3.1/Example3_1.sce new file mode 100755 index 000000000..fe120a9c0 --- /dev/null +++ b/52/CH3/EX3.1/Example3_1.sce @@ -0,0 +1,15 @@ +//Example 3.1
+//Program to Compute the DFT of a Sequence x[n]=[1,1,0,0]
+//and IDFT of a Sequence Y[k]=[1,0,1,0]
+clear;
+clc ;
+close ;
+x = [1,1,0,0];
+//DFT Computation
+X = fft (x , -1);
+Y = [1,0,1,0];
+//IDFT Computation
+y = fft (Y , 1);
+//Display sequence X[k] and y[n] in command window
+disp(X,"X[k]=");
+disp(y,"y[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.11/Example3_11.sce b/52/CH3/EX3.11/Example3_11.sce new file mode 100755 index 000000000..d70d6fd65 --- /dev/null +++ b/52/CH3/EX3.11/Example3_11.sce @@ -0,0 +1,15 @@ +//Example 3.11
+//Program to Compute the 8-point DFT of the following sequences
+//x1[n]=[1,0,0,0,0,1,1,1]
+//x2[n]=[0,0,1,1,1,1,0,0]
+clear;
+clc ;
+close ;
+x1=[1,0,0,0,0,1,1,1];
+x2=[0,0,1,1,1,1,0,0];
+//DFT Computation
+X1 = fft (x1 , -1);
+X2 = fft (x2 , -1);
+//Display sequences X1[k] and X2[k] in command window
+disp(X1,"X1[k]=");
+disp(X2,"X2[k]=");
\ No newline at end of file diff --git a/52/CH3/EX3.13/Example3_13.sce b/52/CH3/EX3.13/Example3_13.sce new file mode 100755 index 000000000..028771f05 --- /dev/null +++ b/52/CH3/EX3.13/Example3_13.sce @@ -0,0 +1,26 @@ +//Example 3.13
+//Program to Compute circular convolution of following sequences
+//x1[n]=[1,-1,-2,3,-1]
+//x2[n]=[1,2,3]
+clear;
+clc ;
+close ;
+x1=[1,-1,-2,3,-1];
+x2=[1,2,3];
+//Loop for zero padding the smaller sequence out of the two
+n1=length(x1);
+n2=length(x2);
+n3=n2-n1;
+if (n3>=0) then
+ x1=[x1,zeros(1,n3)];
+else
+ x2=[x2,zeros(1,-n3)];
+end
+//DFT Computation
+X1=fft(x1,-1);
+X2=fft(x2,-1);
+Y=X1.*X2;
+//IDFT Computation
+y=fft(Y,1);
+//Display sequence y[n] in command window
+disp(y,"y[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.14/Example3_14.sce b/52/CH3/EX3.14/Example3_14.sce new file mode 100755 index 000000000..e7c14b9da --- /dev/null +++ b/52/CH3/EX3.14/Example3_14.sce @@ -0,0 +1,17 @@ +//Example 3.14
+//Program to Compute circular convolution of following sequences
+//x1[n]=[1,2,2,1]
+//x2[n]=[1,2,3,1]
+clear;
+clc ;
+close ;
+x1=[1,2,2,1];
+x2=[1,2,3,1];
+//DFT Computation
+X1=fft(x1,-1);
+X2=fft(x2,-1);
+Y=X1.*X2;
+//IDFT Computation
+y=fft(Y,1);
+//Display sequence y[n] in command window
+disp(y,"y[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.15/Example3_15.sce b/52/CH3/EX3.15/Example3_15.sce new file mode 100755 index 000000000..3862ebe46 --- /dev/null +++ b/52/CH3/EX3.15/Example3_15.sce @@ -0,0 +1,17 @@ +//Example 3.15
+//Program to Compute x3[n] where X3[k]=X1[k].X2[k]
+//x1[n]=[1,2,3,4]
+//x2[n]=[1,1,2,2]
+clear;
+clc ;
+close ;
+x1=[1,2,3,4];
+x2=[1,1,2,2];
+//DFT Computation
+X1=fft(x1,-1);
+X2=fft(x2,-1);
+X3=X1.*X2;
+//IDFT Computation
+x3=fft(X3,1);
+//Display sequence x3[n] in command window
+disp(x3,"x3[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.16/Example3_16.sce b/52/CH3/EX3.16/Example3_16.sce new file mode 100755 index 000000000..985fb2113 --- /dev/null +++ b/52/CH3/EX3.16/Example3_16.sce @@ -0,0 +1,17 @@ +//Example 3.16
+//Program to Compute circular convolution of following sequences
+//x1[n]=[1,1,2,1]
+//x2[n]=[1,2,3,4]
+clear;
+clc ;
+close ;
+x1=[1,1,2,1];
+x2=[1,2,3,4];
+//DFT Computation
+X1=fft(x1,-1);
+X2=fft(x2,-1);
+X3=X1.*X2;
+//IDFT Computation
+x3=fft(X3,1);
+//Display sequence x3[n] in command window
+disp(x3,"x3[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.17/Example3_17.sce b/52/CH3/EX3.17/Example3_17.sce new file mode 100755 index 000000000..e2d971910 --- /dev/null +++ b/52/CH3/EX3.17/Example3_17.sce @@ -0,0 +1,17 @@ +//Example 3.17
+//Program to Compute y[n] where Y[k]=X1[k].X2[k]
+//x1[n]=[0,1,2,3,4]
+//x2[n]=[0,1,0,0,0]
+clear;
+clc ;
+close ;
+x1=[0,1,2,3,4];
+x2=[0,1,0,0,0];
+//DFT Computation
+X1=fft(x1,-1);
+X2=fft(x2,-1);
+Y=X1.*X2;
+//IDFT Computation
+y=round(fft(Y,1));
+//Display sequence y[n] in command window
+disp(y,"y[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.18/Example3_18.sce b/52/CH3/EX3.18/Example3_18.sce new file mode 100755 index 000000000..fc399e533 --- /dev/null +++ b/52/CH3/EX3.18/Example3_18.sce @@ -0,0 +1,36 @@ +//Example 3.18
+//Program to Compute output responce of following sequences
+//x[n]=[1,2,3,1]
+//h[n]=[1,1,1]
+//(1)Linear Convolution
+//(2)Circular Convolution
+//(3)Circular Convolution with zero padding
+clear;
+clc ;
+close ;
+x=[1,2,3,1];
+h=[1,1,1];
+//(1)Linear Convolution Computation
+ylinear=convol (x,h);
+//Display Linear Convoluted Sequence y[n] in command window
+disp(ylinear,"ylinear[n]=");
+//(2)Circular Convolution Computation
+//Now zero padding in h[n] sequence to make length of x[n] and h[n] equal
+h1=[h,zeros(1,1)];
+//Now Performing Circular Convolution by DFT method
+X=fft(x,-1);
+H=fft(h1,-1);
+Y=X.*H;
+ycircular=fft(Y,1);
+//Display Circular Convoluted Sequence y[n] in command window
+disp(ycircular,"ycircular[n]=");
+//(3)Circular Convolution Computation with zero Padding
+x2=[x,zeros(1,2)];
+h2=[h,zeros(1,3)];
+//Now Performing Circular Convolution by DFT method
+X2=fft(x2,-1);
+H2=fft(h2,-1);
+Y2=X2.*H2;
+ycircularp=fft(Y2,1);
+//Display Circular Convoluted Sequence with zero Padding y[n] in command window
+disp(ycircularp,"ycircularp[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.2/Example3_2.jpg b/52/CH3/EX3.2/Example3_2.jpg Binary files differnew file mode 100755 index 000000000..f55657202 --- /dev/null +++ b/52/CH3/EX3.2/Example3_2.jpg diff --git a/52/CH3/EX3.2/Example3_2.sce b/52/CH3/EX3.2/Example3_2.sce new file mode 100755 index 000000000..f927687fd --- /dev/null +++ b/52/CH3/EX3.2/Example3_2.sce @@ -0,0 +1,53 @@ +//Example 3.2
+//Program to Compute the DFT of a Sequence x[n]=1, 0<=n<=2; and 0 otherwise
+//for N=4 and N=8. Plot Magnitude and phase plots of each.
+clear;
+clc ;
+close ;
+//N=4
+x1 = [1,1,1,0];
+//DFT Computation
+X1 = fft (x1 , -1);
+//N=8
+x2 = [1,1,1,0,0,0,0,0];
+//DFT Computation
+X2 = fft (x2 , -1);
+//Display sequence X1[k] and X2[k] in command window
+disp(X1,"X1[k]=");
+disp(X2,"X2[k]=");
+//Plots for N=4
+n1=0:1:3;
+subplot(2,2,1);
+a = gca ();
+a.y_location ="origin";
+a.x_location ="origin";
+plot2d3(n1,abs(X1),2);
+poly1=a.children(1).children (1);
+poly1.thickness=2;
+xtitle('N=4','k','|X1(k)|');
+subplot(2,2,2);
+a = gca ();
+a.y_location ="origin";
+a.x_location ="origin";
+plot2d3(n1,atan(imag(X1),real(X1)),5);
+poly1=a.children(1).children (1);
+poly1.thickness=2;
+xtitle('N=4','k','<X1(k)');
+//Plots for N=8
+n2=0:1:7;
+subplot(2,2,3);
+a = gca ();
+a.y_location ="origin";
+a.x_location ="origin";
+plot2d3(n2,abs(X2),2);
+poly1=a.children(1).children (1);
+poly1.thickness=2;
+xtitle('N=8','k','|X2(k)|');
+subplot(2,2,4);
+a = gca ();
+a.y_location ="origin";
+a.x_location ="origin";
+plot2d3(n2,atan(imag(X2),real(X2)),5);
+poly1=a.children(1).children (1);
+poly1.thickness=2;
+xtitle('N=8','k','<X2(k)');
\ No newline at end of file diff --git a/52/CH3/EX3.20/Example3_20.sce b/52/CH3/EX3.20/Example3_20.sce new file mode 100755 index 000000000..a12978280 --- /dev/null +++ b/52/CH3/EX3.20/Example3_20.sce @@ -0,0 +1,13 @@ +//Example 3.20
+//Program to Compute Linear Convolution of following sequences
+//x[n]=[3,-1,0,1,3,2,0,1,2,1]
+//h[n]=[1,1,1]
+clear;
+clc ;
+close ;
+x=[3,-1,0,1,3,2,0,1,2,1];
+h=[1,1,1];
+// Linear Convolution Computation
+y=convol (x,h);
+//Display Sequence y[n] in command window
+disp(y,"y[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.21/Example3_21.sce b/52/CH3/EX3.21/Example3_21.sce new file mode 100755 index 000000000..a2e9eb126 --- /dev/null +++ b/52/CH3/EX3.21/Example3_21.sce @@ -0,0 +1,13 @@ +//Example 3.21
+//Program to Compute Linear Convolution of following sequences
+//x[n]=[1,2,-1,2,3,-2,-3,-1,1,1,2,-1]
+//h[n]=[1,2]
+clear;
+clc ;
+close ;
+x=[1,2,-1,2,3,-2,-3,-1,1,1,2,-1];
+h=[1,2];
+// Linear Convolution Computation
+y=convol (x,h);
+//Display Sequence y[n] in command window
+disp(y,"y[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.23.a/Example3_23_a.sce b/52/CH3/EX3.23.a/Example3_23_a.sce new file mode 100755 index 000000000..e5105c7e4 --- /dev/null +++ b/52/CH3/EX3.23.a/Example3_23_a.sce @@ -0,0 +1,8 @@ +//Example 3.23 (a)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//N point DFT of delta(n)
+clc;
+syms n k N;
+x=1;
+X=symsum(x*exp(-%i*2*%pi*n*k/N),n,0,0);
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/52/CH3/EX3.23.b/Example3_23_b.sce b/52/CH3/EX3.23.b/Example3_23_b.sce new file mode 100755 index 000000000..ad4c8aef0 --- /dev/null +++ b/52/CH3/EX3.23.b/Example3_23_b.sce @@ -0,0 +1,8 @@ +//Example 3.23 (b)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//N point DFT of delta(n-no)
+clc;
+syms n k N no;
+x=1;
+X=symsum(x*exp(-%i*2*%pi*n*k/N),n,-no,-no);
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/52/CH3/EX3.23.c/Example3_23_c.sce b/52/CH3/EX3.23.c/Example3_23_c.sce new file mode 100755 index 000000000..4db5d7025 --- /dev/null +++ b/52/CH3/EX3.23.c/Example3_23_c.sce @@ -0,0 +1,8 @@ +//Example 3.23 (c)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//N point DFT of a^n
+clc;
+syms a n k N;
+x=a^n;
+X=symsum(x*exp(-%i*2*%pi*n*k/N),n,0,N-1);
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/52/CH3/EX3.23.d/Example3_23_d.sce b/52/CH3/EX3.23.d/Example3_23_d.sce new file mode 100755 index 000000000..bc6373214 --- /dev/null +++ b/52/CH3/EX3.23.d/Example3_23_d.sce @@ -0,0 +1,8 @@ +//Example 3.23 (d)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//N point DFT of x(n)=1, 0<=n<=N/2-1
+clc;
+syms n k N;
+x=1;
+X=symsum(x*exp(-%i*2*%pi*n*k/N),n,0,(N/2)-1);
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/52/CH3/EX3.23.e/Example3_23_e.sce b/52/CH3/EX3.23.e/Example3_23_e.sce new file mode 100755 index 000000000..2a2290156 --- /dev/null +++ b/52/CH3/EX3.23.e/Example3_23_e.sce @@ -0,0 +1,8 @@ +//Example 3.23 (e)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//N point DFT of x(n)=exp(%i*2*%pi*ko*n/N);
+clc;
+syms n k N ko;
+x=exp(%i*2*%pi*ko*n/N);
+X=symsum(x*exp(-%i*2*%pi*n*k/N),n,0,(N/2)-1);
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/52/CH3/EX3.23.f/Example3_23_f.sce b/52/CH3/EX3.23.f/Example3_23_f.sce new file mode 100755 index 000000000..653729e4f --- /dev/null +++ b/52/CH3/EX3.23.f/Example3_23_f.sce @@ -0,0 +1,8 @@ +//Example 3.23 (f)
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//N point DFT of x(n)=1, for n=even and 0, for n=odd
+clc;
+syms n k N ;
+x=1;//x(2n)=1,for all n
+X=symsum(x*exp(-%i*4*%pi*n*k/N),n,0,N/2-1);
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/52/CH3/EX3.24/Example3_24.sce b/52/CH3/EX3.24/Example3_24.sce new file mode 100755 index 000000000..5c911eece --- /dev/null +++ b/52/CH3/EX3.24/Example3_24.sce @@ -0,0 +1,12 @@ +//Example 3.24
+//Program to Compute the DFT of the Sequence x[n]=(-1)^n, for N=4
+clear;
+clc ;
+close ;
+N=4;
+n=0:1:N-1;
+x=(-1)^n;
+//DFT Computation
+X = fft (x,-1);
+//Display Sequence X[k] in command window
+disp(X,"X[k]=");
\ No newline at end of file diff --git a/52/CH3/EX3.25/Example3_25.sce b/52/CH3/EX3.25/Example3_25.sce new file mode 100755 index 000000000..ca3cdfda4 --- /dev/null +++ b/52/CH3/EX3.25/Example3_25.sce @@ -0,0 +1,20 @@ +//Example 3.25
+//Program to Compute the 8-point Circular Convolution of the Sequences
+//x1[n]=[1,1,1,1,0,0,0,0]
+//x2[n]=sin(3*pi*n/8)
+clear;
+clc ;
+close ;
+x1=[1,1,1,1,0,0,0,0];
+n=0:1:7;
+pi=22/7;
+x2=sin(3*pi*n/8);
+//DFT Computation
+X1=fft (x1,-1);
+X2=fft (x2,-1);
+//Circular Convolution using DFT
+Y=X1.*X2;
+//IDFT Computation
+y= fft (Y,1);
+//Display sequence y[n] in command window
+disp(y,"y[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.26/Example3_26.sce b/52/CH3/EX3.26/Example3_26.sce new file mode 100755 index 000000000..5ca6fbe49 --- /dev/null +++ b/52/CH3/EX3.26/Example3_26.sce @@ -0,0 +1,14 @@ +//Example 3.26
+//Program to Compute the Linear Convolution of the following Sequences
+//x[n]=[1,-1,1]
+//h[n]=[2,2,1]
+clear;
+clc ;
+close ;
+x=[1,-1,1];
+h=[2,2,1];
+//Convolution Computation
+y= convol(x,h);
+//Display sequence y[n] in command window
+disp(y,"y[n]=");
+
diff --git a/52/CH3/EX3.27.a/Example3_27_a.sce b/52/CH3/EX3.27.a/Example3_27_a.sce new file mode 100755 index 000000000..ad07aa5ec --- /dev/null +++ b/52/CH3/EX3.27.a/Example3_27_a.sce @@ -0,0 +1,16 @@ +//Example 3.27 (a)
+//Program to Compute the Convolution of the following Sequences
+//x1[n]=[1,1,1]
+//x2[n]=[2,-1,2]
+clear;
+clc ;
+close ;
+x1=[1,1,1];
+x2=[2,-1,2];
+//Convolution Computation
+X1=fft (x1,-1);
+X2=fft (x2,-1);
+Y=X1.*X2;
+y=fft (Y,1);
+//Display Sequence y[n] in command window
+disp(y,"y[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.27.b/Example3_27_b.sce b/52/CH3/EX3.27.b/Example3_27_b.sce new file mode 100755 index 000000000..4008964e1 --- /dev/null +++ b/52/CH3/EX3.27.b/Example3_27_b.sce @@ -0,0 +1,16 @@ +//Example 3.27 (b)
+//Program to Compute the Convolution of the following Sequences
+//x1[n]=[1,1,-1,-1,0]
+//x2[n]=[1,0,-1,0,1]
+clear;
+clc ;
+close ;
+x1=[1,1,-1,-1,0];
+x2=[1,0,-1,0,1];
+//Convolution Computation
+X1=fft (x1,-1);
+X2=fft (x2,-1);
+Y=X1.*X2;
+y= fft (Y,1);
+//Display Sequence y[n] in command window
+disp(y,"y[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.3/Example3_3.sce b/52/CH3/EX3.3/Example3_3.sce new file mode 100755 index 000000000..72af7fbb4 --- /dev/null +++ b/52/CH3/EX3.3/Example3_3.sce @@ -0,0 +1,10 @@ +//Example 3.3
+//Program to Compute the 8-point DFT of the Sequence x[n]=[1,1,1,1,1,1,0,0]
+clear;
+clc ;
+close ;
+x = [1,1,1,1,1,1,0,0];
+//DFT Computation
+X = fft (x , -1);
+//Display sequence X[k] in command window
+disp(X,"X[k]=");
\ No newline at end of file diff --git a/52/CH3/EX3.30/Example3_30.sce b/52/CH3/EX3.30/Example3_30.sce new file mode 100755 index 000000000..4668395e3 --- /dev/null +++ b/52/CH3/EX3.30/Example3_30.sce @@ -0,0 +1,16 @@ +//Example 3.30
+//Program to Calculate N from given data
+//fm=5000Hz
+//df=50Hz
+//t=0.5sec
+clear;
+clc ;
+close ;
+fm=5000 //Hz
+df=50 //Hz
+t=0.5 //sec
+N1=2*fm/df;
+N=2;
+while N<=N1, N=N*2,end
+//Displaying the value of N in command window
+disp(N,"N=");
\ No newline at end of file diff --git a/52/CH3/EX3.32/Example3_32.jpg b/52/CH3/EX3.32/Example3_32.jpg Binary files differnew file mode 100755 index 000000000..a1d06f714 --- /dev/null +++ b/52/CH3/EX3.32/Example3_32.jpg diff --git a/52/CH3/EX3.32/Example3_32.sce b/52/CH3/EX3.32/Example3_32.sce new file mode 100755 index 000000000..1e1905ab1 --- /dev/null +++ b/52/CH3/EX3.32/Example3_32.sce @@ -0,0 +1,22 @@ +//Example 3.32
+//Program to plot the result of the given sequence
+//X[k]=[1,2,2,1,0,2,1,2]
+//y[n]=x[n/2] for n=even,0 for n=odd
+clear;
+clc ;
+close ;
+X=[1,2,2,1,0,2,1,2];
+x = fft (X , 1);
+y=[x(1),0,x(2),0,x(3),0,x(4),0,x(5),0,x(6),0,x(7),0,x(8),0];
+Y = fft (y , -1);
+//Display sequence Y[k] and in command window
+disp(Y,"Y[k]=");
+//Plotting the sequence Y[k]
+k=0:1:15;
+a = gca ();
+a.y_location ="origin";
+a.x_location ="origin";
+plot2d3(k,Y,2);
+poly1=a.children(1).children (1);
+poly1.thickness=2;
+xtitle('Plot of Y(k)','k','Y(k)');
\ No newline at end of file diff --git a/52/CH3/EX3.36/Example3_36.sce b/52/CH3/EX3.36/Example3_36.sce new file mode 100755 index 000000000..1b08ebc44 --- /dev/null +++ b/52/CH3/EX3.36/Example3_36.sce @@ -0,0 +1,12 @@ +//Example 3.36
+//Program to Compute the IDFT of the following Sequence
+//X[k]=[12,-1.5+j2.598,-1.5+j0.866,0,-1.5-j0.866,-1.5-j2.598]
+clear;
+clc ;
+close ;
+j=sqrt(-1);
+X=[12,-1.5+j*2.598,-1.5+j*0.866,0,-1.5-j*0.866,-1.5-j*2.598];
+//IDFT Computation
+x = fft (X , 1);
+//Display Sequence x[n] in command window
+disp(round(x),"x[n]=");
\ No newline at end of file diff --git a/52/CH3/EX3.4/Example3_4.sce b/52/CH3/EX3.4/Example3_4.sce new file mode 100755 index 000000000..9b8af1d59 --- /dev/null +++ b/52/CH3/EX3.4/Example3_4.sce @@ -0,0 +1,11 @@ +//Example 3.4
+//Program to Compute the IDFT of the Sequence X[k]=[5,0,1-j,0,1,0,1+j,0]
+clear;
+clc ;
+close ;
+j=sqrt(-1);
+X = [5,0,1-j,0,1,0,1+j,0]
+//IDFT Computation
+x = fft (X , 1);
+//Display sequences x[n]in command window
+disp(x,"x[n]=");
diff --git a/52/CH3/EX3.7/Example3_7.jpg b/52/CH3/EX3.7/Example3_7.jpg Binary files differnew file mode 100755 index 000000000..c6f68c8e0 --- /dev/null +++ b/52/CH3/EX3.7/Example3_7.jpg diff --git a/52/CH3/EX3.7/Example3_7.sce b/52/CH3/EX3.7/Example3_7.sce new file mode 100755 index 000000000..2873ebb16 --- /dev/null +++ b/52/CH3/EX3.7/Example3_7.sce @@ -0,0 +1,27 @@ +//Example 3.7
+//Program to Compute circular convolution of following sequences
+//x[n]=[1,2,2,1,0]
+//Y[k]=exp(-j*4*pi*k/5).X[k]
+clear;
+clc ;
+close ;
+x=[1,2,2,1,0];
+X=fft(x,-1);
+k=0:1:4;
+j=sqrt(-1);
+pi=22/7;
+H=exp(-j*4*pi*k/5);
+Y=H.*X;
+//IDFT Computation
+y=fft(Y,1);
+//Display sequence y[n] in command window
+disp(round(y),"y[n]=");
+//Plots
+n=0:1:4;
+a = gca ();
+a.y_location ="origin";
+a.x_location ="origin";
+plot2d3(n,round(y),5);
+poly1=a.children(1).children (1);
+poly1.thickness=2;
+xtitle('Plot of sequence y[n]','n','y[n]');
\ No newline at end of file diff --git a/52/CH3/EX3.9/Example3_9.sce b/52/CH3/EX3.9/Example3_9.sce new file mode 100755 index 000000000..04eeebc3f --- /dev/null +++ b/52/CH3/EX3.9/Example3_9.sce @@ -0,0 +1,12 @@ +//Example 3.9
+//Program to remaining samples of the sequence
+//X(0)=12,X(1)=-1+j3,X(2)=3+j4,X(3)=1-j5,X(4)=-2+j2,X(5)=6+j3,X(6)=-2-j3,X(7)=10
+clear;
+clc ;
+close ;
+j=sqrt(-1);
+z=1;
+X(0+z)=12,X(1+z)=-1+j*3,X(2+z)=3+j*4,X(3+z)=1-j*5,X(4+z)=-2+j*2,X(5+z)=6+j*3,X(6+z)=-2-j*3,X(7+z)=10;
+for a=9:1:14 do X(a)=conj(X(16-a)), end;
+//Display the complete sequence X[k] in command window
+disp(X,"X[k]=");
\ No newline at end of file diff --git a/52/CH4/EX4.10/Example4_10.sce b/52/CH4/EX4.10/Example4_10.sce new file mode 100755 index 000000000..9f05bdc5f --- /dev/null +++ b/52/CH4/EX4.10/Example4_10.sce @@ -0,0 +1,10 @@ +//Example 4.10
+//Program to Compute the 4-point DFT of a Sequence x[n]=[0,1,2,3]
+//using DIT-DIF Algorithm.
+clear;
+clc ;
+close ;
+x = [0,1,2,3];
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.11/Example4_11.sce b/52/CH4/EX4.11/Example4_11.sce new file mode 100755 index 000000000..0bbdaa1de --- /dev/null +++ b/52/CH4/EX4.11/Example4_11.sce @@ -0,0 +1,11 @@ +//Example 4.11
+//Program to Compute the IDFT of a Sequence using DIT Algorithm.
+//X[k] = [7,-0.707-j0.707,-j,0.707-j0.707,1,0.707+j0.707,j,-0.707+j0.707]
+clear;
+clc ;
+close ;
+j=sqrt(-1);
+X = [7,-0.707-j*0.707,-j,0.707-j*0.707,1,0.707+j*0.707,j,-0.707+j*0.707];
+//Inverse FFT Computation
+x = fft (X , 1);
+disp(x,'x(n) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.12/Example4_12.sce b/52/CH4/EX4.12/Example4_12.sce new file mode 100755 index 000000000..8e7cc4fa9 --- /dev/null +++ b/52/CH4/EX4.12/Example4_12.sce @@ -0,0 +1,10 @@ +//Example 4.12
+//Program to Compute the 8-point DFT of a Sequence
+//x[n]=[0.5,0.5,0.5,0.5,0,0,0,0] using radix-2 DIT Algorithm.
+clear;
+clc ;
+close ;
+x=[0.5,0.5,0.5,0.5,0,0,0,0];
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.13/Example4_13.sce b/52/CH4/EX4.13/Example4_13.sce new file mode 100755 index 000000000..fbaa5eb00 --- /dev/null +++ b/52/CH4/EX4.13/Example4_13.sce @@ -0,0 +1,10 @@ +//Example 4.13
+//Program to Compute the 8-point DFT of a Sequence
+//x[n]=[0.5,0.5,0.5,0.5,0,0,0,0] using radix-2 DIF Algorithm.
+clear;
+clc ;
+close ;
+x=[0.5,0.5,0.5,0.5,0,0,0,0];
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.14/Example4_14.sce b/52/CH4/EX4.14/Example4_14.sce new file mode 100755 index 000000000..cb82266f0 --- /dev/null +++ b/52/CH4/EX4.14/Example4_14.sce @@ -0,0 +1,10 @@ +//Example 4.14
+//Program to Compute the 4-point DFT of a Sequence x[n]=[1,-1,1,-1]
+//using DIT Algorithm.
+clear;
+clc ;
+close ;
+x=[1,-1,1,-1];
+//FFT Computation
+X =fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.15/Example4_15.sce b/52/CH4/EX4.15/Example4_15.sce new file mode 100755 index 000000000..d6947e6c8 --- /dev/null +++ b/52/CH4/EX4.15/Example4_15.sce @@ -0,0 +1,10 @@ +//Example 4.15
+//Program to Compute the 4-point DFT of a Sequence x[n]=[1,0,0,1]
+//using DIF Algorithm.
+clear;
+clc ;
+close ;
+x=[1,0,0,1];
+//FFT Computation
+X =fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.16.a/Example4_16_a.sce b/52/CH4/EX4.16.a/Example4_16_a.sce new file mode 100755 index 000000000..ec149bf15 --- /dev/null +++ b/52/CH4/EX4.16.a/Example4_16_a.sce @@ -0,0 +1,10 @@ +//Example 4.16 (a)
+//Program to Evaluate and Compare the 8-point DFT of the given Sequence
+//x1[n]=1, -3<=n<=3 using DIT-FFT Algorithm.
+clear;
+clc ;
+close ;
+x1=[1,1,1,1,0,1,1,1];
+//FFT Computation
+X1 = fft (x1 , -1);
+disp(X1,'X1(k) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.16.b/Example4_16_b.sce b/52/CH4/EX4.16.b/Example4_16_b.sce new file mode 100755 index 000000000..00c84c931 --- /dev/null +++ b/52/CH4/EX4.16.b/Example4_16_b.sce @@ -0,0 +1,10 @@ +//Example 4.16 (b)
+//Program to Evaluate and Compare the 8-point DFT of the given Sequence
+//x2[n]=1, 0<=n<=6 using DIT-FFT Algorithm.
+clear;
+clc ;
+close ;
+x2=[1,1,1,1,1,1,1,0];
+//FFT Computation
+X2 = fft (x2 , -1);
+disp(X2,'X2(k) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.17/Example4_17.sce b/52/CH4/EX4.17/Example4_17.sce new file mode 100755 index 000000000..997e2a694 --- /dev/null +++ b/52/CH4/EX4.17/Example4_17.sce @@ -0,0 +1,11 @@ +//Example 4.17
+//Program to find the IDFT of the Sequence using DIF Algorithm.
+//X[k]= [4,1-j2.414,0,1-j0.414,0,1+j0.414,0,1+j2.414]
+clear;
+clc ;
+close ;
+j=sqrt(-1);
+X= [4,1-j*2.414,0,1-j*0.414,0,1+j*0.414,0,1+j*2.414];
+//Inverse FFT Computation
+x = fft (X , 1);
+disp(x,'x(n) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.18/Example4_18.sce b/52/CH4/EX4.18/Example4_18.sce new file mode 100755 index 000000000..14e3c1212 --- /dev/null +++ b/52/CH4/EX4.18/Example4_18.sce @@ -0,0 +1,11 @@ +//Example 4.18
+//Program to find the IDFT of the Sequence X[k]= [10,-2+j2,-2,-2-j2]
+//using DIT Algorithm.
+clear;
+clc ;
+close ;
+j=sqrt(-1);
+X = [10,-2+j*2,-2,-2-j*2];
+//Inverse FFT Computation
+x = fft (X , 1);
+disp(x,'x(n) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.19/Example4_19.sce b/52/CH4/EX4.19/Example4_19.sce new file mode 100755 index 000000000..b76a97815 --- /dev/null +++ b/52/CH4/EX4.19/Example4_19.sce @@ -0,0 +1,9 @@ +//Example 4.19
+//Program to Compute the FFT of given Sequence x[n]=[1,0,0,0,0,0,0,0].
+clear;
+clc ;
+close ;
+x = [1,0,0,0,0,0,0,0];
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.20/Example4_20.sce b/52/CH4/EX4.20/Example4_20.sce new file mode 100755 index 000000000..aa279014e --- /dev/null +++ b/52/CH4/EX4.20/Example4_20.sce @@ -0,0 +1,10 @@ +//Example 4.20
+//Program to Compute the 8-point DFT of given Sequence
+//x[n]=[2,2,2,2,1,1,1,1] using DIT, radix-2,FFT Algorithm.
+clear;
+clc ;
+close ;
+x = [2,2,2,2,1,1,1,1];
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.21/Example4_21.sce b/52/CH4/EX4.21/Example4_21.sce new file mode 100755 index 000000000..d0d9f8469 --- /dev/null +++ b/52/CH4/EX4.21/Example4_21.sce @@ -0,0 +1,10 @@ +//Example 4.21
+//Program to Compute the DFT of given Sequence
+//x[n]=[1,-1,-1,-1,1,1,1,-1] using DIT-FFT Algorithm.
+clear;
+clc ;
+close ;
+x = [1,-1,-1,-1,1,1,1,-1];
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.22/Example4_22.sce b/52/CH4/EX4.22/Example4_22.sce new file mode 100755 index 000000000..c901ac385 --- /dev/null +++ b/52/CH4/EX4.22/Example4_22.sce @@ -0,0 +1,12 @@ +//Example 4.22
+//Program to Compute the DFT of given Sequence
+//x[n]=2^n and N=8 using DIT-FFT Algorithm.
+clear;
+clc ;
+close ;
+N=8;
+n=0:1:N-1;
+x =2^n;
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.23/Example4_23.sce b/52/CH4/EX4.23/Example4_23.sce new file mode 100755 index 000000000..b8a307f1b --- /dev/null +++ b/52/CH4/EX4.23/Example4_23.sce @@ -0,0 +1,13 @@ +//Example 4.23
+//Program to Compute the DFT of given Sequence
+//x[n]=cos(n*pi/2), and N=4 using DIF-FFT Algorithm.
+clear;
+clc ;
+close ;
+N=4;
+pi=22/7;
+n=0:1:N-1;
+x =cos(n*pi/2);
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.24/Example4_24.sce b/52/CH4/EX4.24/Example4_24.sce new file mode 100755 index 000000000..baa8a5dd2 --- /dev/null +++ b/52/CH4/EX4.24/Example4_24.sce @@ -0,0 +1,10 @@ +//Example 4.24
+//Program to Compute the 8-point DFT of given Sequence
+//x[n]=[0,1,2,3,4,5,6,7] using DIF, radix-2,FFT Algorithm.
+clear;
+clc ;
+close ;
+x = [0,1,2,3,4,5,6,7];
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.3/Example4_3.sce b/52/CH4/EX4.3/Example4_3.sce new file mode 100755 index 000000000..9e58c8ec7 --- /dev/null +++ b/52/CH4/EX4.3/Example4_3.sce @@ -0,0 +1,17 @@ +//Example 4.3
+//Program to calculate shortest sequence N such that algorithm B runs //faster than A
+clear;
+clc ;
+close ;
+i=0;
+N=32; //Given
+//Calculation of Twiddle factor exponents for each stage
+while 1==1
+ i=i+1;
+ N=2^i;
+ A=N^2;
+ B=5*N*log2(N);
+ if A>B then break;
+ end;
+end
+disp(N,'SHORTEST SEQUENCE N =');
\ No newline at end of file diff --git a/52/CH4/EX4.4/Example4_4.sce b/52/CH4/EX4.4/Example4_4.sce new file mode 100755 index 000000000..b97272b79 --- /dev/null +++ b/52/CH4/EX4.4/Example4_4.sce @@ -0,0 +1,15 @@ +//Example 4.4
+//Program to calculate Twiddle factor exponents for each stage
+clear;
+clc ;
+close ;
+N=32; //Given
+//Calculation of Twiddle factor exponents for each stage
+for m=1:5
+ disp(m,'Stage: m =');
+ disp(' k =');
+ for t=0:(2^(m-1)-1)
+ k=N*t/2^m;
+ disp(k);
+ end
+end
\ No newline at end of file diff --git a/52/CH4/EX4.6/Example4_6.sce b/52/CH4/EX4.6/Example4_6.sce new file mode 100755 index 000000000..1a4530c8f --- /dev/null +++ b/52/CH4/EX4.6/Example4_6.sce @@ -0,0 +1,10 @@ +//Example 4.6
+//Program to find the DFT of a Sequence x[n]=[1,2,3,4,4,3,2,1]
+//using DIT Algorithm.
+clear;
+clc ;
+close ;
+x = [1,2,3,4,4,3,2,1];
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH4/EX4.8/Example4_8.sce b/52/CH4/EX4.8/Example4_8.sce new file mode 100755 index 000000000..2f3f136d9 --- /dev/null +++ b/52/CH4/EX4.8/Example4_8.sce @@ -0,0 +1,10 @@ +//Example 4.8
+//Program to find the DFT of a Sequence x[n]=[1,2,3,4,4,3,2,1]
+//using DIF Algorithm.
+clear;
+clc ;
+close ;
+x = [1,2,3,4,4,3,2,1];
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
diff --git a/52/CH4/EX4.9/Example4_9.sce b/52/CH4/EX4.9/Example4_9.sce new file mode 100755 index 000000000..c6b4f56d8 --- /dev/null +++ b/52/CH4/EX4.9/Example4_9.sce @@ -0,0 +1,10 @@ +//Example 4.9
+//Program to find the 8-point DFT of a Sequence x[n]=1, 0<=n<=7
+//using DIT,DIF Algorithm.
+clear;
+clc ;
+close ;
+x = [1,1,1,1,1,1,1,1];
+//FFT Computation
+X = fft (x , -1);
+disp(X,'X(z) = ');
\ No newline at end of file diff --git a/52/CH5/EX5.1/Example5_1.sce b/52/CH5/EX5.1/Example5_1.sce new file mode 100755 index 000000000..c62f2ba47 --- /dev/null +++ b/52/CH5/EX5.1/Example5_1.sce @@ -0,0 +1,11 @@ +//Example 5.1
+//To Find out the order of the filter
+clear;
+clc ;
+close ;
+ap=1;//db
+as=30;//db
+op=200;//rad/sec
+os=600;//rad/sec
+N=log(sqrt((10^(0.1*as)-1)/(10^(0.1*ap)-1)))/log(os/op);
+disp(ceil(N),'Order of the filter, N =');
\ No newline at end of file diff --git a/52/CH5/EX5.10/Example5_10.sce b/52/CH5/EX5.10/Example5_10.sce new file mode 100755 index 000000000..1a274a267 --- /dev/null +++ b/52/CH5/EX5.10/Example5_10.sce @@ -0,0 +1,16 @@ +//Example 5.10
+//To Design a H.P.F. with given specifications
+clear;
+clc ;
+close ;
+ap=3;//db
+as=15;//db
+op=500;//rad/sec
+os=1000;//rad/sec
+N=log(sqrt((10^(0.1*as)-1)/(10^(0.1*ap)-1)))/log(os/op);
+disp(ceil(N),'Order of the filter, N =');
+s=%s;
+HS=1/((s+1)*(s^2+s+1));//Transfer Function for N=3
+oc=1000//rad/sec
+HS1=horner(HS,oc/s);
+disp(HS1,'Normalized Transfer Function, H(s) =');
\ No newline at end of file diff --git a/52/CH5/EX5.11/Example5_11.sce b/52/CH5/EX5.11/Example5_11.sce new file mode 100755 index 000000000..7dcefcda2 --- /dev/null +++ b/52/CH5/EX5.11/Example5_11.sce @@ -0,0 +1,16 @@ +//Example 5.11
+//To Design the Filter using Impulse Invarient Method
+clear;
+clc ;
+close ;
+s=%s;
+T=1;
+HS=(2)/(s^2+3*s+2);
+elts=pfss(HS);
+disp(elts,'Factorized HS = ');
+//The poles comes out to be at -2 and -1
+p1=-2;
+p2=-1;
+z=%z;
+HZ=(2/(1-%e^(p2*T)*z^(-1)))-(2/(1-%e^(p1*T)*z^(-1)))
+disp(HZ,'HZ = ');
\ No newline at end of file diff --git a/52/CH5/EX5.12/Example5_12.sce b/52/CH5/EX5.12/Example5_12.sce new file mode 100755 index 000000000..f99a8d680 --- /dev/null +++ b/52/CH5/EX5.12/Example5_12.sce @@ -0,0 +1,13 @@ +//Example 5.12
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//To Design the Filter using Impulse Invarient Method
+clear;
+clc ;
+close ;
+s=%s;
+HS=1/(s^2+sqrt(2)*s+1);
+pp=ilaplace(HS);
+syms n z;
+t=1;
+X= symsum (pp*(z^(-n)),n ,0, %inf );
+disp(X,'Factorized HS = ');
\ No newline at end of file diff --git a/52/CH5/EX5.13/Example5_13.sce b/52/CH5/EX5.13/Example5_13.sce new file mode 100755 index 000000000..4b7da1a21 --- /dev/null +++ b/52/CH5/EX5.13/Example5_13.sce @@ -0,0 +1,13 @@ +//Example 5.13
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//To Design the 3rd Order Butterworth Filter using Impulse Invarient Method
+clear;
+clc ;
+close ;
+s=%s;
+HS=1/((s+1)*(s^2+s+1));
+pp=ilaplace(HS);//Inverse Laplace
+syms n z;
+t=1;
+X= symsum (pp*(z^(-n)),n ,0, %inf );//Z Transform
+disp(X,'H(z)= ');
\ No newline at end of file diff --git a/52/CH5/EX5.15/Example5_15.sce b/52/CH5/EX5.15/Example5_15.sce new file mode 100755 index 000000000..2f8186358 --- /dev/null +++ b/52/CH5/EX5.15/Example5_15.sce @@ -0,0 +1,16 @@ +//Example 5.15
+//To Design the Filter using Impulse Invarient Method
+clear;
+clc ;
+close ;
+s=%s;
+T=0.2;
+HS=10/(s^2+7*s+10);
+elts=pfss(HS);
+disp(elts,'Factorized HS = ');
+//The poles comes out to be at -5 and -2
+p1=-5;
+p2=-2;
+z=%z;
+HZ=T*((-3.33/(1-%e^(p1*T)*z^(-1)))+(3.33/(1-%e^(p2*T)*z^(-1))))
+disp(HZ,'HZ = ');
\ No newline at end of file diff --git a/52/CH5/EX5.16/Example5_16.sce b/52/CH5/EX5.16/Example5_16.sce new file mode 100755 index 000000000..44f853369 --- /dev/null +++ b/52/CH5/EX5.16/Example5_16.sce @@ -0,0 +1,11 @@ +//Example 5.16
+//To Find out Bilinear Transformation of HS=2/((s+1)*(s+2))
+clear;
+clc ;
+close ;
+s=%s;
+z=%z;
+HS=2/((s+1)*(s+2));
+T=1;
+HZ=horner(HS,(2/T)*(z-1)/(z+1));
+disp(HZ,'H(z) =');
\ No newline at end of file diff --git a/52/CH5/EX5.17/Example5_17.sce b/52/CH5/EX5.17/Example5_17.sce new file mode 100755 index 000000000..b03a6a01e --- /dev/null +++ b/52/CH5/EX5.17/Example5_17.sce @@ -0,0 +1,25 @@ +//Example 5.17
+//To Design an H.P.F. monotonic in passband using Bilinear Transform
+clear;
+clc ;
+close ;
+ap=3;//db
+as=10;//db
+fp=1000;//Hz
+fs=350;//Hz
+f=5000;
+T=1/f;
+wp=2*%pi*fp;
+ws=2*%pi*fs;
+op=2/T*tan(wp*T/2);
+os=2/T*tan(ws*T/2);
+N=log(sqrt((10^(0.1*as)-1)/(10^(0.1*ap)-1)))/log(op/os);
+disp(ceil(N),'Order of the filter, N =');
+s=%s;
+HS=1/(s+1)//Transfer Function for N=1
+oc=op//rad/sec
+HS1=horner(HS,oc/s);
+disp(HS1,'Normalized Transfer Function, H(s) =');
+z=%z;
+HZ=horner(HS,(2/T)*(z-1)/(z+1));
+disp(HZ,'H(z) =');
\ No newline at end of file diff --git a/52/CH5/EX5.18/Example5_18.sce b/52/CH5/EX5.18/Example5_18.sce new file mode 100755 index 000000000..12d40dfc1 --- /dev/null +++ b/52/CH5/EX5.18/Example5_18.sce @@ -0,0 +1,11 @@ +//Example 5.18
+//To Find out Bilinear Transformation of H(s)=(s^2+4.525)/(s^2+0.692*s+0.504)
+clear;
+clc ;
+close ;
+s=%s;
+z=%z;
+HS=(s^2+4.525)/(s^2+0.692*s+0.504);
+T=1;
+HZ=horner(HS,(2/T)*(z-1)/(z+1));
+disp(HZ,'H(z) =');
\ No newline at end of file diff --git a/52/CH5/EX5.19/Example5_19.sce b/52/CH5/EX5.19/Example5_19.sce new file mode 100755 index 000000000..adf0d7ff4 --- /dev/null +++ b/52/CH5/EX5.19/Example5_19.sce @@ -0,0 +1,17 @@ +//Example 5.19
+//To Convert a single Pole LPF into BPF
+clear;
+clc ;
+close ;
+s=%s;
+z=%z;
+HZ=(0.5*(1+z^(-1)))/(1-0.302*z^(-2));
+T=1;
+wu=3*%pi/4;
+wl=%pi/4;
+wp=%pi/6;
+k=tan(wp/2)/tan((wu-wl)/2);
+a=cos((wu+wl)/2)/cos((wu-wl)/2);
+transf=-((((k-1)/(k+1))*(z^(-2)))-((2*a*k/(k+1))*(z^(-1)))+1)/(z^(-2)-(2*a*k/(1+k)*z^(-1))+((k-1)/(k+1)));
+HZ1=horner(HZ,transf);
+disp(HZ1,'H(z) of B.P.F =');
\ No newline at end of file diff --git a/52/CH5/EX5.2/Example5_2.sce b/52/CH5/EX5.2/Example5_2.sce new file mode 100755 index 000000000..a6c9a5daf --- /dev/null +++ b/52/CH5/EX5.2/Example5_2.sce @@ -0,0 +1,13 @@ +//Example 5.2
+//To Find out the order of a Low Pass Butterworth Filter
+clear;
+clc ;
+close ;
+ap=3;//db
+as=40;//db
+fp=500;//Hz
+fs=1000;//Hz
+op=2*%pi*fp;
+os=2*%pi*fs;
+N=log(sqrt((10^(0.1*as)-1)/(10^(0.1*ap)-1)))/log(os/op);
+disp(ceil(N),'Order of the filter, N =');
\ No newline at end of file diff --git a/52/CH5/EX5.29/Example5_29.sce b/52/CH5/EX5.29/Example5_29.sce new file mode 100755 index 000000000..2c88c0b92 --- /dev/null +++ b/52/CH5/EX5.29/Example5_29.sce @@ -0,0 +1,35 @@ +//Example 5.29
+//Program to convert given IIR pole-zero Filter into Lattice Ladder Structure.
+clear;
+clc ;
+close ;
+U=1; //Zero Adjust
+a(3+U,0+U)=1;
+a(3+U,1+U)=13/24;
+a(3+U,2+U)=5/8;
+a(3+U,3+U)=1/3;
+a(2+U,0+U)=1; //a(m,0)=1
+a(2+U,3+U)=1/3;
+m=3,k=1;
+a(m-1+U,k+U)=(a(m+U,k+U)-a(m+U,m+U)*a(m+U,m-k+U))/(1-a(m+U,m+U)*a(m+U,m+U));
+m=3,k=2;
+a(m-1+U,k+U)=(a(m+U,k+U)-a(m+U,m+U)*a(m+U,m-k+U))/(1-a(m+U,m+U)*a(m+U,m+U));
+m=2,k=1;
+a(m-1+U,k+U)=(a(m+U,k+U)-a(m+U,m+U)*a(m+U,m-k+U))/(1-a(m+U,m+U)*a(m+U,m+U));
+disp('LATTICE COEFFICIENTS');
+disp(a(1+U,1+U),'k1');
+disp(a(2+U,2+U),'k2');
+disp(a(3+U,3+U),'k3');
+b0=1;
+b1=2;
+b2=2;
+b3=1;
+c3=b3;
+c2=b2-c3*a(3+U,1+U);
+c1=b1-(c2*a(2+U,1+U)+c3*a(3+U,2+U));
+c0=b0-(c1*a(1+U,1+U)+c2*a(2+U,2+U)+c3*a(3+U,3+U));
+disp('LADDER COEFFICIENTS');
+disp(c0,'c0 =');
+disp(c1,'c1 =');
+disp(c2,'c2 =');
+disp(c3,'c3 =');
\ No newline at end of file diff --git a/52/CH5/EX5.4/Example5_4.sce b/52/CH5/EX5.4/Example5_4.sce new file mode 100755 index 000000000..5beb9a6c0 --- /dev/null +++ b/52/CH5/EX5.4/Example5_4.sce @@ -0,0 +1,16 @@ +//Example 5.4
+//To Design an Analog Butterworth Filter
+clear;
+clc ;
+close ;
+ap=2;//db
+as=10;//db
+op=20;//rad/sec
+os=30;//rad/sec
+N=log(sqrt((10^(0.1*as)-1)/(10^(0.1*ap)-1)))/log(os/op);
+disp(ceil(N),'Order of the filter, N =');
+s=%s;
+HS=1/((s^2+0.76537*s+1)*(s^2+1.8477*s+1));//Transfer Function for N=4
+oc=op/(10^(0.1*ap)-1)^(1/(2*ceil(N)));
+HS1=horner(HS,s/oc);
+disp(HS1,'Normalized Transfer Function, H(s) =');
\ No newline at end of file diff --git a/52/CH5/EX5.5/Example5_5.sce b/52/CH5/EX5.5/Example5_5.sce new file mode 100755 index 000000000..211048b0b --- /dev/null +++ b/52/CH5/EX5.5/Example5_5.sce @@ -0,0 +1,18 @@ +//Example 5.5
+//To Design an Analog Butterworth Filter
+clear;
+clc ;
+close ;
+op=0.2*%pi;
+os=0.4*%pi;
+e1=0.9;
+l1=0.2;
+epsilon=sqrt(1/(e1^2)-1);
+lambda=sqrt(1/(l1^2)-1);
+N=log(lambda/epsilon)/log(os/op);
+disp(ceil(N),'Order of the filter, N =');
+s=%s;
+HS=1/((s^2+0.76537*s+1)*(s^2+1.8477*s+1));//Transfer Function for N=4
+oc=op/epsilon^(1/ceil(N));
+HS1=horner(HS,s/oc);
+disp(HS1,'Normalized Transfer Function, H(s) =');
\ No newline at end of file diff --git a/52/CH5/EX5.6/Example5_6.sce b/52/CH5/EX5.6/Example5_6.sce new file mode 100755 index 000000000..0c8717271 --- /dev/null +++ b/52/CH5/EX5.6/Example5_6.sce @@ -0,0 +1,13 @@ +//Example 5.6
+//To Find out the order of the Filter using Chebyshev Approximation
+clear;
+clc ;
+close ;
+ap=3;//db
+as=16;//db
+fp=1000;//Hz
+fs=2000;//Hz
+op=2*%pi*fp;
+os=2*%pi*fs;
+N=acosh(sqrt((10^(0.1*as)-1)/(10^(0.1*ap)-1)))/acosh(os/op);
+disp(ceil(N),'Order of the filter, N =');
\ No newline at end of file diff --git a/52/CH5/EX5.7/Example5_7.sce b/52/CH5/EX5.7/Example5_7.sce new file mode 100755 index 000000000..d8668e5eb --- /dev/null +++ b/52/CH5/EX5.7/Example5_7.sce @@ -0,0 +1,15 @@ +//Example 5.7
+//To Design an analog Chebyshev Filter with Given Specifications
+clear;
+clc ;
+close ;
+os=2;
+op=1;
+ap=3;//db
+as=16;//db
+e1=1/sqrt(2);
+l1=0.1;
+epsilon=sqrt(1/(e1^2)-1);
+lambda=sqrt(1/(l1^2)-1);
+N=acosh(lambda/epsilon)/acosh(os/op);
+disp(ceil(N),'Order of the filter, N =');
\ No newline at end of file diff --git a/52/CH5/EX5.8/Example5_8.sce b/52/CH5/EX5.8/Example5_8.sce new file mode 100755 index 000000000..f4480f66b --- /dev/null +++ b/52/CH5/EX5.8/Example5_8.sce @@ -0,0 +1,11 @@ +//Example 5.8
+//To Find out the order of the poles of the Type 1 Lowpass Chebyshev Filter
+clear;
+clc ;
+close ;
+ap=1;//dB
+as=40;//dB
+op=1000*%pi;
+os=2000*%pi;
+N=acosh(sqrt((10^(0.1*as)-1)/(10^(0.1*ap)-1)))/acosh(os/op);
+disp(ceil(N),'Order of the filter, N =');
\ No newline at end of file diff --git a/52/CH5/EX5.9/Example5_9.sce b/52/CH5/EX5.9/Example5_9.sce new file mode 100755 index 000000000..76cd1d0f8 --- /dev/null +++ b/52/CH5/EX5.9/Example5_9.sce @@ -0,0 +1,11 @@ +//Example 5.9
+//To Design a Chebyshev Filter with Given Specifications
+clear;
+clc ;
+close ;
+ap=2.5;//db
+as=30;//db
+op=20;//rad/sec
+os=50;//rad/sec
+N=acosh(sqrt((10^(0.1*as)-1)/(10^(0.1*ap)-1)))/acosh(os/op);
+disp(ceil(N),'Order of the filter, N =');
\ No newline at end of file diff --git a/52/CH6/EX6.1/Example6_1.sce b/52/CH6/EX6.1/Example6_1.sce new file mode 100755 index 000000000..fdc678d91 --- /dev/null +++ b/52/CH6/EX6.1/Example6_1.sce @@ -0,0 +1,14 @@ +//Example 6.1
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Program to Calculate Group Delay and Phase Delay
+//y(n)=0.25x(n)+x(n-1)+0.25x(n-2)
+clear;
+clc ;
+close ;
+//w=poly(0,"w");
+syms w;
+theeta=-w;
+gd= -diff (theeta,w); //Group Delay
+pd=-theeta/w; //Phase Delay
+disp(gd,'GROUP DELAY =');
+disp(pd,'PHASE DELAY =');
\ No newline at end of file diff --git a/52/CH6/EX6.10/Example6_10.jpg b/52/CH6/EX6.10/Example6_10.jpg Binary files differnew file mode 100755 index 000000000..1d73c59cf --- /dev/null +++ b/52/CH6/EX6.10/Example6_10.jpg diff --git a/52/CH6/EX6.10/Example6_10.sce b/52/CH6/EX6.10/Example6_10.sce new file mode 100755 index 000000000..ca64e9a36 --- /dev/null +++ b/52/CH6/EX6.10/Example6_10.sce @@ -0,0 +1,28 @@ +//Example 6.10
+//Program to Plot Magnitude Responce of given L.P.F. with specifications:
+//N=7,w=pi/4
+//Using Hanning Window
+clear;
+clc ;
+close ;
+N=7;
+alpha=3;
+U=1;
+h_hann=window('hn',N);
+for n=0+U:1:6+U
+if n==4
+hd(n)=0.25;
+else
+hd(n)=(sin(%pi*(n-U-alpha)/4))/(%pi*(n-U-alpha));
+end
+h(n)=hd(n)*h_hann(n);
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of given LPF with N=7');
+xgrid (2);
\ No newline at end of file diff --git a/52/CH6/EX6.11/Example6_11.jpg b/52/CH6/EX6.11/Example6_11.jpg Binary files differnew file mode 100755 index 000000000..298e5757f --- /dev/null +++ b/52/CH6/EX6.11/Example6_11.jpg diff --git a/52/CH6/EX6.11/Example6_11.sce b/52/CH6/EX6.11/Example6_11.sce new file mode 100755 index 000000000..7975abc53 --- /dev/null +++ b/52/CH6/EX6.11/Example6_11.sce @@ -0,0 +1,44 @@ +//Example 6.11
+//Program to Plot Magnitude Responce of given L.P.F. with specifications:
+//wp=20rad/sec, ws=30rad/sec, wsf=100rad/sec
+//as=44.0dB, ap=0.1dB
+//Using Kaiser Window
+clear;
+clc ;
+close ;
+wsf=100//rad/sec
+ws=30;//rad/sec
+wp=20;//rad/sec
+as=44.0//dB
+ap=0.1//dB
+B=ws-wp;
+wc=0.5*(ws+wp);
+wc1=wc*2*%pi/wsf;
+delta1=10^(-0.05*as);
+delta2=(10^(0.05*as)-1)/(10^(0.05*as)+1);
+delta=min(delta1,delta2);
+alphas=-20*log10(delta);
+alpha=0.5842*(alphas-21)^0.4+0.07886*(alphas-21)
+D=(alphas-7.95)/14.36;
+N1=wsf*D/B+1;
+N=ceil(N1);
+U=ceil(N/2);
+win_l=window('kr',N,alpha);
+for n=-floor(N/2)+U:1:floor(N/2)+U
+if n==ceil(N/2);
+hd(n)=0.5;
+else
+hd(n)=(sin(%pi*(n-U)/2))/(%pi*(n-U));
+end
+h(n)=hd(n)*win_l(n);
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of given LPF using Kaiser Window');
+xgrid (2);
+disp(h,"Filter Coefficients,h(n)=");
\ No newline at end of file diff --git a/52/CH6/EX6.12/Example6_12.jpg b/52/CH6/EX6.12/Example6_12.jpg Binary files differnew file mode 100755 index 000000000..e43d4102a --- /dev/null +++ b/52/CH6/EX6.12/Example6_12.jpg diff --git a/52/CH6/EX6.12/Example6_12.sce b/52/CH6/EX6.12/Example6_12.sce new file mode 100755 index 000000000..0d863c555 --- /dev/null +++ b/52/CH6/EX6.12/Example6_12.sce @@ -0,0 +1,50 @@ +//Example 6.12
+//Program to Plot Magnitude Responce of given B.P.F. with specifications:
+//wp1=40pi rad/sec, wp2=60pi rad/sec
+//ws1=20pi rad/sec, ws2=80pi rad/sec
+//as=30dB, ap=0.5dB
+//F=100 Hz
+//Using Kaiser Window
+clear;
+clc ;
+close ;
+wsf=200*%pi;//rad/sec
+ws1=20*%pi;//rad/sec
+ws2=80*%pi;//rad/sec
+wp1=40*%pi;//rad/sec
+wp2=60*%pi;//rad/sec
+as=30//dB
+ap=0.5//dB
+B=min(wp1-ws1,ws2-wp2);
+wc1=wp1-B/2;
+wc2=wp2+B/2;
+wc1=wc1*2*%pi/wsf;
+wc2=wc2*2*%pi/wsf;
+delta1=10^(-0.05*as);
+delta2=(10^(0.05*as)-1)/(10^(0.05*as)+1);
+delta=min(delta1,delta2);
+alphas=-20*log10(delta);
+alpha=0.5842*(alphas-21)^0.4+0.07886*(alphas-21)
+D=(alphas-7.95)/14.36;
+N1=wsf*D/B+1;
+N=ceil(N1);
+U=ceil(N/2);
+win_l=window('kr',N,alpha);
+for n=-floor(N/2)+U:1:floor(N/2)+U
+if n==ceil(N/2);
+hd(n)=0.4;
+else
+hd(n)=(sin(0.7*%pi*(n-U))-sin(0.3*%pi*(n-U)))/(%pi*(n-U));
+end
+h(n)=hd(n)*win_l(n);
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of given LPF using Kaiser Window');
+xgrid (2);
+disp(h,"Filter Coefficients,h(n)=");
\ No newline at end of file diff --git a/52/CH6/EX6.13.a/Example6_13_a.jpg b/52/CH6/EX6.13.a/Example6_13_a.jpg Binary files differnew file mode 100755 index 000000000..8c9b21772 --- /dev/null +++ b/52/CH6/EX6.13.a/Example6_13_a.jpg diff --git a/52/CH6/EX6.13.a/Example6_13_a.sce b/52/CH6/EX6.13.a/Example6_13_a.sce new file mode 100755 index 000000000..da31c195d --- /dev/null +++ b/52/CH6/EX6.13.a/Example6_13_a.sce @@ -0,0 +1,24 @@ +//Example 6.13a
+//Program to Plot Magnitude Responce of ideal differentiator with specifications:
+//N=8,w=pi
+//using Rectangular window
+clear;
+clc ;
+close ;
+N=8;
+alpha=7/2;
+U=1;
+h_Rect=window('re',N);
+for n=0+U:1:7+U
+hd(n)=-(sin(%pi*(n-U-alpha)))/(%pi*(n-U-alpha)*(n-U-alpha));
+h(n)=hd(n)*h_Rect(n);
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of given ideal differentiator using Rectangular Window, N=8');
+xgrid (2)
\ No newline at end of file diff --git a/52/CH6/EX6.13.b/Example6_13_b.jpg b/52/CH6/EX6.13.b/Example6_13_b.jpg Binary files differnew file mode 100755 index 000000000..497af2428 --- /dev/null +++ b/52/CH6/EX6.13.b/Example6_13_b.jpg diff --git a/52/CH6/EX6.13.b/Example6_13_b.sce b/52/CH6/EX6.13.b/Example6_13_b.sce new file mode 100755 index 000000000..02dbf3bd4 --- /dev/null +++ b/52/CH6/EX6.13.b/Example6_13_b.sce @@ -0,0 +1,24 @@ +//Example 6.13b
+//Program to Plot Magnitude Responce of ideal differentiator with specifications:
+//N=8,w=pi
+//using Hamming window
+clear;
+clc ;
+close ;
+N=8;
+alpha=7/2;
+U=1; //Zero Adjust
+h_hamm=window('hm',N);
+for n=0+U:1:7+U
+hd(n)=-(sin(%pi*(n-U-alpha)))/(%pi*(n-U-alpha)*(n-U-alpha));
+h(n)=hd(n)*h_hamm(n);
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of given ideal differentiator using Hamming Window, N=8');
+xgrid (2)
\ No newline at end of file diff --git a/52/CH6/EX6.14.a/Example6_14_a.jpg b/52/CH6/EX6.14.a/Example6_14_a.jpg Binary files differnew file mode 100755 index 000000000..b046a15e0 --- /dev/null +++ b/52/CH6/EX6.14.a/Example6_14_a.jpg diff --git a/52/CH6/EX6.14.a/Example6_14_a.sce b/52/CH6/EX6.14.a/Example6_14_a.sce new file mode 100755 index 000000000..1484a0589 --- /dev/null +++ b/52/CH6/EX6.14.a/Example6_14_a.sce @@ -0,0 +1,26 @@ +//Example 6.14a
+//Program to Plot Magnitude Responce of ideal Hilbert Transformer
+//using Rectangular Window
+//N=11
+clear;
+clc ;
+close ;
+N=11;
+U=6;
+h_Rect=window('re',N);
+for n=-5+U:1:5+U
+if n==6
+hd(n)=0;
+else
+hd(n)=(1-cos(%pi*(n-U)))/(%pi*(n-U));
+end
+h(n)=hd(n)*h_Rect(n);
+end
+[hzm ,fr]= frmag (h,256) ;
+figure
+plot (2*fr ,-hzm);
+a = gca ();
+xlabel ('Frequency w*pi');
+ylabel ('H(exp(j*w))');
+title ('Frequency Response of Hilbert Transformer with N=11 using Rectangular Window');
+xgrid (2);
\ No newline at end of file diff --git a/52/CH6/EX6.14.b/Example6_14_b.jpg b/52/CH6/EX6.14.b/Example6_14_b.jpg Binary files differnew file mode 100755 index 000000000..84a10cb11 --- /dev/null +++ b/52/CH6/EX6.14.b/Example6_14_b.jpg diff --git a/52/CH6/EX6.14.b/Example6_14_b.sce b/52/CH6/EX6.14.b/Example6_14_b.sce new file mode 100755 index 000000000..f256e07fe --- /dev/null +++ b/52/CH6/EX6.14.b/Example6_14_b.sce @@ -0,0 +1,26 @@ +//Example 6.14b
+//Program to Plot Magnitude Responce of ideal Hilbert Transformer
+//using Blackman Window
+//N=11
+clear;
+clc ;
+close ;
+N=11;
+U=6;
+for n=-5+U:1:5+U
+h_balckmann(n) = 0.42+0.5*cos(2*%pi*(n-U)/(N-1))+0.08*cos(4*%pi*(n-U)/(N-1));
+if n==6
+hd(n)=0;
+else
+hd(n)=(1-cos(%pi*(n-U)))/(%pi*(n-U));
+end
+h(n)=hd(n)*h_balckmann(n);
+end
+[hzm ,fr]= frmag (h,256) ;
+figure
+plot (2*fr ,-hzm);
+a = gca ();
+xlabel ('Frequency w*pi');
+ylabel ('H(exp(j*w))');
+title ('Frequency Response of Hilbert Transformer with N=11 using Blackman Window');
+xgrid (2);
\ No newline at end of file diff --git a/52/CH6/EX6.15/Example6_15.sce b/52/CH6/EX6.15/Example6_15.sce new file mode 100755 index 000000000..0b6ff853d --- /dev/null +++ b/52/CH6/EX6.15/Example6_15.sce @@ -0,0 +1,12 @@ +//Example 6.15
+//Program to determine filter coefficients obtained by sampling:
+//N=7,w=pi/2
+clear;
+clc ;
+close ;
+N=7;
+U=1; //Zero Adjust
+for n=0+U:1:N-1+U
+h(n)=(1+2*cos(2*%pi*(n-U-3)/7))/N
+end
+disp(h,"Filter Coefficients,h(n)=")
\ No newline at end of file diff --git a/52/CH6/EX6.16/Example6_16.sce b/52/CH6/EX6.16/Example6_16.sce new file mode 100755 index 000000000..7a47fc3b9 --- /dev/null +++ b/52/CH6/EX6.16/Example6_16.sce @@ -0,0 +1,12 @@ +//Example 6.16
+//Program to determine filter coefficients obtained by sampling:
+//N=15
+clear;
+clc ;
+close ;
+N=15;
+U=1; //Zero Adjust
+for n=0:1:N-1
+h(n+U)=(1+2*cos(2*%pi*(7-n)/N)+2*cos(4*%pi*(7-n)/N)+2*cos(6*%pi*(7-n)/N))/N;
+end
+disp(h,"Filter Coefficients,h(n)=");
diff --git a/52/CH6/EX6.17/Example6_17.sce b/52/CH6/EX6.17/Example6_17.sce new file mode 100755 index 000000000..3d22d5613 --- /dev/null +++ b/52/CH6/EX6.17/Example6_17.sce @@ -0,0 +1,12 @@ +//Example 6.17
+//Program to design bandpass filter with following specifications:
+//N=7, fc1=1000Hz, fc2=3000Hz, F=8000Hz
+clear;
+clc ;
+close ;
+N=7;
+U=1; //Zero Adjust
+for n=0:1:N-1
+h(n+U)=2*(cos(2*%pi*(3-n)/N)+cos(4*%pi*(3-n)/N))/N;
+end
+disp(h,"Filter Coefficients,h(n)=");
diff --git a/52/CH6/EX6.18.a/Example6_18_a.jpg b/52/CH6/EX6.18.a/Example6_18_a.jpg Binary files differnew file mode 100755 index 000000000..3aa5b7330 --- /dev/null +++ b/52/CH6/EX6.18.a/Example6_18_a.jpg diff --git a/52/CH6/EX6.18.a/Example6_18_a.sce b/52/CH6/EX6.18.a/Example6_18_a.sce new file mode 100755 index 000000000..182594604 --- /dev/null +++ b/52/CH6/EX6.18.a/Example6_18_a.sce @@ -0,0 +1,20 @@ +//Example 6.18a
+//Program to design L.P.F. filter with following specifications:
+//N=15, wc=pi/4
+clear;
+clc ;
+close ;
+N=15;
+U=1;
+for n=0+U:1:N-1+U
+h(n)=(1+cos(2*%pi*(7-n)/N))/N;
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure;
+plot (2*fr , hzm_dB );
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of FIR LPF with N=15');
+xgrid (2)
\ No newline at end of file diff --git a/52/CH6/EX6.18.b/Example6_18_b.jpg b/52/CH6/EX6.18.b/Example6_18_b.jpg Binary files differnew file mode 100755 index 000000000..d6ab0a55d --- /dev/null +++ b/52/CH6/EX6.18.b/Example6_18_b.jpg diff --git a/52/CH6/EX6.18.b/Example6_18_b.sce b/52/CH6/EX6.18.b/Example6_18_b.sce new file mode 100755 index 000000000..663ee5aa2 --- /dev/null +++ b/52/CH6/EX6.18.b/Example6_18_b.sce @@ -0,0 +1,20 @@ +//Example 6.18b
+//Program to design L.P.F. filter with following specifications:
+//N=15, wc=pi/4
+clear;
+clc ;
+close ;
+N=15;
+U=1;
+for n=0+U:1:N-1+U
+h(n)=(1+cos(2*%pi*(7-n)/N)+cos(4*%pi*(7-n)/N))/N;
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure;
+plot (2*fr , hzm_dB );
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of FIR LPF with N=11');
+xgrid (2)
\ No newline at end of file diff --git a/52/CH6/EX6.19/Example6_19.jpg b/52/CH6/EX6.19/Example6_19.jpg Binary files differnew file mode 100755 index 000000000..8006ff1bf --- /dev/null +++ b/52/CH6/EX6.19/Example6_19.jpg diff --git a/52/CH6/EX6.19/Example6_19.sce b/52/CH6/EX6.19/Example6_19.sce new file mode 100755 index 000000000..9b8ce9e85 --- /dev/null +++ b/52/CH6/EX6.19/Example6_19.sce @@ -0,0 +1,25 @@ +//Example 6.19
+//Program to Plot Magnitude Responce of given L.P.F. with specifications:
+//N=13,w=pi/6
+clear;
+clc ;
+close ;
+alpha=6;
+U=1;
+for n=0+U:1:12+U
+if n==7
+hd(n)=0.167;
+else
+hd(n)=(sin(%pi*(n-U-alpha)/6))/(%pi*(n-U-alpha));
+end
+end
+[hzm ,fr ]= frmag (hd ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of given LPF with N=13');
+xgrid (2)
+disp(hd,"Filter Coefficients,h(n)=");
\ No newline at end of file diff --git a/52/CH6/EX6.20/Example6_20.jpg b/52/CH6/EX6.20/Example6_20.jpg Binary files differnew file mode 100755 index 000000000..e1aac9a59 --- /dev/null +++ b/52/CH6/EX6.20/Example6_20.jpg diff --git a/52/CH6/EX6.20/Example6_20.sce b/52/CH6/EX6.20/Example6_20.sce new file mode 100755 index 000000000..a981d8666 --- /dev/null +++ b/52/CH6/EX6.20/Example6_20.sce @@ -0,0 +1,30 @@ +//Example 6.20
+//Program to Plot Magnitude Responce of given L.P.F. with specifications:
+//N=13,w=pi/6
+//Using Hamming Window
+clear;
+clc ;
+close ;
+N=13;
+alpha=6;
+U=1;
+h_hamm=window('hm',N);
+for n=0+U:1:12+U
+if n==7
+hd(n)=0.167;
+else
+hd(n)=(sin(%pi*(n-U-alpha)/6))/(%pi*(n-U-alpha));
+end
+h(n)=hd(n)*h_hamm(n);
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of given LPF with N=13');
+xgrid (2)
+disp(h,"Filter Coefficients,h(n)=");
+disp(h,"Filter Coefficients,h(n)=");
\ No newline at end of file diff --git a/52/CH6/EX6.21/Example6_21.jpg b/52/CH6/EX6.21/Example6_21.jpg Binary files differnew file mode 100755 index 000000000..8b895bee2 --- /dev/null +++ b/52/CH6/EX6.21/Example6_21.jpg diff --git a/52/CH6/EX6.21/Example6_21.sce b/52/CH6/EX6.21/Example6_21.sce new file mode 100755 index 000000000..3421ba2e6 --- /dev/null +++ b/52/CH6/EX6.21/Example6_21.sce @@ -0,0 +1,27 @@ +//Example 6.21
+//Program to Plot Magnitude Responce of given L.P.F. with specifications:
+//N=7,fc=1000Hz,F=5000Hz
+clear;
+clc ;
+close ;
+N=7;
+U=4;
+h_Rect=window('re',N);
+for n=-3+U:1:3+U
+if n==4
+hd(n)=0.4;
+else
+hd(n)=(sin(2*%pi*(n-U)/5))/(%pi*(n-U));
+end
+h(n)=hd(n)*h_Rect(n);
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of FIR LPF with N=7');
+xgrid (2)
+disp(h,"Filter Coefficients,h(n)=");
\ No newline at end of file diff --git a/52/CH6/EX6.28/Example6_28.sce b/52/CH6/EX6.28/Example6_28.sce new file mode 100755 index 000000000..8359154d7 --- /dev/null +++ b/52/CH6/EX6.28/Example6_28.sce @@ -0,0 +1,24 @@ +//Example 6.28
+//Program to calculate FIR Filter coefficients for the direct form structure
+//k1=1/2 , k2=1/3 , k3=1/4
+clear;
+clc ;
+close ;
+U=1;
+k1=1/2;
+k2=1/3;
+k3=1/4;
+a(3+U,0+U)=1;
+a(1+U,1+U)=k1;
+a(2+U,2+U)=k2;
+a(3+U,3+U)=k3;
+m=2,k=1;
+a(m+U,k+U)=a(m-1+U,k+U)+a(m+U,m+U)*a(m-1+U,m-k+U);
+m=3,k=1;
+a(m+U,k+U)=a(m-1+U,k+U)+a(m+U,m+U)*a(m-1+U,m-k+U);
+m=3,k=2;
+a(m+U,k+U)=a(m-1+U,k+U)+a(m+U,m+U)*a(m-1+U,m-k+U);
+disp(a(3+U,0+U),'a(3,0)');
+disp(a(3+U,1+U),'a(3,1)');
+disp(a(3+U,2+U),'a(3,2)');
+disp(a(3+U,3+U),'a(3,3)');
\ No newline at end of file diff --git a/52/CH6/EX6.29/Example6_29.sce b/52/CH6/EX6.29/Example6_29.sce new file mode 100755 index 000000000..79543fbdf --- /dev/null +++ b/52/CH6/EX6.29/Example6_29.sce @@ -0,0 +1,21 @@ +//Example 6.29
+//Program to calculate given FIR Filter's Lattice form coefficients.
+clear;
+clc ;
+close ;
+U=1; //Zero Adjust
+a(3+U,0+U)=1;
+a(3+U,1+U)=2/5;
+a(3+U,2+U)=3/4;
+a(3+U,3+U)=1/3;
+a(2+U,0+U)=1; //a(m,0)=1
+a(2+U,3+U)=1/3;
+m=3,k=1;
+a(m-1+U,k+U)=(a(m+U,k+U)-a(m+U,m+U)*a(m+U,m-k+U))/(1-a(m+U,m+U)*a(m+U,m+U));
+m=3,k=2;
+a(m-1+U,k+U)=(a(m+U,k+U)-a(m+U,m+U)*a(m+U,m-k+U))/(1-a(m+U,m+U)*a(m+U,m+U));
+m=2,k=1;
+a(m-1+U,k+U)=(a(m+U,k+U)-a(m+U,m+U)*a(m+U,m-k+U))/(1-a(m+U,m+U)*a(m+U,m+U));
+disp(a(1+U,1+U),'k1');
+disp(a(2+U,2+U),'k2');
+disp(a(3+U,3+U),'k3');
\ No newline at end of file diff --git a/52/CH6/EX6.5/Example6_5.jpg b/52/CH6/EX6.5/Example6_5.jpg Binary files differnew file mode 100755 index 000000000..be6c5c8bc --- /dev/null +++ b/52/CH6/EX6.5/Example6_5.jpg diff --git a/52/CH6/EX6.5/Example6_5.sce b/52/CH6/EX6.5/Example6_5.sce new file mode 100755 index 000000000..a92ebca2f --- /dev/null +++ b/52/CH6/EX6.5/Example6_5.sce @@ -0,0 +1,24 @@ +//Example 6.5
+//Program to Plot Magnitude Responce of ideal L.P.F. with wc=0.5*pi
+//N=11
+clear;
+clc ;
+close ;
+N=11;
+U=6;
+for n=-5+U:1:5+U
+if n==6
+hd(n)=0.5;
+else
+hd(n)=(sin(%pi*(n-U)/2))/(%pi*(n-U));
+end
+end
+[hzm ,fr ]= frmag (hd ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure;
+plot (2*fr , hzm_dB );
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of FIR LPF with N=11');
+xgrid (2);
\ No newline at end of file diff --git a/52/CH6/EX6.6/Example6_6.jpg b/52/CH6/EX6.6/Example6_6.jpg Binary files differnew file mode 100755 index 000000000..a3577499a --- /dev/null +++ b/52/CH6/EX6.6/Example6_6.jpg diff --git a/52/CH6/EX6.6/Example6_6.sce b/52/CH6/EX6.6/Example6_6.sce new file mode 100755 index 000000000..b0a5d9015 --- /dev/null +++ b/52/CH6/EX6.6/Example6_6.sce @@ -0,0 +1,24 @@ +//Example 6.6
+//Program to Plot Magnitude Responce of ideal H.P.F. with wc=0.25*pi
+//N=11
+clear;
+clc ;
+close ;
+N=11;
+U=6;
+for n=-5+U:1:5+U
+if n==6
+hd(n)=0.5;
+else
+hd(n)=(sin(%pi*(n-U))-sin(%pi*(n-U)/4))/(%pi*(n-U));
+end
+end
+[hzm ,fr ]= frmag (hd ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of FIR HPF with N=11');
+xgrid (2);
\ No newline at end of file diff --git a/52/CH6/EX6.7/Example6_7.jpg b/52/CH6/EX6.7/Example6_7.jpg Binary files differnew file mode 100755 index 000000000..fdca17813 --- /dev/null +++ b/52/CH6/EX6.7/Example6_7.jpg diff --git a/52/CH6/EX6.7/Example6_7.sce b/52/CH6/EX6.7/Example6_7.sce new file mode 100755 index 000000000..b855367be --- /dev/null +++ b/52/CH6/EX6.7/Example6_7.sce @@ -0,0 +1,25 @@ +//Example 6.7
+//Program to Plot Magnitude Responce of ideal B.P.F. with
+//wc1=0.25*pi and wc2=0.75*pi
+//N=11
+clear;
+clc ;
+close ;
+N=11;
+U=6;
+for n=-5+U:1:5+U
+if n==6
+hd(n)=0.5;
+else
+hd(n)=(sin(%pi*3*(n-U)/4)-sin(%pi*(n-U)/4))/(%pi*(n-U));
+end
+end
+[hzm ,fr ]= frmag (hd ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of FIR BPF with N=11');
+xgrid (2);
\ No newline at end of file diff --git a/52/CH6/EX6.8/Example6_8.jpg b/52/CH6/EX6.8/Example6_8.jpg Binary files differnew file mode 100755 index 000000000..cd5fd1e4c --- /dev/null +++ b/52/CH6/EX6.8/Example6_8.jpg diff --git a/52/CH6/EX6.8/Example6_8.sce b/52/CH6/EX6.8/Example6_8.sce new file mode 100755 index 000000000..13e506920 --- /dev/null +++ b/52/CH6/EX6.8/Example6_8.sce @@ -0,0 +1,25 @@ +//Example 6.8
+//Program to Plot Magnitude Responce of ideal B.R.F. with
+//wc1=0.33*pi and wc2=0.67*pi
+//N=11
+clear;
+clc ;
+close ;
+N=11;
+U=6;
+for n=-5+U:1:5+U
+if n==6
+hd(n)=0.5;
+else
+hd(n)=(sin(%pi*(n-U))+sin(%pi*(n-U)/3)-sin(%pi*2*(n-U)/3))/(%pi*(n-U));
+end
+end
+[hzm ,fr ]= frmag (hd ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of FIR BRF with N=11');
+xgrid (2);
\ No newline at end of file diff --git a/52/CH6/EX6.9.a/Example6_9_a.jpg b/52/CH6/EX6.9.a/Example6_9_a.jpg Binary files differnew file mode 100755 index 000000000..973fc8a61 --- /dev/null +++ b/52/CH6/EX6.9.a/Example6_9_a.jpg diff --git a/52/CH6/EX6.9.a/Example6_9_a.sce b/52/CH6/EX6.9.a/Example6_9_a.sce new file mode 100755 index 000000000..00f4342f0 --- /dev/null +++ b/52/CH6/EX6.9.a/Example6_9_a.sce @@ -0,0 +1,28 @@ +//Example 6.9a
+//Program to Plot Magnitude Responce of ideal H.P.F.
+//using Hanning Window
+//wc1=0.25*pi
+//N=11
+clear;
+clc ;
+close ;
+N=11;
+U=6;
+h_hann=window('hn',N);
+for n=-5+U:1:5+U
+if n==6
+hd(n)=0.75;
+else
+hd(n)=(sin(%pi*(n-U))-sin(%pi*(n-U)/4))/(%pi*(n-U));
+end
+h(n)=h_hann(n)*hd(n);
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of FIR HRF with N=11 using Hanning Window');
+xgrid (2);
\ No newline at end of file diff --git a/52/CH6/EX6.9.b/Example6_9_b.jpg b/52/CH6/EX6.9.b/Example6_9_b.jpg Binary files differnew file mode 100755 index 000000000..876c36a43 --- /dev/null +++ b/52/CH6/EX6.9.b/Example6_9_b.jpg diff --git a/52/CH6/EX6.9.b/Example6_9_b.sce b/52/CH6/EX6.9.b/Example6_9_b.sce new file mode 100755 index 000000000..b5eb6f916 --- /dev/null +++ b/52/CH6/EX6.9.b/Example6_9_b.sce @@ -0,0 +1,28 @@ +//Example 6.9b
+//Program to Plot Magnitude Responce of ideal H.P.F.
+//using Hamming Window
+//wc1=0.25*pi
+//N=11
+clear;
+clc ;
+close ;
+N=11;
+U=6;
+h_hamm=window('hm',N);
+for n=-5+U:1:5+U
+if n==6
+hd(n)=0.75;
+else
+hd(n)=(sin(%pi*(n-U))-sin(%pi*(n-U)/4))/(%pi*(n-U));
+end
+h(n)=h_hamm(n)*hd(n);
+end
+[hzm ,fr ]= frmag (h ,256) ;
+hzm_dB = 20* log10 (hzm)./ max ( hzm );
+figure
+plot (2*fr , hzm_dB )
+a= gca ();
+xlabel ('Frequency w*pi');
+ylabel ('Magnitude in dB');
+title ('Frequency Response of FIR HRF with N=11 using Hamming Window');
+xgrid (2);
\ No newline at end of file diff --git a/52/CH7/EX7.14/Example7_14.sce b/52/CH7/EX7.14/Example7_14.sce new file mode 100755 index 000000000..07cc102e5 --- /dev/null +++ b/52/CH7/EX7.14/Example7_14.sce @@ -0,0 +1,14 @@ +//Example 7.14
+//To Compare the Varience of Output due to A/D Conversion process
+//y(n)=0.8y(n-1)+x(n)
+clear;
+clc ;
+close ;
+n=8; //Bits
+r=100; //Range
+Q=2*r/(2^n); //Quantization Step Size
+Ve=(Q^2)/12;
+Vo=Ve*(1/(1-0.8^2));
+disp(Q,'QUANTIZATION STEP SIZE =');
+disp(Ve,'VARIANCE OF ERROR SIGNAL =');
+disp(Vo,'VARIANCE OF OUTPUT =');
\ No newline at end of file diff --git a/52/CH7/EX7.2/Example7_2.sce b/52/CH7/EX7.2/Example7_2.sce new file mode 100755 index 000000000..642fa9faa --- /dev/null +++ b/52/CH7/EX7.2/Example7_2.sce @@ -0,0 +1,13 @@ +//Example 7.2
+//To Compute Subtraction
+//(a) 0.25 from 0.5
+clear;
+clc ;
+close ;
+a=0.5;
+b=0.25;
+c=a-b;
+disp(c,'=',b,'-',a,'PART 1');
+//(a) 0.5 from 0.25
+d=b-a;
+disp(d,'=',a,'-',b,'PART 2');
\ No newline at end of file diff --git a/52/CH8/EX8.10/Example8_10.sce b/52/CH8/EX8.10/Example8_10.sce new file mode 100755 index 000000000..34376b327 --- /dev/null +++ b/52/CH8/EX8.10/Example8_10.sce @@ -0,0 +1,12 @@ +//Example 8.10
+//Develop a two band polyphase decomposition for the transfer function
+//H(z)=z^2+z+2/z^2+0.8z+0.6
+clear;
+clc ;
+close ;
+z=%z;
+HZ=(z^2+z+2)/(z^2+0.8*z+0.6);
+HZa=horner(HZ,-z);
+P0=0.5*(HZ+HZa);
+P1=0.5*(HZ-HZa);
+disp(P1/z,'+',P0,'H(z) =')
\ No newline at end of file diff --git a/52/CH8/EX8.9/Example8_9.sce b/52/CH8/EX8.9/Example8_9.sce new file mode 100755 index 000000000..345ad66f8 --- /dev/null +++ b/52/CH8/EX8.9/Example8_9.sce @@ -0,0 +1,16 @@ +//Example 8.9
+//MAXIMA SCILAB TOOLBOX REQUIRED FOR THIS PROGRAM
+//Develop a two component decomposition for the transfer function
+//and determine P0(z) and P1(z)
+clear;
+clc ;
+close ;
+syms z a n;
+HZ=(z)/(z-a);
+hn=a^n;//Inverse Z Transform of HZ
+h2n=a^(2*n);
+P0=symsum(h2n*z^(-n),n,0,%inf);
+h2n1=a^(2*n+1);
+P1=symsum(h2n1*z^(-n),n,0,%inf);
+disp(P0,'P0(Z) = ');
+disp(P1,'P1(Z) = ');
\ No newline at end of file diff --git a/52/CH9/EX9.7.a/Example9_7_a.sce b/52/CH9/EX9.7.a/Example9_7_a.sce new file mode 100755 index 000000000..fa6b75eab --- /dev/null +++ b/52/CH9/EX9.7.a/Example9_7_a.sce @@ -0,0 +1,18 @@ +//Example 9.7 (a)
+//Program To Determine Frequency Resolution of Bartlett,
+//Welch(50% Overlap) and Blackmann-Tukey Methods
+clear;
+clc;
+close;
+//Data
+Q=10;//Quality Factor
+N=1000;//Samples
+//FREQUENCY RESOLUTION CALCULATION
+K=Q;
+rb=0.89*(2*%pi*K/N);
+rw=1.28*(2*%pi*9*Q)/(16*N);
+rbt=0.64*(2*%pi*2*Q)/(3*N);
+//Display the result in command window
+disp(rb,"Resolution of Bartlett Method");
+disp(rw,"Resolution of Welch(50% overlap) Method");
+disp(rbt,"Resolution of Blackmann-Tukey Method");
\ No newline at end of file diff --git a/52/CH9/EX9.7.b/Example9_7_b.sce b/52/CH9/EX9.7.b/Example9_7_b.sce new file mode 100755 index 000000000..12e034cd8 --- /dev/null +++ b/52/CH9/EX9.7.b/Example9_7_b.sce @@ -0,0 +1,17 @@ +//Example 9.7 (b)
+//Program To Determine Record Length of Bartlett,
+//Welch(50% Overlap) and Blackmann-Tukey Methods
+clear;
+clc;
+close;
+//Data
+Q=10;//Quality Factor
+N=1000;//Samples
+//RECORD LENGTH CALCULATION
+lb=N/Q;
+lw=16*N/(9*Q);
+lbt=3*N/(2*Q);
+//Display the result in command window
+disp(lb,"Record Length of Bartlett Method");
+disp(lw,"Record Length of Welch(50% overlap) Method");
+disp(lbt,"Record Length of Blackmann-Tukey Method");
diff --git a/52/CH9/EX9.8.a/Example9_8_a.sce b/52/CH9/EX9.8.a/Example9_8_a.sce new file mode 100755 index 000000000..9a5bf6dd1 --- /dev/null +++ b/52/CH9/EX9.8.a/Example9_8_a.sce @@ -0,0 +1,12 @@ +//Example 9.8 (a)
+//Program To Determine Smallest Record Length of Bartlett Method
+clear;
+clc;
+close;
+//Data
+fr=0.01;//Frequency Resolution
+N=2400; //Samples
+//RECORD LENGTH CALCULATION
+lb=0.89/fr;
+//Display the result in command window
+disp(lb,"Record Length of Bartlett Method");
\ No newline at end of file diff --git a/52/CH9/EX9.8.b/Example9_8_b.sce b/52/CH9/EX9.8.b/Example9_8_b.sce new file mode 100755 index 000000000..9f7552a95 --- /dev/null +++ b/52/CH9/EX9.8.b/Example9_8_b.sce @@ -0,0 +1,13 @@ +//Example 9.8 (b)
+//Program To Determine Quality Factor of Bartlett Method
+clear;
+clc;
+close;
+//Data
+fr=0.01;//Frequency Resolution
+N=2400; //Samples
+lb=0.89/fr;
+//QUALITY FACTOR CALCULATION
+Q=N/lb;
+//Display the result in command window
+disp(Q,"Quality Factor of Bartlett Method");
\ No newline at end of file diff --git a/52/DEPENDENCIES/Program11_5data.sci b/52/DEPENDENCIES/Program11_5data.sci new file mode 100755 index 000000000..7d08e99f8 --- /dev/null +++ b/52/DEPENDENCIES/Program11_5data.sci @@ -0,0 +1,8 @@ +//Program 11.5 data
+//Data for Program 11.5
+A=1;
+B=2;
+C=3;
+X1=4;
+X2=5;
+X3=6;
\ No newline at end of file |
