基于瞬时无功理论的谐波检测的编程问题
我的程序见我的附件，不知道是不是瞬时无功理论编的有问题，出来波形应该是A相里的3*sin(4*w*t)，是个正弦波，我出来的波形不对，请各位帮忙解决一下。 f=50;w=2*pi*f;
m=2;
t=0:0.025:m;
%模拟的三相电流;
ia=sin(w*t)+3*sin(4*w*t);%A相带有四次谐波
ib=sin(w*t-2/3*pi);
ic=sin(w*t+2/3*pi);
ip=sin(w*t).*ia+sin(w*t-2/3*pi).*ib+sin(w*t+2/3*pi).*ic;
iq=cos(w*t).*ia+cos(w*t-2/3*pi).*ib+cos(w*t+2/3*pi).*ic;
b=fir1(11,0.00025,'low');%低通滤波器
ipz=filter(b,1,ip);%ip的直流分量
iqz=filter(b,1,iq);%iq的直流分量
iaf=(2/3)*sin(w*t).*ipz+(2/3)*cos(w*t).*iqz;
ibf=(2/3)*sin(w*t-2/3*pi).*ipz+(2/3)*cos(w*t).*iqz;
icf=(2/3)*sin(w*t+2/3*pi).*ipz+(2/3)*cos(w*t+2/3*pi).*iqz;
iah=ia-iaf;
plot(t,iah)%显示A相谐波电流的波形
[ 本帖最后由 eight 于 2008-4-28 18:29 编辑 ]
页:
[1]