Scilab Home page | Wiki | Bug tracker | Forge | Mailing list archives | ATOMS | File exchange
Please login or create an account
Scilab 5.5.1
Change language to: English - Français - Português - 日本語 -

Please note that the recommended version of Scilab is 6.1.0. This page might be outdated.
See the recommended documentation of this function

mese

maximum entropy spectral estimation

Calling Sequence

`[sm,fr]=mese(x [,npts]);`

Arguments

x

Input sampled data sequence

npts

Optional parameter giving number of points of `fr` and `sm` (default is 256)

sm

Samples of spectral estimate on the frequency grid `fr`

fr

npts equally spaced frequency samples in `[0,.5)`

Description

Calculate the maximum entropy spectral estimate of `x`

Examples

```rand("seed",12345),
y = (0:10);
x1 = sin(2*%pi*y/20);
x2 = sin(3.5*%pi*y/20);
w = .4*(rand(y)-.5*ones(y));
x = x1+2*x2+w;
[sm,fr1] = mese(x,10)
// Define macro which computes the
// Blackman-Tukey periodogram
function [xm, fr]=bt(x)
xsize = max(size(x))
[xf,fr] = frmag(x,256)
xm = xf.*conj(xf)/xsize
endfunction
[xm,fr] = bt(x);
scf();
subplot(3,1,1);
plot(x);
subplot(3,1,2);
plot(fr1,sm)
subplot(3,1,3);
plot(fr,xm)```

Comments

Add a comment:
Please login to comment this page.

 Report an issue << intdec spectral_estimation pspect >>

 Scilab EnterprisesCopyright (c) 2011-2017 (Scilab Enterprises)Copyright (c) 1989-2012 (INRIA)Copyright (c) 1989-2007 (ENPC)with contributors Last updated:Thu Oct 02 14:01:06 CEST 2014