40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
  40m Log  Not logged in ELOG logo
Entry  Fri Jul 17 04:56:37 2015, Ignacio, Update, General, MCL Wiener filtering + FIR to IIR conversion using vectfit 9x
    Reply  Sat Jul 18 06:12:07 2015, Ignacio, Update, General, MCL Wiener filtering + FIR to IIR conversion using vectfit (Update) 15x
Message ID: 11425     Entry time: Sat Jul 18 06:12:07 2015     In reply to: 11424
Author: Ignacio 
Type: Update 
Category: General 
Subject: MCL Wiener filtering + FIR to IIR conversion using vectfit (Update) 

(updateAfter Eric gave me feedback on my previous elog post, I went back and fixed some of the silly stuff I stated.

First of all, I have come to realized that it makes zero sense to plot the ASD's of the mode cleaner against the seismometer noise. These measurements are not only quite different, but elementary, they posess different units. I have focused my attention to the MCL being Wiener filtered with the three siesmometer signals. 

One of the major improvements that I make in the following analysis is,

1) Prefiltering; a band pass filter from 1 to 5 Hz, in order to emphasize subtraction of the bump shown in the figure below.

2) I have used vectfit exclusively in the 1 to ~5 Hz range, in order to model the FIR filter properly, as in, the kind of subtraction that we care about. Limiting myself to the 1 - 5 Hz range has allowed me to play freeley with the number of poles, hence being able to fit the FIIR filter properly with an IIR rational transfer function properly,

The resulting ASD's are shown below, in blue we show the raw MCL output, in blac the Wiener filter (FIR) result, and finally in black, the resultant data being filtered with the calculated IIR Wiener filter.


Now, in the following plots I show the IIR Wiener filters for each of the three seismometers, 

X Seismometer,

For the Y seismometer,

and for the Z seismometer,


The matlab code for this work is attached: code.zip

Attachment 1: Wiener_MCL_seismometers_iir.png  233 kB  | Show | Show all
Attachment 2: seisx_mag.png  64 kB  | Show | Show all
Attachment 3: seisx_mag.png  64 kB  | Show | Show all
Attachment 4: seisx_mag.png  64 kB  | Show | Show all
Attachment 5: seisx_ph.png  53 kB  | Show | Show all
Attachment 6: seisy_mag.png  65 kB  | Show | Show all
Attachment 7: seisy_mag.png  65 kB  | Show | Show all
Attachment 8: seisy_mag.png  65 kB  | Show | Show all
Attachment 9: seisy_ph.png  50 kB  | Show | Show all
Attachment 10: seisz_ph.png  65 kB  | Show | Show all
Attachment 11: seisz_ph.png  65 kB  | Show | Show all
Attachment 12: code.zip  19.560 MB  Uploaded Sat Jul 18 07:16:25 2015
Attachment 13: seisz_mag.png  78 kB  Uploaded Sun Jul 19 17:27:04 2015  | Show | Show all
Attachment 14: seisz_mag.png  78 kB  Uploaded Sun Jul 19 17:28:19 2015  | Show | Show all
Attachment 15: seisz_ph.png  65 kB  Uploaded Sun Jul 19 17:28:29 2015  | Show | Show all
ELOG V3.1.3-