clear all
close all
N=512;
n=1:1:N;
wt=((n-1)/N)*2*pi;
A = xlsread('work_1_1.csv');
ua=A(:,2);
ub=A(:,3);
uc=A(:,4);
uafft=(fft(ua,N))*2/N;
ubfft=(fft(ub,N))*2/N;
ucfft=(fft(uc,N))*2/N;
Maga=abs(uafft(2));Anga=angle(uafft(2));
Magb=abs(ubfft(2));Angb=angle(ubfft(2));
Magc=abs(ucfft(2));Angc=angle(ucfft(2));
U=[1 exp(pi*2/3*i) exp(pi*4/3*i);1 exp(pi*4/3*i) exp(pi*2/3*i);1 1 1]*[uafft(2);ubfft(2);ucfft(2)]/3;
degree=abs(U(2))/abs(U(1));
t=n/N*0.02;
figure(1)
plot(t,ua,t,ub,t,uc),title('三相电压'),xlabel('t/s'),ylabel('U/V')
评论1