半波傅氏算法的matlab实现

Gerneral888

%半波傅氏算法,数据窗10ms
N=12;
h=input('请输入谐波次数(1 or 3)：');
z=1+sqrt(-1)*sqrt(3);            
  B; U- S+ I; R' N9 w$ z/ L' vt=0:0.02/N:0.02;
u= sin(2*pi*50*t)+2*sin(2*3*pi*50*t-pi/6);  
9 h" _! ?1 W) B, L' c6 @  t1=0:0.0001:0.07;   
" O  H$ y( g# Y- t0 G   u1=sin(2*pi*50*t1)+2*sin(2*3*pi*50*t1-pi/6);               
" _& j) @5 V. ~+ A0 @       subplot(211);  plot(t1,u1);title('连续图');
) F: e4 Z# x7 S       subplot(212);  stem(t,u);title('离散图');
                                   i=u./z;  
                              
5 x7 Z& [7 ?& o2 R) B                S=0;C=0;
. R7 I0 K+ U% |: T) y  l        for  k=1:N/2            
' Z# [1 t$ M0 R# }8 T$ K( a# O                             S=S+u(k)*sin(h*(k-1)*2*pi/N);   
                             C=C+u(k)*cos(h*(k-1)*2*pi/N);
. T; L% V) S* c3 N# c% Q                        
        end
                                 Us=S*4/N;
+ c2 b. ]! L0 u0 n) g; t" P1 p                                Uc=4*C/N;
                               U=abs(Us+sqrt(-1)*Uc)/sqrt(2);
% {8 C' k+ b  B! z/ D                             JU=angle(Us+sqrt(-1)*Uc);
7 i' x) W/ c; ~/ o               S=0;C=0;
         for  k=1:N/2
                                S=S+i(k)*sin(h*(k-1)*2*pi/N);
, J2 a5 i- x7 b                              
                                C=C+i(k)*cos(h*(k-1)*2*pi/N);
                           
: m$ V. f4 v+ ~; p0 _# r) p6 b7 W         end
7 u8 f, y3 w) I8 z. S                                Is=S*4/N;
. W( @) F: u" }( V$ t: S                             Ic=4*C/N;
, J+ K9 q4 l1 @/ s/ U5 s                          I=abs(Is+sqrt(-1)*Ic)/sqrt(2);
                         JI=angle(Is+sqrt(-1)*Ic);
( X* o3 y* B5 \9 x# nU                        
I
cos(JU-JI)

lixintmc

谢了，哈哈，