Discussion :

[Bilala92]

Bonjour,
j'aimerai déclarer puis représenter une fonction porte de la sorte : f(t)=1 pour |t|<T/2 et f(t)=0 dans tous les autres cas. T étant la largeur de la porte. Je débute et ne vois vraiment pas comment faire.
Cordialement,

[phryte]

Bonjour,
Voici deux exemples :

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
clear
fs = 10000;
t = -1:1/fs:1;
T=0.1;
x2 = rectpuls(t,T);
plot(t,x2)
grid
 
 
clear
T=0.1;
porte=@(t) abs(t)<T/2;
t=linspace(-1,1,1e5);
figure(2)
plot(t,porte(t))
grid

[Bilala92]

Bonjour, merci beaucoup pour votre aide. Juste une petite vérification, si je calcule le produit de convolution : ma porte * ma porte, cela donne bien une sorte de signal triangle ?
Ensuite j'ai rajouté

Code :

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

1
2
3
y=conv(x2,x2,'same');
 
plot(t,x2,t,y)

mais le graphe est illisible... auriez-vous une petite astuce pour le rendre plus lisible svp ?

[phryte]

Bonjour,

Code :

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

1
2
3
4
5
p=double(porte(t));
y=conv2(p,p,'same')
figure(3)
plot(t,y)
grid

[Bilala92]

Il ne reconnait pas la fonction "porte(t)" et je souhaiterai afficher la porte ainsi que le résultat de la convolution sur le même graphe si possible

[phryte]

Bonjour,
Fais attention aux variables logiques et réelles

[Bilala92]

Bonsoir;
J'ai maintenant

Code :

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

1
2
3
4
5
6
7
8
9
 
clear
fs = 1000;
t = linspace(-5,5,fs);
T=4;
x2 = rectpuls(t,(T/2)); 
y=conv(x2,x2,'same');
figure (1)
plot(t,x2,t,y),axis([-5 5 0 1])

Cependant je n'arrive pas à trouver une échelle qui permettrait de tout voir à l'écran, est-ce possible ?

[phryte]

Bonjour,
Tu peux normaliser :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
clear
fs=1000;
t=linspace(-5,5,fs);
T=4;
x2=rectpuls(t,T/2);
y=conv(x2,x2,'same');
figure(1)
plot(t,max(y)*x2,t,y)
grid
figure(2)
plot(t,x2,t,y/max(y))
grid

[Bilala92]

Merci beaucoup pour votre aide !!