Discussion :

[nonodev]

Bonjour à tous.

Je suis actuellement en stage et mon dernier travail est de modéliser en 3 dimensions le diagramme de rayonnement d'une antenne WIFI avec MATLAB.
Je me permets de préciser que je ne suis pas un mathématicien dans l'âme, veuillez donc m'excuser des possibles stupidités qui vont suivre.

Mon antenne WIFI donc .. L'équation caractéristique est la suivante :
<P> = 1/(360pi) * Io² * L²/lambda² * 1/(4epsilon0*C*r)² * sin²(theta)

Pour faire plus simple, je vous propose de se limiter à l'équation suivante :
f(theta,phi) = theta avec theta € [0,pi] et phi € [0,2pi]
r est considéré comme fixe, disons r=5.
Ceci devrait me donner un cône (une droite qui tourne de 360° autour d'un axe).

Je suis donc en présence de coordonnées sphériques, et mes recherches
m'ont amené au raisonnement suivant :

1. Je pose :
r = 5;
theta = 0:pi/50;pi;
phi = 0:pi/50:2*pi;

2. Transformation des coordonnées sphériques en coordonnées cartésiennes :
> Utilisation de SPH2CART
[x,y,z] = sph2cart(theta,phi,r);

3. Utilisation de matrices régulières (X et Y) :
> Utilisation de MESHGRID
[X,Y] = meshgrid(x,y);

4. Interpolation :
> Utilisation de INTERP2
> Là je suis complètement perdu avec les XI, YI et ZI, mais
> je propose quelque chose comme ça :
[ZI] = interp2(x,y,z,XI,YI);

5. Tracé de la courbe en utilisant MESH


Bien évidemment ça ne marche pas.
Avant toute chose, mon raisonnement est-il correct ?
Si j'ai bien compris, l'interpolation permet de calculer un grand nombre de points de coordonnées x,y,z afin de tracer. Est-ce correct ?
Pourriez-vous me préciser les notions de XI, YI et ZI ?


Merci d'avance,
Arnaud.

[mr_samurai]

Salut,

Essayes comme ceci :

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
 
 
AP = [0 0 0]; % Access Point
 
teta = 0:0.1:2*pi;
phi = 0:0.1:pi;
 
[TETA PHI] = meshgrid(teta,phi);
R = ones(size(TETA));
 
% Directional AP 
PathLoss = sin(TETA).^2 .* abs(sin(PHI));
 
% Plot
[X Y] = sph2cart(TETA,PHI,R);
surf(X,Y,PathLoss)
alpha(0.2)
hold on
scatter3(AP(1),AP(2),AP(3),'*')

[nonodev]

Ok, merci pour le programme.
Ceci dit, je ne comprends pas où est est spécifié l'équation.
Si par exemple, je veux tracer f(theta,phi) = 3 * sin²(theta) ??

Merci d'avance