FFT及周期图法算功率谱密度

t=0:0.0001:0.1;          %时间间隔为0.0001，说明采样频率为10000Hz
x=square(2*pi*1000*t);   %产生基频为1000Hz的方波信号
n=randn(size(t));        %白噪声
f=x;                   %在信号中加入白噪声
figure(1); 
subplot(2,1,1);
plot(t,f);                 %画出原始信号的波形图
ylabel('幅值(V)');
xlabel('时间(s)');
title('原始信号');
y=fft(f,1024);           %对原始信号进行离散傅里叶变换，参加DFT采样点的个数为1024
subplot(2,1,2);
m=abs(y);
f1=(0:length(y)/2-1)'*10000/length(y);%计算变换后不同点对应的幅值 
plot(f1,m(1:length(y)/2)*2/length(y));
ylabel('幅值的模');xlabel('时间(s)');title('原始信号傅里叶变换');
%用周期图法估计功率谱密度
p=y.*conj(y)/1024;      %计算功率谱密度   
ff=10000*(0:499)/1024;  %计算变换后不同点对应的频率值 
figure(2);plot(ff,p(1:500));ylabel('幅值');xlabel('频率(Hz)');
title('功率谱密度（周期图法）'); 

求FFT变换时，纵坐标和横坐标的问题

clear all;clc
fs=500;
t=(0:1:100)/fs;
s=cos(2*pi*50*t);
n=length(t);
figure(1),subplot(211);plot(t,s)
f=fftshift(fft(s));
subplot(212);plot((-n/2:n/2-1)*fs/n,abs(f)*2/n);