I tried the LSC model (tried attaching but does not accept .mat files) with the H2 solver And it worked like normally. You can see these results in the "LSC_H2_Results.pdf" attachment.
When trying the LSC model with the H2/Hinf Mixed solver it struggled to recover even the H2 results. This is with a gamma going from 1e3 down to 1e-1 (not igsq). On the H2 result it looks like the Closed loop is peaking at about 100 so this should be plenty to at least recover the H2 results. but I am left with the results seen in "LSC_H2HinfMix_Results.pdf"
It is possible that this is a problem with the Zinf connection but I will have to investigate further. I also still need to try the damping loop. Also I got multiple color bars working using This.
I am running into a problem choosing the color scheme for the plots for the paper. There are two main challenges to this first is that in the mixed H2/Hinf case the color conveys two pieces of info: the FOM weighting and the gamma limit on the Hinf norm. I have attached a number of single color schemes in the "PaperALS_OL_BODE_colors_comp.pdf" attachment. I don't know how to handle the gamma limit. Ideally there would be a 2d color space to map like in the pycolormap-2d package. Is it worth it?
The other main plot is the RMS plot which also has two color schemes one for the gain margin and one for the FOM weighting. This is easier as long as the colors don't overlap. Also I haven't gotten the dual color bar to work on the RMS Plot. the scaling always is removed by adding the second color bar. See the example RMS plot in "PaperALS_FOM_RMS_Scatter.pdf"
Note: some of the colorbars are messed up in "PaperALS_OL_BODE_colors_comp.pdf"
Remeasured the frequency noise with the 785 after the visibility improvement. Assuming my calibration to Hz/sqrt(hz) is correct, the measurement is still yielding several orders of magnitude larger than what rio claims it to be. Investigating this today.
For reference, the SR785 measures a quantity is V/sqrt(Hz). To convert to frequency noise I do S / (L/c * B*2*pi) where S is the measurement, L is the path length difference between the two arms and B is the fringe amplitude. V/sqrt(Hz) / (m/(m/s) * V) = Hz/sqrt(Hz).
The right axis (Hz) on this plot has a different frequency dependence than the left axis (Rad). So be careful with the specification data. I agree otherwise that there seems to be an excess. It might be intensity noise that needs to be canceled.
Lee came into the lab today and discovered that the TEC controller's range on the integrator is very small, which is why I was having trouble driving the laser frequency. If you center on this range (~10.2 kOhm) you can then drive a signal into the back of the TEC so that there are constant fringes on the Mach-Zehnder read out. I did this via the moku (15 mV @ .1 Hz, corresponds to roughly 1 fringe). "2.png" is this first attempt, red is the readout PD, blue is the driving signal. "3.png" is then hooking up a PD farther upstream to use as a power reference instead of the driving signal (see also cryolab_frequency_noise_measurement_pdf.pdf, the lighter color laser path in that is used). "4.png" is then using the moku's math channels to do (readout PD)/(power reference) to give the fringes at a constant offset. At this point its very easy to align to maximize these fringes. "7.png" shows the improved visibility. "5.png" shows the moku set up used. There are some shenanigans with the input ranges when using an oscilloscope in multi-instrument mode where you either have it on 400mV setting and its saturated, or its on 4V setting and it auto attenuates the signal by 20dB. So the digital filter box is undoing that 20dB attenuation. For unequal powers, I believe the formula is V =(2*sqrt(34*99)/(99+34))=87%, so not perfect but a vast improvement.
I drag and drop the link at the bottom of the post during editing and it should put it inline in the text. - looks like you have to be in the ckeditor for that to work well.