声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 6940|回复: 4

[FFT] 汉宁窗函数的fft变换问题(附程序)??

[复制链接]
发表于 2007-8-30 21:49 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
%%%程序代码如下:
fs=1000;
t=0:1/fs:0.6;
N=length(t);
f=fs*(0:N-1)./N;
n=0:N-1;
w=0.5+0.5*cos(pi*n/N); %Hanning窗
figure(1);
plot(t,w);
title('窗函数');
xlabel('time');
Y=fft(w,N);
PY = Y.* conj(Y)/N;   %%求幅值
figure(2);
plot(f,PY);         
title('窗函数FFT变换');
xlabel('Frequence');
%%%程序结束

问题:上述程序运行后,得到%Hanning窗的频谱特性,但是发现在0赫兹附近有峰外,在结尾处还有峰,不知道这个分析对不对?结尾的峰值是什么原因?还是应该仅对一半时间显示FFT频谱结果?我猜想是fft系数的对称性,但是发现他们的幅值不一样高啊??请教大家啊!

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

 楼主| 发表于 2007-8-31 07:58 | 显示全部楼层
大家有什么意见啊?
发表于 2007-8-31 09:09 | 显示全部楼层
首先你程序里的窗只是一半的hanning窗
发表于 2007-8-31 09:12 | 显示全部楼层
本帖最后由 wdhd 于 2016-9-7 15:34 编辑
原帖由 caichengtao 于 2007-8-30 21:49 发表
%%%程序代码如下:
fs=1000;
t=0:1/fs:0.6;
N=length(t);
f=fs*(0:N-1)./N;
n=0:N-1;
w=0.5+0.5*cos(pi*n/N); %Hanning窗
figure(1);
plot(t,w);
title('窗函数');
xlabel('time');
Y=fft(w,N);
PY ...
问题:上述程序运行后,得到%Hanning窗的频谱特性,但是发现在0赫兹附近有峰外,在结尾处还有峰,不知道这个分析对不对?结尾的峰值是什么原因?还是应该仅对一半时间显示FFT频谱结果?我猜想是fft系数的对称性,但是发现他们的幅值不一样高啊??请教大家啊!

1,楼主给的程序只是一半海宁窗,应改为
w=0.5-0.5*cos(2*pi*n/N);
2,结尾峰值的原因是因对称,应该仅对一半频率(不是时间)显示FFT频谱结果。
3,最后一条谱 线PY(N)对称于PY(2),这两条一样高。

评分

1

查看全部评分

 楼主| 发表于 2007-8-31 14:06 | 显示全部楼层
感谢songzy41 !!!通过你的回答,对信号处理有新的认识,知道了为什么采样定理要2倍了啊!振动论坛!!!每天的好朋友!!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-27 00:33 , Processed in 0.072399 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表