lentickle/frequencycouplingplots.m at master · tobin/lentickle · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
% frequency coupling plots

setupLentickle;

f_numpoints = 200;
f_upperLimit = 5000;
f_lowerLimit = 20;

olcl = 'cl'; %closed loop or open loop

f = logspace(log10(f_lowerLimit),log10(f_upperLimit),f_numpoints).';

inPower = 8;

offsets = [.1 10];%[ 1 -20 -10 -7 7 10 20 ];
offsets = offsets*1e-12;

Nplot = length(offsets);

clear sensors
sensors{Nplot} = [];

for kk = 1:Nplot
    switch kk
        case 1
            sensors{kk} = 'asq';
        otherwise
            sensors{kk} = 'omc';
    end
end

clear result TF
TF{Nplot} = [];

for kk = 1:Nplot
    result(kk) = getEligoResults(f,inPower,offsets(kk),sensors{kk});
    CMclg = pickleTF(result(kk),'CM','CM',olcl);
    calFreq_Phase = 1i*f;
    switch sensors{kk}
        case 'asq'
            calASQ_DARMm = 1./((pickleTF(result(kk),'EX','AS_Q',olcl)-pickleTF(result(kk),'EY','AS_Q',olcl)));
            TF{kk} = [f,calTF(pickleTF(result(kk),'PM','AS_Q',olcl),calASQ_DARMm,CMclg.*calFreq_Phase)];
        case 'omc'
            calOMC_DARMm = 1./((pickleTF(result(kk),'EX','OMC_PD',olcl)-pickleTF(result(kk),'EY','OMC_PD',olcl)));
            TF{kk} = [f,calTF(pickleTF(result(kk),'PM','OMC_PD',olcl),calOMC_DARMm,CMclg.*calFreq_Phase)];
    end
end

figure(45)
SRSbode(TF{:})

% build legend

clear TFleg
TFleg{Nplot} = [];
for kk = 1:Nplot
    switch sensors{kk}
        case 'asq'
            TFleg{kk} = 'RF';
        case 'omc'
            TFleg{kk} = 'DC';
    end

    TFleg{kk} = [TFleg{kk} ' ' num2str(offsets(kk)*1e12) 'pm'];
end

legend(TFleg{:})
title('Frequency Noise Coupling')
ylabel('Magnitude (m/Hz)')