matlab 一维傅里叶变换

gaohaikuo1

1. data=load('bbb.txt');
   
2. t=data(:,1);%设置步进与时间区间
   
3. y=data(:,2);%根据已知写出信号的表达式
   
4. subplot(121);
   
5. plot(t,y);
   
6. 
   
7. fs=1/(t(2)-t(1));
   
8. Y=fft(y,512);%512点的FFT
   
9. P=Y.*conj(Y)/512;%求功率
   
10. f=(0:255)*fs/256%由上面的分析可知，频谱关于奈奎斯特频率对称，所以取其中一半
    
11. subplot(122);
    
12. plot(f,P(1:256))%功率随频率的变化，即功率谱图,绘制出一半
    
13. ylim([0,1e-19]);

复制代码