www.gusucode.com > signal 案例源码程序 matlab代码 > signal/PseudospectrumOfRealSignalExample.m
%% Pseudospectrum of Real Signal % Generate a real signal that consists of the sum of two sinusoids embedded % in white Gaussian noise of unit variance. The signal is sampled at 100 Hz % for 1 second. The sinusoids have frequencies of 25 Hz and 35 Hz. The % lower-frequency sinusoid has twice the amplitude of the other. %% fs = 100; t = 0:1/fs:1-1/fs; s = 2*sin(2*pi*25*t)+sin(2*pi*35*t)+randn(1,100); %% % Use the eigenvector method to compute the pseudospectrum of the signal % between 0 and the Nyquist frequency. Specify a signal subspace dimension % of 2 and a DFT length of 512. peig(s,2,512,fs,'half') %% % It is not possible to resolve the two sinusoids because the signal is % real. Repeat the computation using a signal subspace of dimension 4. peig(s,4,512,fs,'half')