MATLAB - Fast Fourier Transform FFT

[edomar]

Ho trovato in rete un codice che funziona bene:

Codice: Seleziona tutto

Fs = 1000; % Sampling frequency                     

frq_res = 0.01; %desired frequency resolution

t=0 : 1/Fs : 1/frq_res-1/Fs; %time

X = 220*sin(2 * pi * 50 * t);

x1 = 30*sin(2*pi*100*t + 30*(pi/180));
x2 = 10*sin(2*pi*200*t + 50*(pi/180));
x3 = 05*sin(2*pi*300*t + 90*(pi/180));

%% adding the harmonics
Xt = X + x1 + x2 + x3;

%Transformation
Y=fft(Xt); %FFT

df=Fs/length(Y); %frequency resolution

f=(0:1:length(Y)/2)*df; %frequency axis

subplot(2,1,1);
M=abs(Y)/length(Xt)*2; %amplitude spectrum

M_rounded = int16(M(1:size(f, 2))); %Limit the frequency range
ind = find(M_rounded ~= 0);

stem(f(ind), M(ind), 'LineWidth', 0.5);
xlim([0 350]);
grid on; 

xlabel('Frequency (Hz)')
ylabel('Magnitude');

subplot(2,1,2);
P=angle(Y)*180/pi; %phase spectrum (in deg.)
stem(f(ind), P(ind), 'LineWidth', 0.5);
xlim([0 350]);
ylim([-100 100]);
grid on;

xlabel('Frequency (Hz)');
ylabel('Phase (degree)');



Il problema è che come in tutti gli esempi che si trovano, la frequenza di campionamento e, quindi, l'asse dei tempi sono definiti appositamente. Ho provato allora a definire l'asse dei tempi come nel mio programma:

Codice: Seleziona tutto


theta = 0:0.001:2*pi;
omega = 2*pi*50;
t = theta/omega;

Ts = t(34) - t(33);
Fs = 1/Ts;



E il risultato rimane ancora lo stesso:

Poi vado a sostituire nel mio script, cambiando solo Xt con la corrente i1 che voglio analizzare e non funziona più niente...

A questo punto non può che essere un problema di campionamento, di scelta della frequenza di campionamento, ecc ecc... Non riesco a raccapezzarmi

[xyz]

Ripeto, non ha senso scrivere una cosa del genere:

Codice: Seleziona tutto

Ts = t(34) - t(33);

Il tempo di campionamento è un dato fisso e deve essere una costante, ad esempio:

Codice: Seleziona tutto

Ts = 0.001

o imponi la frequenza di campionamento (il reciproco del tempo), ad esempio:

Codice: Seleziona tutto

Fs = 1000

Ricavi l'asse del tempo in base a Ts o Fs, dipende quale valore definisci del codice, come è correttamente fatto nell'esempio che hai trovato funzionante il quale definisce costante la frequenza di campionamento Fs.

[EcoTan]

la corrente i1 che voglio analizzare

non potrebbe darsi che coi tuoi campioni la FFT sia davvero quella?

[edomar]

No, guarda l'andamento della i1 nel mio primo post... Sto provando a modificare il codice a monte, partendo dalla definizione del tempo di campionamento... Vediamo cosa viene fuori

[edomar]

Ripeto, non ha senso scrivere una cosa del genere:

Codice: Seleziona tutto

Ts = t(34) - t(33);

....

Ho provato a fare come hai detto, imponendo a monte la frequenza di campionamento:

Codice: Seleziona tutto

Fs = 10000;

Ts = 1/Fs;

t = 0 : 1/Fs : 2*pi/(2*pi*50) - 1/Fs;

theta = 2*pi*50*t;


Il risultato però rimane invariato... Anche facendo distorcere ulteriormente la i1 il risultato rimane quello

[xyz]

Nel tuo codice esistono parti oscure come questa:

Codice: Seleziona tutto

2*pi/(2*pi*50)

non si capisce il motivo di un calcolo del genere visto che è uguale a 1/50, dovrebbe indicare il tempo massimo.

[edomar]

Nel tuo codice esistono parti oscure come questa:

Codice: Seleziona tutto

2*pi/(2*pi*50)

non si capisce il motivo di un calcolo del genere visto che è uguale a 1/50, dovrebbe indicare il tempo massimo.

Sì esatto, è il tempo massimo... in effetti è inutile scriverlo in questo modo, bastava scrivere 1/50
Però il problema rimane..

[EcoTan]

distorcere ulteriormente la i1

Sarebbe più chiaro fornire, oltre al grafico, anche i valori dei campioni.

[edomar]

Potrei copiare il link del medesimo quesito posto nel forum di Matlab (anche in quel caso senza soluzione) in cui è possibile allegare i file.mat... Ma non so se questa operazione sia in linea con le regole del forum

Intanto io continuo a fare mille tentativi ma non riesco a pervenire a una soluzione...

Alla fine io vorrei trovare l'ampiezza della fondamentale e il THD... Non esistono altri modi? Anche se sarebbe comunque bello trovare gli spettri in ampiezza e fase..