convolution numérique sous matlab

[invite611f7ca3]

Bonjour à tous,

je galere un peu pour faire un petit algorithme qui fait la convolution de deux signaux echantillonné avec matlab.

Je part de deux signaux porte et veut faire la convolution mais je galere je connais la théorie du "bordel" mais alors la mise en pratique s'avere trop compliqué pour moi

je vous met ici mon petit prog si quelqu'un peut m'aider :

x=ones(1,10);
h=ones(1,10);

Nx=length(x);
Nh=length(h);

y=zeros(1,Nx+Nh);

for i=1:Nx+Nh
for k=1:Nh
y(i)=y(i)+h(k)*x(i-k);

end;
end;
plot(y);

??? Subscript indices must either be real positive integers or logicals.

J'espere que qu'elqu'un pourras m'aider je sais pas ou si quelqu'un peut me dire ou je peu trouver de l'aide cela me serais bien util.
merci d'avance
-----

[invite22c0acb2]

Hi
Déja, y'a une erreur dans la dimension de y, il doir etre de la somme moins un
En plus y'a deja la fonction "conv", pour plus d'aide tape 'help conv' tu sauras comment l'utiliser.

[invite611f7ca3]

ouais je sais qu'elle doit surement exister cette fonction mais j'aimerais faire sans.

[invite611f7ca3]

apres quelques heures de prise de tete j'arrive enfin a avancer mais la j'ai un dernier petit probleme mais je sais pas ou il se trouve o_O

x=ones(1,3);
h=ones(1,3);

Nx=length(x);
Nh=length(h);

y=zeros(1,Nx+Nh-1);

for i=1:Nx+Nh-1
for k=1:Nh
if (i-k)>0
if (i-k)<=Nx
y(i)=y(i)+h(k)*x(i-k);
end;
end;
end;
end;
figure(1)
plot(y,'o');

[A voir en vidéo sur Futura]

[invitee6636d48]

svp si vous a fait la solution de la convoltion sous matlab,pouvez vous me l'envoyez car je galaire grave,merci d'avence
#######################

Pas d'adresse email d'après la charte

[soft001]

Bonjour,
C’est un peu retard mais c’est la première fois que je visite ce Forum.
Avec le travail collectif rien n’est possible.
Tient la solution :

clear all
x=[1 2 3 5];
h=[4 5 6 7];
conv_matlab=conv(x,h);
Nx=length(x);
Nh=length(h);
min_lenght=min(Nx,Nh);
max_lenght=max(Nx,Nh);
mon_code=zeros(1,Nx+Nh-1);
mon_code(1)=h(1)*x(1);
for k=2:Nx+Nh-1
for j=1:Nx+Nh
if(k+1-j>0 & k+1-j<=min_lenght &j<=max_lenght)
mon_code(k)=mon_code(k)+h(j)*x (k+1-j);
end
end
end
fprintf('Résultat de la fonction "conv"\n')
conv_matlab
fprintf('\nRésultat de mon code\n')
mon_code