半波傅氏算法的matlab实现

Gerneral888

%半波傅氏算法,数据窗10ms
N=12;
h=input('请输入谐波次数(1 or 3)：');
z=1+sqrt(-1)*sqrt(3);            
+ ?7 {% g4 W0 M& ]; lt=0:0.02/N:0.02;
0 v# X4 _# m4 j0 X" c u= sin(2*pi*50*t)+2*sin(2*3*pi*50*t-pi/6);  
, s  P+ P% }5 K: m) H  [  t1=0:0.0001:0.07;   
+ H) a& R; T5 d9 S- g# g. B   u1=sin(2*pi*50*t1)+2*sin(2*3*pi*50*t1-pi/6);               
: E, L7 R6 }1 P& y1 k       subplot(211);  plot(t1,u1);title('连续图');
       subplot(212);  stem(t,u);title('离散图');
3 H4 @  H1 e& ?+ h$ H5 M2 C                                   i=u./z;  
  y# y1 `/ r8 b' D) Z                              
                S=0;C=0;
        for  k=1:N/2            
# H; ?: z# ^+ L% H  ^. C5 L                             S=S+u(k)*sin(h*(k-1)*2*pi/N);   
                             C=C+u(k)*cos(h*(k-1)*2*pi/N);
8 u3 s$ {) w4 Q' k6 P* D                        
        end
                                 Us=S*4/N;
                                Uc=4*C/N;
                               U=abs(Us+sqrt(-1)*Uc)/sqrt(2);
                             JU=angle(Us+sqrt(-1)*Uc);
               S=0;C=0;
         for  k=1:N/2
- U2 \% ^* C, V2 }( `7 e! P                                S=S+i(k)*sin(h*(k-1)*2*pi/N);
                              
                                C=C+i(k)*cos(h*(k-1)*2*pi/N);
& D$ j3 o  h: Z                           
         end
+ S8 ]& u3 @5 F& x/ A2 v: P                                Is=S*4/N;
                             Ic=4*C/N;
                          I=abs(Is+sqrt(-1)*Ic)/sqrt(2);
# q4 a# N: {/ M* B% F7 k                         JI=angle(Is+sqrt(-1)*Ic);
U                        
I
cos(JU-JI)

lixintmc

谢了，哈哈，