37 lines
530 B
Matlab
Executable file
37 lines
530 B
Matlab
Executable file
clear;
|
|
|
|
f0 = 1100;
|
|
N = 90;
|
|
|
|
Fe1 = 10000;
|
|
Te1 = 1/Fe1;
|
|
T1 = [0:Te1:(N-1)*Te1];
|
|
Y1 = cos(2*pi*f0 * T1);
|
|
|
|
% beau cosinus
|
|
% plot(T1, Y1);
|
|
|
|
N_prime = 2^10;
|
|
% fft discrète + zero padding
|
|
% semilogy([0:N_prime-1]/N_prime*Fe1, abs(fft(Y1, N_prime)));
|
|
|
|
% chelou
|
|
% periodogram(Y1);
|
|
|
|
% periodogram maison
|
|
% plot( [0:N_prime-1]/N_prime*Fe1, abs(fft(Y1, N_prime)).^2/N_prime );
|
|
|
|
% broken
|
|
% plot(xcorr(Y1));
|
|
|
|
|
|
|
|
Fe2 = 1000;
|
|
Te2 = 1/Fe2;
|
|
T2 = [0:Te2:(N-1)*Te2];
|
|
Y2 = cos(2*pi*f0 * T2);
|
|
|
|
% plot(T2, Y2);
|
|
% semilogy([0:N-1]/N*Fe2, abs(fft(Y2)));
|
|
|