diff options
Diffstat (limited to '758/CH6')
31 files changed, 429 insertions, 0 deletions
diff --git a/758/CH6/EX6.1/Ex_6_1.sce b/758/CH6/EX6.1/Ex_6_1.sce new file mode 100755 index 000000000..653fb1c54 --- /dev/null +++ b/758/CH6/EX6.1/Ex_6_1.sce @@ -0,0 +1,27 @@ +//Example 6.1
+clc;clear;close;
+x1=[1 1 2 2];
+x2=[1 2 3 4];
+ylength=length(x1);
+//Calculation of linear convolution
+z=convol(x1,x2);
+//Calculation of circular convolution
+for n=1:ylength
+ y(n)=0;
+ for k=1:ylength,
+ l=n-k+1;
+ if l <= 0 then
+ l=l+ylength;
+ end
+ y(n)=y(n)+(x1(k)*x2(l));
+ end
+end
+//Calculation of circular convolution using DFT and IDFT
+X1=dft(x1,-1);
+X2=dft(x2,-1);
+Y1=X1.*X2;
+y1=dft(Y1,1);
+y1=clean(y1);
+disp(z,'Linear Convolution sequence is z(n): ');
+disp(y,'Circular Convolution sequence is y(n): ');
+disp(y1,'Circular Convolution sequence calculated using DFT-IDFT method is y(n): ');
\ No newline at end of file diff --git a/758/CH6/EX6.10/Ex_6_10.jpg b/758/CH6/EX6.10/Ex_6_10.jpg Binary files differnew file mode 100755 index 000000000..35e314ade --- /dev/null +++ b/758/CH6/EX6.10/Ex_6_10.jpg diff --git a/758/CH6/EX6.10/Ex_6_10.sce b/758/CH6/EX6.10/Ex_6_10.sce new file mode 100755 index 000000000..0719a285e --- /dev/null +++ b/758/CH6/EX6.10/Ex_6_10.sce @@ -0,0 +1,19 @@ +//Example 6.10
+
+clc;clear;close;
+x=[1 1 2 2 3 3];
+//Calculation of DFT
+X=dft(x,-1);
+X=clean(X);
+disp(x,'Given Sequence is x(n): ');
+disp(X,'DFT of the Sequence is X(k): ');
+subplot(3,1,1);
+plot2d3(x);
+title('Given Seqence x[n]:');ylabel('Amplitude-->');xlabel('n-->');
+subplot(3,1,2);
+plot2d3(abs(X));
+title('Magnitude Spectrum |X(k)|');xlabel('k-->');
+subplot(3,1,3);
+plot2d3(atan(X));
+title('Phase Spectrum /_X(k)');xlabel('k-->');
+
diff --git a/758/CH6/EX6.11/Ex_6_11.sce b/758/CH6/EX6.11/Ex_6_11.sce new file mode 100755 index 000000000..6af305e3f --- /dev/null +++ b/758/CH6/EX6.11/Ex_6_11.sce @@ -0,0 +1,12 @@ +//Example 6.11
+
+clc;clear;close;
+N=8;A=1/4;
+n=0:N-1;
+x=A^n;
+//Calculation of DFT
+X=dft(x,-1);
+X=clean(X);
+disp(x,'Given Sequence is x(n): ');
+disp(N,'N=')
+disp(X,'N-point DFT of the Sequence is X(k): ');
\ No newline at end of file diff --git a/758/CH6/EX6.12/Ex_6_12.sce b/758/CH6/EX6.12/Ex_6_12.sce new file mode 100755 index 000000000..0f088a5ed --- /dev/null +++ b/758/CH6/EX6.12/Ex_6_12.sce @@ -0,0 +1,11 @@ +//Example 6.12
+
+clc;clear;close;
+N=4;
+n=0:N-1;
+x=cos(%pi/4*n);
+//Calculation of DFT
+X=dft(x,-1);
+X=clean(X);
+disp(x,'Given Sequence is x(n): ');
+disp(X,'DFT of the Sequence is X(k): ');
\ No newline at end of file diff --git a/758/CH6/EX6.13/Ex_6_13.sce b/758/CH6/EX6.13/Ex_6_13.sce new file mode 100755 index 000000000..41780add6 --- /dev/null +++ b/758/CH6/EX6.13/Ex_6_13.sce @@ -0,0 +1,8 @@ +//Example 6.13
+clc;clear;close;
+X=[1 2 3 4];
+//Calculation of IDFT
+x=dft(X,1);
+x=clean(x);
+disp(X,'DFT of the Sequence is X(k): ');
+disp(x,'Sequence is x(n): ');
diff --git a/758/CH6/EX6.14/Ex_6_14.sce b/758/CH6/EX6.14/Ex_6_14.sce new file mode 100755 index 000000000..c34b8b74e --- /dev/null +++ b/758/CH6/EX6.14/Ex_6_14.sce @@ -0,0 +1,8 @@ +//Example 6.14
+clc;clear;close;
+X=[3 2+%i 1 2-%i];
+//Calculation of IDFT
+x=dft(X,1);
+x=clean(x);
+disp(X,'DFT of the Sequence is X(k): ');
+disp(x,'Sequence is x(n): ');
diff --git a/758/CH6/EX6.15/Ex_6_15.sce b/758/CH6/EX6.15/Ex_6_15.sce new file mode 100755 index 000000000..7ddff5c60 --- /dev/null +++ b/758/CH6/EX6.15/Ex_6_15.sce @@ -0,0 +1,7 @@ +//Example 6.15
+
+clc;clear;
+x=[1 2 3 4 4 3 2 1];
+X=clean(fft(x));
+disp(x,'x(n)=');
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/758/CH6/EX6.16/Ex_6_16.sce b/758/CH6/EX6.16/Ex_6_16.sce new file mode 100755 index 000000000..bf839237c --- /dev/null +++ b/758/CH6/EX6.16/Ex_6_16.sce @@ -0,0 +1,7 @@ +//Example 6.16
+
+clc;clear;
+x=[0 1 2 3 4 5 6 7];
+X=clean(fft(x));
+disp(x,'x(n)=');
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/758/CH6/EX6.17/Ex_6_17.sce b/758/CH6/EX6.17/Ex_6_17.sce new file mode 100755 index 000000000..b2da68756 --- /dev/null +++ b/758/CH6/EX6.17/Ex_6_17.sce @@ -0,0 +1,8 @@ +//Example 6.17
+
+clc;clear;
+n=0:7;
+x=2^n;
+X=clean(fft(x));
+disp(x,'x(n)=');
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/758/CH6/EX6.18/Ex_6_18.sce b/758/CH6/EX6.18/Ex_6_18.sce new file mode 100755 index 000000000..454552e07 --- /dev/null +++ b/758/CH6/EX6.18/Ex_6_18.sce @@ -0,0 +1,7 @@ +//Example 6.18
+
+clc;clear;
+x=[0 1 2 3];
+X=clean(fft(x));
+disp(x,'x(n)=');
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/758/CH6/EX6.19/Ex_6_19.sce b/758/CH6/EX6.19/Ex_6_19.sce new file mode 100755 index 000000000..3649391aa --- /dev/null +++ b/758/CH6/EX6.19/Ex_6_19.sce @@ -0,0 +1,7 @@ +//Example 6.19
+
+clc;clear;
+x=[1 2 3 4 4 3 2 1];
+X=clean(fft(x));
+disp(x,'x(n)=');
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/758/CH6/EX6.2/Ex_6_2.sce b/758/CH6/EX6.2/Ex_6_2.sce new file mode 100755 index 000000000..09659865d --- /dev/null +++ b/758/CH6/EX6.2/Ex_6_2.sce @@ -0,0 +1,21 @@ +//Example 6.2
+
+clc;clear;close;
+x=[1 2];
+h=[1 2 4];
+
+//Calculation of linear convolution
+y=convol(x,h);
+disp(x,'Input Sequence is x(n): ');
+disp(h,'Impulse respnose of FIR filter h(n): ');
+disp(y,'Output sequence is y(n): ');
+subplot(3,1,1);
+plot2d3(x);
+title('Input Seqence x[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,2);
+plot2d3(h);
+title('Impulse Response h[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,3);
+plot2d3(y);
+title('Output Seqence y[n]=x[n]*h[n] :');ylabel('Amplitude-->');xlabel('n-->')
+
diff --git a/758/CH6/EX6.20/Ex_6_20.sce b/758/CH6/EX6.20/Ex_6_20.sce new file mode 100755 index 000000000..b3268a041 --- /dev/null +++ b/758/CH6/EX6.20/Ex_6_20.sce @@ -0,0 +1,8 @@ +//Example 6.20
+
+clc;clear;
+n=0:7;
+x=2^n;
+X=clean(fft(x));
+disp(x,'x(n)=');
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/758/CH6/EX6.21/Ex_6_21.sce b/758/CH6/EX6.21/Ex_6_21.sce new file mode 100755 index 000000000..feee5f4b4 --- /dev/null +++ b/758/CH6/EX6.21/Ex_6_21.sce @@ -0,0 +1,8 @@ +//Example 6.21
+
+clc;clear;
+n=0:7;
+x=n+1;
+X=clean(fft(x));
+disp(x,'x(n)=');
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/758/CH6/EX6.22/Ex_6_22.sce b/758/CH6/EX6.22/Ex_6_22.sce new file mode 100755 index 000000000..c5d0f9bef --- /dev/null +++ b/758/CH6/EX6.22/Ex_6_22.sce @@ -0,0 +1,8 @@ +//Example 6.21
+
+clc;clear;
+n=0:3;
+x=cos(n*%pi/2);
+X=clean(fft(x));
+disp(x,'x(n)=');
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/758/CH6/EX6.23/Ex_6_23.sce b/758/CH6/EX6.23/Ex_6_23.sce new file mode 100755 index 000000000..c3c50b5f8 --- /dev/null +++ b/758/CH6/EX6.23/Ex_6_23.sce @@ -0,0 +1,7 @@ +//Example 6.23
+
+clc;clear;
+X=[6 -2+2*%i -2 -2-2*%i];
+x=clean(ifft(X));
+disp(X,'X(k)=');
+disp(x,'x(n)=');
diff --git a/758/CH6/EX6.24/Ex_6_24.sce b/758/CH6/EX6.24/Ex_6_24.sce new file mode 100755 index 000000000..426a75413 --- /dev/null +++ b/758/CH6/EX6.24/Ex_6_24.sce @@ -0,0 +1,7 @@ +//Example 6.24
+
+clc;clear;
+X=[20 -5.828-2.414*%i 0 -0.172-0.414*%i 0 -0.172+0.414*%i 0 -5.828+2.414*%i];
+x=round(clean(ifft(X)));
+disp(X,'X(k)=');
+disp(x,'x(n)=');
diff --git a/758/CH6/EX6.25/Ex_6_25.sce b/758/CH6/EX6.25/Ex_6_25.sce new file mode 100755 index 000000000..acfb9b9fe --- /dev/null +++ b/758/CH6/EX6.25/Ex_6_25.sce @@ -0,0 +1,7 @@ +//Example 6.25
+
+clc;clear;
+X=[255 48.63+166.05*%i -51+102*%i -78.63+46.05*%i -85 -78.63-46.05*%i -51-102*%i 48.63-166.05*%i];
+x=round(clean(ifft(X)));
+disp(X,'X(k)=');
+disp(x,'x(n)=');
diff --git a/758/CH6/EX6.26/Ex_6_26.sce b/758/CH6/EX6.26/Ex_6_26.sce new file mode 100755 index 000000000..0173d8879 --- /dev/null +++ b/758/CH6/EX6.26/Ex_6_26.sce @@ -0,0 +1,7 @@ +//Example 6.26
+
+clc;clear;
+X=[36 -4+9.656*%i -4+4*%i -4+1.656*%i -4 -4-1.656*%i -4-4*%i -4-9.656*%i ];
+x=round(clean(ifft(X)));
+disp(X,'X(k)=');
+disp(x,'x(n)=');
diff --git a/758/CH6/EX6.27/Ex_6_27.sce b/758/CH6/EX6.27/Ex_6_27.sce new file mode 100755 index 000000000..b3ff94078 --- /dev/null +++ b/758/CH6/EX6.27/Ex_6_27.sce @@ -0,0 +1,9 @@ +//Example 6.27
+
+clc;clear;
+t=0:0.0025:0.0175;
+f=50;
+x=sin(2*%pi*f*t);
+X=clean(fft(x));
+disp(x,'x(n)=');
+disp(X,'X(k)=');
\ No newline at end of file diff --git a/758/CH6/EX6.3/Ex_6_3.sce b/758/CH6/EX6.3/Ex_6_3.sce new file mode 100755 index 000000000..cb432309a --- /dev/null +++ b/758/CH6/EX6.3/Ex_6_3.sce @@ -0,0 +1,21 @@ +//Example 6.3
+
+clc;clear;close;
+x=[1 1 1];
+h=[1 1 1];
+
+//Calculation of linear convolution
+y=convol(x,h);
+disp(x,'First Sequence is x(n): ');
+disp(h,'Second Sequence is h(n): ');
+disp(y,'Output sequence is y(n): ');
+subplot(3,1,1);
+plot2d3(x);
+title('First Seqence x[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,2);
+plot2d3(h);
+title('Second Seqence h[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,3);
+plot2d3(y);
+title('Convolution Seqence y[n]=x[n]*h[n] :');ylabel('Amplitude-->');xlabel('n-->')
+
diff --git a/758/CH6/EX6.34/Ex_6_34.sce b/758/CH6/EX6.34/Ex_6_34.sce new file mode 100755 index 000000000..c89d456e7 --- /dev/null +++ b/758/CH6/EX6.34/Ex_6_34.sce @@ -0,0 +1,27 @@ +//Example 6.34
+
+clc;clear;close;
+h=[2 2 1];
+x=[3 0 -2 0 2 1 0 -2 -1 0];
+M=length(h); //length of impulse response
+L=2^M; //length of FFT/IFFT operation
+N=L-M+1;
+xl=length(x);
+K=ceil(xl/N); //number of iterations
+h=[h zeros(1,L-M)];
+x=[x x(1:K*N-xl)];
+H=fft(h);
+y=zeros(1,M-1);
+for k=0:K-1
+ xk=[x(k*N+1:(k+1)*N) zeros(1,M-1)];
+ Xk=fft(xk);
+ Yk=H.*Xk;
+ yk=ifft(Yk);
+ yk=clean(yk);
+ y=[y(1:k*N) y(k*N+1:k*N+M-1)+yk(1:M-1) yk(M:L)];
+ disp(k+1,'Segment =');
+ disp(xk,'xk(n)=');
+ disp(yk,'yk(n)=');
+end
+y=y(1:xl+M-1);
+disp(y,'Output Sequence is y(n): ');
diff --git a/758/CH6/EX6.35/Ex_6_35.sce b/758/CH6/EX6.35/Ex_6_35.sce new file mode 100755 index 000000000..11e55b4a0 --- /dev/null +++ b/758/CH6/EX6.35/Ex_6_35.sce @@ -0,0 +1,25 @@ +//Example 6.35
+
+clc;clear;close;
+h=[2 2 1];
+x=[3 0 -2 0 2 1 0 -2 -1 0];
+M=length(h); //length of impulse response
+L=2^M; //length of FFT/IFFT operation
+N=L-M+1;
+xl=length(x);
+K=ceil(xl/N); //number of iterations
+h=[h zeros(1,L-M)];
+x=[zeros(1,M-1) x x(1:K*N-xl)];
+H=fft(h);
+for k=0:K-1
+ xk=x(k*N+1:(k+1)*N+M-1);
+ Xk=fft(xk);
+ Yk=H.*Xk;
+ yk=ifft(Yk);
+ yk=clean(yk);
+ y=[y(1:k*N) yk(M:L)];
+ disp(k+1,'Segment =');
+ disp(xk,'xk(n)=');
+ disp(yk,'yk(n)=');
+end
+disp(y,'Output Sequence is y(n): ');
diff --git a/758/CH6/EX6.36/Ex_6_36.sce b/758/CH6/EX6.36/Ex_6_36.sce new file mode 100755 index 000000000..1cc3bb2fe --- /dev/null +++ b/758/CH6/EX6.36/Ex_6_36.sce @@ -0,0 +1,28 @@ +//Example 6.36
+
+clc;clear;close;
+x=[1 0 0 1];
+h=[4 3 2 1];
+ylength=length(x)+length(h)-1;
+xlength=length(x);
+x=[zeros(1,length(h)-1) x zeros(1,length(h)-1)];
+y=0;
+//Calculation of cross correlation
+for n=1:ylength;
+ y(n)=x*[zeros(1,n-1) h zeros(1,ylength-n)]'; //this instruction performs cross correlation of x & h
+end
+
+disp(x,'First Sequence is x(n): ');
+disp(h,'Second Sequence is h(n): ');
+disp(y,'Correlation Sequence y[n] is');
+figure;
+subplot(3,1,1);
+plot2d3(x);
+title('First Seqence x[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,2);
+plot2d3(h);
+title('Second Seqence h[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,3);
+plot2d3(y);
+title('Correlation Seqence y[n]:');ylabel('Amplitude-->');xlabel('n-->')
+
diff --git a/758/CH6/EX6.37/Ex_6_37.sce b/758/CH6/EX6.37/Ex_6_37.sce new file mode 100755 index 000000000..42711368e --- /dev/null +++ b/758/CH6/EX6.37/Ex_6_37.sce @@ -0,0 +1,34 @@ +//Example 6.37
+
+clc;clear;close;
+x=[1 0 0 1];
+h=[4 3 2 1];
+ylength=length(x);
+y=0;
+//Calculation of circular correlation
+for n=1:ylength,
+ y(n)=0;
+ for k=1:ylength,
+ l=k-n+1;
+ if l <= 0 then
+ l=l+ylength;
+ end
+ y(n)=y(n)+(x(k)*h(l));
+ end
+ y(n)=y(n)/4;
+end
+
+disp(x,'First Sequence is x(n): ');
+disp(h,'Second Sequence is h(n): ');
+disp(y,'Correlation Sequence y[n] is');
+figure;
+subplot(3,1,1);
+plot2d3(x);
+title('First Seqence x[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,2);
+plot2d3(h);
+title('Second Seqence h[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,3);
+plot2d3(y);
+title('Correlation Seqence y[n]:');ylabel('Amplitude-->');xlabel('n-->')
+
diff --git a/758/CH6/EX6.4/Ex_6_4.sce b/758/CH6/EX6.4/Ex_6_4.sce new file mode 100755 index 000000000..95a51d54e --- /dev/null +++ b/758/CH6/EX6.4/Ex_6_4.sce @@ -0,0 +1,25 @@ +//Example 6.4
+
+clc;clear;close;
+a=0.5;
+n=1:50;
+x=ones(1,50);
+h=a^n;
+
+//Calculation of linear convolution
+for i=1:50
+ y(1,i)=sum(h(1:i));
+end
+disp('First Sequence is x(n)=u(n) ');
+disp(a,'Second Sequence is h(n)=a^n*u(n) where a= ');
+disp(y,'Output sequence is y(n): ');
+subplot(3,1,1);
+plot2d3(x);
+title('First Seqence x[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,2);
+plot2d3(h);
+title('Second Seqence h[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,3);
+plot2d3(y);
+title('Convolution Seqence y[n]=x[n]*h[n] :');ylabel('Amplitude-->');xlabel('n-->')
+
diff --git a/758/CH6/EX6.5/Ex_6_5.sce b/758/CH6/EX6.5/Ex_6_5.sce new file mode 100755 index 000000000..807ed3286 --- /dev/null +++ b/758/CH6/EX6.5/Ex_6_5.sce @@ -0,0 +1,22 @@ +//Example 6.5
+clc;clear;close;
+x=[1 2 3];xmin=0;nx=xmin:length(x)+xmin-1;
+h=[1 2 -2 -1];hmin=-1;nh=length(h)+hmin-1;
+
+//Calculation of linear convolution
+y=convol(x,h);
+ymin=xmin+hmin;ny=ymin:length(y)+ymin-1;
+
+disp(x,'First Sequence is x(n): ');
+disp(h,'Second Sequence is h(n): ');
+disp(y,'Output sequence is y(n): ');
+subplot(3,1,1);
+plot2d3(nx,x);
+title('First Seqence x[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,2);
+plot2d3(nh,h);
+title('Second Seqence h[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,3);
+plot2d3(ny,y);
+title('Convolution Seqence y[n]=x[n]*h[n] :');ylabel('Amplitude-->');xlabel('n-->')
+
diff --git a/758/CH6/EX6.6/Ex_6_6.sce b/758/CH6/EX6.6/Ex_6_6.sce new file mode 100755 index 000000000..73dd97420 --- /dev/null +++ b/758/CH6/EX6.6/Ex_6_6.sce @@ -0,0 +1,23 @@ +//Example 6.6
+
+clc;clear;close;
+x=[1 1 0 1 1];xmin=-2;nx=xmin:length(x)+xmin-1;
+h=[1 -2 -3 4];hmin=-3;nh=length(h)+hmin-1;
+
+//Calculation of linear convolution
+y=convol(x,h);
+ymin=xmin+hmin;ny=ymin:length(y)+ymin-1;
+
+disp(x,'First Sequence is x(n): ');
+disp(h,'Second Sequence is h(n): ');
+disp(y,'Output sequence is y(n): ');
+subplot(3,1,1);
+plot2d3(nx,x);
+title('First Seqence x[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,2);
+plot2d3(nh,h);
+title('Second Seqence h[n]:');ylabel('Amplitude-->');xlabel('n-->')
+subplot(3,1,3);
+plot2d3(ny,y);
+title('Convolution Seqence y[n]=x[n]*h[n] :');ylabel('Amplitude-->');xlabel('n-->')
+
diff --git a/758/CH6/EX6.8/Ex_6_8.sce b/758/CH6/EX6.8/Ex_6_8.sce new file mode 100755 index 000000000..4c03f0861 --- /dev/null +++ b/758/CH6/EX6.8/Ex_6_8.sce @@ -0,0 +1,10 @@ +//Example 6.8
+
+clc;clear;close;
+L=3;A=1/4;
+x=A*ones(1,L);
+//Calculation of DFT
+X=dft(x,-1);
+X=clean(X);
+disp(x,'Given Sequence is x(n): ');
+disp(X,'DFT of the Sequence is X(k): ');
\ No newline at end of file diff --git a/758/CH6/EX6.9/Ex_6_9.sce b/758/CH6/EX6.9/Ex_6_9.sce new file mode 100755 index 000000000..978ddc838 --- /dev/null +++ b/758/CH6/EX6.9/Ex_6_9.sce @@ -0,0 +1,11 @@ +//Example 6.9
+
+clc;clear;close;
+L=3;A=1/5;
+n=-1:1;
+x=A*ones(1,L);
+//Calculation of DFT
+X=dft(x,-1);
+X=clean(X);
+disp(x,'Given Sequence is x(n): ');
+disp(X,'DFT of the Sequence is X(k): ');
\ No newline at end of file |
