tdehaeze/nass-micro-station-measurements
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Update tfestimate command

Browse Source
This commit is contained in:
Thomas Dehaeze 2019-04-18 09:37:33 +02:00
parent 030dae7ab1
commit b3d626cd0c
3 changed files with 34 additions and 34 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

BIN
huddle-test-geophones/figs/tf_geophones.png
View File

Binary file not shown.
Side by Side Swipe Overlay

Before

Width:  |  Height:  |  Size: 124 KiB

After

Width:  |  Height:  |  Size: 132 KiB

64
huddle-test-geophones/index.html
View File

@ -3,7 +3,7 @@
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" lang="en" xml:lang="en">
<head>
<!-- 2019-04-18 jeu. 09:20 -->
<!-- 2019-04-18 jeu. 09:37 -->
<meta http-equiv="Content-Type" content="text/html;charset=utf-8" />
<meta name="viewport" content="width=device-width, initial-scale=1" />
<title>SpeedGoat</title>
@ -254,22 +254,22 @@ for the JavaScript code in this tag.
<h2>Table of Contents</h2>
<div id="text-table-of-contents">
<ul>
<li><a href="#org4a1a539">1. Setup</a></li>
<li><a href="#org474b62b">2. Signal Processing</a>
<li><a href="#org1a9211f">1. Setup</a></li>
<li><a href="#org85eb6b7">2. Signal Processing</a>
<ul>
<li><a href="#org336ec3e">2.1. Load data</a></li>
<li><a href="#org2cf9057">2.2. Time Domain Data</a></li>
<li><a href="#orgb221001">2.3. Compute PSD</a></li>
<li><a href="#org655ba9d">2.4. Take into account sensibility of Geophone</a></li>
<li><a href="#orgd36c537">2.5. Transfer function between the two geophones</a></li>
<li><a href="#orgc69f0d0">2.1. Load data</a></li>
<li><a href="#org8e9f3dd">2.2. Time Domain Data</a></li>
<li><a href="#orgade2d36">2.3. Compute PSD</a></li>
<li><a href="#orgbee0cb3">2.4. Take into account sensibility of Geophone</a></li>
<li><a href="#org92a562d">2.5. Transfer function between the two geophones</a></li>
</ul>
</li>
</ul>
</div>
</div>
<div id="outline-container-org4a1a539" class="outline-2">
<h2 id="org4a1a539"><span class="section-number-2">1</span> Setup</h2>
<div id="outline-container-org1a9211f" class="outline-2">
<h2 id="org1a9211f"><span class="section-number-2">1</span> Setup</h2>
<div class="outline-text-2" id="text-1">
<p>
Two L22 geophones are used.
@ -283,14 +283,14 @@ The voltage amplifiers include a low pass filter with a cut-off frequency at 1kH
</p>
<div id="orga919daa" class="figure">
<div id="orgb098801" class="figure">
<p><img src="./figs/setup.jpg" alt="setup.jpg" width="500px" />
</p>
<p><span class="figure-number">Figure 1: </span>Setup</p>
</div>
<div id="orgb58846e" class="figure">
<div id="org9c92645" class="figure">
<p><img src="./figs/geophones.jpg" alt="geophones.jpg" width="500px" />
</p>
<p><span class="figure-number">Figure 2: </span>Geophones</p>
@ -298,12 +298,12 @@ The voltage amplifiers include a low pass filter with a cut-off frequency at 1kH
</div>
</div>
<div id="outline-container-org474b62b" class="outline-2">
<h2 id="org474b62b"><span class="section-number-2">2</span> Signal Processing</h2>
<div id="outline-container-org85eb6b7" class="outline-2">
<h2 id="org85eb6b7"><span class="section-number-2">2</span> Signal Processing</h2>
<div class="outline-text-2" id="text-2">
</div>
<div id="outline-container-org336ec3e" class="outline-3">
<h3 id="org336ec3e"><span class="section-number-3">2.1</span> Load data</h3>
<div id="outline-container-orgc69f0d0" class="outline-3">
<h3 id="orgc69f0d0"><span class="section-number-3">2.1</span> Load data</h3>
<div class="outline-text-3" id="text-2-1">
<div class="org-src-container">
<pre class="src src-matlab">load<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-string">'mat/data_001.mat', 't', 'x1', 'x2'</span><span class="org-rainbow-delimiters-depth-1">)</span>;
@ -313,8 +313,8 @@ dt = t<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-high
</div>
</div>
<div id="outline-container-org2cf9057" class="outline-3">
<h3 id="org2cf9057"><span class="section-number-3">2.2</span> Time Domain Data</h3>
<div id="outline-container-org8e9f3dd" class="outline-3">
<h3 id="org8e9f3dd"><span class="section-number-3">2.2</span> Time Domain Data</h3>
<div class="outline-text-3" id="text-2-2">
<div class="org-src-container">
<pre class="src src-matlab"><span class="org-type">figure</span>;
@ -329,7 +329,7 @@ xlim<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-rainbo
</div>
<div id="orgc7e063b" class="figure">
<div id="org8d315fb" class="figure">
<p><img src="figs/data_time_domain.png" alt="data_time_domain.png" />
</p>
<p><span class="figure-number">Figure 3: </span>Time domain Data</p>
@ -349,7 +349,7 @@ xlim<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-rainbo
</div>
<div id="orgf38ff94" class="figure">
<div id="org15ef996" class="figure">
<p><img src="figs/data_time_domain_zoom.png" alt="data_time_domain_zoom.png" />
</p>
<p><span class="figure-number">Figure 4: </span>Time domain Data - Zoom</p>
@ -357,8 +357,8 @@ xlim<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-rainbo
</div>
</div>
<div id="outline-container-orgb221001" class="outline-3">
<h3 id="orgb221001"><span class="section-number-3">2.3</span> Compute PSD</h3>
<div id="outline-container-orgade2d36" class="outline-3">
<h3 id="orgade2d36"><span class="section-number-3">2.3</span> Compute PSD</h3>
<div class="outline-text-3" id="text-2-3">
<div class="org-src-container">
<pre class="src src-matlab"><span class="org-rainbow-delimiters-depth-1">[</span>pxx1, f1<span class="org-rainbow-delimiters-depth-1">]</span> = pwelch<span class="org-rainbow-delimiters-depth-1">(</span>x1, hanning<span class="org-rainbow-delimiters-depth-2">(</span>ceil<span class="org-rainbow-delimiters-depth-3">(</span><span class="org-highlight-numbers-number">1</span><span class="org-type">/</span>dt<span class="org-rainbow-delimiters-depth-3">)</span><span class="org-rainbow-delimiters-depth-2">)</span>, <span class="org-highlight-numbers-number">0</span>, <span class="org-rainbow-delimiters-depth-2">[]</span>, <span class="org-highlight-numbers-number">1</span><span class="org-type">/</span>dt<span class="org-rainbow-delimiters-depth-1">)</span>;
@ -368,8 +368,8 @@ xlim<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-rainbo
</div>
</div>
<div id="outline-container-org655ba9d" class="outline-3">
<h3 id="org655ba9d"><span class="section-number-3">2.4</span> Take into account sensibility of Geophone</h3>
<div id="outline-container-orgbee0cb3" class="outline-3">
<h3 id="orgbee0cb3"><span class="section-number-3">2.4</span> Take into account sensibility of Geophone</h3>
<div class="outline-text-3" id="text-2-4">
<p>
The Geophone used are L22.
@ -389,7 +389,7 @@ ylabel<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-stri
</div>
<div id="org8d58394" class="figure">
<div id="org434e518" class="figure">
<p><img src="figs/geophone_sensibility.png" alt="geophone_sensibility.png" />
</p>
<p><span class="figure-number">Figure 5: </span>Sensibility of the Geophone</p>
@ -427,18 +427,18 @@ xlim<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-rainbo
</div>
<div id="org9606210" class="figure">
<div id="org16b2e10" class="figure">
<p><img src="figs/psd_velocity.png" alt="psd_velocity.png" />
</p>
<p><span class="figure-number">Figure 6: </span>Spectral density of the velocity</p>
</div>
</div>
</div>
<div id="outline-container-orgd36c537" class="outline-3">
<h3 id="orgd36c537"><span class="section-number-3">2.5</span> Transfer function between the two geophones</h3>
<div id="outline-container-org92a562d" class="outline-3">
<h3 id="org92a562d"><span class="section-number-3">2.5</span> Transfer function between the two geophones</h3>
<div class="outline-text-3" id="text-2-5">
<div class="org-src-container">
<pre class="src src-matlab"><span class="org-rainbow-delimiters-depth-1">[</span>T12, f12<span class="org-rainbow-delimiters-depth-1">]</span> = tfestimate<span class="org-rainbow-delimiters-depth-1">(</span>x1, x2, hanning<span class="org-rainbow-delimiters-depth-2">(</span><span class="org-highlight-numbers-number">1</span><span class="org-type">/</span>dt<span class="org-rainbow-delimiters-depth-2">)</span>, <span class="org-highlight-numbers-number">0</span>, <span class="org-rainbow-delimiters-depth-2">[]</span>, <span class="org-highlight-numbers-number">1</span><span class="org-type">/</span>dt<span class="org-rainbow-delimiters-depth-1">)</span>;
<pre class="src src-matlab"><span class="org-rainbow-delimiters-depth-1">[</span>T12, f12<span class="org-rainbow-delimiters-depth-1">]</span> = tfestimate<span class="org-rainbow-delimiters-depth-1">(</span>x1, x2, hanning<span class="org-rainbow-delimiters-depth-2">(</span>ceil<span class="org-rainbow-delimiters-depth-3">(</span>length<span class="org-rainbow-delimiters-depth-4">(</span>x1<span class="org-rainbow-delimiters-depth-4">)</span><span class="org-type">/</span><span class="org-highlight-numbers-number">100</span><span class="org-rainbow-delimiters-depth-3">)</span><span class="org-rainbow-delimiters-depth-2">)</span>, <span class="org-rainbow-delimiters-depth-2">[]</span>, <span class="org-rainbow-delimiters-depth-2">[]</span>, <span class="org-highlight-numbers-number">1</span><span class="org-type">/</span>dt<span class="org-rainbow-delimiters-depth-1">)</span>;
</pre>
</div>
@ -458,12 +458,12 @@ yticks<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-rain
xlabel<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-string">'Frequency </span><span class="org-string"><span class="org-rainbow-delimiters-depth-2">[</span></span><span class="org-string">Hz</span><span class="org-string"><span class="org-rainbow-delimiters-depth-2">]</span></span><span class="org-string">'</span><span class="org-string"><span class="org-rainbow-delimiters-depth-1">)</span></span><span class="org-string">; ylabel</span><span class="org-string"><span class="org-rainbow-delimiters-depth-1">(</span></span><span class="org-string">'Phase'</span><span class="org-rainbow-delimiters-depth-1">)</span>;
linkaxes<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-rainbow-delimiters-depth-2">[</span>ax1,ax2<span class="org-rainbow-delimiters-depth-2">]</span>,<span class="org-string">'x'</span><span class="org-rainbow-delimiters-depth-1">)</span>;
xlim<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-rainbow-delimiters-depth-2">[</span><span class="org-highlight-numbers-number">2</span>, <span class="org-highlight-numbers-number">500</span><span class="org-rainbow-delimiters-depth-2">]</span><span class="org-rainbow-delimiters-depth-1">)</span>;
xlim<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-rainbow-delimiters-depth-2">[</span><span class="org-highlight-numbers-number">1</span>, <span class="org-highlight-numbers-number">500</span><span class="org-rainbow-delimiters-depth-2">]</span><span class="org-rainbow-delimiters-depth-1">)</span>;
</pre>
</div>
<div id="org9571503" class="figure">
<div id="org628446d" class="figure">
<p><img src="figs/tf_geophones.png" alt="tf_geophones.png" />
</p>
<p><span class="figure-number">Figure 7: </span>Estimated transfer function between the two geophones</p>
@ -474,7 +474,7 @@ xlim<span class="org-rainbow-delimiters-depth-1">(</span><span class="org-rainbo
</div>
<div id="postamble" class="status">
<p class="author">Author: Thomas Dehaeze</p>
<p class="date">Created: 2019-04-18 jeu. 09:20</p>
<p class="date">Created: 2019-04-18 jeu. 09:37</p>
<p class="validation"><a href="http://validator.w3.org/check?uri=referer">Validate</a></p>
</div>
</body>

4
huddle-test-geophones/index.org
View File

@ -162,7 +162,7 @@ The cut-off frequency is set at 1kHz.
[[file:figs/psd_velocity.png]]
** Transfer function between the two geophones
#+begin_src matlab :results none
[T12, f12] = tfestimate(x1, x2, hanning(1/dt), 0, [], 1/dt);
[T12, f12] = tfestimate(x1, x2, hanning(ceil(length(x1)/100)), [], [], 1/dt);
#+end_src
#+begin_src matlab :results none
@ -181,7 +181,7 @@ The cut-off frequency is set at 1kHz.
xlabel('Frequency [Hz]'); ylabel('Phase');
linkaxes([ax1,ax2],'x');
xlim([2, 500]);
xlim([1, 500]);
#+end_src
#+NAME: fig:tf_geophones