# 0.11 Lab 7c - power spectrum estimation  (Page 2/2)

 Page 2 / 2

Use your Pgram function to estimate the power spectrumof $y\left(n\right)$ , ${P}_{yy}\left({\omega }_{k}\right)$ . Plot ${P}_{yy}\left({\omega }_{k}\right)$ vs. ${\omega }_{k}$ .

Next, estimate the power spectrum of $y$ , ${P}_{yy}\left({\omega }_{k}\right)$ using $1/4$ of the samples of $y$ . Do this only using samples ${y}_{0}$ , ${y}_{1}$ , $\cdots$ , ${y}_{255}$ . Plot ${P}_{yy}\left({\omega }_{k}\right)$ vs. ${\omega }_{k}$ .

## Inlab report

1. Hand in your labeled plots and your Pgram code.
2. Compare the two plots. The first plot uses 4 times as many samples as the second one. Is the firstone a better estimation than the second one? Does the first give you a smoother estimation?
3. Judging from the results, when the number of samples of a discrete random variable becomes larger,will the estimated power spectrum be smoother?

## Averaged periodogram

The periodogram is a simple method, but it does not yield very good results.To get a better estimate of the power spectrum, we will introduce Bartlett's method, also known as the averaged periodogram . This method has three steps. Suppose we have a length-N sequence $x\left(n\right)$ .

1. Subdivide $x\left(n\right)$ into $K$ nonoverlapping segments of length $M$ . Denote the ${i}^{th}$ segment as ${x}_{i}\left(n\right)$ .
${x}_{i}\left(n\right)=x\left(n+iM\right),\phantom{\rule{0.277778em}{0ex}}\phantom{\rule{0.277778em}{0ex}}i=0,1,\cdots ,K-1,\phantom{\rule{0.277778em}{0ex}}\phantom{\rule{0.277778em}{0ex}}n=0,1,\cdots ,M-1$
2. For each segment ${x}_{i}\left(n\right)$ , compute its periodogram
${P}_{{x}_{i}{x}_{i}}^{\left(i\right)}\left({\omega }_{k}\right)=\frac{1}{M}{\left|\sum _{n=0}^{M-1},{x}_{i},\left(n\right),{e}^{-j{\omega }_{k}n}\right|}^{2}$
${\omega }_{k}=\frac{2\pi k}{M}$
where $k=0,1,\cdots ,M-1$ and $i=0,1,\cdots ,K-1$ .
3. Average the periodograms over all $K$ segments to obtain the averaged periodogram, ${P}_{xx}^{A}\left(k\right)$ ,
${P}_{xx}^{A}\left({\omega }_{k}\right)=\frac{1}{K}\sum _{i=0}^{K-1}{P}_{{x}_{i}{x}_{i}}^{\left(i\right)}\left({\omega }_{k}\right)$
${\omega }_{k}=\frac{2\pi k}{M}$
where $k=0,1,\cdots ,M-1$ .

Write a Matlab function called AvPgram to calculate the averaged periodogram, using the above steps.The syntax for this function should be

[P,w] = AvPgram(x, K)

where $x$ is a discrete random sequence of length $N$ and the $K$ is the number of nonoverlapping segments. The outputs of this command are $P$ , the samples of the averaged periodogram, and $w$ , the corresponding frequencies of the samples. Both $P$ and $w$ should be vectors of length M where N=KM. You may use your Pgram function.

The command A=reshape(x,M,K) will orient length $M$ segments of the vector $x$ into $K$ columns of the matrix $A$ .

Use your Matlab function AvPgram to estimate the power spectrum of $y\left(n\right)$ which was generated in the previous section. Use all 1024 samples of $y\left(n\right)$ , and let $K=16$ . Plot $P$ vs. $w$ .

## Inlab report

1. Submit your plot and your AvPgram code.
2. Compare the power spectrum that you estimated using the averaged periodogram withthe one you calculated in the previous section using the standard periodogram. What differences do you observe? Which do you prefer?

## Power spectrum and lti systems

Consider a linear time-invariant system with frequency response $H\left({e}^{j\omega }\right)$ , where ${S}_{xx}\left(\omega \right)$ is the power spectrum of the input signal, and ${S}_{yy}\left(\omega \right)$ is the power spectrum of the output signal.It can be shown that these quantities are related by

${S}_{yy}\left(\omega \right)=\mid H\left({e}^{j\omega }\right){\mid }^{2}{S}_{xx}\left(\omega \right)\phantom{\rule{4pt}{0ex}}.$

In the "Periodogram" section, the sequence $y\left(n\right)$ was generated byfiltering an i.i.d. Gaussian (mean=0, variance=1) sequence $x\left(n\right)$ , using the filter in [link] . By hand, calculate the power spectrum ${S}_{xx}\left(\omega \right)$ of $x\left(n\right)$ , the frequency response of the filter, $H\left({e}^{j\omega }\right)$ , and the power spectrum ${S}_{yy}\left(\omega \right)$ of $y\left(n\right)$ .

In computing ${S}_{xx}\left(\omega \right)$ , use the fact that $\mid ab{\mid }^{2}=a{b}^{*}\phantom{\rule{4pt}{0ex}}.$
Plot ${S}_{yy}\left(\omega \right)$ vs. $\omega$ , and compare it with the plots from the "Periodogram" and "Averaged Periodogram" sections. What do you observe?

## Inlab report

2. Submit your analytical calculations for ${S}_{xx}\left(\omega \right)$ , $H\left({e}^{j\omega }\right)$ , ${S}_{yy}\left(\omega \right)$ .
3. Compare the theoretical power spectrum of ${S}_{yy}\left(\omega \right)$ with the power spectrum estimated using the standard periodogram. What can you conclude?
4. Compare the theoretical power spectrum of ${S}_{yy}\left(\omega \right)$ with the power spectrum estimated using the averaged periodogram. What can you conclude?

## Estimate the power spectrum of a speech signal

Download the file speech.au for this section. For help on the following Matlab topics select the corresponding link: how to load and play audio signals and specgram function .

The methods used in the last two sections can only be applied to stationary random processes. However, most signals in nature are not stationary. For a nonstationary random process,one way to analyze it is to subdivide the signal into segments (which may be overlapping)and treat each segment as a stationary process. Then we can calculate the power spectrum of each segment. This yieldswhat we call a spectrogram .

While it is debatable whether or not a speech signal is actually random, in many applications it is necessary to model it as being so.In this section, you are going to use the Matlab command specgram to calculate the spectrogram of a speech signal. Read the help for the specgram function . Find out what the command does, and how to calculate and draw aspectrogram.

Draw the spectrogram of the speech signal in speech.au . When using the specgram command with no output arguments, the absolute value of the spectrogram will be plotted.Therefore you can use

speech=auread('speech.au');

to read the speech signal and use

specgram(speech);

to draw the spectrogram.

## Inlab report

1. Hand in your spectrogram plot.
2. Describe the information that the spectrogram is giving you.

how can chip be made from sand
are nano particles real
yeah
Joseph
Hello, if I study Physics teacher in bachelor, can I study Nanotechnology in master?
no can't
Lohitha
where we get a research paper on Nano chemistry....?
nanopartical of organic/inorganic / physical chemistry , pdf / thesis / review
Ali
what are the products of Nano chemistry?
There are lots of products of nano chemistry... Like nano coatings.....carbon fiber.. And lots of others..
learn
Even nanotechnology is pretty much all about chemistry... Its the chemistry on quantum or atomic level
learn
da
no nanotechnology is also a part of physics and maths it requires angle formulas and some pressure regarding concepts
Bhagvanji
hey
Giriraj
Preparation and Applications of Nanomaterial for Drug Delivery
revolt
da
Application of nanotechnology in medicine
has a lot of application modern world
Kamaluddeen
yes
narayan
what is variations in raman spectra for nanomaterials
ya I also want to know the raman spectra
Bhagvanji
I only see partial conversation and what's the question here!
what about nanotechnology for water purification
please someone correct me if I'm wrong but I think one can use nanoparticles, specially silver nanoparticles for water treatment.
Damian
yes that's correct
Professor
I think
Professor
Nasa has use it in the 60's, copper as water purification in the moon travel.
Alexandre
nanocopper obvius
Alexandre
what is the stm
is there industrial application of fullrenes. What is the method to prepare fullrene on large scale.?
Rafiq
industrial application...? mmm I think on the medical side as drug carrier, but you should go deeper on your research, I may be wrong
Damian
How we are making nano material?
what is a peer
What is meant by 'nano scale'?
What is STMs full form?
LITNING
scanning tunneling microscope
Sahil
how nano science is used for hydrophobicity
Santosh
Do u think that Graphene and Fullrene fiber can be used to make Air Plane body structure the lightest and strongest. Rafiq
Rafiq
what is differents between GO and RGO?
Mahi
what is simplest way to understand the applications of nano robots used to detect the cancer affected cell of human body.? How this robot is carried to required site of body cell.? what will be the carrier material and how can be detected that correct delivery of drug is done Rafiq
Rafiq
if virus is killing to make ARTIFICIAL DNA OF GRAPHENE FOR KILLED THE VIRUS .THIS IS OUR ASSUMPTION
Anam
analytical skills graphene is prepared to kill any type viruses .
Anam
Any one who tell me about Preparation and application of Nanomaterial for drug Delivery
Hafiz
what is Nano technology ?
write examples of Nano molecule?
Bob
The nanotechnology is as new science, to scale nanometric
brayan
nanotechnology is the study, desing, synthesis, manipulation and application of materials and functional systems through control of matter at nanoscale
Damian
Got questions? Join the online conversation and get instant answers!