Discussion :

[cedric973]

Bonjour,

Je réalise un petit programme calculant le diagramme de rayonnement d'une antenne à ouverture circulaire, je vous met le code çi-dessous :

Code :

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

1
2
3
4
5
6
7
[theta]=meshgrid(0:1:180);
a=2;
fr=10*10^10;
lambda=3*10^8/(fr);
u = (2*pi*a*sin(theta)/lambda)
 
fonction=2*(besselj(1,(u)))/u;

Cependant je reçois en retour :

Code :

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

1
2
Warning: Matrix is singular to working precision
> In ouverturecirculaire at 7

J'ai regardé sur un ancien post, pour ce dernier il s'agissait d'un pb avec la fonction "inv" mais moi je ne l'utilise pas et je ne comprends pas trop mon erreur, si quelqu'un pouvait me donner un coup de main svp???

Merci d'avance.

[Jerome Briot]

Il me semble que le problème vient du fait que theta contient des valeurs identiques

D'ailleurs, cette variable ne devrait-elle pas être un vecteur plutôt qu'une matrice ?

Code :

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

1
2
3
4
5
6
7
8
theta = linspace(0,180,181);
 
a = 2;
fr = 10*10^10;
lambda = 3*10^8/(fr);
u = (2*pi*a*sin(theta)/lambda);
 
fonction = 2*(besselj(1,(u)))/u;

[AF_2.8]

Le problème vient du fait que ta matrice est singulière (donc mal conditionnée) et donc ceci pause problème pour l'inverser. Ce n'est pas un problème de matlab, mais un problème de calcul matriciel.
Normalement une matrice singulière n'est pas inversible, il te faut alors calculer son pseudo-inverse.

[cedric973]

Merci dut, bien en fait ma "fonction" est en fonction de theta et me renvoie une valeur comprise entre 0 et 1 pour chaque valeurs de theta entre 0 et 180.
Dans le cas de ce code, je n'obtiens qu'une valeur de "fonction"...

andnicam, merci également de ta réponse, mais je ne comprends pas très bien : de quelle matrice parles tu qui est mal conditionnée? et ou procède t'on à l'inversion de cette dernière?

[AF_2.8]

U n'est pas un scalaire, mais une matrice, et tu divises par u (la division revient à multiplier par la matrice inversée).
Pour faire la division par U matlab inverse U et là il te dit que ta matrice est mal conditionnée. D'ou ton message d'erreur de matlab. Celà n'a aucun rapport avec le fait de diviser par zéro qui produirait un autre message.