besoin de commentaires pour un prog matlab

[invite43361057]

Bonjour,

J'ai le programme suivant qui calcule le produit de convolution de deux fonctions, le problème c'est que je ne comprends pas toutes les lignes du programme et comment il marche exactement.

function y=convol()
close all;
t=0:0.1i;
x=cos(t);
g=sin(t);
X = length(x);
G = length(g);
K = X+G-1;
y(1:K)=0;
x(X+1:K)=0;
for k=1:K
for l=k-G:k-1
if l>=0
y(k)=y(k)+x(l+1)*g(k-l);
axis([0 100 -10 10])
subplot(2,1,1)
plot(k,y(k));
hold on
end
end
end
subplot(2,1,2)
plot(conv(x,g))

C'est surtout au niveau du principe de calcul que je vois pas.
Est ce que quelqu'un peut m'aider?
Merci d'avance


Ps: c'est assez urgent svp
-----

[invite86822278]

Bonjour,

function y=convol()
close all;
t=0:0.1i;
x=cos(t);
g=sin(t);

Constitution des vecteurs - initialisations

X = length(x);
G = length(g);
K = X+G-1;

tailles des vecteurs

y(1:K)=0;
x(X+1:K)=0;

On complete x avec des 0 on initialise y a 0 sur K elements

for k=1:K
for l=k-G:k-1
if l>=0
y(k)=y(k)+x(l+1)*g(k-l);
axis([0 100 -10 10])
subplot(2,1,1)
plot(k,y(k));
hold on
end
end
end

A y qui doit contenir le produit de convolution on associe la somme correspond a la definition discrete d'un produit de convolution :
En terme de fonction :


En terme de suite discrete :



ou (equivalent)



En considerant dans la somme en i que les valeurs non definis sont nulles.

Ce qui permet de retomber sur ta somme entre k-G et k-1

subplot(2,1,2)
plot(conv(x,g))

surement pour comparer ? je ne connais pas conv, mais le nom semble explicite.

[invite43361057]

conv calcule le produit de convolution de 2 fonctions

quand je compare les deux tracés, il y a trop de points avec le 1°

où dois je mettre le 1°plot pour que ce soit pareil?

[invite86822278]

essaye de decaler et de remplacer ta somme par
y(k)=y(k)+x(l)*g(k-l+1);

Dans ce cas, les bornes de ta somme interne doivent etre de k-G+1 a k.
Je n'ai rien teste, n'ayant pas matmab sous la main, mais a premiere vue, si tu dis qu'il y a juste un decalage, ca pourrait etre ca.

[A voir en vidéo sur Futura]

[invite43361057]

en fait c'est bon j'ai trouver, je posterai dans peu de temps le prog fini pour ceux qui veulent