Rezime

Spectre (transformée de Fourier) d'un signal bruité

Spectrum (Fourier transform) of a signal with noise

Auteur/author : Christophe Dang Ngoc Chan (cdang)

Réalisé avec/made with Scilab

clear;clf;

// ******************
// variable ajustable

p=0.1; // pas

// ******************
// variables de base

x=[0:p:20*%pi]'; // x du sinus
N=size(x);
n=N(1,1); // nombre de points
s=1/p*(1:n); // fréquences

// nombres aleatoires, loi uniforme

y2=0.2*(rand(n,1)-0.5);

// ******************
// Calcul

y1=sin(x);
y3=y1+y2;
z1=abs(fft(y1))^2;
z2=abs(fft(y2))^2;
z3=abs(fft(y3))^2;

for i=1:n
  if z1(i)>500 then z1(i)=500; end
  if z3(i)>500 then z3(i)=500; end
end

// ******************
// affichage

subplot(3,2,1)
plot2d(x,y1,style=2,rect=[0,-1.1,60,1.1])
xtitle('signal sinusoïdal','Temps','Intensité')

subplot(3,2,2)
plot2d(s,z1,style=2,rect=[0,0,0.03,50])
xtitle('TF signal sinusoïdal','Fréq.','Puissance')

subplot(3,2,3)
plot2d(x,y2,style=2,rect=[0,-1.1,60,1.1])
xtitle('bruit','Temps','Intensité')

subplot(3,2,4)
plot2d(s,z2,style=2,rect=[0,0,0.03,50])
xtitle('TF bruit','Fréq.','Puissance')

subplot(3,2,5)
plot2d(x,y3,style=2,rect=[0,-1.1,60,1.1])
xtitle('signal sinusoïdal bruité','Temps','Intensité')

subplot(3,2,6)
plot2d(s,z3,style=2,rect=[0,0,0.03,50])
xtitle('TF signal sinusoïdal bruité','Fréq.','Puissance')

Kondisyon itlizasyon

Moi, en tant que détenteur des droits d’auteur sur cette œuvre, je la publie sous les licences suivantes :

Ou bay pèmisyon pou kopye, distribye ak/oswa modifye dokiman sa a dapre kondisyon GNU Free Documentation License, Vèsyon. 1.2 oswa nenpòt vèsyon ki pibliye pita pa Free Software Foundation; san okenn Seksyon ki pa ka chanje, pa gen okenn Tèks ak premye paj kouvèti, epi pa gen okenn Tèks Denyè-Kouvèti. Gen yon kopi lisans lan nan seksyon ki rele GNU Free Documentation License.http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue

	Fichye sa a distribiye selon lisans Creative Commons.
	Ou gen dwa: pou pataje – pou kopye, distribiye epi difize ouvraj la pou remiks – pou adapte ouvraj la Selon kondisyon sa a yo: Rekonèt avèk non otè – Ou dwe bay enfòmasyon apwopriye sou otè a, bay yon lyen pou lisans lan, epi endike si te gen chanjman ki fèt. Ou ka fè sa nan nenpòt fason rezonab, men okenn fason ki sijere ke otè a soutni oswa apwouve itilizasyon ou an. Enfòmasyon sa yo ka pataje – Si ou fè rekopye, transfòme, oswa baze sou materyèl sa, ou dwe distribye kontribisyon ou anba menm oswa yon lisans konpatib avèk orijinal la.
	Ce bandeau de licence a été ajouté à ce fichier dans le cadre de la procédure de mise à jour des licences des images sous GFDL.http://creativecommons.org/licenses/by-sa/3.0/CC BY-SA 3.0Creative Commons Attribution-Share Alike 3.0truetrue

Chwazi lisans ou vle a.