이 페이지의 최신 내용은 아직 번역되지 않았습니다. 최신 내용은 영문으로 볼 수 있습니다.

ellip

타원 필터 설계

설명

예제

[b,a] = ellip(n,Rp,Rs,Wp)는 정규화된 통과대역 경계 주파수 Wp를 갖는 n차 저역통과 디지털 타원 필터의 전달 함수 계수를 반환합니다. 결과로 생성되는 필터는 Rp(단위: 데시벨)의 피크 간 통과대역 리플과 피크 통과대역 값을 기준으로 한 Rs(단위: 데시벨)의 저지대역 감쇠량을 갖습니다.

예제

[b,a] = ellip(n,Rp,Rs,Wp,ftype)ftype의 값과 Wp의 요소 개수에 따라 저역통과, 고역통과, 대역통과 또는 대역저지 타원 필터를 설계합니다. 결과로 생성되는 대역통과 설계와 대역저지 설계는 차수가 2n입니다.

참고:   전달 함수를 구성하는 데 영향을 미치는 수치적 문제에 대한 자세한 내용은 제한 사항 항목을 참조하십시오.

예제

[z,p,k] = ellip(___)은 저역통과, 고역통과, 대역통과 또는 대역저지 디지털 타원 필터를 설계하고 이 필터의 영점, 극점, 이득을 반환합니다. 이 구문은 위에 열거한 구문에 있는 어떤 입력 인수도 포함할 수 있습니다.

예제

[A,B,C,D] = ellip(___)은 저역통과, 고역통과, 대역통과 또는 대역저지 디지털 타원 필터를 설계하고 이 필터의 상태공간 표현을 지정하는 행렬을 반환합니다.

예제

[___] = ellip(___,'s')는 통과대역 경계 각주파수 Wp, Rp(단위: 데시벨)의 통과대역 리플, Rs(단위: 데시벨)의 저지대역 감쇠량을 갖는 저역통과, 고역통과, 대역통과 또는 대역저지 아날로그 타원 필터를 설계합니다.

예제

모두 축소

5dB의 통과대역 리플 및 40dB의 저지대역 감쇠량과, 1000Hz로 샘플링된 데이터에 대해 0.6π rad/sample에 해당하는 수치인 300Hz의 통과대역 경계 주파수를 갖는 6차 저역통과 타원 필터를 설계합니다. 필터의 크기 응답과 위상 응답을 플로팅합니다. 이를 사용하여 1000개 샘플로 구성된 랜덤 신호를 필터링합니다.

[b,a] = ellip(6,5,40,0.6);
freqz(b,a)

dataIn = randn(1000,1);
dataOut = filter(b,a,dataIn);

0.2π rad/sample과 0.6π rad/sample의 정규화된 경계 주파수, 5dB의 통과대역 리플, 50dB의 저지대역 감쇠량을 갖는 6차 타원 대역저지 필터를 설계합니다. 필터의 크기 응답과 위상 응답을 플로팅합니다. 이를 사용하여 랜덤 데이터를 필터링합니다.

[b,a] = ellip(3,5,50,[0.2 0.6],'stop');
freqz(b,a)

dataIn = randn(1000,1);
dataOut = filter(b,a,dataIn);

1000Hz로 샘플링된 데이터에 대해 0.6π rad/sample에 해당하는 수치인 300Hz의 통과대역 경계 주파수를 갖는 6차 고역통과 타원 필터를 설계합니다. 통과대역 리플을 3dB로 지정하고 50dB의 저지대역 감쇠량을 50dB로 지정합니다. 크기 응답과 위상 응답을 플로팅합니다. fvtool에서 사용할 수 있도록 영점, 극점, 이득을 2차섹션형(SOS)으로 변환합니다.

[z,p,k] = ellip(6,3,50,300/500,'high');
sos = zp2sos(z,p,k);
fvtool(sos,'Analysis','freq')

500Hz의 저역 통과대역 주파수와 560Hz의 고역 통과대역 주파수를 갖는 20차 타원 대역통과 필터를 설계합니다. 통과대역 리플을 3dB로, 저지대역 감쇠량을 40dB로, 샘플 레이트를 1500Hz로 지정합니다. 상태공간 표현을 사용합니다. designfilt를 사용하여 동일한 필터를 설계합니다.

[A,B,C,D] = ellip(10,3,40,[500 560]/750);
d = designfilt('bandpassiir','FilterOrder',20, ...
    'PassbandFrequency1',500,'PassbandFrequency2',560, ...
    'PassbandRipple',3, ...
    'StopbandAttenuation1',40,'StopbandAttenuation2',40, ...
    'SampleRate',1500);

상태공간 표현을 2차섹션형(SOS)으로 변환합니다. fvtool을 사용하여 주파수 응답을 시각화합니다.

sos = ss2sos(A,B,C,D);
fvt = fvtool(sos,d,'Fs',1500);
legend(fvt,'ellip','designfilt')

2GHz의 차단 주파수를 갖는 5차 아날로그 버터워스 저역통과 필터를 설계합니다. 2π를 곱하여 주파수를 초당 라디안 값으로 변환합니다. 4096개 점에서 필터의 주파수 응답을 계산합니다.

n = 5;
f = 2e9;

[zb,pb,kb] = butter(n,2*pi*f,'s');
[bb,ab] = zp2tf(zb,pb,kb);
[hb,wb] = freqs(bb,ab,4096);

동일한 경계 주파수와 3dB의 통과대역 리플을 갖는 5차 체비쇼프 유형 I 필터를 설계합니다. 필터의 주파수 응답을 계산합니다.

[z1,p1,k1] = cheby1(n,3,2*pi*f,'s');
[b1,a1] = zp2tf(z1,p1,k1);
[h1,w1] = freqs(b1,a1,4096);

동일한 경계 주파수와 30dB의 저지대역 감쇠량을 갖는 5차 체비쇼프 유형 II 필터를 설계합니다. 필터의 주파수 응답을 계산합니다.

[z2,p2,k2] = cheby2(n,30,2*pi*f,'s');
[b2,a2] = zp2tf(z2,p2,k2);
[h2,w2] = freqs(b2,a2,4096);

동일한 경계 주파수, 3dB의 통과대역 리플, 30dB의 저지대역 감쇠량을 갖는 5차 타원 필터를 설계합니다. 필터의 주파수 응답을 계산합니다.

[ze,pe,ke] = ellip(n,3,30,2*pi*f,'s');
[be,ae] = zp2tf(ze,pe,ke);
[he,we] = freqs(be,ae,4096);

감쇠량(단위: 데시벨)을 플로팅합니다. 주파수를 기가헤르츠 단위로 표현합니다. 필터를 비교합니다.

plot(wb/(2e9*pi),mag2db(abs(hb)))
hold on
plot(w1/(2e9*pi),mag2db(abs(h1)))
plot(w2/(2e9*pi),mag2db(abs(h2)))
plot(we/(2e9*pi),mag2db(abs(he)))
axis([0 4 -40 5])
grid
xlabel('Frequency (GHz)')
ylabel('Attenuation (dB)')
legend('butter','cheby1','cheby2','ellip')

버터워스 필터와 체비쇼프 유형 II 필터는 평탄한 통과대역과 넓은 천이 대역을 가집니다. 체비쇼프 유형 I 필터와 타원 필터는 더 빨리 롤오프되지만 통과대역 리플을 가집니다. 체비쇼프 유형 II 설계 함수에 대한 주파수 입력값은 통과대역의 끝값이 아니라 저지대역의 시작값을 설정합니다.

입력 인수

모두 축소

필터 차수로, 정수 스칼라로 지정됩니다.

데이터형: double

피크 간 통과대역 리플로, 양의 스칼라로 지정됩니다(단위: 데시벨).

사양 ℓ이 선형 단위인 경우 Rp = 40 log10((1+ℓ)/(1–ℓ))을 사용하여 데시벨로 변환할 수 있습니다.

데이터형: double

피크 통과대역 값을 기준으로 한 저지대역 감쇠량으로, 데시벨을 단위로 하는 양의 스칼라로 지정됩니다.

사양 ℓ이 선형 단위인 경우 Rs = –20log10ℓ을 사용하여 데시벨로 변환할 수 있습니다.

데이터형: double

통과대역 경계 주파수로, 스칼라나 요소를 2개 가진 벡터로 지정됩니다. 통과대역 경계 주파수는 필터의 크기 응답이 –Rp(단위:데시벨)인 주파수입니다. 통과대역 리플 Rp의 값이 작고 저지대역 감쇠량 Rs의 값이 클수록 천이 대역이 더 넓어집니다.

  • Wp가 스칼라인 경우, ellip은 경계 주파수 Wp를 갖는 저역통과 필터나 고역통과 필터를 설계합니다.

    Wp가 요소를 2개 가진 벡터 [w1 w2](여기서 w1 < w2)이면 ellip은 하한 경계 주파수 w1과 상한 경계 주파수 w2를 갖는 대역통과 필터나 대역저지 필터를 설계합니다.

  • 디지털 필터의 경우, 통과대역 경계 주파수는 0과 1 사이에 있어야 합니다. 여기서 1은 샘플 레이트의 절반, 즉 π rad/sample인 나이퀴스트 레이트에 해당합니다.

    아날로그 필터의 경우, 통과대역 경계 주파수는 초당 라디안으로 표현되어야 하고 모든 양의 값을 받을 수 있습니다.

데이터형: double

필터 유형으로, 다음 중 하나로 지정됩니다.

  • 'low'는 통과대역 경계 주파수 Wp를 갖는 저역통과 필터를 지정합니다. 'low'는 스칼라 Wp의 디폴트 값입니다.

  • 'high'는 통과대역 경계 주파수 Wp를 갖는 고역통과 필터를 지정합니다.

  • 'bandpass'Wp가 요소를 2개 가진 벡터인 경우 차수가 2n인 대역통과 필터를 지정합니다. 'bandpass'Wp가 2개의 요소를 가지는 경우 디폴트 값입니다.

  • 'stop'Wp가 요소를 2개 가진 벡터인 경우 차수가 2n인 대역저지 필터를 지정합니다.

출력 인수

모두 축소

필터의 전달 함수 계수로, 저역통과 필터와 고역통과 필터에 대해서는 길이가 n + 1인 행 벡터로 반환되고 대역통과 필터와 대역저지 필터에 대해서는 길이가 2n + 1인 행 벡터로 반환됩니다.

  • 디지털 필터의 경우, 전달 함수는 다음과 같이 ba로 표현됩니다.

    H(z)=B(z)A(z)=b(1)+b(2)z1++b(n+1)zna(1)+a(2)z1++a(n+1)zn.

  • 아날로그 필터의 경우, 전달 함수는 다음과 같이 ba로 표현됩니다.

    H(s)=B(s)A(s)=b(1)sn+b(2)sn1++b(n+1)a(1)sn+a(2)sn1++a(n+1).

데이터형: double

필터의 영점, 극점, 이득으로, 길이가 n(대역통과 설계와 대역저지 설계의 경우 2n임)인 두 개의 열 벡터와 하나의 스칼라로 반환됩니다.

  • 디지털 필터의 경우, 전달 함수는 다음과 같이 z, p, k로 표현됩니다.

    H(z)=k(1z(1)z1)(1z(2)z1)(1z(n)z1)(1p(1)z1)(1p(2)z1)(1p(n)z1).

  • 아날로그 필터의 경우, 전달 함수는 다음과 같이 z, p, k로 표현됩니다.

    H(s)=k(sz(1))(sz(2))(sz(n))(sp(1))(sp(2))(sp(n)).

데이터형: double

필터의 상태공간 표현으로, 행렬로 반환됩니다. 저역통과 설계와 고역통과 설계에 대해 m = n이고 대역통과 필터와 대역저지 필터에 대해 m = 2n이면 A는m × m,B는m × 1,C는 1 × m, 그리고D는 1 × 1입니다.

  • 디지털 필터의 경우, 상태공간 행렬은 상태 벡터 x, 입력값 u, 출력값 y와 다음 관계를 가집니다.

    x(k+1)=Ax(k)+Bu(k)y(k)=Cx(k)+Du(k).

  • 아날로그 필터의 경우, 상태공간 행렬은 상태 벡터 x, 입력값 u, 출력값 y와 다음 관계를 가집니다.

    x˙=Ax+Buy=Cx+Du.

데이터형: double

세부 정보

모두 축소

제한 사항

전달 함수 구문의 수치적 불안정성

일반적으로, [z,p,k] 구문을 사용하여 IIR 필터를 설계하십시오. 그런 다음, zp2sos[z,p,k] 출력값을 사용하여 필터를 분석하거나 구현할 수 있습니다. [b,a] 구문을 사용하여 필터를 설계하면 수치적 문제가 발생할 수 있습니다. 이러한 문제는 반올림 오차로 인해 발생하며, n이 4 정도로 낮을 때 발생할 수 있습니다. 다음 예제에서는 이러한 한계를 보여줍니다.

n = 6; 
Rp = 0.1;
Rs = 80;
Wn = [2.5e6 29e6]/500e6; 
ftype = 'bandpass'; 

% Transfer Function design 
[b,a] = ellip(n,Rp,Rs,Wn,ftype);            % This filter is unstable 

% Zero-Pole-Gain design 
[z,p,k] = ellip(n,Rp,Rs,Wn,ftype);
sos = zp2sos(z,p,k);

% Plot and compare the results
hfvt = fvtool(b,a,sos,'FrequencyScale','log'); 
legend(hfvt,'TF Design','ZPK Design')

알고리즘

타원 필터는 버터워스 필터나 체비쇼프 필터보다 더욱 가파른 롤오프 특성을 제공하지만, 통과대역과 저지대역 모두에 등리플이 존재합니다. 일반적으로, 타원 필터는 모든 필터 유형의 가장 낮은 차수로 주어진 성능 사양을 충족합니다.

ellip은 다음 5단계의 알고리즘을 사용합니다.

  1. 함수 ellipap를 사용하여 저역통과 아날로그 프로토타입 극점, 영점, 이득을 구합니다.

  2. 극점, 영점, 이득을 상태공간 형식으로 변환합니다.

  3. 필요한 경우, 상태공간 변환을 사용하여 원하는 주파수 제약 조건을 갖는 대역통과 필터, 고역통과 필터 또는 대역저지 필터로 저역통과 필터를 변환합니다.

  4. 디지털 필터 설계의 경우, bilinear를 사용하여 주파수 사전굽힘(Prewarping)을 사용한 쌍선형 변환을 통해 아날로그 필터를 디지털 필터로 변환합니다. 주파수를 세심하게 조정하여 아날로그 필터와 디지털 필터가 Wp 또는 w1w2에서 동일한 주파수 응답 크기를 가지도록 할 수 있습니다.

  5. 필요한 경우, 상태공간 필터를 다시 전달 함수나 영점-극점-이득 형식으로 변환합니다.

확장 기능

R2006a 이전에 개발됨