半波傅氏算法的matlab实现

Gerneral888

%半波傅氏算法,数据窗10ms
N=12;
5 t, L3 W0 a' P/ m h=input('请输入谐波次数(1 or 3)：');
/ a; |, p! L3 ^* f5 S/ Cz=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;   
   u1=sin(2*pi*50*t1)+2*sin(2*3*pi*50*t1-pi/6);               
  S. z+ h! k) Y       subplot(211);  plot(t1,u1);title('连续图');
       subplot(212);  stem(t,u);title('离散图');
                                   i=u./z;  
+ o0 ?3 e5 U( A4 a: U6 i                              
* w  s( H' H& ?& r                S=0;C=0;
        for  k=1:N/2            
6 X, p: Z" @1 O* W                             S=S+u(k)*sin(h*(k-1)*2*pi/N);   
( B7 [; _2 P) M' v9 M                             C=C+u(k)*cos(h*(k-1)*2*pi/N);
                        
        end
                                 Us=S*4/N;
                                Uc=4*C/N;
                               U=abs(Us+sqrt(-1)*Uc)/sqrt(2);
" z, Q6 [% W" O                             JU=angle(Us+sqrt(-1)*Uc);
& o6 E* T! |8 X               S=0;C=0;
9 O1 W5 G) `: v1 |+ ^% P5 J/ N8 b$ k         for  k=1:N/2
                                S=S+i(k)*sin(h*(k-1)*2*pi/N);
7 Y$ k3 W% d- u4 {$ K. ?                              
# ?  N! U5 V( C+ G4 t1 I- T+ d" d                                C=C+i(k)*cos(h*(k-1)*2*pi/N);
                           
+ m' }$ [5 L- G: N1 l         end
                                Is=S*4/N;
                             Ic=4*C/N;
                          I=abs(Is+sqrt(-1)*Ic)/sqrt(2);
  U1 R# ^3 W; h- R3 S* m: M+ Q                         JI=angle(Is+sqrt(-1)*Ic);
1 K5 i6 N; ?# X9 a% M$ p3 }& T7 x8 HU                        
I
cos(JU-JI)

lixintmc

谢了，哈哈，