文档

fwind2

二维FIR滤波器采用二维窗口法

语法

h = fwind2(Hd, win)
h = fwind2(f1, f2, Hd, win)

描述

使用fwind2利用窗口法设计了二维FIR滤波器。fwind2使用二维窗口规范设计一个基于期望频率响应的二维FIR滤波器高清fwind2工作与二维窗口;使用fwind1使用一维窗口。

h = fwind2(Hd, win)生成二维FIR滤波器h使用所需频率响应的傅里叶反变换高清靠窗做乘法赢得高清是一个矩阵,它包含在笛卡尔平面上等间距点处所需的频率响应。fwind2返回h作为一个计算分子,这是合适的形式使用filter2h大小是一样的吗赢得

为了得到准确的结果,使用返回的频率点freqspace创建高清.(参见条目freqspace为更多的信息。)

h = fwind2(f1, f2, Hd, win)让您指定所需的频率响应高清在任意频率(f1f2)在x -y相互重合。频率向量f1f2应该在-1.0到1.0之间,其中1.0对应一半的采样频率,或π弧度。h大小是一样的吗赢得

类的支持金宝app

输入矩阵高清可以是一流的或任何整数类。所有其他输入fwind2必须是一流的.所有的输出都是class的

例子

全部折叠

这个例子展示了如何使用二维窗口法设计一个近似圆对称的二维带通滤波器。

创建频率范围向量f1f2使用freqspace.这些向量的长度是21。

(f1、f2) = freqspace(21日“meshgrid”);

计算每个位置到中心频率的距离。

r =√f1。^ 2 + f2。^ 2);

创建一个矩阵高清它包含所需的带通响应。在本例中,所需的通频带在0.1和0.5之间(归一化频率,其中1.0对应采样频率的一半,或弧度)。

高清= 1 (21);高清((r < 0.1) | (r > 0.5)) = 0;

显示理想的带通响应。

colormap (parula(64)网(f1, f2,高清)

创建一个2-D高斯窗口fspecial.规范化的窗口。

赢得= fspecial (“高斯”21岁,2);Win = Win ./ max(Win (:));

显示窗口。

网格(赢)

使用二维窗口,设计最佳产生所需频率响应的滤波器

h = fwind2(高清,赢得);

显示此滤波器的实际频率响应。

freqz2 (h)

算法

fwind2计算h用傅里叶反变换和乘以二维窗口赢得

h d n 1 n 2 1 2 π 2 π π π π H d ω 1 ω 2 e j ω 1 n 1 e j ω 2 n 2 d ω 1 d ω 2

h n 1 n 2 h d n 1 n 2 w n 1 n 2

参考文献

[1] Lim, Jae S,二维信号和图像处理, Englewood Cliffs, NJ, Prentice Hall, 1990,第202-213页。

之前介绍过的R2006a

这个话题有用吗?