Main Content

polyscale

Scale roots of polynomial

Syntax

b = polyscale(a,alpha)

Description

b = polyscale(a,alpha) scales the roots of a polynomial in the z-plane, where a is a vector containing the polynomial coefficients and alpha is the scaling factor.

If alpha is a real value in the range [0 1], then the roots of a are radially scaled toward the origin in the z-plane. Complex values for alpha allow arbitrary changes to the root locations.

Examples

collapse all

Express the solutions to the equation x7=1 as the roots of a polynomial. Plot the roots in the complex plane.

pp = [1 0 0 0 0 0 0 -1];
zplane(pp,1)

Scale the roots of p in and out of the unit circle. Plot the results.

hold on

for sc = [1:-0.2:0.2 1.2 1.4];
    b = polyscale(pp,sc);
    plot(roots(b),'o')
end

axis([-1 1 -1 1]*1.5)

hold off

Load a speech signal sampled at Fs=7418Hz. The file contains a recording of a female voice saying the word "MATLAB®."

load mtlb

Model a 100-sample section of the signal using a 12th-order autoregressive polynomial. Perform bandwidth expansion of the signal by scaling the roots of the autoregressive polynomial by 0.85.

Ao = lpc(mtlb(1000:1100),12);
Ax = polyscale(Ao,0.85);

Plot the zeros, poles, and frequency responses of the models.

subplot(2,2,1)
zplane(1,Ao)
title('Original')

subplot(2,2,3)
zplane(1,Ax)
title('Flattened')

subplot(1,2,2)
[ho,w] = freqz(1,Ao);
[hx,w] = freqz(1,Ax);
plot(w/pi,abs([ho hx]))
legend('Original','Flattened')

Tips

By reducing the radius of the roots in an autoregressive polynomial, the bandwidth of the spectral peaks in the frequency response is expanded (flattened). This operation is often referred to as bandwidth expansion.

Version History

Introduced before R2006a

See Also

|