spectrum.m

consider the power and spectrum

function spectrum
ts = 0.001;
fs = 1/ts;
t = [0:ts:10];
x = cos(2*pi*47*t) + cos(2*pi*219*t);
p = spower(x);
% pause
psd = spectrum(x);
p
specplot(psd, fs)