1 | function induced_amplitude6D |
---|
2 | %function induced_amplitude |
---|
3 | % |
---|
4 | % Plot induced amplitudes in achromat and straight section |
---|
5 | % transported back to the beginning of the lattice |
---|
6 | % NB: use files computed w/ solamor2.lat |
---|
7 | % |
---|
8 | % Written by Laurent S. Nadolski, SOLEIL 04/2004 |
---|
9 | |
---|
10 | hold_state = ishold; |
---|
11 | hold on |
---|
12 | |
---|
13 | %% Doesnot work for Pascale |
---|
14 | %PATH=[getenv('HOME') '/matlab/soleil/'] |
---|
15 | PATH='/home/nadolski/matlab/soleil/'; |
---|
16 | |
---|
17 | %% amplitude induite dans les achromates |
---|
18 | file = [PATH 'amp_ind_achromat.out']; |
---|
19 | [header data] = hdrload(file); |
---|
20 | x_achromat = sqrt(data(:,8).*data(:,4)).*abs(data(:,1))*1e3; |
---|
21 | |
---|
22 | %% amplitude induite dans sections moyennes |
---|
23 | file = [PATH 'amp_ind_SDM.out']; |
---|
24 | [header data] = hdrload(file); |
---|
25 | x_SDM = sqrt(data(:,8).*data(:,4)).*abs(data(:,1))*1e3; |
---|
26 | |
---|
27 | %% amplitude induite dans sections courtes |
---|
28 | file = [PATH 'amp_ind_SDC.out']; |
---|
29 | [header data] = hdrload(file); |
---|
30 | x_SDC = sqrt(data(:,8).*data(:,4)).*abs(data(:,1))*1e3; |
---|
31 | |
---|
32 | %% ampltudes induites dans sections longues |
---|
33 | file = [PATH 'amp_ind_SDL.out']; |
---|
34 | [header data] = hdrload(file); |
---|
35 | x_SDL = sqrt(data(:,8).*data(:,4)).*abs(data(:,1))*1e3; |
---|
36 | delta = data(:,1)*100; |
---|
37 | |
---|
38 | h=plot(delta,2*x_SDL,'k-.',delta,2*x_SDC,'k--',delta,2*x_achromat,'k-'); |
---|
39 | set(h,'LineWidth',3) |
---|
40 | legend('SDL/SDM','SDC', 'Achromat', 'Physical') |
---|