Масштабируйте корни полинома
b = polyscale(a,alpha)
b = polyscale(a,alpha) масштабирует корни полинома в z-плоскости, где a является вектором, содержащим полиномиальные коэффициенты, и alpha является масштабным коэффициентом.
Если alpha является действительным значением в области значений [0 1], то корни a радиально масштабируются к источнику в z-плоскости. Комплексные числа для alpha позволяют произвольные изменения в корневых местоположениях.
Выразите решения уравнения как корни полинома. Постройте корни в комплексной плоскости.
pp = [1 0 0 0 0 0 0 -1]; zplane(pp,1)
Масштабируйте корни p в и из модульного круга. Постройте график результатов.
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
Загрузите речевой сигнал, выбранный в . Файл содержит запись розеточной речи, говоря слово "MATLAB®".
load mtlbСмоделируйте раздел с 100 выборками сигнала с помощью 12-го порядка авторегрессивный полином.
Ao = lpc(mtlb(1000:1100),12); Ax = polyscale(Ao,0.85);
Выполните расширение пропускной способности сигнала путем масштабирования корней авторегрессивного полинома 0,85. Постройте нули, полюса и частотные характеристики моделей.
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')
Путем сокращения радиуса корней в авторегрессивном полиноме пропускная способность спектрального peaks в частотной характеристике расширена (сглаженная). Эта операция часто упоминается как расширение пропускной способности.