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: 11424     Entry time: Fri Jul 17 04:56:37 2015     Reply to this: 11425
Author: Ignacio 
Type: Update 
Category: General 
Subject: MCL Wiener filtering + FIR to IIR conversion using vectfit 

 

We took data for the mode cleaner a while ago, June 30th I believe. This data contained signals from the six accelerometers and the three seismometers. In here I have only focused on the seimometer signals as witnesses in order to construct Wiener filters for each of the three seismometer signals (x,y,z) and for the combined seismometers signal. The following plot showing the ASD's shows the results,

 Wiener filtering works beautifully for the seismometers. Note that subtraction is best when we use all three seismometers as the witnesses in the Wiener filter calculation, as can be clearly seen in the first plot above.

Now, I used vectfit to conver the Wiener FIR filters for each seismometer to their IIR versions. The following are the bode plots for the IIR filters,

For the x-direction seismometer,

 

For the y-direction seismometer

,

 

And for the z-direction seismometer,

 The IIR filters were computed using 5 zeros and 5 poles using vectfit. That was the maximum number of poles that I could use wihtout running into trouble with matrices being almost singular in Matlab. I still need to figure out how to deal with this issue in more detail as fitting the y-seismometer was a bit problematic. I think having a greater number of poles will make the fitting a bit easier.

Attachment 1: Wiener_MCL_seismometers.png  282 kB  | Hide | Hide all
Wiener_MCL_seismometers.png
Attachment 2: seisx_mag.png  40 kB  | Hide | Hide all
seisx_mag.png
Attachment 3: seisx_mag.png  40 kB  | Hide | Hide all
seisx_mag.png
Attachment 4: seisx_mag.png  40 kB  | Hide | Hide all
seisx_mag.png
Attachment 5: seisx_phase.png  34 kB  | Hide | Hide all
seisx_phase.png
Attachment 6: seisy_mag.png  42 kB  | Hide | Hide all
seisy_mag.png
Attachment 7: seisy_phase.png  36 kB  | Hide | Hide all
seisy_phase.png
Attachment 8: seisz_mag.png  40 kB  | Hide | Hide all
seisz_mag.png
Attachment 9: seisz_phase.png  33 kB  | Hide | Hide all
seisz_phase.png
ELOG V3.1.3-