spectrum.m

来自「基于半导体激光器的matlab程序产生混沌」· M 代码 · 共 10 行

function [fw,pw]=spectrum(fs,y)

%yt=sqrt(y);
nfft=2^(nextpow2(length(y)));                         %采样点个数
nup=ceil((nfft+1)/2);  %the number of unique points 'NumUniquePts'=nup
yw=fft(y,nfft);
fw=fs*(0:nup-1)/nfft;  %nup=(nfft/2)+1
yw=yw(1:nup);   
pw=yw.*conj(yw)/length(y);
pw=10*log10(pw);