在matlab中,通过FFT得到幅值-频率谱图后,怎样求得基频,

imagesc(Q)

n要取最接近数据长度的2的整数次方,命令是2^nextpow2(N),其中N是实际数据长度,因为这样的n可以使fft更快.比如,N=1021,执行n=2^nextpow2(N)后,n=2^10=102

一个长度为n的信号A其fft之后得到的信号长度也是n得到的fft数据除了0频率分量幅度为时域值的n陪外其他亲率的幅度值是n/2倍所以要获得某个非零频率的幅度,需要将fft之后的数据除以n乘以2b_S=

fft是一维傅里叶变换,即将时域信号转换位频域fftshift是针对频域的,将FFT的DC分量移到频谱中心即对频域的图像,（假设用一条水平线和一条垂直线将频谱图分成四块）对这四块进行对角线的交换与反对

1、NFFT=2^nextpow2(N)*2这里是否乘以2,应该怎么选择,貌似乘以3结果也一样,此处只是确定做FFT的点数,由于是单频信号,所以,在最后的结果上面会出现两个对称的尖峰,和你的FFT点数

%下面的程序里Pn存的就是基波相位如果求的是谐波相位,稍微修改即可x=load('data.dat');%load数据fs=10000;%采样频率,自己根据实际情况设置N=length(x);%x是待

fft值大小反应信号在这个频率点的幅值大小或功率大小再问：得到fft值后怎么算出它的振幅？再答：如果对傅里叶变换的公式不熟悉,用matlab看看吧symstAf;x1=A*sin(2*pi*f*t);

做FFT的函数是fft,注意结果要取模（abs）以便于显示,再fftshift一下将0频移到中间.取最大值的函数是max.

fftshift的作用正是让正半轴部分和负半轴部分的图像分别关于各自的中心对称.因为直接用fft得出的数据与频率不是对应的,fftshift可以纠正过来以下是Matlab的帮助文件中对fftshift

fft(x)就是fft(x,x.length);N默认为x的长度.还有就是*Ts就是把时域变成数字域.

大于x的长度时就在x原数据后面补零,而补零不影响信号的分辨率,只是将频谱平滑处理,因为时域补零相当于频域插值；大于x的长度时就直接截取x数据中相应长度的元素.再问：小于x的长度的长度是怎样处理的呢？比

fft是内建函数,不是matlab写的,看不到源代码的下面是我写的一个fft,可以用functionxn=myfft(x)N=length(x);M=log2(N);xtmp=zeros(1,N);v

一维快速傅里叶变换,低频成分在矩阵的两边,要将低频成分放到中间的话用fftshift(fft())

matlanb会自动补零的

x=[114.37594.75126105.75108.625146.75123.5101.25126.25103.875111.875147.5112.375109.375150.5108.8751

FFT的频谱本身就是对称的,一般画频谱图时取前半部分就行.再问：怎么取？谢谢再答：假如你想画功率谱，设频率w，频谱F。则plot(w(1:Lt/2),F(1:Lt/2);%%Lt是转换成波谱后的信号长

有关matlab中fft（傅里叶变换）的问题请高手解惑爱是包容而不是放纵爱是关怀而不是宠爱爱是相互交融而不是单相思爱是百味而不全是甜蜜

表示对数据向量x进行16点的快速傅里叶变换,不管x长度是大于还是小于16

最后加一句y=real(y);real为求复数实部的函数这样x和y就一样了,因为反变换后虚部系数都为零

最大值所在的Bin：max_index=max(ans)该Bin代表的是基频的频率：Ffund=max_index*FsamplingRate/NumberOfSamples