summaryrefslogtreecommitdiff
path: root/macros/pyulear.sci
diff options
context:
space:
mode:
Diffstat (limited to 'macros/pyulear.sci')
-rw-r--r--macros/pyulear.sci2
1 files changed, 0 insertions, 2 deletions
diff --git a/macros/pyulear.sci b/macros/pyulear.sci
index 55278d8..1c031de 100644
--- a/macros/pyulear.sci
+++ b/macros/pyulear.sci
@@ -1,5 +1,4 @@
function [psd,f_out] = pyulear(x,poles,freq,Fs,range,method,plot_type)
-
//Calculates a Yule-Walker autoregressive (all-pole) model of the data "x" and computes the power spectrum of the model.
//Calling Sequence
//[psd,f_out] = pyulear(x,poles,freq,Fs,range,method,plot_type)
@@ -12,7 +11,6 @@ function [psd,f_out] = pyulear(x,poles,freq,Fs,range,method,plot_type)
// range: 'half', 'onesided' : frequency range of the spectrum is from zero up to but not including sample_f/2. Power from negative frequencies is added to the positive side of the spectrum. 'whole', 'twosided' : frequency range of the spectrum is -sample_f/2 to sample_f/2, with negative frequencies stored in "wrap around" order after the positive frequencies; e.g. frequencies for a 10-point 'twosided' spectrum are 0 0.1 0.2 0.3 0.4 0.5 -0.4 -0.3 -0.2 -0.1 'shift', 'centerdc' : same as 'whole' but with the first half of the spectrum swapped with second half to put the zero-frequency value in the middle. (See "help fftshift". If "freq" is vector, 'shift' is ignored. If model coefficients "ar_coeffs" are real, the default range is 'half', otherwise default range is 'whole'.
// method: 'fft': use FFT to calculate power spectral density. 'poly': calculate spectral density as a polynomial of 1/z N.B. this argument is ignored if the "freq" argument is a vector. The default is 'poly' unless the "freq" argument is an integer power of 2.
// plot_type: 'plot', 'semilogx', 'semilogy', 'loglog', 'squared' or 'db' specifies the type of plot. The default is 'plot', which means linear-linear axes. 'squared' is the same as 'plot'. 'dB' plots "10*log10(psd)". This argument is ignored and a spectrum is not plotted if the caller requires a returned value.
-
//Description
//This function is being called from Octave.
//This function is a wrapper for aryule and ar_psd.