Discussion :

[cedric06000]

bonjour,
j'ai un signal carré non périodique (figure pièce jointe) échantillonné à une fréquence 44,1 khz.
comment je peux faire pour enlever les "0" qui se trouve entre les "1"?
ça veut dire avoir un signal continu entre le premier "1" et le dernier "1".
sachant que si j'ai un "1" entre 500 "0" à droite et 500 "0" à gauche, je considère que c'est un "0" et pas un "1".
est ce que je doit faire un filtre ou calculer une somme ou calculer une moyenne?
merci.
je travail sur scilab.

[cedric06000]

comment je calcule le nombre des "1" pour chaque 1000 ech de signal

[phryte]

Bonjour.
Un exemple pour éliminer un 0 entre deux 1 :

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
N=100;
x=ones(1,N).*floor(2*rand(1,N))
%plot(1:N,x,'.','markersize',15)
plot(1:N,x)
axis([0 N 0 1.5])
grid
for k=2:N-1
  if x(k-1) ==1 & x(k) ==0 & x(k+1)==1
      x(k)=1;
  end
end
figure(2)
plot(1:N,x)
grid
axis([0 N 0 1.5])

[cedric06000]

bonjour phryte,
est ce que vous pouvez faire une fenetre de 1000 échantillon (pour un signal de 50000 échantillon) qui se décale de debut de signal jusqu'à la fin en calculant a chaque 1000 échantillon le nombre des "1" et des "0"?

[cedric06000]

quelq'un d'autre à une idée?

[phryte]

le nombre des "1" et des "0"?

Voici un exemple basique pour 10 et 500 :

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
clear
N=500;
x=ones(1,N).*floor(2*rand(1,N));
plot(1:N,x)
axis([0 N 0 1.5])
grid
H=10;
pas=1;
for k=1:H:N
    nbzero(pas)=0;
    nbun(pas)=0;
     for j=k:k+H-1
if x(j) == 0
    nbzero(pas)=nbzero(pas)+1;
else
    nbun(pas)=nbun(pas)+1;
end
 
end
pas=pas+1;
end
nbun
nbzero

[cedric06000]

exactement ce que je veux
merci

[Caro-Line]

AXIS est une commande MATLAB qui n'existe pas sous Scilab.
Je te donne l'aide de MATLAB dans l'utilisation qu'en fait phryte, peut-être connais-tu un équivalent ? (enfin sinon c'est juste pour faire plus joli, tu peux essayer en mettant cette ligne en commentaire) :

AXIS Control axis scaling and appearance.
AXIS([XMIN XMAX YMIN YMAX]) sets scaling for the x- and y-axes
on the current plot.

[phryte]

Pardon, je me croyais sous matlab.
Il faut :

Code :

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

1
2
mtlb_axis([0 N 0 1.5])
xgrid()

[cedric06000]

probleme résolut.
merci à vous

[Jerome Briot]

Voici un exemple basique pour 10 et 500 :

Ou directement :

Code :

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

1
2
3
4
5
6
7
8
9
10
N=500;
m = 10;
x=ones(1,N).*floor(2*rand(1,N));
 
n = N/m;
 
x = reshape(x,m,n);
 
nbun = sum(x);
nbzero = m-nbun;

[phryte]

Bonjour.

Code :

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

1
2
3
4
5
6
n = N/m;
 
x = reshape(x,m,n);
 
nbun = sum(x);
nbzero = m-nbun;

C'est parfait Dut. Je m'en veux de ne pas y avoir pensé.