Scilab 5.4.0
Please note that the recommended version of Scilab is 2025.0.0. This page might be outdated.
See the recommended documentation of this function
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*%asn(1,m)-n*%k(m1) //チェック sample=0:0.01:20; //等高線に正の実軸をマップ... z=alpha*%asn(sample/omegac,m)+Beta*ones(sample); plot(sample,ell1mag(epsilon,m1,z))
参照
- buttmag — バターワースフィルタの応答
Report an issue | ||
<< convol | filters | eqfir >> |