luboq.m

这是一个关于hht变换很有用的工具箱

function  q=luboq(nx,ny,m)
 
% The function LUBOQ returns a 2-D filter of nx by ny.
% 
% Both nx and ny should be even. 
% m should be between 0 and 1; 
% The smaller m is the wider filter is generated.
%
% Calling sequence-
% q=luboq(nx,ny,m)
% 
% Input-
%	nx	- number of filter points in 1st dimension
%	ny	- number of filter points in 2nd dimension
%	m	- smoothing window
% Output-
%	q	- filter
 
% Z.Shen (JHU)			1997 Initial

%----- Construct a 2-D filter
q1=fir1(nx,m);
q2=fir1(ny,m);
[Q1,Q2]=meshgrid(q1,q2);
q=Q1.*Q2;
q=q/(max(max(q)));