freq_shear.m

来自「该工具箱为时频分析中的一类重要变换,离散时频分布的工具箱,其中提供基本算法的程序」· M 代码 · 共 30 行

function y = freq_shear(x, c)
% freq_shear -- chirp multiplication
%
%  Usage
%    y = freq_shear(x, c)
%
%  Inputs
%    x     signal vector, must have an odd length (to have a center point)
%    c     chirp rate
%
%  Outputs
%    y     the frequency sheared signal
%
% If x has length N and bandwidth B, then y will have length N and 
% bandwidth c*N*B/pi.  One must be careful to avoid aliasing.

% Copyright (C) -- see DiscreteTFDs/Copyright

error(nargchk(2, 2, nargin));

x = x(:);
N = length(x);
if (rem(N,2)==0)
  error('signal length must be odd')
end
M = (N-1)/2;
n = (-M:M)';

y = x .* exp(j*c/2*n.^2);