0
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 1 <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 2 "http://www.w3.org/TR/1999/REC-html401-19991224/loose.dtd">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 3
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 4 <html lang="en">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 5 <head>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 6 <meta name="generator" content=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 7 "HTML Tidy for Mac OS X (vers 1st December 2004), see www.w3.org">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 8 <meta http-equiv="Content-Type" content=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 9 "text/html; charset=us-ascii">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 10
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 11 <title>IFO/Temperature Example - Pre-processing (LTPDA Toolbox)</title>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 12 <link rel="stylesheet" href="docstyle.css" type="text/css">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 13 <meta name="generator" content="DocBook XSL Stylesheets V1.52.2">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 14 <meta name="description" content=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 15 "Presents an overview of the features, system requirements, and starting the toolbox.">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 16 </head>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 17
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 18 <body>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 19 <a name="top_of_page" id="top_of_page"></a>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 20
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 21 <p style="font-size:1px;"> </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 22
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 23 <table class="nav" summary="Navigation aid" border="0" width=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 24 "100%" cellpadding="0" cellspacing="0">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 25 <tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 26 <td valign="baseline"><b>LTPDA Toolbox</b></td><td><a href="../helptoc.html">contents</a></td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 27
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 28 <td valign="baseline" align="right"><a href=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 29 "ltpda_training_topic_2_8.html"><img src="b_prev.gif" border="0" align=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 30 "bottom" alt="Split and join AOs"></a> <a href=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 31 "ltpda_training_topic_3.html"><img src="b_next.gif" border="0" align=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 32 "bottom" alt="Topic 3 - Spectral Analysis"></a></td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 33 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 34 </table>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 35
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 36 <h1 class="title"><a name="f3-12899" id="f3-12899"></a>IFO/Temperature Example - Pre-processing</h1>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 37 <hr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 38
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 39 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 40 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 41 Now we return to the IFO/Temperature example that was started in Topic 1.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 42 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 43 <h2>Loading and checking the calibrated data sets from topic2</h2>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 44 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 45 In the last topic you should have saved your calibrated data files as<br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 46 <tt>ifo_temp_example/ifo_disp.xml</tt> and <br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 47 <tt>ifo_temp_example/temp_kelvin.xml</tt>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 48 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 49 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 50 Now load each file into an AO, and simplify the name of the objects by assigning them the name of the Matlab workspace variable containing them:<br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 51 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 52 <div class="fragment"><pre>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 53 ifo = ao(<span class="string">'ifo_temp_example/ifo_disp.xml'</span>);
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 54 temp = ao(<span class="string">'ifo_temp_example/temp_kelvin.xml'</span>);
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 55 ifo.setName;
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 56 temp.setName;
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 57 </pre></div>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 58 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 59 Let's see what kind of pre-processing we have to apply to our data prior to further analysis.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 60 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 61 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 62 You can have a look at the data by for example displaying the AOs on the terminal and by
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 63 plotting them, of course. Since the two data series have different Y units, we should plot
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 64 them on subplots. To do that with <tt>iplot</tt>, pass the key 'arrangement' in a plist. For
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 65 example:
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 66 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 67 <div class="fragment"><pre>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 68 pl = plist(<span class="string">'arrangement'</span>, <span class="string">'subplots'</span>);
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 69 </pre></div>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 70 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 71 If you plot the two time-series you should see something like the following:<br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 72 <img src="images/ltpda_training_1/topic2/ifotempraw.png" alt="ifotempdata" border="1">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 73 <br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 74 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 75 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 76 Some points to note:
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 77 <ol>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 78 <li>The two data streams:
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 79 <ul>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 80 <li> do not have the same sampling frequency.</li>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 81 <li> are not of the same length (nsecs).</li>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 82 </ul></li>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 83 <li>The interferometer data has a small transient at the start</li>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 84 </ol>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 85 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 86 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 87 To have a closer look at the samples we plot markers at each sample and only plot a zoomed-in section. You can
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 88 do this by creating a parameter list with the following key/value pairs (remember the 'key' properties for <tt>plist</tt> entries is not case sensitive):
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 89 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 90 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 91 <table cellspacing="0" class="body" cellpadding="2" border="0" width="50%">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 92 <colgroup>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 93 <col width="35%"/>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 94 <col width="65%"/>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 95 </colgroup>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 96 <thead>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 97 <tr valign="top">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 98 <th class="categorylist">Key</th>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 99 <th class="categorylist">Value</th>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 100 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 101 </thead>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 102 <tbody>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 103 <!-- Key 'arrangement' -->
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 104 <tr valign="top">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 105 <td bgcolor="#f3f4f5">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 106 <p><tt>ARRANGEMENT</tt></p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 107 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 108 <td bgcolor="#f3f4f5">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 109 <span class="string">'subplots'</span>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 110 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 111 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 112 <!-- Key 'arrangement' -->
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 113 <tr valign="top">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 114 <td bgcolor="#f3f4f5">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 115 <p><tt>LINESTYLES</tt></p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 116 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 117 <td bgcolor="#f3f4f5">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 118 {<span class="string">'none'</span>,<span class="string">'none'</span>}
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 119 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 120 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 121 <!-- Key 'markers' -->
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 122 <tr valign="top">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 123 <td bgcolor="#f3f4f5">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 124 <p><tt>MARKERS</tt></p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 125 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 126 <td bgcolor="#f3f4f5">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 127 {<span class="string">'+'</span>,<span class="string">'+'</span>}
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 128 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 129 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 130 <!-- Key 'xranges' -->
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 131 <tr valign="top">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 132 <td bgcolor="#f3f4f5">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 133 <p><tt>XRANGES</tt></p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 134 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 135 <td bgcolor="#f3f4f5">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 136 {<span class="string">'all'</span>, [200 210]}
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 137 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 138 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 139 <!-- Key 'yranges' -->
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 140 <tr valign="top">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 141 <td bgcolor="#f3f4f5">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 142 <p><tt>YRANGES</tt></p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 143 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 144 <td bgcolor="#f3f4f5">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 145 {[2e-7 3e-7], [200 350]}
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 146 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 147 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 148 </tbody>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 149 </table>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 150 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 151 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 152 <table cellspacing="0" class="note" summary="Note" cellpadding="5" border="1">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 153 <tr width="90%">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 154 <td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 155 Notice the use of the keyword 'all' in the value for the 'XRANGES' parameter. Many
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 156 of the <tt>iplot</tt> options support this keyword, which tells <tt>iplot</tt> to use
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 157 the same value for all plots and subplots. For the 'YRANGES' we specify different values
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 158 for each subplot.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 159 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 160 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 161 </table>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 162 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 163 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 164 <table cellspacing="0" class="note" summary="Note" cellpadding="5" border="1">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 165 <tr width="90%">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 166 <td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 167 Please store your parameter lists in 2 different variables. We can reuse them for
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 168 plotting our results later. If you are working on a pipeline insteaad of a script, you
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 169 can use two <tt>plist</tt> constructor blocks and pass these as an input to an
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 170 <tt>iplot</tt> block.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 171 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 172 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 173 </table>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 174 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 175 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 176 Passing such a parameter list to <tt>iplot</tt> together with the two AO time-series should yield
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 177 a plot something like:
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 178 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 179 <br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 180 <img src="images/ltpda_training_1/topic2/samples.png" alt="samples" border="1">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 181 <br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 182 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 183 From this plot you may be able to see that the temperature data is <bb>unevenly sampled</bb>.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 184 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 185 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 186 To confirm this, enter the following (standard) MATLAB commands on the terminal:
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 187 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 188 <div class="fragment"><pre>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 189 dt = diff(temp.x);
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 190 min(dt)
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 191 max(dt)
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 192 </pre></div>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 193 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 194 <table cellspacing="0" class="note" summary="Note" cellpadding="5" border="1">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 195 <tr width="90%">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 196 <td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 197 Don't forget the semicolon at the end of the <tt>diff</tt> calculation; this is a long data series
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 198 and will be printed to the terminal if you do forget.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 199 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 200 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 201 </table>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 202 You see that the minimum and maximum difference in the time-stamps of the data is different, showing
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 203 that the data are not evenly sampled.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 204 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 205 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 206 Before we proceed with the later analysis of this data, we need to
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 207 <ul>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 208 <li>Fix the uneven sampling of the temperature data</li>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 209 <li>Resample both data streams to the same rate</li>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 210 <li>Resample both data streams on to the same timing grid</li>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 211 <li>Select the segment of interferometer data that matches the temperature data</li>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 212 </ul>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 213 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 214 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 215 Each of these steps can, in principle, be done by hand. However, LTPDA provides
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 216 a 'data fixer' method called <tt>ao/consolidate</tt> which attempts to automate this
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 217 process. The call to <tt>consolidate</tt> is shown below:
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 218 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 219 <div class="fragment"><pre>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 220 [temp_fixed ifo_fixed] = consolidate(temp, ifo, plist(<span class="string">'fs'</span>,1));
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 221 </pre></div>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 222 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 223 We tell <tt>consolidate</tt> that we want to have our data resampled to 1 Hz by specifying the
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 224 parameter key 'fs'.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 225 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 226 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 227 Now we can inspect the time-series of these data. The result should look something like
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 228 the figure below:
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 229 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 230 <br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 231 <img src="images/ltpda_training_1/topic2/ifo_temp_consolidated.png" alt="consolidated" border="1">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 232 <br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 233 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 234 <table cellspacing="0" class="note" summary="Note" cellpadding="5" border="1">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 235 <tr width="90%">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 236 <td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 237 Note that the time origin above the plots has now changed from zero to 13.105 which
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 238 was the time of the first sample in the temperature measurement.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 239 </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 240 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 241 </table>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 242 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 243 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 244 If we also plot the zoomed-in view again, we should see something like:
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 245 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 246 <img src="images/ltpda_training_1/topic2/samples_conso.png" alt="consolidated samples" border="1">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 247 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 248 As you can see <tt>consolidate</tt> solved all our issues with these two data streams.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 249 They now start at the same time and are evenly sampled at the same sampling frequency.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 250 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 251 <p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 252 In the next topic, we will look at the
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 253 spectral content and coherence of the data before and after the pre-processing. For now,
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 254 finish by saving the consolidated data ready for the next topic.
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 255 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 256 <div class="fragment"><pre>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 257 save(temp_fixed,<span class="string">'ifo_temp_example/temp_fixed.xml'</span>);
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 258 save(ifo_fixed,<span class="string">'ifo_temp_example/ifo_fixed.xml'</span>);
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 259 </pre></div>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 260
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 261
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 262
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 263 </p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 264
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 265 <br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 266 <br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 267 <table class="nav" summary="Navigation aid" border="0" width=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 268 "100%" cellpadding="0" cellspacing="0">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 269 <tr valign="top">
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 270 <td align="left" width="20"><a href="ltpda_training_topic_2_8.html"><img src=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 271 "b_prev.gif" border="0" align="bottom" alt=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 272 "Split and join AOs"></a> </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 273
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 274 <td align="left">Split and join AOs</td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 275
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 276 <td> </td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 277
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 278 <td align="right">Topic 3 - Spectral Analysis</td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 279
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 280 <td align="right" width="20"><a href=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 281 "ltpda_training_topic_3.html"><img src="b_next.gif" border="0" align=
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 282 "bottom" alt="Topic 3 - Spectral Analysis"></a></td>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 283 </tr>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 284 </table><br>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 285
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 286 <p class="copy">©LTP Team</p>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 287 </body>
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
+ − 288 </html>