wvd.m

时频分布算法中的短时fourier变换和维格纳变换的程序。

function [Wz, Z] = WVD(z, delta_n)
N = length(z);
L = (N-1)/2;
m = -L : L;
k = [-L/2 : L/2]';
In_M = k*m;
W_M = exp(-j*4*pi/N*In_M);
Z = zeros(N, (N-1)/delta_n);
for n = -L : delta_n : L
    m = max(-L-n, -L+n) : min(L-n, L+n);
    Z(m+L+1, n+L+1) = z(n+m+L+1) .* conj(z(n-m+L+1));
end
Wz = W_M * Z;