半波傅氏算法的matlab实现

Gerneral888

%半波傅氏算法,数据窗10ms
N=12;
% P8 z# p# N2 o$ F5 ^$ G h=input('请输入谐波次数(1 or 3)：');
z=1+sqrt(-1)*sqrt(3);            
  M8 e' j3 f8 [' a0 \/ {5 _: C# ~! m% A/ n) \t=0:0.02/N:0.02;
! F! i$ `$ Z: x3 ~ 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);               
       subplot(211);  plot(t1,u1);title('连续图');
3 D5 j$ y: S& W, T, `, z1 v: G- w* ~       subplot(212);  stem(t,u);title('离散图');
                                   i=u./z;  
                              
$ R3 ^6 \0 w- C( U: l, X, o7 [                S=0;C=0;
& l0 n+ u/ P. V  ^; m+ W( L        for  k=1:N/2            
+ r% q4 S0 ?1 [" {; U3 Q                             S=S+u(k)*sin(h*(k-1)*2*pi/N);   
                             C=C+u(k)*cos(h*(k-1)*2*pi/N);
                        
7 ~3 c  T/ I# @3 Q1 a1 T        end
                                 Us=S*4/N;
# B( e. U. W! ]3 `6 `- u8 K3 P                                Uc=4*C/N;
! V" n/ [6 X9 p' g* l: [% c7 C$ [                               U=abs(Us+sqrt(-1)*Uc)/sqrt(2);
$ m8 s  j" T, l- @' d) T0 x, [: a5 j                             JU=angle(Us+sqrt(-1)*Uc);
               S=0;C=0;
) C. [1 G1 x/ O# J# Y7 o         for  k=1:N/2
- X) r  }7 |' u: K3 n% h! A                                S=S+i(k)*sin(h*(k-1)*2*pi/N);
                              
                                C=C+i(k)*cos(h*(k-1)*2*pi/N);
$ j, R' K) [8 h5 C7 ]                           
         end
0 W6 f' V% r8 E, k+ p9 }                                Is=S*4/N;
                             Ic=4*C/N;
                          I=abs(Is+sqrt(-1)*Ic)/sqrt(2);
                         JI=angle(Is+sqrt(-1)*Ic);
U                        
I
- X9 p0 n6 o- n9 s5 Ocos(JU-JI)

lixintmc

谢了，哈哈，