基于瞬时无功理论的谐波检测的编程问题

开心女孩

我的程序见我的附件，不知道是不是瞬时无功理论编的有问题，出来波形应该是A相里的3*sin(4*w*t)，是个正弦波，我出来的波形不对，请各位帮忙解决一下。

1. f=50;
   
2. w=2*pi*f;
   
3. m=2;
   
4. t=0:0.025:m;
   
5. %模拟的三相电流;
   
6. ia=sin(w*t)+3*sin(4*w*t);%A相带有四次谐波
   
7. ib=sin(w*t-2/3*pi);
   
8. ic=sin(w*t+2/3*pi);
   
9. ip=sin(w*t).*ia+sin(w*t-2/3*pi).*ib+sin(w*t+2/3*pi).*ic;
   
10. iq=cos(w*t).*ia+cos(w*t-2/3*pi).*ib+cos(w*t+2/3*pi).*ic;
    
11. b=fir1(11,0.00025,'low');%低通滤波器
    
12. ipz=filter(b,1,ip);%ip的直流分量
    
13. iqz=filter(b,1,iq);%iq的直流分量
    
14. iaf=(2/3)*sin(w*t).*ipz+(2/3)*cos(w*t).*iqz;
    
15. ibf=(2/3)*sin(w*t-2/3*pi).*ipz+(2/3)*cos(w*t).*iqz;
    
16. icf=(2/3)*sin(w*t+2/3*pi).*ipz+(2/3)*cos(w*t+2/3*pi).*iqz;
    
17. iah=ia-iaf;
    
18. plot(t,iah)%显示A相谐波电流的波形

复制代码

[ 本帖最后由 eight 于 2008-4-28 18:29 编辑 ]