%partie enregistrement
% paramètres de l'enregistrement d'un son et de son écoute
Fe = 22050; % fréquence d'échantillonnage
n_bits = 16; % nombre de bits
n_channels = 1; % nombre de canaux
d = 2; % durée d'enregistrement
% création d'une variable objet associée à l'enregistrement sonore
x = audiorecorder(Fe,n_bits,n_channels);
% enregistrement d'un son pendant d secondes
recordblocking(x,d);
% écoute du son enregistré
play(x);
% conversion de l'objet x en un vecteur de valeurs
s = getaudiodata(x);
% création d'une variable objet associée à un vecteur d'échantillons s
y = audioplayer(s,Fe);
% écoute du son y
play(y);
% Visualisation de la partie isolée
temps = (0:length(s)-1) / Fe;
figure;
plot(temps*1000, s);
% la partie isolation du signal + visualisation
debut = 7570;
fin = 8992;
u = s(debut:fin);
tu = (0:length(u)-1) / Fe;
figure;
plot(tu*1000, u);
% la partie sauvgardage
save('signal_2.mat', 'u', 'Fe');
%partie upload et application sur enregistrement
load("signal_2.mat");
whos u Fe
t=(0:length(u)-1)*1000/Fe;
y = audioplayer(u,Fe);
% écoute du son y
play(y)
figure('Name','a');
plot(t,u);
d_ordre = 20;
[Den,Pe]=iden(u,d_ordre);
%partie analyse du voyelle
abs_poles=abs(poles);
figure('Name','b');
freqz(1,Den,1024,Fe);
figure('Name','c');
impz(1,Den);
poles= roots(Den);
figure('Name','d');
zplane(1,Den);
%partie FFT
U_fft = fft(u);
U_fft = U_fft(1:floor(N/2)); % partie positive
f = (0:length(U_fft)-1) * Fe / N;
A_mag = abs(U_fft);
figure('Name','e');
plot(f, 20*log10(A_mag));
xlabel('Fréquence (Hz)');
ylabel('Magnitude (dB)');
title('Spectre du voyelle FT');
grid on;
xlim([0 Fe/2]);
%partie synthese
Delta_n= round(T*Fe);
A=sqrt(Pe*Delta_n);
e=zeros(N, 1);
i=1:Delta_n:N;
e(i)=A;
syn = filter(1,Den,e);
figure('Name','f');
plot(t,syn);
plat(syn)