半波傅氏算法的matlab实现

Gerneral888

%半波傅氏算法,数据窗10ms
3 W2 V" p( s+ F' f  k7 ?N=12;
h=input('请输入谐波次数(1 or 3)：');
# \5 u8 @1 L. L2 a; e: y3 m9 Sz=1+sqrt(-1)*sqrt(3);            
8 w) t  n4 [% d/ X7 I0 bt=0:0.02/N:0.02;
6 [/ d1 F' E+ ?7 }! f8 u 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);               
$ @$ z' m# a( Z* T6 Z% R       subplot(211);  plot(t1,u1);title('连续图');
       subplot(212);  stem(t,u);title('离散图');
  m/ [# x8 T: K) f/ B: S5 r                                   i=u./z;  
                              
                S=0;C=0;
# @4 H8 }: h/ W' Z        for  k=1:N/2            
                             S=S+u(k)*sin(h*(k-1)*2*pi/N);   
                             C=C+u(k)*cos(h*(k-1)*2*pi/N);
                        
& P  E0 _5 g4 L, a        end
7 P/ Y0 B) J. g" F                                 Us=S*4/N;
/ U  B, k( ^. q+ h, F                                Uc=4*C/N;
( s2 M( n) P  y5 h: @5 X                               U=abs(Us+sqrt(-1)*Uc)/sqrt(2);
                             JU=angle(Us+sqrt(-1)*Uc);
/ r: s) I7 y" {8 {) C               S=0;C=0;
         for  k=1:N/2
% X& z, Z9 i& t' }& [$ d                                S=S+i(k)*sin(h*(k-1)*2*pi/N);
4 ?/ y0 a* J  k3 p: W* L8 Z                              
                                C=C+i(k)*cos(h*(k-1)*2*pi/N);
                           
         end
* A9 }2 J6 I& p! V                                Is=S*4/N;
                             Ic=4*C/N;
                          I=abs(Is+sqrt(-1)*Ic)/sqrt(2);
                         JI=angle(Is+sqrt(-1)*Ic);
U                        
I
cos(JU-JI)

lixintmc

谢了，哈哈，