半波傅氏算法的matlab实现

Gerneral888

%半波傅氏算法,数据窗10ms
N=12;
/ ?0 o  O2 G% T/ A2 K- }+ k4 D h=input('请输入谐波次数(1 or 3)：');
z=1+sqrt(-1)*sqrt(3);            
t=0:0.02/N:0.02;
u= sin(2*pi*50*t)+2*sin(2*3*pi*50*t-pi/6);  
  t1=0:0.0001:0.07;   
/ L* H) _1 s% s. }7 }8 c   u1=sin(2*pi*50*t1)+2*sin(2*3*pi*50*t1-pi/6);               
       subplot(211);  plot(t1,u1);title('连续图');
& o3 U4 H# w- T) F8 @       subplot(212);  stem(t,u);title('离散图');
                                   i=u./z;  
                              
                S=0;C=0;
9 A- K/ Y/ T6 X  d$ M        for  k=1:N/2            
                             S=S+u(k)*sin(h*(k-1)*2*pi/N);   
9 A" c  \% v7 t2 ?2 y/ d                             C=C+u(k)*cos(h*(k-1)*2*pi/N);
  {$ k  d4 b) u0 _7 \6 E' t                        
        end
                                 Us=S*4/N;
' I2 m  n( `% ~  i3 B                                Uc=4*C/N;
                               U=abs(Us+sqrt(-1)*Uc)/sqrt(2);
7 X4 z/ R4 F$ v, [& K                             JU=angle(Us+sqrt(-1)*Uc);
               S=0;C=0;
         for  k=1:N/2
' J; F- u% x+ T                                S=S+i(k)*sin(h*(k-1)*2*pi/N);
                              
                                C=C+i(k)*cos(h*(k-1)*2*pi/N);
. ^# [3 {. q$ G! g* C. T7 s7 F                           
         end
                                Is=S*4/N;
                             Ic=4*C/N;
                          I=abs(Is+sqrt(-1)*Ic)/sqrt(2);
# `) Q% ~2 V4 z6 z) n: ^                         JI=angle(Is+sqrt(-1)*Ic);
$ U" _6 `' i) i' xU                        
I
! O. W; i2 ?% W8 v& N& U& Ccos(JU-JI)

lixintmc

谢了，哈哈，