From 1a6d30f676605accf50e0bfb3490fb5013fbac74 Mon Sep 17 00:00:00 2001 From: Thomas Dehaeze Date: Tue, 28 Jan 2020 15:07:45 +0100 Subject: [PATCH] Update tangled file --- ground-motion/data/ground_meas_id31.zip | Bin 2958065 -> 2958304 bytes ground-motion/matlab/ground_meas_id31.m | 33 ++++++++++++++++++++++-- 2 files changed, 31 insertions(+), 2 deletions(-) diff --git a/ground-motion/data/ground_meas_id31.zip b/ground-motion/data/ground_meas_id31.zip index ad8f547961d77cbe8f819c2cc8f153f65152d3c8..cf6ac2e7280f35c7d9ef362c7a6f82a198cb74a1 100644 GIT binary patch delta 2208 zcmY+``6JVh1IO{rn0sc9921*i%jE8}q;hUW&X8M#awRvBA$RFRIdZjJF;~ge%ALYm zuH?vaD`&z|LO$=lzkPpryncE81JBJj4d9uc25?L~IG!yY63-sb5f5$P4D}^%g=?D! z8YgyK;sk--iGn}~5E67X&_Cq5yNjQPYmkeVySnOeKWi&)5Cp5E;`~=gI4=myz61e* z{x_|*erEFo_@75EuFEj9VGc16pl4fULe*A_Y*@nW01?%8UtFR()jwRxqw zG5XE^Vy|8;pzjwzqdI%#OxPU%=v*}I0jsm`>WKD*)uSX#gynNoeV|zk12-KwIKE}EFI~Bcc(MvNBthT;ODBf%*jr436 zk~@r@BM)^4iwiLBwT8*VTW2w5vIBmCWG=NnSzacU*6`Hua=;ye{sE`Y?-1CQgTel} zGAhejO7@|K!tn&L)7&{wz7mmcq2TVVk5NCYzxCH=9z3k;Efo41Wf8Dq@sQoIUz-T! z4)QR91a_f27Yo*y9E_a)LWrI;?W{Rma88?ZuY$C<@^5prB8I@^otNTY`lZmc;zR5G4x|A4xI`2MM^JH^2 zS=_$(q(+qR(%jyBE;FeY-acn<4rYTy4hY+f8NFDXG>gEvA91| z#2vy^531iawr<1Z$MA|CYf~ydqC$$}vi4({aEX6f676k-RfS^KnCEve6#otKD`lIE z!c!({k5=(}s;^nSLe=aC*4I(YaJ&U(+RJ*x)-YM%^Z7Oa6U-?WPvR`DYj`C4*$>wp)wj4kt+v&mtpLb0)F}^u>{oLmp zO}yP4XUO@wjB6#fOpnx+Ptq_*tj)IjUMpGh*fEmey`q@M+gEft=2F)B;KX-t_RbDV zd6^+sm&H$>&Id+iPXiZ}ds?3q?|ape1Rb>QNsBEPC>y<6D7h6cnLBZX96h>|*B<6u zvq5wa)^!kOjQdRe>KKHxB8_N0hS&AV;|P;8BFy|$wb(Yd2UfE#T~o!}be8`#crj8- z=eMQMNZsnHA|X3XBZ7GHi}x+14^_|T!39la^hGU>N=L0qPrXa^cEj7ci0H)kSz3-& zVt1xP`)~L9t+xmrx2i&t)hZTEZRgvX^LpjZqJ&!R)2mO}P5Nh7SboT;U`Q-8Bm$m} ziCbWyyLz7sT87js=QbbYz11X&GOJna<8a^4gf?&f6Om=~IsPbCw zIR4;eRtG^uI}W4xA@SiPyr{*Ez?_5VN#V|xOX4$%+-$Qr4^NSKZ2~u>RcGBp)=&7_ zPTb%nW?sVE1+XW~2V{wlX_{J@JldMsb|T2uAG!9!>H!IhJcUr?Nf0O!A1BACw7>A#4v-LX(A(IcH^+NscSH4XOl*ivRO zTE{lvIY>&e?~c$7A;i*L>6o3(K#Nvq&!;L|-)i!Np7rL6_k*3#FCiZGquTvXuT>bq ziuLSq0e)eghxP%O2$ZwcnefflI@7SHA1bIxWGu|v>{~nZu=1$Z?$@f7B*+O@!O{;N zYn4#3kA|wFm}3j{1k+IGL;fV0AWi@mfE&OA;05pj{(dlk000gU1RwxH0AYX#KooEU zAO;WzNB|@ONPraJC_ow@13&@L09gPAAP0~KC;+ek9N-v0u~{YvPpyUhe?sJpp_CN= xwn9z~N(q^RLMt6*69~@HMJu^O3ikwa-u)Jo$|*-FiRFAoE5V>PXeAH`^dD#n0KWhL delta 1967 zcmY+_`9IT-1HkcZ#$uTx5ewx`9NU*^z^n-F1Q^7ZQsp;#jcEeU9Xt zpYDRUMztPb)NhKt`8b37e0Pbs0JlHMRGQzlERC^vP>k+|@w?m0G=^8J2(1oWmqa~o zz}zVVGJzHss29O0Hrr?y*DOb_&Y6Z|Yfy)b`CRv%LP7So5T^uNzAmeeFxt$jDbuhT((t9H+!LDAH=^?1j@R9ExUvY-wR&6pHNp>( zeM_`1Rx+$L@5sE#$|eWTPM9k)+YEb=xb|`EYh^nBmQsT8_CyCV<=q)0RUzv)#&0KW zNV#`3TZx73asJT!;1##-w`n3~7eU`mH4s(7Uhbn@G;@jBt@gibs#E=MJN>$>D5WUx zZ`ntUmwPz&OU2r$edZG{?s}ofrX=S!f$*_Ft$BSbd~1glo4C-5o1Rwxb++zl(w zLKlOhFh<&WNsM$ILRAG@=`NKk8{btxj&~|4|D|Ta-~w?js;XbUk{R&IJI6I`Kex~y zjbamF4KJHAC{~dtmQ(L|-;{B6=$%AY%N<43SL38tw>}yjrLEXV?CrsziLP#QH|j!S zBg&rQU$AIc8}RY*T&Y%g37Vwe{15`o9LPBlK0feV!|dwRpB=XjJD1)lRj%QwCf7D> zcdm%j8CbIwYB9u_1`n>1_fP9rhlo&%jQh^det7o(@^^ZP3(<)fT!e;@LX&Aj=g1dp z?&SK-Tty~Xn4^SwmFIF0PueS3?N&o2%(R+LRR1+kC zcCrs<1!h^$#p%+};W+e<=Mo`Ep-qaRwT&$10n8W79@7d=Mf;fMo#_?`efD~qc&r#6 zr8eej&{-*i-1=9lwHj^+++N11 z0hRiqeY~uzYScp!fp*PE5Kc1Sipx&l@G9g9mt;Z>CxchnOPLAvQBeQz&dlOSFv~VZ5Py)V z{{XT-xauNJXjTtfoMVw+yRZL+)lxsTc^j@;gh&*_lv-A}MIGsi>d4OKnUJUvWINFo z?xe#O^c}fH!&ZD=UL~V*tO2A6pG&w+xuxhnL2(GUv41q;RrS5JV+&wuOj)lAXI^n% zBdx536B2s$Nyn1dsGe^n3oq>+_IU97hD&85PqDcL$PJM7(wVOR1MR$yM*etnr zBYD`QgGs2W_o{*w6dESH-V{{sPSrVLUh6dUYo68FOaeN2Ab0QkV50Yy+!bZNr9kXr zp_-zgk%rY+H(Kj$$=Q(Grl{)Yk$AG6=db*@`n|(lo#cDpbo>eoK5NN~-TZ19rYstq zpY;7}T;yjsSNFn}*C9M-K$k3dZ@0$PKUGh0r=#_jMymOL4za|tWMW}}BLGQ&6aWr5 z3XldM05SksfE+*`a15XTPy{FejsuhdDgaf08bBSO0YCyY0Vn_(fB|3uS^#YT4xj_j z1)Kor0ra2A#2Rp$5dVKv_8%C${y*zvb1-;ywk#HpgrH*B1T3B|*7RMD{o;q52D=D@ RS7!HP@dyb98V>@2{tHgTeZ>F( diff --git a/ground-motion/matlab/ground_meas_id31.m b/ground-motion/matlab/ground_meas_id31.m index 87cd5b2..e39fa04 100644 --- a/ground-motion/matlab/ground_meas_id31.m +++ b/ground-motion/matlab/ground_meas_id31.m @@ -8,7 +8,7 @@ s = zpk('s'); data = load('mat/data_028.mat', 'data'); data = data.data; -% Time domain plots +% Time domain plots of the measured voltage figure; hold on; @@ -82,6 +82,35 @@ psd_gv = px_dc./abs(squeeze(freqresp(G0*S, f, 'Hz'))).^2; psd_gm = psd_gv./(2*pi*f).^2; +% Time domain plots of the ground motion +% We can inverse the dynamics of the geophone to convert the measured voltage into the estimated ground motion. + + +est_vel = lsim(inv(G0*S)*(s/2/pi)/(1+s/2/pi), data(:, 1), data(:, 3)); % Estimated velocity above 1Hz +est_vel = est_vel - mean(est_vel(data(:,3)>10)); % The mean value of the velocity if removed +est_dis = lsim(1/(1+s/2/pi), est_vel, data(:, 3)); % The velocity is integrated above 1Hz + +figure; +hold on; +plot(data(:, 3), est_vel); +hold off; +xlabel('Time [s]'); ylabel('Velocity [m/s]'); +xlim([10, 100]); + + + +% #+NAME: fig:time_domain_velocity +% #+CAPTION: Time domain velocity ([[./figs/time_domain_velocity.png][png]], [[./figs/time_domain_velocity.pdf][pdf]]) +% [[file:figs/time_domain_velocity.png]] + + +figure; +hold on; +plot(data(:, 3), est_dis); +hold off; +xlabel('Time [s]'); ylabel('Displacement [m]'); +xlim([10, 100]); + % Computation of the ASD of the velocity and displacement % The ASD of the measured velocity is shown on figure [[fig:ground_motion_id31_asd_velocity]]. @@ -135,7 +164,7 @@ xlim([0.1, 500]); ylim([1e-13, 1e3]); % Save % We save the PSD of the ground motion for further analysis. -save('./mat/psd_gm', 'f', 'psd_gm'); +save('./mat/psd_gm.mat', 'f', 'psd_gm', 'psd_gv'); % Load the measurement data % First we load the measurement data.