THD calculation using MATLAB
显示 更早的评论
Hello,
I'm a novice in using MATLAB. I want to calculate THD value from my experimental data using MATLAB. I have 10,001 data points, sampling frequency is 100kHz and fundamental frequency is 50 Hz (I'm not in USA). Below is my code. Could you please check that my code is correct or not ?
>> data = importdata('sample_1.xlsx'); % 10,001 data points, sampling rate = 100kHz
>> spect = fft(data)/10001; % frequency of each element = 0, 10, 20, 30,...Hz (100kHz/10,001 data)
>> harmonic = abs(spect(11:5:201)); % here, fundamental = 50 Hz so, frequency 100, 150, 200,....,2000 Hz correspond to 11th, 16th, 21st,....,201st element of 'spect' and 'harmonic' has 39 elements
>> harmonic_sqr = harmonic.^2; % square each element
>> THD = sqrt(sum(harmonic_sqr(1:39)))/abs(spect(6))*100 % 50 Hz corresponds to 6th element of 'spect'
回答(0 个)
类别
在 帮助中心 和 File Exchange 中查找有关 Spectral Measurements 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!