Discussion :

[geogeo4209]

Bonjour à tous !

Je suis nouveau sur ce forum et un novice en Matlab, je sollicite votre aide :

Pour une partie d'un projet sur lequel je travaille, je dois réaliser une application Matlab qui permette d'enregistrer un sinus balayé (on demande la fréquence min, max et la durée de l'enregistrement à l'utilisateur) au format .wav.

Exemple : sur une durée de 10 seconde, je veux faire l'enregistrement d'un sinus qui monte en fréquence de 10 à 1000 Hz.

Etant novice en Matlab, j'aimerais un peu d'aide...

Merci, cdt

Geoffrey'

[Jerome Briot]

Montre nous ce que tu as déjà fait (même si c'est faux)

[geogeo4209]

c'est bon j'ai trouvé
voici en pièce jointe ce que j'ai fais.

[Tanki1001]

Bonjour, ce message arrive un peu tard. Nous avons fait un TD sur Octave (même programmation que MatLab). Je dois faire un sinus Balayé en Java, et je m'aide donc de ce code en Octave, s'il sert aujourd'hui, voilà un sinus balayé très correct (accompagné de sa FFT):

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
% TRAITEMENT DU SIGNAL
% TD - Applications en vibrations des structures
% Emmanuel FOLTETE
 
% ACTIVITE 2 - SIGNAUX
% 2.1 Sinus balayé
 
clear all
close all
 
% Paramètres d'échantillonnage
SamplingFrequency = 40000 ;
SignalDuration = 1 ;
 
NbSamples = SignalDuration*SamplingFrequency;
FrequencyResolution = SamplingFrequency/NbSamples,
Time = linspace(0,SignalDuration,NbSamples);
 
% Sinus balayé
SignalAmplitude = 100 ;
StartFrequency =  100 ;
StopFrequency =   500 ;
 
L2 = NbSamples * StopFrequency / (StopFrequency - StartFrequency);
w = StopFrequency * pi / (SamplingFrequency * L2);
k = L2-NbSamples:L2-1;
Signal = SignalAmplitude*sin(w * k.^2);
 
soundsc(Signal,SamplingFrequency);
 
figure(1);
plot(Time,Signal);
axis([Time(1),Time(end),-SignalAmplitude*1.1,SignalAmplitude*1.1]);
xlabel('Temps (s)');
ylabel('Signal (V)')
title('Sinus balayé','FontWeight','Bold')
 
% Spectre unilatéral
Spectrum = fft(Signal)/NbSamples;
Spectrum = 2*Spectrum(1:NbSamples/2+1);
Frequency = SamplingFrequency/2*linspace(0,1,NbSamples/2+1);
figure(2)
semilogy(Frequency,abs(Spectrum))
xlabel('Fréquence (Hz)')
ylabel('Amplitude (V)')
title('Sinus balayé','FontWeight','Bold')

Bien cordialement, Pierre.
Ce code a été écrit par mon professeur : Emmanuel FOLTETE, professeur en étude des vibrations mécaniques.