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
| CpA=inputdlg('Entrez la valeur de CpA');
>> CpA=str2num(CpA{1});
>> CpB=inputdlg('Entrez la valeur de CpB');
>> CpB=str2num(CpB{1});
>> for i=1:1:50
tmpA = abs(CCpA-CpA);
[idxA idxA] = min(tmpA);
PlusProcheA = CCpA(idxA); %Valeur la plus proche de CpA
indice1=find(CCpA==PlusProcheA);
tmpB = abs(CCpB-CpB);
[idxB idxB] = min(tmpB);
PlusProcheB = CCpB(idxB); %Valeur la plus proche de CpB
indice2=find(CCpB==PlusProcheB);
CCpA(indice1)=9999;
CCpB(indice2)=9999;
m(i,:)=[PlusProcheA PlusProcheB indice1 indice2];
end
A=m(:,3);
B=m(:,4);
C=intersect(A,B);
disp(['L angle Alpha vaut : ' num2str(Alpha(C))]);
disp(['L angle Beta vaut : ' num2str(Beta(C))]);
L angle Alpha vaut : -32.1012
L angle Beta vaut : -2.6944 |
