Discussion :

[saimonpresshill]

Bonjour à tous,
Je souhaite tracer avec un double axe en abscisse une courbe de réflectivité.
Je trace mon signal en fonction de la longueur d'onde en microns. Jusque là pas de problème. Seulement je souhaiterais aussi ajouter un axe horizontal en nombre d'onde (i.e. en cm-1).
Pour ce faire, je défini un axe en microns puis le correspondant en cm-1. Ensuite je trace le tout en me référant à l'aide de Matlab : Using Multiple X- and Y-Axes.
Cela donne le code suivant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
 
clear all;close all;
 
%longueur d'onde en microns
lambda1 = 5; 
lambda2 = 20; 
 
%nb points
N = 100;
 
%definition des axes
Axe_microns = linspace(lambda1,lambda2,N);
Axe_cm = 1e4./Axe_microns;
 
% signal a afficher
R = rand(1,100);
 
% affichage
figure
hold on
hl1 = line(Axe_microns,R,'Color','r');
ax1 = gca;
ax2 = axes('Position',get(ax1,'Position'),...
           'XAxisLocation','top',...
           'YAxisLocation','left',...
           'Color','none',...
           'XColor','k','YColor','k','Xdir','reverse');
hl2 = line(Axe_cm,R,'Color','k','Parent',ax2);
 
% Ajustement des Ticks
xlimits = get(ax1,'XLim');
ylimits = get(ax1,'YLim');
xinc = (xlimits(2)-xlimits(1))/5;
yinc = (ylimits(2)-ylimits(1))/5;
set(ax1,'XTick',[xlimits(1):xinc:xlimits(2)],...
        'YTick',[ylimits(1):yinc:ylimits(2)])
 
xlimits = get(ax2,'XLim');
ylimits = get(ax2,'YLim');
xinc = (xlimits(2)-xlimits(1))/5;
yinc = (ylimits(2)-ylimits(1))/5;
set(ax2,'XTick',[xlimits(1):xinc:xlimits(2)],...
        'YTick',[ylimits(1):yinc:ylimits(2)])

Et là c'est clair que ça marche pas. 8 microns n'ont jamais fait 1700cm-1.
Le problème vient du fait que les deux axes sont inversement proportionnels et du coup le pas en cm-1 n'est pas régulier, et à partir de là je ne sais plus quoi faire.
Comment faire coïncider correctement mes deux axes ??
Merci d'avance.

[Jerome Briot]

En fait, il faut modifier le premier Axes avant de créer le second

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
figure
 
plot(rand(1,5),rand(1,5));
 
% Modification des propriétés du premier Axes
set(gca,'xlim',[0 1])
 
% Construction du second Axes
xt = get(gca,'xtick');
axes('Position',get(gca,'Position'),...
           'XAxisLocation','top',...
           'YAxisLocation','left',...
           'Color','none',...
           'XColor','k','YColor','k', ...
           'ytick',[],...
           'xticklabel',num2str(1./xt(:),'%.2f'));

[saimonpresshill]

Parfait, merci !