ell1mag

楕円フィルタの振幅

呼び出し手順

[v]=ell1mag(eps,m1,z)

引数

eps: 通過域のリップル=1/(1+eps^2)
m1: 阻止域のリップル=1/(1+(eps^2)/m1)
z: 複素平面上の標本ベクトルの値
v: 標本点における楕円フィルタの値

説明

楕円フィルタの二乗振幅を得る際に使用される関数. 通常は,m1=eps*eps/(a*a-1)です. s=%sn(z,m1)として, v=real(ones(z)./(ones(z)+eps*eps*s.*s))を返します.

例

deff('[alpha,BeTa]=alpha_beta(n,m,m1)',...
'if 2*int(n/2)==n then, BeTa=K1; else, BeTa=0;end;'+...
'alpha=%k(1-m1)/%k(1-m);')
epsilon=0.1;A=10;  //リップルパラメータ
m1=(epsilon*epsilon)/(A*A-1);n=5;omegac=6;
m=find_freq(epsilon,A,n);omegar = omegac/sqrt(m)
%k(1-m1)*%k(m)/(%k(m1)*%k(1-m))-n   //チェック...
[alpha,Beta]=alpha_beta(n,m,m1)
alpha*delip(1,sqrt(m))-n*%k(m1)      //チェック
sample=0:0.01:20;
//等高線に正の実軸をマップ...
z=alpha*delip(sample/omegac,sqrt(m))+Beta*ones(sample);
plot(sample,ell1mag(epsilon,m1,z))

参照

buttmag — バターワースフィルタの応答

Report an issue
<< convol	Filters	eqfir >>

Copyright (c) 2022-2024 (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:
Mon Feb 12 23:12:40 CET 2018