Discussion :

[code_source]

Bonjour,

Je reviens encore vers vous car je n'ai toujours pas de réponse concernant mon calcul d'intégrale d'une matrice.

Je vous explique le problème: J'ai une matrice de dimension (23x30) et je souhaite calculer son intégrale double. Tout d'abord, j'ai définit une grille (avec la fonction meshgrid) pour avoir les bornes de l'intégrale. J'ai essayé la fonction dblquad, mais j'ai des erreurs car il s'agit d'une matrice et non d'une fonction.

Si quelqu'un peut m'apporter des explications, ça m'aidera beaucoup
Merci

[Kangourou]

Salut,

ben, c'est deux problèmes différents, non ?

Pour calculer une intégrale double, l'entrée doit être une fonction de deux variables qui renvoie une valeur scalaire.

Quand on dispose d'une matrice, on peut imaginer que les valeurs de la matrice correspondent à des valeurs de la fonction après un échantillonnage selon une grille régulière. Le plus simple est de calculer la somme des éléments de la matrice, et de multiplier par la résolution en x et y. Cela correspond à une discrétisation de l'intégrale double.

Est-ce plus clair ?

[code_source]

Salut,

ben, c'est deux problèmes différents, non ?

Pour calculer une intégrale double, l'entrée doit être une fonction de deux variables qui renvoie une valeur scalaire.

Quand on dispose d'une matrice, on peut imaginer que les valeurs de la matrice correspondent à des valeurs de la fonction après un échantillonnage selon une grille régulière. Le plus simple est de calculer la somme des éléments de la matrice, et de multiplier par la résolution en x et y. Cela correspond à une discrétisation de l'intégrale double.

Est-ce plus clair ?

Merci beaucoup pour ta réponse (Kangourou), je commençais vraiment à désespérer.
Si j'ai bien compris, pour calculer l'intégrale double d'une matrice, il faut d'abord définir une grille régulière de même taille que la matrice:

Code :

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

1
2
3
xedges=linspace(0,1344,30);
yedges=linspace(0,1024,23);
[Xbox,Ybox]=meshgrid(xedges,yedges);

Ensuite, je calcule la somme des éléments de la matrice que je multiplie par la résolution en x et y:

Code :

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

res=sum(MaMatrice)*diff(xedges)*diff(yedges);

La somme d'une matrice est bien un vecteur? (on somme sur chaque colonne). Le résultat n'est donc pas une seule valeur, mais un vecteur. Est ce que c'est juste.

Maintenant si je veux calculer l'intégrale de cette matrice, mais sur chaque case de la grille, je fais comment??

[FR119492]

Salut!

J'ai une matrice de dimension (23x30) et je souhaite calculer son intégrale double.

Cette phrase n'a aucun sens pour deux raisons:

1. "matrice" est un terme d'algèbre linéaire; à première vue, il semble que tu as plutôt un tableau de valeurs.
2. On calcule l'intégrale d'une fonction et non d'un ensemble de valeurs discrètes.

Jean-Marc Blanc

[code_source]

Salut,

ben, c'est deux problèmes différents, non ?

Pour calculer une intégrale double, l'entrée doit être une fonction de deux variables qui renvoie une valeur scalaire.

Quand on dispose d'une matrice, on peut imaginer que les valeurs de la matrice correspondent à des valeurs de la fonction après un échantillonnage selon une grille régulière. Le plus simple est de calculer la somme des éléments de la matrice, et de multiplier par la résolution en x et y. Cela correspond à une discrétisation de l'intégrale double.

Est-ce plus clair ?

Merci (Kangourou), le code marche bien, j'ai bien sommé les éléments de la matrice et multiplié par la résolution en x et y.

Code :

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

sum(sum(FLUX.*dS));