Discussion :

[membreComplexe12]

Bonjour tous,

j'ai fais un petit programme d'une façon très bête pour effectuer des sinusoides interrompues.
Je voudrais réécrire ce programme d'une façon beaucoup plus intelligente avec des boucles ou autre mais je n'y arrive pas, j'espere que vous pourrez m'aider

1°) ce que je veux faire
-je veux faire une premiere sinusoide de frequence 0,2 pendant un temps "t", (pour cela pas de problème)
-ensuite je veux faire une courbe comprenant: une sinusoide de frequence 1 Hz, une partie de courbe nulle, une sinusoide de frequence 1hz, une partie de courbe nulle...
-de même que la précédente mais avec 5Hz

2°) contrainte à respecter
je veux pour les courbe 2 et 3 que les temps où démarre la sinusoide soient les meme temps que pour la courbe 1

3°) mon programme
le programme ci dessous permet de faire ce que je veux mais il n'est pas du tout optimisé car je ne vois pas comment faire...

4°) voici mon programme

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
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
 
wa=2*pi*0.2;
ta1=[2.5:0.1:7.5];
ta2=[7.5:0.1:12.5];
ta3=[12.5:0.1:17.5];
ta4=[17.5:0.1:22.5];
a1=30*sin(wa*ta1);
a2=30*sin(wa*ta2);
a3=30*sin(wa*ta3);
a4=30*sin(wa*ta4);
ta=[ta1,ta2,ta3,ta4];
a=[a1,a2,a3,a4];
figure(1)
plot(ta,a)
%------------------------------------------
wb=2*pi*1;
tb1=[2.5:0.1:3.5];
tb2=[3.5:0.1:7.5];
tb3=[7.5:0.1:8.5];
tb4=[8.5:0.1:12.5];
tb5=[12.5:0.1:13.5];
tb6=[13.5:0.1:17.5];
tb7=[17.5:0.1:18.5];
tb8=[18.5:0.1:22.5];
b1=30*sin(wb*tb1);
b2=zeros(1,length(tb2));
b3=30*sin(wb*tb3);
b4=zeros(1,length(tb4));
b5=30*sin(wb*tb5);
b6=zeros(1,length(tb6));
b7=30*sin(wb*tb7);
b8=zeros(1,length(tb8));
tb=[tb1,tb2,tb3,tb4,tb5,tb6,tb7,tb8];
b=[b1,b2,b3,b4,b5,b6,b7,b8];
figure(2)
plot(tb,b)
%------------------------------------------
wc=2*pi*5;
tc1=[2.5:0.01:2.7];
tc2=[2.7:0.01:7.5];
tc3=[7.5:0.01:7.7];
tc4=[7.7:0.01:13.5];
tc5=[13.5:0.01:13.7];
tc6=[13.7:0.01:18.5];
tc7=[18.5:0.01:18.7];
tc8=[18.7:0.01:22.5];
c1=30*sin(wc*tc1);
c2=zeros(1,length(tc2));
c3=30*sin(wc*tc3);
c4=zeros(1,length(tc4));
c5=30*sin(wc*tc5);
c6=zeros(1,length(tc6));
c7=30*sin(wc*tc7);
c8=zeros(1,length(tc8));
tc=[tc1,tc2,tc3,tc4,tc5,tc6,tc7,tc8];
c=[c1,c2,c3,c4,c5,c6,c7,c8];
figure(3)
plot(tc,c)

[Jerome Briot]

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
wa=2*pi*0.2;
ta1=[2.5:0.1:7.5];
ta2=[7.5:0.1:12.5];
ta3=[12.5:0.1:17.5];
ta4=[17.5:0.1:22.5];
a1=30*sin(wa*ta1);
a2=30*sin(wa*ta2);
a3=30*sin(wa*ta3);
a4=30*sin(wa*ta4);
ta=[ta1,ta2,ta3,ta4];
a=[a1,a2,a3,a4];
figure(1)
plot(ta,a)

est équivalent à :

Code :

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

1
2
3
4
5
6
wa = 2*pi*0.2;
ta = 2.5:0.1:22.5;
a = 30*sin(wa*ta);
 
figure
plot(ta,a)

A quelques répétitions près... je ne comprends pas vraiment ce que tu cherches à faire

[membreComplexe12]

est équivalent à :

Code :

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

1
2
3
4
5
wa = 2*pi*0.2;
ta = 2.5:0.1:22.5;
a = 30*sin(wa*ta);
figure
plot(ta,a)

A quelques répétitions près... je ne comprends pas vraiment ce que tu cherches à faire

oui, merci DUT. En fait pour ce cas là j'avais remarqué mais pour les autres je ne sais pas comment faire
(en fait pour le "a" j'ai gardé la même structure que les autres juste pour expliquer ce que je cherche à faire)

[tubaas]

salut
je te propose ça

Code :

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

1
2
3
4
5
6
7
8
9
10
11
f=[0.2;1;5];
N=4; % nombre de périodes
dt=0.01;
ti=2.5;
t=ti:dt:ti+N/f(1)-dt;
a=zeros(3,length(t));
a(1,:)=sin(2*pi*f(1)*t);
a(2,:)=repmat([sin(2*pi*f(2)*t(1:1/(dt*f(2)))) zeros(1,1/(dt*f(1))-1/(dt*f(2)))],1,N);
a(3,:)=repmat([sin(2*pi*f(3)*t(1:1/(dt*f(3)))) zeros(1,1/(dt*f(1))-1/(dt*f(3)))],1,N);
figure;
plot(t,a);

mais tu vois qu'il faut un nombre entier de périodes.
Sinon tu peux toujours placer tes "impulsions" où il faut mais ce peut-être un peu plus lourd

[membreComplexe12]

merci beaucoup d'avoir pris le temps de repondre.

pourrais tu m'expliquer ce que fait repmat car j'ai regardé sur le net et je n'ai pas vraiment saisi

[tubaas]

repmat(A,i,j) réplique une matrice A i fois verticalement et j fois horizontalement.
Ici je m'en sers car je construis le signal sur une période du premier signal et je le duplique pour le bon nombre de périodes. C'est à cause de ça que je dois avoir un nombre de périodes entier.

Edit : une autre méthode. Je ne sais pas si c'est très efficace mais c'était rigolo.

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
17
18
function creaSignal
f=0.2;
df=[1;5;18];
dt=0.01;
ti=2.5;
t=ti:dt:ti+1366*dt;
a=arrayfun(@(x,y)partialPeriod(x,y),repmat(t,length(df),1),repmat(df,1,length(t)));
figure;
plot(t,a);
 
    function a=partialPeriod(x,y)
        if mod(floor((x-ti)*y*f),y)==0
            a=sin(2*pi*f*y*x);
        else
            a=0;
        end
    end
end

ici tu vois que j'introduis un vecteur df qui correspond à un multiple par rapport à la fréquence de base. Il doit être entier.

[membreComplexe12]

j'ai un peu de mal avec cette ligne:

Code :

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

a=arrayfun(@(x,y)partialPeriod(x,y),repmat(t,3,1),repmat(df,1,length(t)));

si je comprends bien arrayfun est une fonction et entre parenthèse on a les arguments? c'est quoi ce truc bizarre :

Code :

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

@(x,y)partialPeriod(x,y)

merci beaucoup tubaas pour ton aide!!!

[tubaas]

arrayfun applique une fonction à chaque élément d'un tableau (ou ici de deux tableaux)
dans la fonction appelée x va correspondre à un élément du tableau dans lequel j'ai le temps, l'élément y à mon élément multiplicateur.

c'est équivalent à

Code :

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

1
2
3
4
5
6
7
8
a=zeros(length(df),length(t));
for i=1:length(df)
    for j=1:length(t)
        if mod(floor((t(j)-ti)*df(i)*f),df(i))==0
            a(i,j)=sin(2*pi*f*df(i)*t(j));
        end
    end
end

edit : le truc bizarre c'est parce que la Matlab va automatiquement passer deux argument à arrayfun car j'ai deux tableaux, je les récupère avant d'appeler la fonction et je les repasse à la fonction. cette notation ne sert à rien ici on peut se contenter de @partialPeriod

[membreComplexe12]

merci beaucoup de ton aide

(je vais conserver la premiere solution c'est plus facile pour moi)