请问用fft如何实现离散wigner分布

butler

仿真信号为x＝fmlin(128,-0.5,0.5)，要用fft实现信号的离散wigner分布，离散形式采用Claasen，我求出了时延m的矩阵，老师说将m进行fft就行了，在时频域上画出的图是两条交叉的直线，但我对这个过程不是很了解。因为是新手，请高手不吝赐教，帮忙编个程序或者推荐一些相关的资料，谢谢！！

zhlong

wigner程序本版就有的

butler

能把链接给出来吗？我比较菜，找不出来，谢谢！！！

zhlong

http://forum.vibunion.com/forum/thread-40903-1-30.html

butler

sig=fmlin(128);
sig =real(sig);        
siglen = length(sig);           %  获取信号的长度   
fs=1;
N = 128;
N = min(N, siglen);
for n = 1:siglen,
    m = min([n-1,siglen-n,round(N/2)-1]);  
    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.5,0.5,N)*fs;
t = (0: siglen-1)/fs;                                          
contour(t,f,abs(WVD));
axis([min(t) max(t) min(f) max(f)]);
ylabel('f/ Hz');
xlabel('t / Sec');


请问 m = min([n-1,siglen-n,round(N/2)-1]);和Temp = rem(N + tau,N) + 1;这两句的作用是什么呢？谢谢！！

shenlong1hao

非常感谢！