Compute the DFS coefficients

Assignment Help MATLAB Programming

Reference no: EM131724778

Lab - Fourier Series and Fourier Transform

Theory:
By sampling the discrete in time Fourier transform in the frequency domain we can obtain the discrete Fourier series (DFS). The DFS can then be extended to finite- duration sequences which results in a new transform, the discrete Fourier transform (DFT). The numerical computation of the DFT for long sequences can be prohibitively time consuming. As a result, algorithms were developed to compute the DFT efficiently and is collectively known as the fast Fourier transform (FFT).

Prelab Assignment:

Lab:

1. Compute the DFS coefficients of the following periodic sequences using the DFS definition and then verify your answers using MATLAB.

a. x₁^~(n) = {4,1, -1, 1} N = 4
b. x₂^~(n) = {2, 0, 0, 0, -1, 0, 0, 0}, N = 8
c. x₃^~(n) = {1, 0, -1, -1, 0}, N = 5
d. x₄^~(n) = {0, 0, 2 j, 0, 2j, 0), N = 6
e. x₅^~(n) = (3, 2, 1), N = 3

2. Let X(e^jω) be the DTFT of a finite lenght sequence

a. Let

y₁(n) = IDFS^10-point[X(e^j0), X(e^j2^Π/10), X(e^j4^Π/10),......,X(e^j18^Π/10) ]

Determine y₁(n) using the frequency sampling theorem. Verify your answer using MATLAB.

b. Let

y₂(n) = IDFS^200-point[X(e^j0), X(e^j2^Π/200), X(e^j4^Π/200),......,X(e^j398^Π/200) ]

Determine y₂(n) using the frequency sampling theorem. Verify your answer using MATLAB.

c. Comment on your results in parts (a) and (b).

3. Consider the following finite-length sequence:

x(n) = { sinc²{(n - 50)/2}, 0 ≤ n ≤ 100;

{ 0, else.

a. Determine the DFT X(k) of x(n). Plot (using the stem function) its magnitude and phase.

b. Plot the magnitude and phase of the DTFT X(e_jω) of x(n) using MATLAB.

c. Verify that the above DFT is the sampled version of X(e^jω). It might be helpful to combine the above two plots in one graph using the hold function.

d. Is it possible to reconstruct the DTFT X(e^jω) from the DFT X(k)? If possible, give the necessary interpolation formula for reconstruction. If not possible, state why this reconstruction cannot be done.

4. Let a finite-length sequence be given by:

x(n) = { 2e^-0.9|n|, -5 ≤ n ≤ 5;

{ 0, otherwise.

Plot the DTFT X(e^jω) of the above sequence using DFT as a computation tool. Choose the length N of the DFT so that this plot appears to be a smooth graph.

5. Plot the DTFT magnitude and angle of each of the following sequences using the DFT as a computation tool. Make an educated guess about the length N so that your plots are meaningful.

a. x₁ (n) = (0.6)^|n|[u(n + 10) - u(n - 10)].
b. x₂(n) = n(0.9)ⁿ, 0 ≤ n ≤ 20.
c. x₃(n) = cos(0.5Πn) + j sin(0.5Πn), 0 ≤ n ≤ 50.
d. x(n) = {1, 2, 3, 4, 3, 2, 1}.
↑
e. x(n) = (-1,-2,-3,0,3,2,1).
↑

6. A 512-point DFT X(k) of a real-valued sequence x(n) has the DFT values:

X(0) = 20 + jα; X(5) = 20 + j30; X(k1) = -10 + j15; x(152) = 17 + j23;

X(k₂) = 20 - j30; X(k₃) = 17 - j23; X(480) = -10 - j15; x(256) = 30 + jβ;

And all other values are known to be zero.
a. Determine the real-valued coefficients α and β.
b. Determine the values of the integers k₁, k₂ and k₃.
c. Determine the energy of the signal x(n).
d. Express the sequence x(n) in a closed form.

7. Compute the N-point circular convolution for the following sequences. Plot their samples.

a. x₁(n) = sin(Πn/3)R₆(n),   x₂(n) = cos(Πn/4)R₈(n);   N = 10
b. x₁(n) = cos(2Πn/3)R_N(n),   x₂(n) = sin(2Πn/4)R_N(n);   N = 32
c. x₁(n) = (0.8)ⁿR_N(n),   x₂(n) = (-0.8)ⁿR_N(n); N = 10
d. x₁(n) = nR_N(n),   x₂(n) = (N - n)R_N(n); N = 10
e. x₁(n) = (0.8)ⁿR₂₀x₂(n),   x₂(n) = u(n) -u(n-40); N = 10

8. Let x(n) be a uniformly distributed random number between [-1,1] for 0 ≤ n ≤ 10⁶.

Let

h(n) = sin(0.4Πn), 0 ≤ n ≤ 100

a. Consider the conv function, determine the output sequence y(n) = x(n)*h(n)

b. Consider the overlap and save method of block convolution along with the FFT algorithm to implement high speed block convolution. Using this approach, determine y(n) with FFT sizes of 1024, 2048 and 4096.

c. Compare these approaches in terms of the convolution results and their execution times.

Reference no: EM131724778

Questions Cloud

Discuss an attacker interposes during key exchange : Man-in-the-Middle Attack: An attacker interposes during key exchange, acting as the client to the server and as the server to the client

What is the probability of answering exactly : Given a 12?-question ?true-false quiz, what is the probability of answering exactly 8 of the questions correctly by? guessing?

What will be your approximate house payment each month : Discuss, realistically, your professional aspirations and how you will attain them.

Briefly compare the data for the three countries you selecte : Briefly compare the data for the three countries you selected. Identify the types of additional information needed to make a final decision.

Compute the DFS coefficients : Fourier Series Compute the DFS coefficients of the following periodic sequences using the DFS definition and then verify your answers using MATLAB.

Discuss health insurance portability and accountability act : such as the Health Insurance Portability and Accountability Act (HIPAA) in the U.S. and similar laws in the U.K. encourage

Explain how benefits are communicated to employees : Explain how benefits are communicated to employees. Analyze the environmental factors that have influenced the compensation and benefits of the organization.

What is the definition of socioeconomic status : Based on the readings and presentation, what is the definition of socioeconomic status? Do you know someone who is from low SES families

Description of the cultural consideration : When you begin work with couples and families, you likely will conclude after only a few cases that all couples and families are different in at least a few way

User Account

All Pages