Bonjour,
Voila je fais un TIPE sur la compression JPEG et jessaye dimplementer la DCT avec Maple
seulement jai un petit probleme car je ne retrouve pas les meme coeff dans ma matrice transformée que ceux donné en exemple sur internet ... Pourriez vous jeter un oeil a ma procedure et me dire ce qui cloche ?
Code:
dct := proc (M) 
N := coldim(M); P := matrix(N, N); 
for i to N do 
for j to N do 
if i = 1 then ci := 1/sqrt(2) else ci := 1 end if; if j = 1 then cj := 1/sqrt(2) else cj := 1 end if; 
P[i, j] := floor(2*ci*cj*(sum(sum(evalf(1/2*M[x, y]*cos(1/2*(2*x-1)*i*Pi/N)*cos*(2*y-1)*j*Pi/N), x = 1 .. N), y = 1 .. N))/N) 
end do end do; evalm(P) end proc
je vous remercie davance