求叠加了高斯噪声的正弦信号的功率谱密度

clear all
t=0:0.001:0.6;
x=0.4*sin(2*pi*50*t)+0.4*sin(2*pi*320*t);
y=x+randn(size(t));
subplot(211);
plot(t(1:100),y(1:100))
title('0均值的随机信号');
xlabel('时间（秒）')
Nf=length(t);
F=fft(y,Nf);
Pyy=abs(F).^2/Nf;
f=1000*(0:(Nf-1)/2)/Nf;//序列间隔为fs/Nf
size(Pyy)
subplot(212);
plot(f,Pyy(1:(Nf-1)/2+1));
title('信号的功率谱密度');
xlabel('频率(Hz)')