fft_dat.m

使用matlab变成

function[hmr_i]=FFT_DAT(ch1);

for i=1:1000
    a(i)=ch1(i+250);
end

%FFT计算频谱；
a1=abs(fft(a))/100;
a2=a1/a1(2)*100;
a2(1)=0;

%计算2-25次谐波的总畸变率
hmr=0;
for i=3:26
    hmr=hmr+a2(i)^2;
end
hmr_i=sqrt(hmr);

hmr_5=a2(4);%特定次谐波FFT分析

%绘制原始数据波形图；
figure;
subplot(1,1,1);
plot(ch1);grid on;axis([0,2500,-1.5,1.5]);
title(' 电 流 波 形  ');
xlabel('\fontsize{14}\bfωt\rightarrow'),                      
ylabel('\fontsize{14}\bfi\rightarrow'),