Distribution du champ E.M pour une structure de 20 paires de couches

[invite08b82e0f]

bonsoir à tous
alors mon problème est que je veux dessiner les 20 paires de couches derrière ma courbe de distribution comme celle la:

-----

[Jack]

Avec les infos que tu donnes, je dirais, une feuille de papier, des feutres de couleur et une règle.

[invite08b82e0f]

lol ,

vous avez raison,
alors, j'ai utilisé ce code :

x = 0:1:10;
y = 1;
 
figure
 
nx = numel(x);
 
b = bar(repmat(y,1,nx), 1);
 
bb = get(b,'children');
c = ones(nx,1);
c(2:2:end) = 2;
set(bb, 'facevertexcdata', c)
 
axis tight
 
colormap([0.71 0.87 0.90 ; 0.19 0.52 0.61])

afin de dessiner les barres souhaité
et je veux fusionner les 2 figures ensembles.

premièrement matlab m'affiche toujours une erreur :" Attempt to execute SCRIPT bar as a function:
I:\Rihab Work\simulation MATLAB\Intensity\bar.m

Error in barr (line 8)
b = bar(repmat(y,1,nx), 1);"

je ne sais pas comment le corriger.
et deuxièmement je veux fusionner ces deux résultats sur une même figure.
merci d'avance.

[Jack]

C'est mieux

Un peu de patience et tu devrais obtenir des réponses.

[A voir en vidéo sur Futura]

[invite08b82e0f]

concernant la première partie du question je l'ai corrigé: en fait: mon script est nommé "bar" , alors le script apparaît avant la fonction sur le chemin,et les commandes seront utilisées à la place de la fonction. donc j'ai juste changé le nom de fichier, et ça fonctionne.
sinon j'attend la réponse sur ma deuxième question.
et merci

[invite896757ff]

Non mais là, tu n'es pas dans de la retouche d'image; on ne "fusionne" pas.
Ton code doit générer une seule image contenant les deux types d'information.
Alors soit tu génère alternativement soit successivement

[Antoane]

Bonjour,

% Prepare the Figure
figure; hold on
 
% plot bars
x = 0:1:10;
y = 1;
 
nx = numel(x);
 
b = bar(repmat(y,1,nx), 1);
 
bb = get(b,'children');
c = ones(nx,1);
c(2:2:end) = 2;
set(bb, 'facevertexcdata', c)
 
axis tight
 
colormap([0.71 0.87 0.90 ; 0.19 0.52 0.61])

% Add |E|²
x1 = linspace(-1, 13, 1e3) ;
E = 0.5*sin(x1*5*pi).*exp(-(x1-6).^2/10) +0.5;
plot(x1, E)