Discussion :

[matlabconv]

Bonjour
J'ai cherché à faire la convolution discrète entre deux matrices B=ones(3,3) et A=magic(3), j'ai tapé le code suivant :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
A = magic(3);
B=ones(3,3);
conv2(A,B)
 
ans =
 
     8     9    15     7     6
    11    17    30    19    13
    15    30    45    30    15
     7    21    30    23     9
     4    13    15    11     2

J'ai calculé la convolution discrète par sa formule ordinaire et je n'ai pas obtenu la même réponse, où est le problème?

Merci

[Jerome Briot]

Tu trouveras l'algorithme utilisé par la fonctionCONV2 dans la page de la documentation correspondante :

Code :

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

doc conv2

[matlabconv]

Je ne suis pas convaincu, je sait que conv(A,B)=FFT^-1(FFT(A).FFT(B)) en utilisant zero padding,
mais si je calcul conv(A,B)=sum_i=1,N(sum_j=1,N) A(i)B(j-i) (i=(i1,i2), j=(j1,j2)) (sans passe par FFT), je n obtient pas la meme reponse, c est ici que je me deande ou est le probleme
Merci

[matlabconv]

j´ai tape le code en fortran de la convolution discrete (sans FFT) en fortran j´ai obtenu seulement la premiere colone est identique a celle de FFT, svp pouvez vous me clarifier ce point
voila mon code

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
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
 
 
program gh
integer i,j, i1,i2,j1,j2
double precision A(0:2,0:2), B(0:2,0:2),C(0:2,0:2)
A(0,0)=0.8147
A(0,1)=0.9134
A(0,2)=0.2785
A(1,0)=0.9058
A(1,1)=0.6324
A(1,2)=0.5469
A(2,0)=0.1270
A(2,1)=0.0975
A(2,2)=0.9575
 
B(0,0)=0.9649
B(0,1)=0.9572
B(0,2)=0.1419 
B(1,0)=0.1576
B(1,1)=0.4854
B(1,2)=0.4218 
B(2,0)=0.9706
B(2,1)=0.8003
B(2,2)=0.9157
 
 
 
 
 
do i1=0,2
   do j1=0,2
   C(i1,j1)=0.d0
      do i2=0,2
	     do j2=0,2
		    C(i1,j1)=C(i1,j1)+A(i1-i2,j1-j2)*B(i2,j2)
	     enddo
      enddo
	 print*, i1,j1,C(i1,j1)
	enddo
enddo
 
end

[Invité]

Bonsoir,

Ne connaissant pas le Fortran, je ne peux que poser quelques questions:
- Pour i1,j1,i2,j2 dans {0, 1, 2}, tu auras plusieurs fois i1-i2 et j1-j2 en dehors des indices {0, 1, 2}, et donc dans ces cas, que retourne

Code :

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

A(i1-i2,j1-j2)

- De plus je vois que tu déclares C matrice 3x3 (initialisée avec des zéros?), donc est-ce que conv2(A,B,'same') répond à tes attentes?

[matlabconv]

Alors comment definir le produit de convolution discret

[Invité]

Ta définition est bonne, cependant il te faut faire attention aux indices qui posent problème.

[matlabconv]

Quést ce quíl faut modifier alors dans le programme

[Invité]

Par exemple un test if i1-i2>=0 && j1-j2>=0 avant la ligne C(i1,j1)=C(i1,j1)+A(i1-i2,j1-j2)*B(i2,j2)