Tracer une fonction avec conditions sur variable

Bonjour,

je souhaites tout simplement tracer la fonction suivant avec MATLAB

Code:

---

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

Bonsoir,

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

Code:

---

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

---

Pourrait convenir?

Bonne soirée,
Duf

Cela suppose que a est positif, non?

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

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

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:

---

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

---

Maintenant je suis d'accords!:ccool:

Oui a est forcement positif

j'ai essaye ca comme code

Code:

---

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!

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

Code:

---

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:

---

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...

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:

---

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

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.

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

Merci beaucoup

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

Code:

---

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);

---

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:

---

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

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

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

Code:

---

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

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:

---

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

---

Et pourquoi ne pas faire directement :

Code:

---

1 2 3 4

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

---

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

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:

---

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