Reference no: EM132274291
Questions -
Questions related to filter design, signal processing sequences, ODE, Fourier transform, detail descriptions is attached.
Exercise 1 - Fix a constant λ > 0. Find a solution to the delay differential equation
y'(x) = - λy(x - 1), y(0) = 1,
where we are looking for a function y = y(x) that satisfies this equation.
Does your solution decay over time? At what rate? Do you think the decay rate is somehow related to the delay in the differential equation? (Here the delay is 1, since we have a x -1 in the DE.) Test your idea by solving a related equation y'(x) = - λy(x - x0).
Exercise 2 - A sequence a = (..., a-2, a-1, a0, a1, a2, ...) is called anti-symmetric if $$a_{-n} = -a_n \mbox{ for all $n$. }$$ For instance, the sequence a = (..., 0, 0, 1, - 3, 5, 0, - 5, 3, - 1, 0, 0, ...) is anti-symmetric, when we label the middle 0 as a0.
- Show that the Fourier transform a^(θ) of an anti-symmetric sequence is purely imaginary (its real part is zero.)
- Show that a^(1 - θ) = -a^(θ) for all θ if a is anti-symmetric.
- Show that the convolution at a * a is symmetric whenever a is anti-symmetric.
Hint: If you are finding this hard for general sequences a, try to solve it for that one example
a = (..., 0, 0, 1, - 3, 5, 0, - 5, 3, -1, 0, 0, ...).
Exercise 3 - Here's a little trick to keep convolution from making our sequences "too long."
Recall if you convolve two sequence a, b of length N, the result is of length 2N-1.
To "fix" this, we can define circular convolution as c = a ⊗ b, using the formula $$c_k = \sum_{j=0}^{N-1} a_{k-j}b_j, \mbox{ with $k-j$ computed modulo $N$.}$$
So, for instance, the convolution of length-3 sequences is given as
(a0, a1, a2) ⊗ (b0, b1, b2) = (a0b0 + a1b2 + a2b1, a0b1 + a1b0 + a2b2, a0b2 + a1b1 + a2b0).
(Notice each component on the RHS has 3 terms, and the indices in each component add up to the same number, mod 3.)
1. Compute the circular convolution: (1, 1, 0, 0) ⊗ (1, 2, 3, 4). The result should have length 4.
2. Compute the regular convolution: (1, 1, 0, 0) * (1, 2, 3, 4). The result should have length 7.
3. Compare your two answers. What parts are the same, what are different?
4. Compute the DFT of the three sequences a = (1, 1, 0, 0), b = (1, 2, 3, 4), c = (1, 1, 0, 0) ⊗ (1, 2, 3, 4). Either by hand, or using Python. (Each result will be a length-4 sequence, with complex values.)
5. Verify your results in #4, by showing that the DFT of c is equal to the product of the DFT of a times the DFT of b.
DFT = Discrete Fourier transform. It takes vectors of length N and produces vectors of length N, using that formula we saw in class.
Exercise 4 - In the notes to lecture 18, we took the 2D Fourier transform of the identity matrix. As an image, this took a line pointing at 45 degrees down, and transformed it into a line at 45 degrees up.
Write some code to explore this idea further. Take a big array of zeros, and then fill in a line at some angle. Take the 2D FFT and display it -- it should look like a line, going at a different direction.
Once you get the lines working, see if you can describe the relation between the angle (or slope) of the input and output lines.
Can you prove this relationship between the angles?
Exercise 5 - This is an exploration exercise.
- Take a geometric figure, place into a NxN array to create an image.
- It's a good idea to take N to be a power of 2, like 128 or 256.
- Take the 2D FFT of the image, and save the absolute values as a NxN array.
- Plot the original image and absolute value of the FFT of the image.
- You might also try fftshift on the image (see Lecture 18 again).
- Can you see any connection between the geometric figure and its FFT?
For instance, we saw in Exercise 4 that the FFT of a line was a line.
What is the FFT of a square? A circle? A triangle? A polygon?
Try the outline of a circle, and a solid (filled in) circle. Also for the other shapes.
Try creating a pattern in the image.
- Parallel lines
- Grid lines in x and y direction
- Rectangular array of dots
- Hexagonal array of dots. What does the FFT of these things look like?
Note: There is no right answer here, I just want you to explore some shapes and patterns.
Historical note: In the early 1900s when scientists were trying to discover the atomic structure of crystals, they could essentially "see" the FFT of the arrangement of atoms in the crystal, but not see the arrangement itself. So it was an important challenge to be able to guess the arrangement only by looking at the FFT of the image. (Actually, they could only look at the absolute value of the FFT.)
The structure of DNA, a double helix, was discovered by Watson and Crick in 1953. Their work depended heavily upon the X-ray images produced by Rosalind Franklin, which is essentially the FFT images as above. Watson and Crick received the Nobel prize in 1962; Franklin died of cancer in 1958.
Attachment:- Assignment Files.rar