spectralanalysis.nb

analysis of spectral in mathematica

(************** Content-type: application/mathematica **************
                     CreatedBy='Mathematica 5.0'

                    Mathematica-Compatible Notebook

This notebook can be used with any Mathematica-compatible
application, such as Mathematica, MathReader or Publicon. The data
for the notebook starts with the line containing stars above.

To get the notebook into a Mathematica-compatible application, do
one of the following:

* Save the data starting with the line of stars above into a file
  with a name ending in .nb, then open the file inside the
  application;

* Copy the data starting with the line of stars above to the
  clipboard, then use the Paste menu command inside the application.

Data for notebooks contains only printable 7-bit ASCII and can be
sent directly in email or through ftp in text mode.  Newlines can be
CR, LF or CRLF (Unix, Macintosh or MS-DOS style).

NOTE: If you modify the data for this notebook not in a Mathematica-
compatible application, you must delete the line below containing
the word CacheID, otherwise Mathematica-compatible applications may
try to use invalid cache data.

For more information on notebooks and Mathematica-compatible 
applications, contact Wolfram Research:
  web: http://www.wolfram.com
  email: [e-mail]
  phone: [phone] (U.S.)

Notebook reader applications are available free of charge from 
Wolfram Research.
*******************************************************************)

(*CacheID: 232*)


(*NotebookFileLineBreakTest
NotebookFileLineBreakTest*)
(*NotebookOptionsPosition[    135597,       2951]*)
(*NotebookOutlinePosition[    136241,       2973]*)
(*  CellTagsIndexPosition[    136197,       2969]*)
(*WindowFrame->Normal*)



Notebook[{
Cell[BoxData[{
    StyleBox[\(ClearAll;\),
      FontSize->18], "\[IndentingNewLine]", 
    StyleBox[\(ampSpectrum[x_List] := 
        Block[{len, y}, \[IndentingNewLine]len = 
            1 + Floor[Length[x]/2]; \[IndentingNewLine]y = 
            Take[InverseFourier[x], \[IndentingNewLine]{1, 
                len}\[IndentingNewLine]]; \[IndentingNewLine]y = 
            Abs[y]; \[IndentingNewLine]y /= Max[y]; \[IndentingNewLine]Return[
            y];\[IndentingNewLine]]\),
      FontSize->18], "\[IndentingNewLine]", 
    StyleBox[\(SpectrumPlot[x_List] := \[IndentingNewLine]ListPlot[
            ampSpectrum[x], PlotJoined \[Rule] True, 
            PlotRange \[Rule] {0, 1}];\),
      FontSize->18], "\[IndentingNewLine]", 
    StyleBox[\(ListSpectrumPlot[x_List] := \[IndentingNewLine]ListPlot[
          ampSpectrum[x], PlotRange \[Rule] {0, 1.02}, 
          PlotStyle -> PointSize[0.018]\[IndentingNewLine]]\),
      FontSize->18]}], "Input"],

Cell[CellGroupData[{

Cell[BoxData[{
    StyleBox[\(sig[j_] := 2\ Mod[Floor[j/32.0], 2] - 1\),
      FontSize->18], "\[IndentingNewLine]", 
    StyleBox[\(x = Table[sig[j], {j, 256}];\),
      FontSize->18], "\[IndentingNewLine]", 
    StyleBox[\(ListPlot[x, PlotJoined \[Rule] True]\),
      FontSize->18], "\[IndentingNewLine]", 
    StyleBox[\(SpectrumPlot[x];\),
      FontSize->18]}], "Input"],

Cell[GraphicsData["PostScript", "\<\
%!
%%Creator: Mathematica
%%AspectRatio: .61803 
MathPictureStart
/Mabs {
Mgmatrix idtransform
Mtmatrix dtransform
} bind def
/Mabsadd { Mabs
3 -1 roll add
3 1 roll add
exch } bind def
%% Graphics
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10  scalefont  setfont
% Scaling calculations
0.0238095 0.00372024 0.309017 0.294302 [
[.20982 .29652 -6 -9 ]
[.20982 .29652 6 0 ]
[.39583 .29652 -9 -9 ]
[.39583 .29652 9 0 ]
[.58185 .29652 -9 -9 ]
[.58185 .29652 9 0 ]
[.76786 .29652 -9 -9 ]
[.76786 .29652 9 0 ]
[.95387 .29652 -9 -9 ]
[.95387 .29652 9 0 ]
[.01131 .01472 -12 -4.5 ]
[.01131 .01472 0 4.5 ]
[.01131 .16187 -24 -4.5 ]
[.01131 .16187 0 4.5 ]
[.01131 .45617 -18 -4.5 ]
[.01131 .45617 0 4.5 ]
[.01131 .60332 -6 -4.5 ]
[.01131 .60332 0 4.5 ]
[ 0 0 0 0 ]
[ 1 .61803 0 0 ]
] MathScale
% Start of Graphics
1 setlinecap
1 setlinejoin
newpath
0 g
.25 Mabswid
[ ] 0 setdash
.20982 .30902 m
.20982 .31527 L
s
[(50)] .20982 .29652 0 1 Mshowa
.39583 .30902 m
.39583 .31527 L
s
[(100)] .39583 .29652 0 1 Mshowa
.58185 .30902 m
.58185 .31527 L
s
[(150)] .58185 .29652 0 1 Mshowa
.76786 .30902 m
.76786 .31527 L
s
[(200)] .76786 .29652 0 1 Mshowa
.95387 .30902 m
.95387 .31527 L
s
[(250)] .95387 .29652 0 1 Mshowa
.125 Mabswid
.06101 .30902 m
.06101 .31277 L
s
.09821 .30902 m
.09821 .31277 L
s
.13542 .30902 m
.13542 .31277 L
s
.17262 .30902 m
.17262 .31277 L
s
.24702 .30902 m
.24702 .31277 L
s
.28423 .30902 m
.28423 .31277 L
s
.32143 .30902 m
.32143 .31277 L
s
.35863 .30902 m
.35863 .31277 L
s
.43304 .30902 m
.43304 .31277 L
s
.47024 .30902 m
.47024 .31277 L
s
.50744 .30902 m
.50744 .31277 L
s
.54464 .30902 m
.54464 .31277 L
s
.61905 .30902 m
.61905 .31277 L
s
.65625 .30902 m
.65625 .31277 L
s
.69345 .30902 m
.69345 .31277 L
s
.73065 .30902 m
.73065 .31277 L
s
.80506 .30902 m
.80506 .31277 L
s
.84226 .30902 m
.84226 .31277 L
s
.87946 .30902 m
.87946 .31277 L
s
.91667 .30902 m
.91667 .31277 L
s
.99107 .30902 m
.99107 .31277 L
s
.25 Mabswid
0 .30902 m
1 .30902 L
s
.02381 .01472 m
.03006 .01472 L
s
[(-1)] .01131 .01472 1 0 Mshowa
.02381 .16187 m
.03006 .16187 L
s
[(-0.5)] .01131 .16187 1 0 Mshowa
.02381 .45617 m
.03006 .45617 L
s
[(0.5)] .01131 .45617 1 0 Mshowa
.02381 .60332 m
.03006 .60332 L
s
[(1)] .01131 .60332 1 0 Mshowa
.125 Mabswid
.02381 .04415 m
.02756 .04415 L
s
.02381 .07358 m
.02756 .07358 L
s
.02381 .10301 m
.02756 .10301 L
s
.02381 .13244 m
.02756 .13244 L
s
.02381 .1913 m