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); //コンテキストを回復
| Report an issue | ||
| << cspect | Spectral estimation | mese >> |