MATLAB CSV读入与分析

高铁是我国工程技术的名片，电力系统为高铁提供动力。某高铁牵引变电站由220 kV高压电缆供电，为研究高铁的电能质量问题需要统计线芯电流谐波含量，A相的线芯电流存储在通道1_[2019年10月27日 14时13分56秒].cvs文件中，B相的线芯电流存储在通道2_[2019年10月27日 14时13分56秒].cvs文件中，C相的线芯电流存储在通道3_[2019年10月27日 14时13分56秒].cvs文件中，采样率为20 kS/s，电力系统的频率为50 Hz。题目要求同学们在一张图上画出三个周期的电流波形，ABC三相用黄绿红色分别表示，线宽为2 pt，画出图例。分别使用快速傅里叶变换做出A、B、C三相的频谱图（3张），横坐标为谐波次数，纵坐标为谐波电流含量百分比（谐波电流/基波电流）。（20分）

报告要求：三个周期的电流选取为第（学号后两位）个周波；谐波分析仅需分析三个周波中的第一个周波。

%2022.6.15

%by 秦舞阳

n=input('请输入你的学号后两位');
filename1='D:\QQ\FileRecv\MATLAB\通道1_[2019年10月27日 14时13分51秒].csv';
data1 = csvread(filename1);

filename2='D:\QQ\FileRecv\MATLAB\通道2_[2019年10月27日 14时13分51秒].csv';
data2 = csvread(filename2);

filename3='D:\QQ\FileRecv\MATLAB\通道3_[2019年10月27日 14时13分51秒].csv';
data3 = csvread(filename3);

figure(1)
plot(data1(1200*(n-1):1200*n,1),data1(1200*(n-1):1200*n,2),'y','linewidth',2)
hold on
plot(data2(1200*(n-1):1200*n,1),data2(1200*(n-1):1200*n,2),'g','linewidth',2)
hold on
plot(data3(1200*(n-1):1200*n,1),data3(1200*(n-1):1200*n,2),'r','linewidth',2)
hold on
legend('A','B','C')
xlabel('采样点')
ylabel('电流')

figure(2)
Y = fft(data1(:,2));
plot(data1(n*400:n*400+400,1),Y(n*400:n*400+400)); 
title('Single-Sided Amplitude Spectrum of A')
xlabel('f (Hz)')
ylabel('|P1(f)|')
figure(2)
Y = fft(data2(:,2));
plot(data2(n*400:n*400+400,1),Y(n*400:n*400+400))
title('Single-Sided Amplitude Spectrum of B')
xlabel('f (Hz)')
ylabel('|P1(f)|')

figure(3)
Y = fft(data3(:,2));
plot(data3(n*400:n*400+400,1),Y(n*400:n*400+400))
title('Single-Sided Amplitude Spectrum of A')
xlabel('f (Hz)')
ylabel('|P1(f)|')