减少频谱泄露的rife-vincent窗
这个窗有人知道吗?比汉宁窗还要好的窗,有人知道吗? 附我自己编的程序,估计那里需要技巧处理.大家可参考<基于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());
subplot(212)
plot(f-fs/2,abs()); %经rife-vecent处理
运行结果是不能减少频谱泄露的，不知道那里出错了
[ 本帖最后由 kevin19821 于 2007-9-17 19:54 编辑 ]
页:
[1]