半波傅氏算法的matlab实现

Gerneral888

%半波傅氏算法,数据窗10ms
N=12;
h=input('请输入谐波次数(1 or 3)：');
$ d% i" l& L1 k8 {2 p: _& z% w& Q1 Yz=1+sqrt(-1)*sqrt(3);            
, o4 u; ^8 s: r) vt=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;   
2 ^  |& H( T; L8 w   u1=sin(2*pi*50*t1)+2*sin(2*3*pi*50*t1-pi/6);               
5 c- l  A6 d2 Z0 `; i7 u2 }       subplot(211);  plot(t1,u1);title('连续图');
3 z' K  k& Z6 c2 `       subplot(212);  stem(t,u);title('离散图');
                                   i=u./z;  
                              
                S=0;C=0;
- U" J3 k& W  {0 I! W- `0 f        for  k=1:N/2            
6 \  V, \# f  K* D  y; {. B5 j                             S=S+u(k)*sin(h*(k-1)*2*pi/N);   
% ?+ v' f' q. [! Z( R/ e                             C=C+u(k)*cos(h*(k-1)*2*pi/N);
* B0 k8 m" i$ V; Z/ E                        
1 N  e. P# ~3 h! K) Y) y        end
                                 Us=S*4/N;
                                Uc=4*C/N;
                               U=abs(Us+sqrt(-1)*Uc)/sqrt(2);
7 `; ?0 Z/ k6 j. H/ l( K                             JU=angle(Us+sqrt(-1)*Uc);
' K% Y; a, Q) A# e8 @               S=0;C=0;
) _, ~; ~$ p$ N         for  k=1:N/2
                                S=S+i(k)*sin(h*(k-1)*2*pi/N);
! r8 K% P0 M( c) \                              
4 R# O8 c' X4 [$ h7 T& ?$ x( b                                C=C+i(k)*cos(h*(k-1)*2*pi/N);
                           
         end
& i( d7 ^7 d% |& e- W; Q& m! o                                Is=S*4/N;
6 g, ^! }" k* x- F' b                             Ic=4*C/N;
9 e. F1 ?2 k1 t0 L% D8 z2 P: O# f                          I=abs(Is+sqrt(-1)*Ic)/sqrt(2);
                         JI=angle(Is+sqrt(-1)*Ic);
U                        
: w% ^, u3 L$ v" D2 T, GI
cos(JU-JI)

lixintmc

谢了，哈哈，