减少频谱泄露的rife-vincent窗

kevin19821

这个窗有人知道吗?比汉宁窗还要好的窗,有人知道吗?

kevin19821

附我自己编的程序,估计那里需要技巧处理.
大家可参考<基于FFT算法的间谐波分析>  中国电机工程学报
程序如下:
clc;
clear all;
fs=2000;
dt=1./fs;
N=512;
% r=ones(1,N)
% Wd=r;
n=0:1:N-1;
t=n*dt;
a0=0.5;a1=0.46;a2=0.46;a3=0.5;
Wd=a0-a1*cos(2*pi*n/(N-1))+a2*cos(4*pi*n/(N-1))-a3*cos(6*pi*n/(N-1));
%  Wd=hann(N);
% Wd=hamming(N);
f=n*fs./N;
t=n*dt;
% s=2*cos(2*30*t);
% % s=2*sin(2*pi*50*t)+38*sin(2*pi*30.5*t);
s=3*sin(2*pi*50*t)+2*sin(2*pi*25*t)+3*sin(2*pi*75*t);
% wvtool(s)
sc=s.*Wd;
figure(1)
subplot(211)
plot(s)
subplot(212)
plot(sc)
s1=fft(s);
s2=fft(sc);
figure(2)
subplot(211)
plot(f-fs/2,abs([s1(N/2+1:N) s1(1:N/2)]));
subplot(212)
plot(f-fs/2,abs([s2(N/2+1:N) s2(1:N/2)])); %经rife-vecent处理
运行结果是不能减少频谱泄露的，不知道那里出错了

[ 本帖最后由 kevin19821 于 2007-9-17 19:54 编辑 ]