Change language to:
English - Français - 日本語 - Русский

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

# mese

maximum entropy spectral estimation

### Syntax

`[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)```
 Report an issue << czt Spectral estimation pspect >>

 Copyright (c) 2022-2023 (Dassault Systèmes)Copyright (c) 2017-2022 (ESI Group)Copyright (c) 2011-2017 (Scilab Enterprises)Copyright (c) 1989-2012 (INRIA)Copyright (c) 1989-2007 (ENPC)with contributors Last updated:Tue Feb 25 08:52:30 CET 2020