L = membrane
L = membrane(k)
L = membrane(k, n)
membrane(...)
| Paramètre | Description |
|---|---|
| k | un entier positif : l'indice de la fonction propre (par défaut 1). |
| n | un entier positif : le paramètre de grille ; le résultat est (2*n+1)-par-(2*n+1) (par défaut 15). |
| Paramètre | Description |
|---|---|
| L | une matrice (2*n+1)-par-(2*n+1) : la k-ième fonction propre échantillonnée sur la grille. |
membrane renvoie la k-ième fonction propre de la membrane en forme de L, échantillonnée sur une grille carrée (2*n+1)-par-(2*n+1). Le domaine est le carré privé d'un quadrant ; la valeur est nulle sur la frontière et en dehors de ce domaine.
La première fonction propre (k = 1) est le mode fondamental, une unique bosse positive ; les indices supérieurs ajoutent des lignes nodales et donc des valeurs négatives. Le résultat est normalisé de sorte que son amplitude maximale vaut 1 avec un pic positif.
Appelée sans argument de sortie, membrane trace la surface.
La fonction propre est obtenue par une discrétisation en différences finies du domaine, résolue par un calcul de valeurs propres en mode shift-invert (eigs) pour les plus petites valeurs propres.
L = membrane(1, 10);
size(L)L = membrane(1, 15);
[min(L(:)), max(L(:))]f = figure();
surf(membrane(1, 50), 'EdgeColor', 'none');
axis('off');| Version | Description |
|---|---|
| 2.0.0 | version initiale |