Discussion :

[Leonpolou]

Bonjour,

je souhaites tout simplement tracer la fonction suivant avec MATLAB

Code :

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

1
2
3
 
f(t)=asin(wt)+C si sin(wt)>0
f(t)=C sinon

En gros c'est une sinusoïde qui est filtre avec un pont de diode a laquelle on rajoute un offset C en sortie pour décaler le signal.

Comment fait ton pour instaurer de telles conditions sous matlab?

Merci
Jeremy

[duf42]

Bonsoir,

Je peux pas tester mais est-ce que quelquechose du genre:

Code :

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

1
2
 
f(t)=max(asin(wt)+C, C);

Pourrait convenir?

Bonne soirée,
Duf

[Jean Dumoncel]

Cela suppose que a est positif, non?

Sinon pourquoi ne pas faire une bonne vieille boucle avec une condition?

[duf42]

Ah oui tout à fait, j'avais pas pensé à toutes les conditions

Cela dit, il doit quand même y avoir un moyen d'éviter de faire une boucle (en mettant une condition sur le signe de A je pense car sauf autre erreur de ma part, la valeur limite est toujours C)

[EDIT] Donc plutôt un truc du genre:

Code :

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

1
2
3
4
5
6
 
if (a>0)
f(t)=max(a*sin(w*t)+C, C);
else
f(t)=min(a*sin(w*t)+C, C);
end

[Jean Dumoncel]

Maintenant je suis d'accords!

[Leonpolou]

Oui a est forcement positif

j'ai essaye ca comme code

Code :

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

1
2
3
4
5
6
7
8
9
 
for t=0:0.01:10
   if sin(o*t)>0
      y=a*sin(o*t)+c;
   else
      y=c;
   end
end
plot(t,y)

mais la figure ne m'affiche aucune courbe!

[Jean Dumoncel]

Si a est forcément positive, alors utilise la première solution de duf42 :

Code :

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

f(t)=max(asin(wt)+C, C);

Si tu n'es pas sûr que a soit toujours positive alors utilise la deuxième solution toujours proposé par duf42 :

Code :

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

1
2
3
4
5
if (a>0)
f(t)=max(a*sin(w*t)+C, C);
else
f(t)=min(a*sin(w*t)+C, C);
end

Le mieux serait que tu comprennes ces codes afin que tu puisse les refaire tout seul...

[Leonpolou]

j'ai finalement utilise la premiere methode de Duf, malheureusement j'ai toujours le meme probleme, j'ai les valeurs mais pas de courbes

Code :

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

1
2
3
4
5
6
7
8
9
 
a=3;
o=30;
c=-0.5;
for t=0:0.1:10
    y=max(a*sin(o*t)+c,c)
end
t=0:0.1:10;
plot(t,y);

Et puis rien sur la figure
[EDIT]
Oui il a effectivement une courbe, mais en fait c'est une constante y=c
Pourtant au vu des valeurs y n'est pas une constante

[Jean Dumoncel]

Bon, alors duf42 te propose une méthode sans boucle, et toi tu en rajoute une!

Mais bon, à la limite cela te permet de comprendre le principe.
Il y a une erreur de programmation dans ton code : dans la boucle for, à chaque itération, tu viens écraser la valeur de y, donc au final tu auras "y = à la dernière valeur".

Pour éviter cela il faut placer les différentes valeurs de y dans un tableau. Je te laisse chercher un peu pour faire cela.

Une deuxième méthode consiste à ne pas utiliser de boucle : le calcul de y se fait une ligne.

[Leonpolou]

Oui ca y'est j'ai trouve, effectivement j'ai créer un tableau (n,1)

Merci beaucoup

[Jean Dumoncel]

Et maintenant
le même code, mais sans boucle :

Code :

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

1
2
3
4
5
6
a=3;
o=30;
c=-0.5;
t=0:0.1:10;
y=max(a*sin(o*t)+c,c)
plot(t,y);

[Leonpolou]

Re bonjour,


alors je suis en train de ploter une differente fonction finalement seulement j'ai un probleme avec cette indexation dans un tableau (car cette fois ci la boucle est indispensable)

par exemple

Code :

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

1
2
3
4
 
for t=0:0.1:10
X(t)=1;
end

ca me renvoi un message d'erreur!
ou est-ce que j'ai fais une erreur?

Merci
jeremy

[Jerome Briot]

Voir la : Subscript indices must either be real positive integers or logicals.

[Leonpolou]

Merci pourtant j'ai bien essaye de remplir ces conditions avec ca

Code :

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

1
2
3
4
5
6
 
X=[0:0.1:10];
for a=0:0.1:10
b=a*10+1;
X(b)=2;
end

Je comprends pas pourquoi ca ne marche toujours pas

[Jean Dumoncel]

Salut,

dans ton message d'erreur il y a écrit cela :
??? Subscript indices must either be real positive integers or logicals.

Tu peux forcer b à être un entier :

Code :

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

1
2
3
4
5
X=[0:0.1:10];
for a=0:0.1:10
b=a*10+1
X(round(b))=2;
end

[Jean Dumoncel]

Et pourquoi ne pas faire directement :

Code :

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

1
2
3
4
X=[0:0.1:10];
for a=1:101
X(a)=2;
end

et sans boucle :
X(1:101)=2;

[Leonpolou]

Merci beaucoup, ca marche comme ca, seulement je crois que j'ai un probleme dans l'algo,

je veux une sinusoide d'amplitue a qui croit deux fois plus vite qu'elle ne decroit(o1=02/3 avec o la pulsation propre) avec toujours un offset -c

Code :

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

1
2
3
4
5
 
si cos(o1*t)>0
y=(max(a*sin(o1*t)-c,-c)
sinon
y=(max(a*sin(o1*t)-c,-c)

J'obtien une droite affine partant de -c en x=0

Quelqu'un aurait une petite idee? ( Peut etre que cela est plus une question mathematique, dans ce cas la desole)

jeremy