powerspectrum.m

来自「是“MATLAB混合编程与工程应用”一书的源码」· M 代码 · 共 18 行

function [m]=PowerSpectrum()
%global x;
%Original Signal 
global y; 
%Noisy Signal
if (isempty(y)) 
m =[];
return;
end
f=1000*(0:256)/512;
y = fft(y,512);
Pyy=y.* conj(y)/512;
subplot(212)
plot(f,Pyy(1:257))
title('Frequency content of y')
xlabel('frequency (Hz)')
m=Pyy;