cheby1

R2026b

Chebyshev Type I filter design

Syntax

[b,a] = cheby1(n,Rp,Wp)

[b,a] = cheby1(n,Rp,Wp,fType)

[z,p,k] = cheby1(___)

[A,B,C,D] = cheby1(___)

[___] = cheby1(___,"s")

[B,A] = cheby1(n,Rp,Wp,"ctf")

[___] = cheby1(n,Rp,Wp,fType,"ctf")

[___,gS] = cheby1(___)

Description

[b,a] = cheby1(n,Rp,Wp) designs an nth-order lowpass digital Chebyshev Type I filter with normalized passband edge frequency Wp and Rp decibels of peak-to-peak passband ripple. The cheby1 function returns the numerator and denominator coefficients of the filter transfer function.

[b,a] = cheby1(n,Rp,Wp,fType) designs a lowpass, highpass, bandpass, or bandstop digital Chebyshev Type I filter, depending on the value of fType and the number of elements of Wp. The resulting bandpass and bandstop designs are of order 2n.

Note

You might encounter numerical instabilities when designing IIR filters with transfer functions for orders as low as 4. See Transfer Functions and CTF for more information about numerical issues that affect forming the transfer function.

[z,p,k] = cheby1(___) designs a digital Chebyshev Type I filter and returns its zeros, poles, and gain. This syntax can include any of the input arguments in previous syntaxes.

[A,B,C,D] = cheby1(___) designs a digital Chebyshev Type I filter and returns the matrices that specify its state-space representation.

[___] = cheby1(___,"s") designs an analog Chebyshev Type I filter using any of the input or output arguments in previous syntaxes.

[B,A] = cheby1(n,Rp,Wp,"ctf") designs a lowpass digital Chebyshev Type I filter using second-order Cascaded Transfer Functions (CTF). The function returns matrices that list the denominator and numerator polynomial coefficients of the filter transfer function, represented as a cascade of filter sections. This approach generates IIR filters with improved numerical stability compared to single-section transfer functions. (since R2024b)

[___] = cheby1(n,Rp,Wp,fType,"ctf") designs a lowpass, highpass, bandpass, or bandstop digital Chebyshev Type I filter, and returns the filter representation using the CTF format. The resulting design sections are of order 2 (lowpass and highpass filters) or 4 (bandpass and bandstop filters). (since R2024b)

[___,gS] = cheby1(___) also returns the overall gain of the system. You must specify "ctf" to return gS. (since R2024b)

Examples

Lowpass Chebyshev Type I Transfer Function

Design a 6th-order lowpass Chebyshev Type I filter with 10 dB of passband ripple and a passband edge frequency of 300 Hz, which, for data sampled at 1000 Hz, corresponds to $$0.6\pi$$ rad/sample. Plot its magnitude and phase responses. Use it to filter a 1000-sample random signal.

fc = 300;
fs = 1000;

[b,a] = cheby1(6,10,fc/(fs/2));

freqz(b,a,[],fs)

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

Bandstop Chebyshev Type I Filter

Design a 6th-order Chebyshev Type I bandstop filter with normalized edge frequencies of $$0.2\pi$$ and $$0.6\pi$$ rad/sample and 5 dB of passband ripple. Plot its magnitude and phase responses. Use it to filter random data.

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

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

Highpass Chebyshev Type I Filter

Design a 9th-order highpass Chebyshev Type I filter with 0.5 dB of passband ripple and a passband edge frequency of 300 Hz, which, for data sampled at 1000 Hz, corresponds to $$0.6\pi$$ rad/sample. Convert the zeros, poles, and gain to second-order sections. Plot the magnitude and phase responses of the filter.

[z,p,k] = cheby1(9,0.5,300/500,'high');
sos = zp2sos(z,p,k);
freqz(sos)

Bandpass Chebyshev Type I Filter

Design a 20th-order Chebyshev Type I bandpass filter with a lower passband frequency of 400 Hz and a higher passband frequency of 560 Hz. Specify a passband ripple of 3 dB and a sample rate of 1500 Hz. Use the state-space representation.

fs = 1500;

[A,B,C,D] = cheby1(10,3,[400 560]/(fs/2));

Convert the state-space representation to second-order sections and visualize the filter response.

sos = ss2sos(A,B,C,D);
freqz(sos,[],fs)

Design an identical filter using designfilt. Visualize the filter response.

d = designfilt("bandpassiir",FilterOrder=20, ...
    PassbandFrequency1=400,PassbandFrequency2=560, ...
    PassbandRipple=3,SampleRate=fs);
freqz(d,[],fs)

Comparison of Analog IIR Lowpass Filters

Design a fifth-order analog Butterworth lowpass filter with a cutoff frequency of 2 GHz. Multiply by $$2\pi$$ to convert the frequency to radians per second. Compute the frequency response of the filter at 4096 points.

n = 5;
wn = 2*pi*2e9;
w = 2*pi*1e9*logspace(-2,1,4096)';

[zb,pb,kb] = butter(n,wn,"s");
[bb,ab] = zp2tf(zb,pb,kb);
[hb,wb] = freqs(bb,ab,w);
gdb = -diff(unwrap(angle(hb)))./diff(wb);

Design a fifth-order Chebyshev Type I filter with a passband edge frequency of 2 GHz and 3 dB of passband ripple. Compute its frequency response.

wp = wn;

[z1,p1,k1] = cheby1(n,3,wp,"s");
[b1,a1] = zp2tf(z1,p1,k1);
[h1,w1] = freqs(b1,a1,w);
gd1 = -diff(unwrap(angle(h1)))./diff(w1);

Design a fifth-order Chebyshev Type II filter with a stopband edge frequency of 2.5 GHz and 30 dB of stopband attenuation. Compute its frequency response.

ws = 2*pi*2.5e9;

[z2,p2,k2] = cheby2(n,30,ws,"s");
[b2,a2] = zp2tf(z2,p2,k2);
[h2,w2] = freqs(b2,a2,w);
gd2 = -diff(unwrap(angle(h2)))./diff(w2);

Design a fifth-order elliptic filter with the same passband and stopband edge frequencies, 3 dB of passband ripple, and 30 dB of stopband attenuation. Compute its frequency response.

[ze,pe,ke] = ellip(n,3,30,wp,"s");
[be,ae] = zp2tf(ze,pe,ke);
[he,we] = freqs(be,ae,w);
gde = -diff(unwrap(angle(he)))./diff(we);

Design a fifth-order Bessel filter with the same edge frequency. Compute its frequency response.

[zf,pf,kf] = besself(n,wn);
[bf,af] = zp2tf(zf,pf,kf);
[hf,wf] = freqs(bf,af,w);
gdf = -diff(unwrap(angle(hf)))./diff(wf);

Plot the attenuation in decibels. Express the frequency in gigahertz. Compare the filters.

fGHz = [wb w1 w2 we wf]/(2e9*pi);
plot(fGHz,mag2db(abs([hb h1 h2 he hf])))
axis([0 5 -45 5])
grid on
xlabel("Frequency (GHz)")
ylabel("Attenuation (dB)")
legend(["butter" "cheby1" "cheby2" "ellip" "besself"], ...
    Location="southwest")

The Butterworth and Chebyshev Type II filters have flat passbands and wide transition bands. The Chebyshev Type I and elliptic filters roll off faster but have passband ripple. The frequency input to the Chebyshev Type II design function sets the beginning of the stopband rather than the end of the passband. Elliptic filters offer steeper rolloff characteristics than Butterworth and Chebyshev filters, but they are equiripple in both the passband and the stopband. Of these four classical filter types, elliptic filters usually meet a given set of filter performance specifications with the lowest filter order.

Plot the group delay in samples. Express the frequency in gigahertz and the group delay in nanoseconds. Compare the filters. The Bessel filter has approximately constant group delay along the passband.

gdns = [gdb gd1 gd2 gde gdf]*1e9;
gdns(gdns<0) = NaN;
loglog(fGHz(2:end,:),gdns)
grid on
xlabel("Frequency (GHz)")
ylabel("Group delay (ns)")
legend(["butter" "cheby1" "cheby2" "ellip" "besself"], ...
    Location="southwest")

Highpass Chebyshev Type I Filter in Cascade

Design a ninth-order highpass Chebyshev Type I filter with a cutoff frequency of 300 Hz and sampling rate of 1000 Hz. The passband ripple is 5 dB. Return the coefficients of the filter system as a cascade of second-order sections.

Wn = 300/(1000/2);
[B,A] = cheby1(9,5,Wn,"high","ctf")

B = 5×3

    0.1336   -0.1336         0
    0.1336   -0.2671    0.1336
    0.1336   -0.2671    0.1336
    0.1336   -0.2671    0.1336
    0.1336   -0.2671    0.1336

A = 5×3

    1.0000    0.9022         0
    1.0000    1.6117    0.8336
    1.0000    1.1993    0.8788
    1.0000    0.8296    0.9291
    1.0000    0.6335    0.9767

Plot the magnitude response of the filter.

filterAnalyzer(B,A)

Input Arguments

n — Filter order

integer scalar up to 500

Filter order, specified as an integer scalar less than or equal to 500. For bandpass and bandstop designs, n represents one-half the filter order.

Data Types: double

Rp — Peak-to-peak passband ripple in dB

positive scalar

Peak-to-peak passband ripple, specified as a positive scalar expressed in decibels.

If your specification, ℓ, is in linear units, you can convert it to decibels using Rp = 40 log₁₀((1+ℓ)/(1–ℓ)).

Data Types: double

Wp — Passband edge frequency

scalar | two-element vector

Passband edge frequency, specified as a scalar or a two-element vector. The passband edge frequency is the frequency at which the magnitude response of the filter is –Rp decibels. Smaller values of passband ripple, Rp, result in wider transition bands.

• If Wp is a scalar, then cheby1 designs a lowpass or highpass filter with edge frequency Wp.

  If Wp is the two-element vector [w1 w2], where w1 < w2, then cheby1 designs a bandpass or bandstop filter with lower edge frequency w1 and higher edge frequency w2.

• For digital filters, the passband edge frequencies must lie between 0 and 1, where 1 corresponds to the Nyquist rate—half the sample rate or π rad/sample.

  For analog filters, the passband edge frequencies must be expressed in radians per second and can take on any positive value.

Data Types: double

fType — Filter type

"low" | "bandpass" | "high" | "stop"

Filter type, specified as one of the following:

• "low" specifies a lowpass filter with passband edge frequency Wp. "low" is the default for scalar Wp.

• "high" specifies a highpass filter with passband edge frequency Wp.

• "bandpass" specifies a bandpass filter of order 2n if Wp is a two-element vector. "bandpass" is the default when Wp has two elements.

• "stop" specifies a bandstop filter of order 2n if Wp is a two-element vector.

Output Arguments

b,a — Transfer function coefficients

row vectors

Transfer function coefficients of the filter, returned as row vectors. Given the filter order n, the function returns b and a as vectors with r elements, where r = n+1 for lowpass and highpass filters and r = 2*n+1 for bandpass and bandstop filters.

The transfer function is expressed in terms of b = [b₁ b₂ ⋯ b_r] and a = [a₁ a₂ ⋯ a_r] as one of these:

• $$H(z)=\frac{b_1+b_2{z}^{-1}+\cdots +b_r{z}^{-(r-1)}}{a_1+a_2{z}^{-1}+\cdots +a_r{z}^{-(r-1)}}$$ for digital filters.

• $$H(s)=\frac{b_1{s}^{r-1}+b_2{s}^{r-2}+\cdots +b_r}{a_1{s}^{r-1}+a_2{s}^{r-2}+\cdots +a_r}$$ for analog filters.

Data Types: double

z,p,k — Zeros, poles, and gain

column vectors, scalar

Zeros, poles, and gain of the filter, returned as two column vectors and a scalar. Given the filter order n, the function returns z and p as vectors with r elements, where r = n for lowpass and highpass filters and r = 2*n for bandpass and bandstop filters.

The transfer function is expressed in terms of z = [z₁ z₂ ⋯ z_r], p = [p₁ p₂ ⋯ p_r], and k as one of these:

• $$H(z)=k\frac{(1-z_1{z}^{-1})(1-z_2{z}^{-1})\cdots (1-z_r{z}^{-1})}{(1-p_1{z}^{-1})(1-p_2{z}^{-1})\cdots (1-p_r{z}^{-1})}$$ for digital filters.

• $$H(s)=k\frac{(s-z_1)(s-z_2)\cdots (s-z_r)}{(s-p_1)(s-p_2)\cdots (s-p_r)}$$ for analog filters.

Data Types: double

A,B,C,D — State-space representation

matrices

State-space representation of the filter, returned as matrices. If r = n for lowpass and highpass designs and r = 2n for bandpass and bandstop filters, then A is r × r, B is r × 1, C is 1 × r, and D is 1 × 1.

The state-space matrices relate the state vector x, the input u, and the output y through one of these equation systems.

• For digital filters:

  $$\begin{array}{c}x(k+1)=\text{A}x(k)+\text{B}u(k)\\ y(k)=\text{C}x(k)+\text{D}u(k).\end{array}$$

• For analog filters:

  $$\begin{array}{l}\dot{x}=\text{A}x+\text{B}u\\ y=\text{C}x+\text{D}u.\end{array}$$

Data Types: double

B,A — Cascaded transfer function (CTF) coefficients

row vectors | matrices

Cascaded transfer function (CTF) coefficients, returned as a row vectors or matrices. B and A list the numerator and denominator coefficients of the cascaded transfer function, respectively.

The sizes for B and A are L-by-(m+1) and L-by-(n+1), respectively. The function returns the first column of A as 1, thus A(1) = 1 when A is a row vector.

• L represents the number of CTF sections.

• m represents the order of the CTF numerators.

• n represents the order of the CTF denominators.

The cheby1 function returns the CTF coefficients with these order specifications:

• m = n = 2 for lowpass and highpass filters.

• m = n = 4 for bandpass and bandstop filters.

Note

To customize the CTF coefficient computation, such as setting a different order in the CTF coefficients or customizing the gain scaling, specify to return z,p,k and then use zp2ctf to obtain B,A.

For more information about the cascaded transfer function format and coefficient matrices, see Return Digital Filters in CTF Format.

gS — Overall system gain

real-valued scalar

Since R2024b

Overall system gain, returned as a real-valued scalar.

• If you specify to return gS, the cheby1 function normalizes the numerator coefficients so that the first column of B is 1 and returns the overall system gain in gS.

• If you do not specify to return gS, the cheby1 function uniformly distributes the system gain across all system sections using the scaleFilterSections function.

More About

Transfer Functions and CTF

Numerical Instability of Transfer Function Syntax

In general, use cascaded transfer functions ("ctf" syntaxes) to design IIR digital filters. If you design the filter using transfer functions (any of the [b,a] syntaxes), you might encounter numerical instabilities. These instabilities are due to round-off errors and can occur for an order n as low as 4. This example illustrates this limitation.

n = 6;
Rp = 0.6;
Fs = 200e6;
Wn = [0.5e6 6e6]/(Fs/2);
ftype = "bandpass";

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

% CTF design
[B,A] = cheby1(n,Rp,Wn,ftype,"ctf");

% Compare frequency responses
[hTF,f] = freqz(b,a,8192,Fs);
hCTF = freqz(B,A,8192,Fs);

semilogx(f/1e6,db(hTF),".-",f/1e6,db(hCTF))
grid on
legend(["TF Design" "CTF Design"])
xlabel("Frequency (MHz)")
ylabel("Magnitude (dB)")

Algorithms

Chebyshev Type I filters are equiripple in the passband and monotonic in the stopband. Type I filters roll off faster than Type II filters, but at the expense of greater deviation from unity in the passband.

The cheby1 function uses this algorithm:

1. It finds the lowpass analog prototype poles, zeros, and gain using the function cheb1ap.

2. It converts the poles, zeros, and gain into state-space form.

3. When designing digital filters, it prewarps the specified digital frequency Wp to the corresponding analog frequencies.

4. If required, it uses a state-space transformation to convert the lowpass filter to a highpass, bandpass, or bandstop filter with the desired frequency constraints.

5. When designing digital filters, it uses bilinear to convert the analog filter into a digital filter through a bilinear transformation.

6. It converts the state-space filter back to the transfer-function or zero-pole-gain form, as required.

Extended Capabilities

C/C++ Code Generation

Generate C and C++ code using MATLAB® Coder™.

Usage notes and limitations:

All inputs must be constants. Expressions or variables are allowed if their values do not change.

GPU Code Generation

Generate CUDA® code for NVIDIA® GPUs using GPU Coder™.

Refer to the usage notes and limitations in the C/C++ Code Generation section. The same usage notes and limitations apply to GPU code generation.

Thread-Based Environment

Run code in the background using MATLAB® backgroundPool or accelerate code with Parallel Computing Toolbox™ ThreadPool.

The cheby1 function fully supports thread-based environments. For more information, see Run MATLAB Functions in Thread-Based Environment.

Version History

Introduced before R2006a

R2026b: Use thread-based environments

You can use cheby1 in thread-based environments. For more information, see Run MATLAB Functions in Thread-Based Environment.

R2024b: Design digital filters using cascaded transfer functions

The cheby1 function supports outputs in the cascaded transfer function (CTF) format.