Documentation

This is machine translation

Translated by Microsoft
Mouseover text to see original. Click the button below to return to the English verison of the page.

Note: This page has been translated by MathWorks. Please click here
To view all translated materals including this page, select Japan from the country navigator on the bottom of this page.

ftrans2

2-D FIR filter using frequency transformation

Syntax

h = ftrans2(b,t)
h = ftrans2(b)

Description

h = ftrans2(b,t) produces the two-dimensional FIR filter h that corresponds to the one-dimensional FIR filter b using the transform t. b must be a one-dimensional, Type I (even symmetric, odd-length) filter such as can be returned by fir1, fir2, or firpm in the Signal Processing Toolbox software. The transform matrix t contains coefficients that define the frequency transformation to use.

example

h = ftrans2(b) uses the McClellan transform matrix t.

t = [1 2 1; 2 -4 2; 1 2 1]/8;

Examples

collapse all

Use ftrans2 to design an approximately circularly symmetric two-dimensional bandpass filter with passband between 0.1 and 0.6 (normalized frequency, where 1.0 corresponds to half the sampling frequency, or π radians). Since ftrans2 transforms a one-dimensional FIR filter to create a two-dimensional filter, first design a one-dimensional FIR bandpass filter using the Signal Processing Toolbox function firpm.

colormap(jet(64))
b = firpm(10,[0 0.05 0.15 0.55 0.65 1],[0 0 1 1 0 0]);
[H,w] = freqz(b,1,128,'whole');
plot(w/pi-1,fftshift(abs(H)))

Use ftrans2 with the default McClellan transformation to create the desired approximately circularly symmetric filter.

h = ftrans2(b);
freqz2(h)

Input Arguments

collapse all

1-D FIR filter, specified as a numeric matrix. b must be a 1-D Type I (even symmetric, odd-length) filter such as can be returned by fir1, fir2, or firpm in the Signal Processing Toolbox software,

Example:

Data Types: double

The transform matrix, specified as a numeric matrix. t contains coefficients that define the frequency transformation to use.

Example:

Data Types: double

Output Arguments

collapse all

2-D FIR filter, returned as a numeric matrix. ftrans2 returns h as a computational molecule, which is the appropriate form to use with filter2. If t is m-by-n and b has length Q, then h is size ((m-1)*(Q-1)/2+1)-by-((n-1)*(Q-1)/2+1).

Algorithms

The transformation below defines the frequency response of the two-dimensional filter returned by ftrans2.

H(ω1,ω2)=B(ω)|cosω=T(ω1,ω2),

where B(ω) is the Fourier transform of the one-dimensional filter b:

B(ω)=n=NNb(n)ejωn

and T(ω1,ω2) is the Fourier transform of the transformation matrix t:

T(ω1,ω2)=n2n1t(n1,n2)ejω1n1ejω2n2.

The returned filter h is the inverse Fourier transform of H(ω1,ω2):

h(n1,n2)=1(2π)2ππππH(ω1,ω2)ejω1n1ejω2n2dω1dω2.

References

[1] Lim, Jae S., Two-Dimensional Signal and Image Processing, Englewood Cliffs, NJ, Prentice Hall, 1990, pp. 218-237.

See Also

| | | |

Introduced before R2006a

Was this topic helpful?