二次系统课程设计(4)

参考文献

[1]《配电网分析及应用》何正友中国科学出版社

[2]《广域同步相量测量技术及其应用》 A.G. Phadke著毕天姝译中国电力出版社

[3]《配电系统及其自动化技术》陈堂赵祖康中国电力出版社 [4]《城市电网规划与改造》陈章潮唐德光中国电力出版社 [5]《配电系统自动化及其发展》王明俊于尔铿中国电力出版社

[6]《配电网络规划与设计》 E.Lakervi等著，范明天等译中国电力出版社

附录

去噪算法程序：

A1=10; A2=0.3; P1=30; P2=-30; F1=50; Fs=256; N=256; figure

t=0:0.0001:0.1;

u=A1*cos(2*pi*F1*t+pi*P1/180)+rand(1,length(t)); plot(t,u)

title('电压原始信号'); figure

t=0:0.0001:0.1;

i=A2*cos(2*pi*F1*t)+(rand(1,length(t)))/20; plot(t,i)

title('电流原始信号'); figure

t1=[0:1/Fs:N/Fs];

u1=A1*cos(2*pi*F1*t1+pi*P1/180)+rand(1,length(t1)); Y1 = fft(u1,N); Ayy1 = (abs(Y1)); Ayy1=Ayy1/(N/2); Ayy1(1)=Ayy1(1)/2; F=([1:N]-1)*Fs/N;

plot(F(1:N/2),Ayy1(1:N/2)); title('电压幅度-频率曲线图'); Ayy1(51) figure

t1=[0:1/Fs:N/Fs];

i1=A2*cos(2*pi*F1*t1)+(rand(1,length(t1)))/20; Y2 = fft(i1,N); Ayy2 = (abs(Y2)); Ayy2=Ayy2/(N/2); Ayy2(1)=Ayy2(1)/2; F=([1:N]-1)*Fs/N;

plot(F(1:N/2),Ayy2(1:N/2)); title('电流幅度-频率曲线图'); Ayy2(51)

figure

Pyy1=[1:N/2]; for i=1:N/2

Pyy1(i)=phase(Y1(i)); Pyy1(i)=Pyy1(i)*180/pi; end;

plot(F(1:N/2),Pyy1(1:N/2)); title('电压相位-频率曲线图'); Pyy1(51) figure

Pyy2=[1:N/2]; for i=1:N/2

Pyy2(i)=phase(Y2(i)); Pyy2(i)=Pyy2(i)*180/pi; end;

plot(F(1:N/2),Pyy2(1:N/2)); title('电流相位-频率曲线图') Pyy2(51) figure

t=0:0.0001:0.1;

u3=Ayy1(51)*cos(2*pi*F1*t+Pyy1(51)) plot(t,u3)

title('电压基波信号'); figure

t=0:0.0001:0.1;

u4=Ayy2(51)*cos(2*pi*F1*t+Pyy2(51)) plot(t,u4)

title('电流基波信号'); S=Ayy1(51)*Ayy2(51)/2

P=Ayy1(51)*Ayy2(51)*cos((Pyy1(51)-Pyy2(51))*pi/180)/2 Q= Ayy1(51)*Ayy2(51)*sin((Pyy1(51)-Pyy2(51))*pi/180)/2