1 | function run13decembre10 |
---|
2 | |
---|
3 | %% DEFINITION DES VALEURS DE COURANT S |
---|
4 | |
---|
5 | % alphaby10 sur optique MAHER lin 1 - optimisation de l'ouverture dynamique |
---|
6 | IS= [24.48 -71.77 -217.19 256.27 -133.51 144.71 -117.02 133.12 -204.43 231.92 24.48]; |
---|
7 | |
---|
8 | %% ACTION |
---|
9 | |
---|
10 | % appliquer_Q(IQ10) |
---|
11 | % appliquer_S(IS10) |
---|
12 | |
---|
13 | %% function sauvegarde |
---|
14 | function IQ_depart = sauvegarde_Q |
---|
15 | for iQ = 1:10 |
---|
16 | Name = ['Q' num2str(iQ)]; |
---|
17 | A = getam((Name),'Model'); |
---|
18 | IQ(iQ) = A(1); |
---|
19 | end |
---|
20 | |
---|
21 | IQ_depart = IQ; |
---|
22 | |
---|
23 | function IS_depart = sauvegarde_S |
---|
24 | for iS = 1:11 |
---|
25 | Name = ['S' num2str(iS)]; |
---|
26 | A = getam((Name),'Model'); |
---|
27 | IS(iS) = A(1); |
---|
28 | end |
---|
29 | |
---|
30 | IS_depart = IS; |
---|
31 | |
---|
32 | %% fonction appliquer |
---|
33 | function appliquer_Q(IQ) |
---|
34 | for iQ = 1:10 |
---|
35 | iQ |
---|
36 | Name = ['Q' num2str(iQ)]; |
---|
37 | A = IQ(iQ); |
---|
38 | %s = getfamilydata((Name)); |
---|
39 | %B = A * ones(length(s.DeviceList),1); |
---|
40 | setsp((Name),A); |
---|
41 | end |
---|
42 | |
---|
43 | function appliquer_S(IS) |
---|
44 | for iS = 1:11 |
---|
45 | iS |
---|
46 | Name = ['S' num2str(iS)]; |
---|
47 | A = IS(iS); |
---|
48 | %s = getfamilydata((Name)); |
---|
49 | %B = A * ones(length(s.DeviceList),1); |
---|
50 | setsp((Name),A); |
---|
51 | end |
---|
52 | |
---|
53 | %% Vï¿œrification des tendances de variation en courant |
---|
54 | % figure(100); |
---|
55 | % |
---|
56 | % hold on ; plot(IQ10./IQ10,'Color',nxtcolor) |
---|
57 | % hold on ; plot(IQ20./IQ10,'Color',nxtcolor) |
---|
58 | % hold on ; plot(IQ30./IQ10,'Color',nxtcolor) |
---|
59 | % hold on ; plot(IQ40./IQ10,'Color',nxtcolor) |
---|
60 | % hold on ; plot(IQ50./IQ10,'Color',nxtcolor) |
---|
61 | % hold on ; plot(IQ60./IQ10,'Color',nxtcolor) |
---|
62 | % |
---|
63 | % figure(200); |
---|
64 | % |
---|
65 | % hold on ; plot(IS10./IS10,'Color',nxtcolor) |
---|
66 | % hold on ; plot(IS20./IS10,'Color',nxtcolor) |
---|
67 | % hold on ; plot(IS30./IS10,'Color',nxtcolor) |
---|
68 | % hold on ; plot(IS40./IS10,'Color',nxtcolor) |
---|
69 | % hold on ; plot(IS50./IS10,'Color',nxtcolor) |
---|
70 | % hold on ; plot(IS60./IS10,'Color',nxtcolor) |
---|
71 | |
---|
72 | % sensibilite alphas |
---|
73 | |
---|
74 | ******** Summary for 'alphaby10_AMOR_new_mod_nov10_lin_auto_0' ******** |
---|
75 | Quadrupole change for quadFamList change of 0.100000 A |
---|
76 | Q1 (1A): Delta MCF1 = -2.35e-06 Delta MCF2 = -2.63e-05 |
---|
77 | Q2 (1A): Delta MCF1 = -6.25e-06 Delta MCF2 = -6.54e-05 |
---|
78 | Q3 (1A): Delta MCF1 = -1.29e-06 Delta MCF2 = -1.57e-05 |
---|
79 | Q4 (1A): Delta MCF1 = -1.90e-06 Delta MCF2 = +1.06e-05 |
---|
80 | Q5 (1A): Delta MCF1 = -9.29e-06 Delta MCF2 = +5.12e-05 |
---|
81 | Q6 (1A): Delta MCF1 = -2.37e-07 Delta MCF2 = -3.11e-06 |
---|
82 | Q7 (1A): Delta MCF1 = -3.77e-06 Delta MCF2 = -3.59e-05 |
---|
83 | Q8 (1A): Delta MCF1 = -2.22e-06 Delta MCF2 = -2.22e-05 |
---|
84 | Q9 (1A): Delta MCF1 = -1.63e-06 Delta MCF2 = +2.03e-05 |
---|
85 | Q10 (1A): Delta MCF1 = -9.00e-06 Delta MCF2 = +9.52e-05 |
---|
86 | ******** Summary for 'alphaby10_AMOR_new_mod_nov10_lin_auto_0' ******** |
---|
87 | Sextupole change for sextuFamList change of 0.100000 A |
---|
88 | S1 (1A) : Delta MCF1 = +4.48e-17 Delta MCF2 = -4.25e-06 |
---|
89 | S2 (1A) : Delta MCF1 = -1.51e-16 Delta MCF2 = -4.40e-07 |
---|
90 | S3 (1A) : Delta MCF1 = -3.12e-16 Delta MCF2 = -3.65e-06 |
---|
91 | S4 (1A) : Delta MCF1 = -1.00e-16 Delta MCF2 = -7.45e-06 |
---|
92 | S5 (1A) : Delta MCF1 = -1.24e-16 Delta MCF2 = -1.03e-07 |
---|
93 | S6 (1A) : Delta MCF1 = +1.33e-16 Delta MCF2 = -5.48e-07 |
---|
94 | S7 (1A) : Delta MCF1 = -3.58e-17 Delta MCF2 = -1.88e-07 |
---|
95 | S8 (1A) : Delta MCF1 = -3.31e-16 Delta MCF2 = -1.05e-06 |
---|
96 | S9 (1A) : Delta MCF1 = +1.64e-16 Delta MCF2 = -3.65e-06 |
---|
97 | S10 (1A) : Delta MCF1 = +5.27e-16 Delta MCF2 = -1.69e-05 |
---|
98 | |
---|
99 | % pour comparaison, les memes resultats pour l'optique low alpha maher |
---|
100 | ******** Summary for 'alphaby10_nouveau_modele_dec08_opt_lin_1' ******** |
---|
101 | Quadrupole change for quadFamList change of 0.100000 A |
---|
102 | Q1 (1A): Delta MCF1 = -1.88e-06 Delta MCF2 = +1.80e-04 |
---|
103 | Q2 (1A): Delta MCF1 = -4.10e-06 Delta MCF2 = +3.95e-04 |
---|
104 | Q3 (1A): Delta MCF1 = -7.99e-07 Delta MCF2 = +7.47e-05 |
---|
105 | Q4 (1A): Delta MCF1 = -2.84e-06 Delta MCF2 = -1.33e-04 |
---|
106 | Q5 (1A): Delta MCF1 = -1.37e-05 Delta MCF2 = -6.46e-04 |
---|
107 | Q6 (1A): Delta MCF1 = -1.24e-08 Delta MCF2 = -2.99e-07 |
---|
108 | Q7 (1A): Delta MCF1 = -8.39e-07 Delta MCF2 = -2.13e-05 |
---|
109 | Q8 (1A): Delta MCF1 = +9.54e-07 Delta MCF2 = +2.41e-05 |
---|
110 | Q9 (1A): Delta MCF1 = -2.43e-06 Delta MCF2 = +1.27e-04 |
---|
111 | Q10 (1A): Delta MCF1 = -1.34e-05 Delta MCF2 = +6.16e-04 |
---|
112 | ******** Summary for 'alphaby10_nouveau_modele_dec08_opt_lin_1' ******** |
---|
113 | Sextupole change for sextuFamList change of 0.100000 A |
---|
114 | S1 (1A) : Delta MCF1 = -1.56e-16 Delta MCF2 = +2.47e-06 |
---|
115 | S2 (1A) : Delta MCF1 = -3.28e-16 Delta MCF2 = -2.20e-07 |
---|
116 | S3 (1A) : Delta MCF1 = +1.89e-16 Delta MCF2 = -7.01e-06 |
---|
117 | S4 (1A) : Delta MCF1 = -3.52e-16 Delta MCF2 = -1.40e-05 |
---|
118 | S5 (1A) : Delta MCF1 = -1.65e-17 Delta MCF2 = -4.73e-09 |
---|
119 | S6 (1A) : Delta MCF1 = -1.59e-16 Delta MCF2 = -8.27e-08 |
---|
120 | S7 (1A) : Delta MCF1 = -1.31e-16 Delta MCF2 = -9.16e-09 |
---|
121 | S8 (1A) : Delta MCF1 = +1.45e-16 Delta MCF2 = -1.62e-07 |
---|
122 | S9 (1A) : Delta MCF1 = -1.00e-16 Delta MCF2 = -6.60e-06 |
---|
123 | S10 (1A) : Delta MCF1 = -3.81e-16 Delta MCF2 = -3.04e-05 |
---|