半波傅氏算法的matlab实现

Gerneral888

%半波傅氏算法,数据窗10ms
& O! D1 t% L5 k% @% j+ m- VN=12;
h=input('请输入谐波次数(1 or 3)：');
z=1+sqrt(-1)*sqrt(3);            
- o* i; C+ Q* A  m0 Et=0:0.02/N:0.02;
6 t. i+ Y! D! P2 C% k4 ^ u= sin(2*pi*50*t)+2*sin(2*3*pi*50*t-pi/6);  
, W8 \, x' C" x: {3 c0 K. e  t1=0:0.0001:0.07;   
4 b( U' l2 L3 F; Q/ a' I   u1=sin(2*pi*50*t1)+2*sin(2*3*pi*50*t1-pi/6);               
       subplot(211);  plot(t1,u1);title('连续图');
9 i: @  Q( J- h1 Q5 I; g: K       subplot(212);  stem(t,u);title('离散图');
                                   i=u./z;  
+ t, c+ M9 g8 r% d9 M  t- Q! Y5 X                              
9 K" [4 d+ N& k% _# |/ b                S=0;C=0;
$ {$ ~4 b* }' f+ k6 @        for  k=1:N/2            
5 D+ u5 r  y% I- N# `/ n: \                             S=S+u(k)*sin(h*(k-1)*2*pi/N);   
" V9 M9 I  _) f$ b7 W; q                             C=C+u(k)*cos(h*(k-1)*2*pi/N);
                        
& w4 f: p' h1 W8 t. B        end
. h- H% k7 q8 X2 H+ t4 W                                 Us=S*4/N;
4 `- s' B' A3 S- i# e                                Uc=4*C/N;
# i$ |- C1 @8 _- S                               U=abs(Us+sqrt(-1)*Uc)/sqrt(2);
# b2 W" E! C. y1 ~3 L  O                             JU=angle(Us+sqrt(-1)*Uc);
               S=0;C=0;
         for  k=1:N/2
. z$ _: a+ P/ P+ |8 J2 O/ P* x3 E                                S=S+i(k)*sin(h*(k-1)*2*pi/N);
# s. g, E4 r+ H! k- ^7 _* G                              
4 @. K% }! x' x* V& U# R                                C=C+i(k)*cos(h*(k-1)*2*pi/N);
8 F" ]" ~& F2 L3 s                           
         end
                                Is=S*4/N;
                             Ic=4*C/N;
" F% r7 \& R4 t; ^( I                          I=abs(Is+sqrt(-1)*Ic)/sqrt(2);
% ~4 P0 n+ y$ U0 x. l; l                         JI=angle(Is+sqrt(-1)*Ic);
U                        
! k9 B( f2 j" q: \I
2 K2 F+ N) g3 ]# _) e% `1 ^9 kcos(JU-JI)

lixintmc

谢了，哈哈，