Spectrum Analysis

Download as pptx, pdf, or txt
Download as pptx, pdf, or txt
You are on page 1of 8

SPECTRUM ANALYSIS

Analogy
Time domain

Space domain

Time t (second)

Distance x (meter)

Angular frequency (rad/s)

Wavenumber k (rad/meter) or radial average


wavenumber r

Frequency f (cycle/s)

Wavenumber k (cycle/meter)

Fourier Transform

Fourier analysis, a methodology that maps functions of space (or time) into the functions of
wavenumber (or frequency).
The fourier transform is, in general, a complex function with real and imaginary parts, that is,
It also can be written as
Where
and
The functions , , and ln are called the amplitude, phase spectrum, and power-density
spectrum respectively

Discrete Fourier Transform


The
discrete Fourier transform is periodic with a period inversely proportional to the
sample interval
= sampling wavenumber (rad/m), or
In cycle/m
= Nyquist wavenumber (rad/m), or
In cycle/m.
Because the discrete Fourier transform repeats itself each , all unique information lies
between . Hence, the Nyquist wavenumber is the largest wavenumber at our disposal.

Effect of Depth (Spector and Grant, 1970); (Blakely,


1995)
Te
ensemble average depth enters only into the factor
For value of which are , therefore
And logarithm of this factor approximates a straight line whose slope is . The
term is invariably the domination factor in the power spectrum. Hence, we can
write this with

is called the power-density spectrum, and is the slope of linear equation above.
Finally, we have

Matlab code
function DGRAV
clf
clc;
clear all;
CBA=xlsread('data.xlsx');
x=CBA(:,1);
y=CBA(:,2);
dx=100; %interval sampling
ks=1/dx; % sampling wavenumber
xi=min(x):dx:max(x);
yi=spline(x,y,xi);

You might also like