butler 发表于 2008-4-15 21:08

关于离散wigner分布的问题

sig=fmlin(128,0,0.25);               %sig为仿真信号            
M = length(sig);         %获取信号的长度   
fs=1;                           %fs为采样频率
N = 128;                        %取采样点数
N = min(N, M);
WVD=zeros(N,M);
for n = 1:M
    m = min();
    tau = -m : m;
    Temp = rem(N + tau,N) + 1;
    WVD(Temp,n) = sig(n + tau) .* conj(sig(n - tau));    %计算x*(t-1/2τ) x(t+1/2τ)
end
WVD = fft(WVD)/N;
f = linspace(-0.25,0.25,N)*fs;
t = (0: M-1)/fs;                                          
contour(t,f,abs(WVD));
axis();
ylabel('f');
xlabel('t / Sec');

上面是对一线性调频信号的离散wigner分布的程序,
仿真信号的频率范围为,为什么图像的频率会在[-0.25,0]呢?
在语句f = linspace(-0.25,0.25,N)*fs;中,若改为f = linspace(-0.5,0.5,N)*fs,为什么所得的等高线图是一样的呢?
程序中不对的地方请高手指正,谢谢
页: [1]
查看完整版本: 关于离散wigner分布的问题