clear ; clc ; close ; N = 7 ; U = 4 ; a = window ( ' re ' ,N ) ; for n = - 3 + U : 1 : 3 + U if ( n == 4 ) h ( n ) = 0.4 ; else h ( n ) = ( sin ( 2 * %pi * ( n - U ) / 5 ) / %pi * ( n - U ) ) ; end h1 ( n ) = h ( n ) * a ( n ) ; end [ av,ba ] = frmag ( h1, 256 ) ; pu = 20 * log10 ( av ) / max ( av ) ; plot ( 2 * ba,pu ) ; xgrid ( 2 ) ;
Comments
Post a Comment