Spektrálna analýza Parametric a non-parametrické metódy

[seifer]

Zdravím všetky ima začiatočník v MATLAB a potřebuju trochu pomoct
Chcem porovnať parametrické a non-parametrické metódy
Potrebujem kód signál a potom kód, ako urobiť psd tohto signálu s parametrické alebo non-parametrická metóda .. a uvidíte, ktorý spôsob je lepší ... já jen príklad .. nádej som sám jasné
ďakujem všetkým vopred!

 

[asjohnas]

V MATLAB môžete použiť spectrum.music alebo psd funkcie.
Potrebujete typu v MATLAB
doc spectrum.music
alebo
doc psd

 

[seifer]

Avšak aj tento súbor našiel ... ale jak jsem řekl ima začiatočník v MATLAB ... a zdá trochu zložité understant .. to je aj dôvod, prečo potrebujete niekoho, aby mi sprievodca ..
Potrebujem signál ... i, ktoré možno ľahko kedykoľvek zmeniť N (vzorky) a frekvencia a potom sa toolbox vidieť psd ... i need príklad takhle!
thx za Vašu odpoveď!

 

[asjohnas]

Napríklad môžete vykonávať niektoré jednoduché, ako je tento:
len = 1000;
fd = 10;
t = 0:1 / fd: len/fd-1/fd;
s = sin (t * 3 * 2 * pi) randn (1, len);
psd (s, len, fd)

 

[seifer]

môže u dať rovnaký príklad s SALPE toolbox ... som pripravený danej kódy:
[y, t] = sg_cissoid (beta, f, Dim, <OSF>, <SNR>)
%
% Parametre: - beta: Vektorový komplexné amplitúdy
% - F: Vektorové s frekvenciou hodnôt (-0,5 <f <0,5)
% - Dim: Počet vzoriek
% - OSF: převzorkování koeficient (predvolené: OSF = 1)
% - SNR: odstupom signálu od šumu v [dB] (predvolené: SNR = inf)
%
% Výstup: - y: vektor sa vzorky signálu
% - T: Časové osi
%
% Príklad: sg_cissoid ([1 0.3 0.5 0.7], [0.1 0.2 0.25 0.48], 100, 5, 30);

a vhodné psd
Volajte: [PSD, f] = psd_welch (y, <Win>, <N_Overlap>, <N_FFT>, <fs>)
%
% Parametre: - y: vektor sa vzorky signálu
% - Win: Signál je rozdelený do K_seg prekrývajúce
% Segmenty o dĺžke rovnajúcej sa dĺžke
% Na okná a okná každý segment
% (Predvolené: 8 segmentov s Hamming okne)
% - N_Overlap: Počet prekrývaniu vzoriek zo segmentu
% Segmentu (predvolené: 50% prekrytie)
% - N_FFT: Počet bodov FFT (predvolené: N_FFT = 256)
% - Fs: Vzorkovacia frekvencia v Hz (predvolené: fs = 1 Hz)
%
% Výstup: - PSD: spektrálna hustota výkonu
% - F: Frekvencia osi
%
% Príklad: y = sg_cissoid ([1 0.5 0.25], [0.2 0.3 0.4], 100, 1, 10);
% Psd_welch (y, hanna (50), 49, 1024);