找回密码
 注册
关于网站域名变更的通知
查看: 289|回复: 2
打印 上一主题 下一主题

[仿真讨论] 关于fft与加窗fft

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2022-9-20 16:18 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
用matlab仿真fft算法时,分别采用直接fft和加窗fft仿真,为什么加窗后的结果还不如不加窗的呢?# P) i) F$ L* N9 p" \5 h
程序如下:( W+ o- f) y( [$ T
直接fft程序:
5 N; J- Z+ }. n  zclear;! w1 O- |3 t* ]! z8 ]
fs=12800;
; e9 x8 m* @9 \8 w) @; V  ], K6 }- XN=256;" x. {3 n6 r) D$ e
t = 0:1/fs:0.02-1/fs;: D7 m/ X5 M4 w* v8 _
y=cos(2*pi*50*t+10*pi/180)+cos(2*pi*300*t+20*pi/180); %模拟输入信号
. M/ y! \+ E/ N- c5 w9 c( DY = fft(y,256);6 L# x: ^- m. \
Pyy =abs(Y(1:256))*2/256; %幅值9 H5 Q: Z0 U' N% v1 G7 b0 r8 Y* q* C  f
pha=angle(Y(1:256))/pi*180; %相位' s/ F- i. K/ }. G0 p7 G/ `
' V8 W4 x. b! X; u3 K: Q: h
加窗fft程序:
2 V5 R& y: M7 n" [clear;
8 G3 u3 s6 f8 [' Y/ nfs=12800;9 z$ Y: @3 b" p+ L- K
N=256;
3 M" Y6 v( K0 a  M9 b" B7 N- kt=0:1/fs:0.02-1/fs;
; o; J% n, ?' i' y: Wx=cos(2*pi*50*t+10*pi/180)+cos(2*pi*300*t+20*pi/180); %模拟输入信号
' J) h- n+ `, p# |& L9 _; E- I& L2 In=[0:N-1];
5 o0 {. X0 ~- ]7 K+ Rw=0.5-0.5*cos(2*pi*n/N); % 汉宁5 @, E; u2 c7 T
y=x.*w;
* F( B- T4 d6 x8 jY = fft(y,N);$ s0 D) s) S9 `) M  B* [
Pyy =abs(Y(1:N))*2*2/N; %幅值
9 @2 b/ D& s5 _2 Opha=angle(Y(1:N))/pi*180; %相位
5 u& N7 D( V/ g, {1 Z1 W
8 c' y( k% G1 L5 f加窗的频谱泄漏为何更严重?, Y# r0 A) U3 B: L% w( e3 o* ]

该用户从未签到

2#
发表于 2022-9-20 16:51 | 只看该作者
我没有用你的程序仿真,单一般这样的问题造成的原因是,窗函数不正确,窗的类型没选对,比如衰减度并不符合你原来信号频谱的要求。* F2 B, D5 i8 P3 T1 i5 l1 t

该用户从未签到

3#
发表于 2022-9-20 17:09 | 只看该作者
你可以试试用MATLAB自配的窗函数试试,用凯撒窗通过调BETA值可以调整衰减度,多试几次对比一下,慢慢就能看出结果了。$ }  A% e3 \. L4 l9 F8 b6 D: n- K6 s3 i
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-6-28 13:19 , Processed in 0.093750 second(s), 23 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

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