Discussion :

[Erzvinh]

Bonjour,
je vous expose mon problème, je n'arrive pas à tracer une fonction de type créneau qui je devrai utiliser pour un produit de convolution.

Code :

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

1
2
3
4
5
6
7
8
 
function y=f_conv(x)
if((x>1) || (x<-1))
	y=0;
else
	y=1;
endif
endfunction

la fonction plot(x,f_conv(x)) ne marche pas mais si je tape f_conv(1) ou autre il me donne la bonne réponse....je pense que cela vient du fait que le resultat est un scalaire....mais n'ayant pas du tout l'habitude de ce logiciel je ne vois pas comment faire .
Merci d'avance

[rostomus]

Le probleme est le fait que x est un vecteur et non pas une seule valeur
vous pouvez faire une boucle, mais c'est deconseillé.
je vous propose ca:

Code :

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

1
2
3
4
 
function y=f_conv(x)
y=zeros(size(x));
y( (x >= -1) & (x <= 1))=1;

[Erzvinh]

Merci beaucoup pour votre réponse précise et rapide!!
cela marche parfaitement.

[tug83]

tu n'as plus qu'à mettre le sujet en résolu!
bye

[Erzvinh]

Bonjour,
j'ai encore un autre souci du meme ordre....j'aimerai definir une fonction de ce type

Code :

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

1
2
3
4
5
 
function y=g_conv(x)
y=zeros(size(x));
y( (x >= -1) & (x <= 1))=1-abs(x);	
endfunction

je me suis inspiré de la fonction précedente mais impossible de la tracer ....pourtant les resultats retourner en sortie sont bons....le fonctionement devrait pourtant etre identique à la précédente non?

Merci d'avance pour vos réponses.

[rostomus]

Bonjour,
la taille de y( (x>= -1) & (x<= 1)) est differente de celui de abs(x), donc la solution est:

Code :

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

1
2
3
4
5
 
function y=g_conv(x)
y=zeros(size(x));
y( (x >= -1) & (x <= 1))=1-abs(x( (x >= -1) & (x <= 1) ) );	
endfunction

ou bien :

Code :

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

1
2
3
4
5
6
 
function y=g_conv(x)
y=zeros(size(x));
ind=find((x >= -1) & (x <= 1));
y(ind)=1-abs(x( ind) );	
endfunction

[Erzvinh]

Merci encore pour votre réponse!
Quand on a la solution cela parait évident ....
Je passe le sujet en résolu!

[Jerome Briot]

Code :

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

1
2
3
4
5
6
 
function y=g_conv(x)
y=zeros(size(x));
ind=find((x >= -1) & (x <= 1));
y(ind)=1-abs(x( ind) );	
endfunction

Bonjour,

n'utilisez FIND que lorsque cela est vraiment nécessaire :

Code :

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

1
2
3
4
5
function y=g_conv(x)
y=zeros(size(x));
ind=(x >= -1) & (x <= 1);
y(ind)=1-abs(x( ind) );	
%endfunction

Question: qu'est ce que endfunction ?

On peut aussi encore optimiser ce code comme ceci :

Code :

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

1
2
3
4
function x=g_conv(x)
ind=(x >= -1) & (x <= 1);
x=x.*ind;
x(ind)=1-abs(x( ind));

[rostomus]

Bonjour,

ind=(x >= -1) & (x <= 1);

absolument, ca m'a echappé

pour le endfunction, je ne sais rien, et je pense qu'elle n'est même pas definie