matlab filter函数源码,求MATLAB中filter2函数的源代码求MATLAB中二维 爱问知识人

若需要函数体，我可以传文件给你。下面是其代码function y = filter2(b,x,shape)

%FILTER2 Two-dimensional digital filter。

% Y = FILTER2(B,X) filters the data in X with the 2-D FIR

% filter in the matrix B。

The result, Y, is computed

% using 2-D correlation and is the same size as X。

%

% Y = FILTER2(B,X,'shape') returns Y computed via 2-D

% correlation with size specified by 'shape':

% 'same' - (default) returns the central part of the

% correlation that is the same size as X。

% 'valid' - returns only those parts of the correlation

% that are computed without the zero-padded

% edges, size(Y) size(X)。

%

% FILTER2 uses CONV2 to do most of the work。 2-D correlation

% is related to 2-D convolution by a 180 degree rotation of the

% filter matrix。

%

% Class support for inputs B,X:

% float: double, single

%

% See also FILTER, CONV2。

% Copyright 1984-2004 The MathWorks, Inc。

% $Revision: 5。13。4。2 $ $Date: 2004/03/09 16:16:19 $error(nargchk(2,3,nargin));

if nargin mx*nx)

% The filter is bigger than the input。

This is a nontypical

% case, and it may be counterproductive to check the

% separability of the stencil。

y = conv2(x,stencil,shape);

else

separable = false;

if all(isfinite(stencil(:)))

% Check rank (separability) of stencil

[u,s,v] = svd(stencil);

s = diag(s);

tol = length(stencil) * eps(max(s));

rank = sum(s > tol);

separable = (rank ==1);

end

if separable

% Separable stencil

hcol = u(:,1) * sqrt(s(1));

hrow = conj(v(:,1)) * sqrt(s(1));

if (all(all((round(stencil) == stencil))) & all(all((round(x) == x))))

% Output should be integer

y = round(conv2(hcol, hrow, x, shape));

else

y = conv2(hcol, hrow, x, shape);

end

else

% Nonseparable stencil

y = conv2(x,stencil,shape);

end

end

end。

全部