声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2167|回复: 4

[FFT] czt变换与fft变换,幅值不等问题求教

[复制链接]
发表于 2017-11-6 19:56 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 信号小虾米 于 2017-11-6 20:07 编辑

    我对一个随机信号进行fft变换以及czt变换后发现同频(20HZ——50HZ)的幅值是不一样的,而且差距还挺大,想问问各位大神,这是什么原因造成的。下面是我写的程序以及得到的频谱图,希望各位大神帮我看看,感激万分!!!程序:
data=importdata('68.txt')
figure(1)
plot(data)
xlabel('采样点数');ylabel('接触压力/N')


fs=1000
N=1024
n=0:N-1
y1=fft(data,N)
mag1=abs(y1)
mag2=mag1*2/N
f=n*fs/N
figure(2)
subplot(1,2,1)
plot(f(1:N/2),mag2(1:N/2))
axis([20,50,0,20])
xlabel('频率/Hz');ylabel('振幅')
title('fft变换频谱')


f1=20;f2=50;M=1024;
w=exp(-j*2*pi*(f2-f1)/(fs*M))
a=exp(j*2*pi*f1/fs)
y=czt(data,M,w,a)
mag3=abs(y)*2/N
h=0:1:M-1
f0=(f2-f1)/M*h+20
subplot(1,2,2)
plot(f0,mag3)
xlabel('频率/Hz');ylabel('振幅')
title('CZT变换后的细化频谱')

频谱图

频谱图
回复
分享到:

使用道具 举报

发表于 2017-11-7 14:25 | 显示全部楼层
幅值差这么多?
 楼主| 发表于 2017-11-8 14:43 | 显示全部楼层

对啊,查了差不多一倍,不知道是什么原因呢?
发表于 2017-11-9 12:58 | 显示全部楼层

回帖奖励 +1 点体能

是不是得加窗啊?
发表于 2017-11-17 17:11 | 显示全部楼层
关注,学习
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-4-25 18:23 , Processed in 0.102548 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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