陷波滤波器matlab程序

陷波器是一种特殊的带阻滤波器,其阻带在理想情况下只有一个频率点,因此也被称为点阻滤波器。这种滤波器主要用于消除某个特定频率的干扰,例如,在各种测量仪器和数据采集系统中用于消除电源干扰的工频陷波器。

function g = dftfilt（f H classout）
%DFTFILT Performs frequency domain filtering.
%   g = DFTFILT（f H CLASSOUT） filters f in the frequency domain
%   using the filter transfer function H. The output g is the
%   filtered image which has the same size as f. 
%
%   Valid values of CLASSOUT are
%
%	‘original‘	The ouput is of the same class as the input. 
%               This is the default if CLASSOUT is not included
%               in the call.
%	‘fltpoint‘	The output is floating point of class single unless
%               both f and H are of class double in which case the
%               output also is of class double.
%
%   DFTFILT automatically pads f to be the same size as H. Both f
%   and H must be real. In addition H must be an uncentered
%   circularly-symmetric filter function.  

%   Copyright 2002-2009 R. C. Gonzalez R. E. Woods and S. L. Eddins
%   From the book Digital Image Processing Using MATLAB 2nd ed.
%   Gatesmark Publishing 2009.
%
%   Book web site: http://www.imageprocessingplace.com
%   Publisher web site: http://www.gatesmark.com/DIPUM2e.htm

% Convert the input to floating point.
[f revertClass] = tofloat（f）;

% Obtain the FFT of the padded input.
F = fft2（f size（H 1） size（H 2））;

% Perform filtering. 
g = ifft2（H.*F）;

% Crop to original size.
g = g（1:size（f 1） 1:size（f 2））; % g is of class single here.

% Convert the output to the same class as the input if so specified.
if nargin == 2 || strcmp（classout ‘original‘）
   g = revertClass（g）;
elseif strcmp（classout ‘fltpoint‘） % g is floating point already.
   return
else
   error（‘Undefined class for the output image.‘）
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1697  2013-10-22 00:41  源码\dftfilt.m

     文件       2191  2013-10-22 00:41  源码\gscale.m

     文件        560  2014-10-27 14:51  源码\main.m

     文件       3938  2014-10-27 11:09  源码\recnotch.m

     文件       1132  2013-10-22 00:41  源码\tofloat.m

     目录          0  2014-11-06 15:51  源码

----------- ---------  ---------- -----  ----

                 9518                    6