Project DSP
Project DSP
Problem ! :
clc;
clear all;
close all;
N = 16;
xn = [1 1 0];
xk = fft(xn, N);
subplot(2,2,1);
n = 0:length(xn)-1;
stem(n, xn);
xlabel('n');
ylabel('x(n)');
title('Original Signal');
k = 0:N-1;
subplot(2,2,2);
stem(k, abs(xk));
xlabel('k');
ylabel('|X(k)|');
title('Magnitude Spectrum');
subplot(2,2,3);
stem(k, angle(xk));
xlabel('k');
ylabel('?X(k)');
title('Phase Spectrum');
Output :
clc;
clear all;
close all;
N = 16;
i = sqrt(-1);
xn = [1 1 0 1 1 1 1 1 0 1 0 0 1 0 1 1];
%xk = fft(xn, N);
xk = zeros(1,N);
for k =0:1:N-1;
for n=0:1:N-1;
xk(k+1)=xk(k+1)+xn(n+1)*exp(-i*2*pi*k*n/N);
end
end
subplot(2,2,1);
n = 0:length(xn)-1;
stem(n, xn);
xlabel('n');
ylabel('x(n)');
title('Original Signal');
k = 0:N-1;
subplot(2,2,2);
stem(k, abs(xk));
xlabel('k');
ylabel('|X(k)|');
title('Magnitude Spectrum');
subplot(2,2,3);
stem(k, angle(xk));
xlabel('k');
ylabel('?X(k)');
title('Phase Spectrum');
Output :
Problem 3 :
clc;
clear all;
close all;
N = 3;
xn = [1 1 0];
xk = fft(xn, N);
subplot(2,2,1);
n = 0:length(xn)-1;
stem(n, xn);
xlabel('n');
ylabel('x(n)');
title('Original Signal');
k = 0:N-1;
subplot(2,2,2);
stem(k, abs(xk));
xlabel('k');
ylabel('|X(k)|');
title('Magnitude Spectrum');
subplot(2,2,3);
stem(k, angle(xk));
xlabel('k');
ylabel('?X(k)');
title('Phase Spectrum');
Output :
Problem 4 :
lc;
clear all;
close all;
N = 16;
xn = zeros(1, N);
for n = 0:N-1
if isprime(n)
xn(n+1)=3;
elseif mod(n,3)==0
xn(n+1)=1;
elseif mod(n,2)==0
xn(n+1)=2;
else
xn(n+1)=0;
end
end
xk = fft(xn, N);
subplot(2,2,1);
n = 0:length(xn)-1;
stem(n, xn);
xlabel('n');
ylabel('x(n)');
title('Original Signal');
k = 0:N-1;
subplot(2,2,2);
stem(k, abs(xk));
xlabel('k');
ylabel('|X(k)|');
title('Magnitude Spectrum');
subplot(2,2,3);
stem(k, angle(xk));
xlabel('k');
ylabel('?X(k)');
title('Phase Spectrum');
Output :