Code covered by the BSD License

# Ingle & Proakis Chapter 3 Solutions

### Edgar Guevara (view profile)

13 Sep 2007 (Updated )

P308.m
```% P3.8
colordef white; clear;  clc;
M = 250;    k = -M:M;   w = (pi/M)*k;   % [0, pi] axis divided into 501 points.

n = -50:50;
x = exp(j*0.1*pi*n).*(stepseq(0,-50,50)-stepseq(20,-50,50));
[xe,xo,m] = evenoddcomplex(x,n);
xr = real(x);               xi = imag(x);
XE = dft(xe,length(m));     XO = dft(xo,length(m));
fXE = idft(XE,length(XE));  fXO = idft(XO,length(XO));
figure(1);
subplot(2,2,1); plot(m,fXE,'k'); axis tight;
xlabel('frequency in pi units'); title('Conjugate Symmetric Part Inverse Fourier Transform'); ylabel('|XE-1|');
subplot(2,2,3); plot(m,fXO,'k'); axis tight;
xlabel('frequency in pi units'); title('Conjugate Antisymmetric Part Inverse Fourier Transform'); ylabel('|XO-1|');
subplot(2,2,2); plot(m,xr,'k');  axis tight;
xlabel('frequency in pi units'); title('Real part of x(n)'); ylabel('Re((x(n))');
subplot(2,2,4); plot(m,xi,'k');    axis tight;
xlabel('frequency in pi units'); title('Imaginary part of x(n)'); ylabel('Im(x(n))');```