13t.m

来自「包含50个数字信号处理方面的matlab程序」· M 代码 · 共 14 行

t=0:0.001:1;
x=sin(2*pi*60*t)+sin(2*pi*200*t);
subplot(2,1,1);stem(t,x,'.');
title('signial x(n)');grid;
y=x+1.5*randn(1,length(t));
Y=fft(y,1024);
p=Y.*conj(Y)/1024;
N=1:1024;n=N/1000*1024;
subplot(2,1,2);plot(n,p);
axis([0,600,0,280]);
title('signial y(n)');grid;
xlabel('频率（单位：Hz)');
ylabel('p');grid