作业帮 > 数学 > 作业

关于离散傅里叶变换后频谱的问题

来源:学生作业帮 编辑:搜搜做题作业网作业帮 分类:数学作业 时间:2024/06/23 20:05:50
关于离散傅里叶变换后频谱的问题
我想分析一组变形数据(共2036个点)的频率,我对一组变形数据(2036个)做了快速傅里叶变换,得出了2048组复数值,如何由这组复数值求频谱图.今儿确定频率,主要是频谱图中的左边周如何确定!万谢
关于离散傅里叶变换后频谱的问题
建议你做fft时做一个频谱矫正搬移,即做fft时,要这样做,fftshift(fft(x,2048)).这样之后,假设时域采样数据的频率为F,那么你得到的2048组复数值所对应的频率值依次为F*n/2048,其中n=1,2.2047.对DFT所得复数值取模值,画出随频率变化的曲线,即为幅频特性(频谱图);取相位角,画出随频率变化的曲线,即为相频特性.这只是具体方法,一般的原理你可以参考数字信号处理的书,或者继续交流也可,
再问: 想问一下,由于我的数据时每天一个数,那么时域采样数据的频率F 应该怎么假设?
假如F=2036的话,画出的频谱图0和2036近处巨大的突出,而中间的值相对很小会是怎么回事?请详解下,谢谢
再答: 关于时域采样频率,按照你的采样情况,采样周期是1天,一般应以秒计,不会有错,但你可以试试用天作为单位时间,这个我没这样做过;但是你说F=2036,我不理解,以秒计的话,频率值为F=1/(24*3600)。另一个问题,你说的两个突出处位于采样频率区间的两端,是0频位置,说明你的信号有较大成分的直流分量。
再问: 不好意思,这么晚才回复您,这是我的代码,我把直流分量去掉还是不行,哪不对么?FFT始终看不出怎么看出aa的周期?不知怎么回事。 另外感觉频率 f不对?帮忙改改,非常感谢
aa=aa-mean(aa);
NFFT=2048;L=2036;
fs=1;
f=fs*(0:NFFT-1)/NFFT;
y=fft(aa,NFFT);
plot(f,abs(y))

这里aa是观测量,采样频率按 1 来算
再答: 你的程序我看不出什么问题。FFT可以看出信号的谱分布,没有周期性;频率怎么不对?一般得出的频谱图并不是书本上的一条一条的线,因为有噪声!
再问: 图画出来了,按那个最高振幅对应的频率算出来的周期正好是我的离散数据段的时间,2000多天,(我的频率是次/天)这是怎么回事,应该去除吧?还有这个频率如何选取啊,太宽了
再答: 你的问题问的比较笼统。。你要明确你到底是要分析什么。你上面提到“最高振幅对应的频率算出来的周期正好是我的离散数据段的时间”,可是为什么要去除呢?这个或为巧合或有理论依据,只要不妨碍你的分析目标,就不用去除呀!你最后的问题,频率如何选取,是要选择什么频率?(很抱歉,前几天在路上,没有及时回复)