function[hmr_i]=FFT_DAT(ch1);

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

%FFT計(jì)算頻譜；
a1=abs(fft(a))/100;
a2=a1/a1(2)*100;
a2(1)=0;

%計(jì)算2-25次諧波的總畸變率
hmr=0;
for i=3:26
    hmr=hmr+a2(i)^2;
end
hmr_i=sqrt(hmr);

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

%繪制原始數(shù)據(jù)波形圖；
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'),