Written by Yuta
Now that the mystery of the PD signal chain and the BS actuation efficiency for Michelson are solved (see alog #10213 and #10127), I corrected what we've found missing/wrong in the noise budget (NB) model.
[Motivation]
We wanted to check the validity of the NB model.
[What are missing/wrong?]
Things that were missing/wrong in the NB model "DRMI_Live.slx" copied from LLO was;
[What files do I use?]
The NB model and functions for our Michelson lives in /ligo/svncommon/NbSVN/aligonoisebudget/trunk/MICH/H1
. They are based on LLO DRMI NB model but corrected the things mentioned above.
The essential files are
run_NB.m
: main script
make_MICH.m
: function called from run_NB.m
./Params/paramNbH1MICHtest.m
: parameter file
MICH_Live.slx
: simulink model
They use the following Optickle files which live in /ligo/svncommon/IscCVS/iscmodeling/LentickleAligo/PRMI/ParamFiles
for optical simulation.
optL1DRM.m
: constructs Optickle model
probesH1DRM_00.m
: puts probes for the Optickle model
paramH1MICH.m
: parameter file
[Result]
Attached plot shows the comparison of the OLTF of the MICH loop from the measurement (see alog #10127) and from the NB model. They agree within ~20%. This error mainly comes from the power measurement error.
Note that the Matlab function "linmod" doesn't work correctly when the Optickle block is put in the NB simulink model. Do something like;
[systm,flexTfs] = linFlexTf('MICH_Live');
systm = prescale(systm,{freq(1) freq(end)});
systm = linFlexTfFold(systm,flexTfs);
to plot transfer functions (Thanks to Chris Wipf!).
Thank you for the excellent work on validating this model!