MATLAB FFT 信号处理

Y=fft(y,1024);
pyy=Y.*conj(Y)/1024;
这两句讲的是对随机信号y的功率谱估计，可以参考数字信号处理中功率谱估计直接法部分，这个是定义式，要说由来要看 胡广书著《数字信号处理理论、算法与实现》的10.4节、11.2.1节，其他相关函数、平稳随机信号的功率谱的知识也可以看看便于理解

f=1000*(0:512)/1024;
这句讲的是构建真实的频率轴，这个频率轴是0Hz ～ 512*(1000/1024)Hz，要理解这个就要理解DFT变换的定义的推导过程，简单的说对时域进行以Ts为间隔的采样相当于对频谱进行fs＝1/Ts的周期延拓，反之对频域进行以fs为间隔的采样相当于对时域波形进行Ts＝1/fs的周期延拓，这个程序中Ts＝0.001，对应fs＝1000，也就是频谱一个周期长度为1000Hz，由于进行的是1024点的fft就相当于频率采样间隔为1000/1024，再来，程序中只想画出正频段部分因此就只去了0～512这几个点，DFT知识可以参考胡广书著《数字信号处理理论、算法与实现》3.5节
求采纳为满意回答。

基本正确吧，其实不用横坐标，直接就能得到傅里叶谱。采样频率直接决定最大频率，所以采样数不变时，改变采样频率肯定改变了频率的振幅。