Scilab Website | Contribute with GitLab | Mailing list archives | ATOMS toolboxes
Scilab Online Help
5.5.2 - Français

Change language to:
English - 日本語 - Português - Русский

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

Aide de Scilab >> Traitement du Signal > spectral_estimation > mese

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

Sample

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
<< intdec 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:
Wed Apr 01 10:21:40 CEST 2015