dat_m.m

使用matlab变成

hold off;
clear;

%导入数据文件：CH1.DAT(单列的数据文件1×2500个数据)
load ch1.dat;load ch2.dat;load ch3.dat;load ch4.dat;load ch5.dat;load ch6.dat;load ch7.dat;load ch8.dat;
load ch9.dat;load ch10.dat;load ch10.dat;load ch11.dat;load ch12.dat;load ch13.dat;load ch14.dat;load ch15.dat;
load ch16.dat;load ch17.dat;load ch18.dat;load ch19.dat;load ch20.dat;load ch21.dat;

hmr_i25(1)=FFT_DAT(ch1);
hmr_i25(2)=FFT_DAT(ch2);
hmr_i25(3)=FFT_DAT(ch3);
hmr_i25(4)=FFT_DAT(ch4);
hmr_i25(5)=FFT_DAT(ch5);
hmr_i25(6)=FFT_DAT(ch6);
hmr_i25(7)=FFT_DAT(ch7);
hmr_i25(8)=FFT_DAT(ch8);
hmr_i25(9)=FFT_DAT(ch9);
hmr_i25(10)=FFT_DAT(ch10);
hmr_i25(11)=FFT_DAT(ch11);
hmr_i25(12)=FFT_DAT(ch12);
hmr_i25(13)=FFT_DAT(ch13);
hmr_i25(14)=FFT_DAT(ch14);
hmr_i25(15)=FFT_DAT(ch15);
hmr_i25(16)=FFT_DAT(ch16);
hmr_i25(17)=FFT_DAT(ch17);
hmr_i25(18)=FFT_DAT(ch18);
hmr_i25(19)=FFT_DAT(ch19);
hmr_i25(20)=FFT_DAT(ch20);
%hmr_i25(21)=FFT_DAT(ch21);

for i=1:20
    hmr(i)=hmr_i25(21-i);
end

%绘制棒形频谱图
plot(hmr);grid on;
%axis([0,21,0,70]);
title('不同输出电流时3次谐波电流THD'),
%title('不同输出电流时谐波电流总THD'),
xlabel('\fontsize{14}\bfI(A)\rightarrow'),                      
ylabel('\fontsize{14}\bfTHD I (%)\rightarrow'),