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
| clear all
close all
Fs = 32768; % Sampling frequency
T = 1/Fs; % Sample time
L = 2048; % Length of signal
t=(0:L-1)*T; % 2048*dt Time vector
f_array = 1:0.1:3200;
SimSig = zeros(length(f_array) , L); % Preallocation
for Index_f = 1 : length (f_array)
Phi = rand; % Random Phase for each frequencies
f = f_array(Index_f); % All frequencies required
SimSig(Index_f,:) = sin(2*pi*f*t+Phi); % Fill Preallocation
end
SimSigFin=sum(SimSig); % Sum of all frequencies
figure('name','Signal','numbertitle','off');
n = 1:length(f_array);
plot(t, SimSig(n,:))
figure('name','Full FFT','numbertitle','off');
FFT_size=2048;
X=fft(SimSigFin,FFT_size);
f=(0:FFT_size/2)/(FFT_size/2)*Fs/2;
plot(f,abs(X(1:FFT_size/2+1)));
title('FFT - Pré');
xlabel('Fréquence(Hz)');
ylabel('Amplitude'); |
Partager