半波傅氏算法的matlab实现

Gerneral888

%半波傅氏算法,数据窗10ms
N=12;
h=input('请输入谐波次数(1 or 3)：');
6 N0 p/ p5 h: x6 V% Wz=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);               
       subplot(211);  plot(t1,u1);title('连续图');
; z( t8 c2 k" f8 `/ w; d& i2 ^1 j$ o       subplot(212);  stem(t,u);title('离散图');
                                   i=u./z;  
                              
+ `. X8 |3 ~  @                S=0;C=0;
        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);
( A0 G; G, j3 _/ F6 u. D# c                        
        end
                                 Us=S*4/N;
: c  H) V3 a! m7 n                                Uc=4*C/N;
" ^- B5 l+ s( {( o% e% y0 `  e                               U=abs(Us+sqrt(-1)*Uc)/sqrt(2);
                             JU=angle(Us+sqrt(-1)*Uc);
               S=0;C=0;
3 z' j$ m5 m: \0 _& C, R         for  k=1:N/2
0 q7 _0 e+ g# J8 \7 H$ H4 `, Q, e                                S=S+i(k)*sin(h*(k-1)*2*pi/N);
                              
6 C/ t* o0 _4 s- g  R                                C=C+i(k)*cos(h*(k-1)*2*pi/N);
3 u6 n5 k2 `! a5 b5 F- i  f                           
         end
7 _) B: k+ j3 X% I                                Is=S*4/N;
                             Ic=4*C/N;
8 B6 H) E  {5 r" p                          I=abs(Is+sqrt(-1)*Ic)/sqrt(2);
$ k: f) n! X2 t' m! j1 k                         JI=angle(Is+sqrt(-1)*Ic);
7 L* [1 N0 H/ m6 e1 r# uU                        
I
cos(JU-JI)

lixintmc

谢了，哈哈，