1 | function [CorFamRes,CorDevRes,CorElemRes,MaxEffRes,BpmFamRes,BpmDevRes]= mec(QuadFamily,QuadDev,Plane,CorBpmResp,CorNumb) |
---|
2 | % MEC |
---|
3 | % function [CorFam,CorDev,MaxEff]= mec(QuadFamily,QuadDev,Plane,CorBpmResp) |
---|
4 | % |
---|
5 | % INPUTS -> VALUES |
---|
6 | % |
---|
7 | % 1. QuadFamily: Family of the quadrupole studied -> 'Q2' |
---|
8 | % 2. QuadDev: [cell elt] of the quadrupole studied -> [1 1] |
---|
9 | % 3. Plane: Plane=1 horizontal only |
---|
10 | % =2 vertical only |
---|
11 | % 4. CorBpmResp: Response matrix of correctors obtained from "respcor(plane)" |
---|
12 | % |
---|
13 | % |
---|
14 | % OUPUTS |
---|
15 | % |
---|
16 | % 1. CorFamRes: Family of the Most Effective Corrector found |
---|
17 | % 2. CorDevRes: [cell elt] of the MEC found |
---|
18 | % 3. MaxEffRes: Efficacity of the MEC in mm/A |
---|
19 | % |
---|
20 | |
---|
21 | if ~iscellstr(QuadFamily) |
---|
22 | error('la famille entree n est pas au format cell') |
---|
23 | end |
---|
24 | |
---|
25 | tailleQDev=size(QuadDev(1,:)); |
---|
26 | tailleCorBpm=size(CorBpmResp); |
---|
27 | |
---|
28 | if isempty(QuadDev) |
---|
29 | error('Veuillez entrer le device du quadripole à étudier au format QuadDev=[cell elt]') |
---|
30 | elseif tailleQDev(2)~=2 |
---|
31 | error('Trop d''elements pour 1 device de QP: QuadDev=[cell elt]') |
---|
32 | end |
---|
33 | |
---|
34 | if Plane<0 || Plane>2 |
---|
35 | error('Veuillez entrer le plan étudié: Horizontal: Plane=1 ; Vertical:Plane=2 ; Both:Plane=0') |
---|
36 | end |
---|
37 | |
---|
38 | if tailleCorBpm(1)~=56 & tailleCorBpm(2)~=120 |
---|
39 | error('Veuillez vérifier la matrice réponse correcteurs/BPM') |
---|
40 | end |
---|
41 | |
---|
42 | |
---|
43 | |
---|
44 | |
---|
45 | |
---|
46 | load('respcor.mat') |
---|
47 | |
---|
48 | if Plane==1 |
---|
49 | BpmFam='BPMx'; |
---|
50 | CorFam='HCOR'; |
---|
51 | %CorBpmResp=efficacy.HPLane; %effic(Cor,Bpm) |
---|
52 | elseif Plane==2 |
---|
53 | BpmFam='BPMz'; |
---|
54 | CorFam='VCOR'; |
---|
55 | %CorBpmResp=efficacy.VPlane; |
---|
56 | else |
---|
57 | error('input must "1" for horizontal plane or "2" for vertical plane'); |
---|
58 | end |
---|
59 | |
---|
60 | |
---|
61 | |
---|
62 | %Most Effective Corrector |
---|
63 | |
---|
64 | %BPM le plus proche du Qp (QuadFamily,[QuadDev]); |
---|
65 | [BpmFam2,BpmDev,BpmPos]=proche2(BpmFam,QuadFamily,QuadDev,0); |
---|
66 | |
---|
67 | indBpm=dev2elem(BpmFam,BpmDev); |
---|
68 | % Faire un tri sur EffCor de maniÚre ᅵ pouvoir choisir par exemple le 2Úme correcteurs le plus |
---|
69 | % efficace si le premier ne convient pas... |
---|
70 | EffCor=CorBpmResp(:,indBpm); % vecteur colonne de longueur nb HCOR ou nb VCOR |
---|
71 | [posEff,indEff]=sort(abs(EffCor),'descend'); |
---|
72 | |
---|
73 | %MaxEff=max(abs(EffCor)); % extraction de l'efficacité max |
---|
74 | MaxEff=posEff(CorNumb); |
---|
75 | %indCor=find(abs(EffCor)==MaxEff); % recherche de l'indice de l'efficacité max dans "EffCor" |
---|
76 | indCor=indEff(CorNumb); |
---|
77 | CorDevList=getlist(CorFam); |
---|
78 | |
---|
79 | %Rï¿œsultats |
---|
80 | CorFamRes=CorFam; |
---|
81 | CorDevRes=CorDevList(indCor,:); |
---|
82 | CorElemRes=indCor; % Recherche du correcteur dᅵlivrant l'efficacitᅵ prᅵcᅵdente |
---|
83 | MaxEffRes=MaxEff; |
---|
84 | BpmFamRes=BpmFam; |
---|
85 | BpmDevRes=BpmDev; |
---|
86 | |
---|
87 | |
---|
88 | |
---|
89 | |
---|
90 | %************************************************************************** |
---|
91 | |
---|
92 | |
---|