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 - Русский
Scilabヘルプ >> Signal Processing > Spectral estimation > czt

czt

チャープz変換アルゴリズム

呼び出し手順

[czx]=czt(x,m,w,phi,a,theta)

引数

x

入力データ

m

cztはz平面のm 個の点で評価されます

w

振幅乗数

phi

位相増分

a

振幅初期値

theta

位相初期値

czx

チャープz変換出力

説明

k=0,1,...,m-1として 点[a*exp(j*theta)][w^kexp(j*k*phi)]における z平面の螺旋上でz変換を計算するチャープz変換アルゴリズム.

//カレントのコンテキストを保存
ax=gca(); [ffr,bds]=(ax.axes_bounds, ax.data_bounds(:)');
//
rect=[-1.2,-1.2*sqrt(2),1.2,1.2*sqrt(2)];
a = .7*exp(%i*%pi/6);
t=2*%pi*(0:179)/179;xsetech([0,0,0.5,1]);
plot2d(sin(t)',cos(t)',[2],"012",' ',rect)
plot2d([0 real(a)]',[0 imag(a)]',[3],"000")
xsegs([-1.0,0;1.0,0],[0,-1.0;0,1.0])
w0=.93*exp(-%i*%pi/15);w=exp(-(0:9)*log(w0));z=a*w;
zr=real(z);zi=imag(z);
plot2d(zr',zi',[5],"000")
xsetech([0.5,0,0.5,1]);
plot2d(sin(t)',cos(t)',[2],"012",' ',rect)
plot2d([0 real(a)]',[0 imag(a)]',[-1],"000")
xsegs([-1.0,0;1.0,0],[0,-1.0;0,1.0])
w0=w0/(.93*.93);w=exp(-(0:9)*log(w0));z=a*w;
zr=real(z);zi=imag(z);
plot2d(zr',zi',[5],"000")
xsetech(ffr,bds); //コンテキストを回復
Scilab Enterprises
Copyright (c) 2011-2017 (Scilab Enterprises)
Copyright (c) 1989-2012 (INRIA)
Copyright (c) 1989-2007 (ENPC)
with contributors
Last updated:
Thu Feb 14 15:02:12 CET 2019