带注解的快速傅里叶变换(fft)matlab程序

Fs=128;%采样频率T=1/Fs;%采样时间L=256;%信号长度t=(0:L-1)*T;%时间x=5+7*cos(2*pi*15*t-30*pi/180)+3*cos(2*pi*40*t-90*pi/180);%cos为底原始信号y=x+randn(size(t));%添加噪声figure;plot(t,y)title('加噪声的信号')xlabel('时间(s)')N=2^nextpow2(L);%采样点数，采样点数越大，分辨的频率越精确，N=L，超出的部分信号补为0Y=fft(y,N)/N*2;%除以N乘以2才是真实幅值，N越大，幅值精度越高f=Fs/N*(0:1:N-1);%频率A=abs(Y);%幅值P=angle(Y);%相值figure;subplot(211);plot(f(1:N/2),A(1:N/2));%函数fft返回值的数据结构具有对称性,因此我们只取前一半title('幅值频谱');xlabel('频率(Hz)');ylabel('幅值');subplot(212);plot(f(1:N/2),P(1:N/2));title('相位谱频');xlabel('频率(Hz)');ylabel('相位');