Why thd() function captures incorrect harmonics?

6 visualizaciones (últimos 30 días)
Ying Wu
Ying Wu el 4 de Oct. de 2022
Comentada: Ying Wu el 7 de Oct. de 2022
Hi, I use thd() to calculate one time-domain singal, and I try three commands but all fails. Can anyone help me to solve this question? Thanks!
load signal.mat
% 1. time-domain signal as inputs
% error: the fundemental freq is captured, but the other five harmonics are not correctly captured)
thd(signal,100,6); xlim([0 0.3])
% 2. psd as inputs (same error)
[pxx,f] = periodogram(signal,rectwin(length(signal)),length(signal),100);
thd(pxx,f,6,'psd'); xlim([0 0.3]);
% 3. power spectural as inputs
[sxx,f] = periodogram(signal,rectwin(length(signal)),length(signal),100,'power');
rbw = enbw(rectwin(length(signal)),100);
thd(sxx,f,rbw,'power'); xlim([0 0.3]);

Respuestas (1)

Geetla Sindhu
Geetla Sindhu el 7 de Oct. de 2022
Editada: Geetla Sindhu el 7 de Oct. de 2022
Hello Ying Wu,
I understand that you are trying to calculate THD of a signal using three approaches.
It would be helpful if you could elaborate more on how you are assuring that the results obtained through all the approaches are wrong. This is because we cannot manually calculate the THD for such a signal.
Thank you.
  1 comentario
Ying Wu
Ying Wu el 7 de Oct. de 2022
@Geetla Sindhu Hi Geetla, sorry about the unclear description of this question. To my knowledge, the definition of thd() in matlab is the ratio of sum of all harmonics amplitude to that of main frequency. In my case, I set the thd() to capture 6 frequencies including one main frequency (marked as "F" in blue in the figure) and 5 harmonics (marked from 2-5 in the orange in the figure). But obviously, the harmonics are not correctly captured like NO.2, 4, 6 are not harmonic peaks, while the real harmonics are not captured. So the calculated thd is also incorrect. Could you please help me to solve this problem?

Iniciar sesión para comentar.

Categorías

Más información sobre Spectral Measurements en Help Center y File Exchange.

Productos


Versión

R2022a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by