Line 23: Line 23:
  
 
Q: Does anyone have any insight in determining the power spectral density of y?  I've read through the links and listened to older lectures from previous semesters, but am still lost.
 
Q: Does anyone have any insight in determining the power spectral density of y?  I've read through the links and listened to older lectures from previous semesters, but am still lost.
 +
A1 : Probably this might help : http://en.wikipedia.org/wiki/Wiener%E2%80%93Khinchin_theorem

Revision as of 21:43, 23 January 2013

Hello Class, before you start doing Lab2, here are some background information you might want to read:
http://en.wikipedia.org/wiki/Window_function
http://en.wikipedia.org/wiki/Spectral_density
http://en.wikipedia.org/wiki/Orthogonal_transformation

Also as Additional Information: http://www.mathworks.com/help/matlab/ref/mtimes.html http://www.mathworks.com/help/matlab/ref/times.html

Q&A Section
Q: Anyone have an idea of how to compute the squared DFT magnitude in MATLAB? Is this the sum of pixels across the window multipled by cos(n)+isin(n) where n is the location between 0-64 in the window?
A: I would just x.^2 to be honest :)

Q: What does window function do?
A: Rememeber how rect goes to sinc from lab1? Everytime you sample an input segment and perform FFT, you are essentially applying rectangle window to the input samples. Therefore to reduce the effect of artifacts, people apply different type of windows before applying orthogonal transforms.

Q: In the second part, why are we adding 127 to each pixel value before displaying them, instead of scaling them to 0 - 255? Because, in my understanding, adding 127, would increase the (0,0) frequency component, and thus giving us the illusion that the image is smoother.
A1: I think it's just for visualization. For me, min(min(y)) = -121.79 and you can't display that. So, shifting by 127 makes pixel values non-negative.
A2: As mentioned from A1. Given X~U(-0.5,0.5), to correctly display the image using 8bit gray scale, you have to make the best out of your quantization levels. Therefore first add 0.5 to the average then multiply by 255 as stated from the problem statement.

Q: Is anyone having trouble using the image() command in Matlab for displaying their 'random image' in part 2? Mine comes up as just a black square every time. I set the colormap to gray(256) and I'm converting to an 8-bit representation using uint8(). I have checked my matrix and the values range from 0 to 255. Any insight?
A1: Did you add 127 to y (y + 127)?

Q: Does anyone have any insight in determining the power spectral density of y? I've read through the links and listened to older lectures from previous semesters, but am still lost. A1 : Probably this might help : http://en.wikipedia.org/wiki/Wiener%E2%80%93Khinchin_theorem

Alumni Liaison

Ph.D. on Applied Mathematics in Aug 2007. Involved on applications of image super-resolution to electron microscopy

Francisco Blanco-Silva