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 Sep 23 17:57:34 2011, Paul, Update, SUS, Oplev filter optimization for 2 poles and 2 zeros optimised2p2z_v1.pngoptimised2p2z_v1_TFs.png
    Reply  Sat Sep 24 01:51:02 2011, rana, Update, SUS, Oplev filter optimization for 2 poles and 2 zeros 
       Reply  Tue Sep 27 08:51:29 2011, Paul, Update, SUS, Oplev filter optimization for 2 poles and 2 zeros 
    Reply  Sat Sep 24 02:09:43 2011, kiwamu, Update, SUS, Re:Oplev filter optimization for 2 poles and 2 zeros 
       Reply  Sat Sep 24 17:45:56 2011, Paul, Update, SUS, Re:Oplev filter optimization for 2 poles and 2 zeros 
Message ID: 5554     Entry time: Tue Sep 27 08:51:29 2011     In reply to: 5536
Author: Paul 
Type: Update 
Category: SUS 
Subject: Oplev filter optimization for 2 poles and 2 zeros 

Quote:

Quote:

I have made a function to optimise the overall gain, pole frequencies and zero frequencies for the oplev filter. The script will optimize any user defined number of poles and zeros in order to minimise the RMS motion below a certain cut off frequency (in this case 20Hz). The overall gain is adjusted so that each trial filter shape always has a UGF of 10 Hz.

I think this is a nice start. Its clear that we don't want to use this feedback law, but the technique can be tweaked to do what we want by just tweaking our cost function.

Let's move the scripts into the SUS/ scripts area and then start putting in weights that do what we want:

1) Limit the gain peaking at the upper UGF to 6 dB.

2) Minimum phase margin of 45 deg.

3) Minimum gain margin of 10 dB.

4) Lower UGF = 0.1 Hz / Upper UGF = 10 Hz.

5) Assume a A2L coupling of 0.003 m/rad and constrain the injected noise at the test mass to be less than the seismic + thermal level.

6) Looser noise contraint above 50 Hz for the non TM loops.

 I moved two matlab scripts into the folder /cvs/cds/rtcds/caltech/c1/scripts/SUS/Oplev_filter_optimization

These are the function 'filter_optimiser_zeros_and_poles.m', and the example script to run the function 'run_filter_optimiser.m'. Type 'help filter_optimiser_zeros_and_poles.m' to get details about the function.

I haven't implemented the new weights yet. I've pasted them into the the file header to remind me/us of the work to be done on the function.

ELOG V3.1.3-