matlab中Y=fft(y,512);

imagesc(Q)

y(:,k)表示取数组y中第k列的所有元素；类似的,y(:,1:8)指取y中1到8列的所有元素.两边没有数的冒号表示所有行的意思

length(y)是总的采样点数fs是采样频率,就是频域横轴的最大值.那么频域上的刻度,就是df=fs/length(y)从0到fs取n个点,那么就有n-1个区间,就是数列的项的问题,你推推看吧.f(

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

你的那个程序对初学者来说比较难理解,我帮你改了一下,概念比较清晰,你可以对照一下Fs=128;%%采样频率N=512;%%采样点数dt=1/Fs;%%时域最小间隔,即时域分辨率t=(0:N-1)*dt

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

x为从-8到8以0.5为步长的列向量,即x=（-8,-7.5……）’y是x的转置,即y是行向量ones(size(y))生成和y的元素个数相同且元素全为1的向量,当然格式也和y相同即也为行向量X=x中

因为sin(2*pi*f*t)=(e^(j*2*pi*f*t)-e^(-j*2*pi*f*t))/2j.所以你的A*sin做FFT之后得到的幅值应该是在+/-(NFFT*f)的地方各有一个幅度为A*N

Y(:,2)表示的是Y的第二行的所有元素,出来的是一个行向量!矩阵的这种描述还有如下形式：Y(k,:),表示的Y的第k列Y(i:j,k:l),表示的事Y的第i-j行,与k-l列的交叉处得数据形成的矩阵

y(:)表示y数组（或向量）中的每一个元素y=y(:);这句话就是把y数组中的每个元素赋给y数组,看起来是一句废话,但是MATLAB中默认数组都是列向量（列数组）,似乎不管原来y是行向量还是列向量,赋

匿名函数

这个是做FFT时的点数,如果x长度比NFFT小,则在x后面补0,直到x长度等于NFFT

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

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

假如你变换后了,那么幅值=abs（FFT（y））；f=2*pi*k/N;(k=0、1、2、3、.、N-1）N是原始数据的个数.再问：那如果用plot这个函数画出频谱图，那横坐标怎么用原来的时间表示成频

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

X为信号,n为变换点数.Y=fft(X)是对信号X进行快速傅里叶变换；Y=fft(X,n)就是对信号X的前n个点进行快速傅里叶变换,如果n大于x的点数,则直接取前n个点,若小于n,则X先进行补零扩展为

将y=x^y转换成x=y^(1/y),然后用绘图命令,即可得到图形

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

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