Scilab Home page | Wiki | Bug tracker | Forge | Mailing list archives | ATOMS | File exchange
Please login or create an account
Change language to: English - Français - Português - Русский

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

Scilab help >> Signal Processing > spectral_estimation > mese

mese

スペクトル推定の最大エントロピ

呼び出し手順

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

引数

x

入力標本化データ系列

npts

fr および smの点の数 (デフォルトは256)を指定するオプションパラメータ

sm

周波数グリッドfrにおける スペクトル推定の標本

fr

[0,.5)の範囲のnpts個の等間隔周波数の標本

説明

xの最大エントロピスペクトル推定値を計算します.

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)
// Blackman-Tukey ペリオドグラムを計算するマクロを定義
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)
Scilab Enterprises
Copyright (c) 2011-2017 (Scilab Enterprises)
Copyright (c) 1989-2012 (INRIA)
Copyright (c) 1989-2007 (ENPC)
with contributors
Last updated:
Mon Oct 01 17:40:28 CEST 2012