Субдискретизация — искажение
В этом примере показано, как постараться не искажать при субдискретизации сигнала. Если дискретное время основополосная спектральная поддержка сигнала не ограничивается интервалом
радианов ширины, прореживающих
результатами в искажении. Искажение является искажением, которое происходит, когда перекрывающиеся копии спектра сигнала добавляются вместе. Чем больше основополосная спектральная поддержка сигнала превышает радианы, тем более серьезный искажение. Продемонстрируйте искажение в сигнале, прореженном два. Основополосная спектральная поддержка сигнала превышает
радианы по ширине.
Создайте сигнал с основополосной спектральной поддержкой, равной
радианам. Используйте fir2 спроектировать сигнал. Постройте спектр сигнала.
F = [0 0.2500 0.5000 0.7500 1.0000]; A = [1.00 0.6667 0.3333 0 0]; Order = 511; B1 = fir2(Order,F,A); [Hx,W] = freqz(B1,1,8192,'whole'); Hx = [Hx(4098:end) ; Hx(1:4097)]; omega = -pi+(2*pi/8192):(2*pi)/8192:pi; plot(omega,abs(Hx)) xlim([-pi pi]) grid title('Magnitude Spectrum') xlabel('Radians/Sample') ylabel('Magnitude')
Вы видите, что основополосная спектральная поддержка сигнала превышает
.
Downsample сигнал фактором 2 и график спектр прореженного сигнала со спектром исходного сигнала.
y = downsample(B1,2,0); [Hy,W] = freqz(y,1,8192,'whole'); Hy = [Hy(4098:end) ; Hy(1:4097)]; hold on plot(omega,abs(Hy),'r','linewidth',2) legend('Original Signal','Downsampled Signal') text(-2.5,0.35,'\downarrow aliasing','HorizontalAlignment','center') text(2.5,0.35,'aliasing \downarrow','HorizontalAlignment','center') hold off
В дополнение к амплитудному масштабированию спектра суперпозиция наложения спектральных копий вызывает искажение исходного спектра для
.
Увеличьте основополосную спектральную поддержку сигнала к
и проредите сигнал 2. Постройте исходный спектр наряду со спектром прореженного сигнала.
F = [0 0.2500 0.5000 0.7500 7/8 1.0000]; A = [1.00 0.7143 0.4286 0.1429 0 0]; Order = 511; B2 = fir2(Order,F,A); [Hx,W] = freqz(B2,1,8192,'whole'); Hx = [Hx(4098:end) ; Hx(1:4097)]; omega = -pi+(2*pi/8192):(2*pi)/8192:pi; plot(omega,abs(Hx)) xlim([-pi pi]) y = downsample(B2,2,0); [Hy,W] = freqz(y,1,8192,'whole'); Hy = [Hy(4098:end) ; Hy(1:4097)]; hold on plot(omega,abs(Hy),'r','linewidth',2) grid legend('Original Signal','Downsampled Signal') xlabel('Radians/Sample') ylabel('Magnitude') hold off
Увеличенная спектральная ширина приводит к более явному искажению в спектре прореженного сигнала, потому что больше энергии сигнала снаружи.
Наконец, создайте сигнал с основополосной спектральной поддержкой, ограниченной. Downsample сигнал фактором 2 и график спектр исходных и прореженных сигналов. Прореженный сигнал является полной полосой, но форма спектра сохраняется, потому что спектральные копии не накладываются. Нет никакого искажения.
F = [0 0.250 0.500 0.7500 1]; A = [1.0000 0.5000 0 0 0]; Order = 511; B3 = fir2(Order,F,A); [Hx,W] = freqz(B3,1,8192,'whole'); Hx = [Hx(4098:end) ; Hx(1:4097)]; omega = -pi+(2*pi/8192):(2*pi)/8192:pi; plot(omega,abs(Hx)) xlim([-pi pi]) y = downsample(B3,2,0); [Hy,W] = freqz(y,1,8192,'whole'); Hy = [Hy(4098:end) ; Hy(1:4097)]; plot(omega,abs(Hx)) hold on plot(omega,abs(Hy),'r','linewidth',2) grid legend('Original Signal','Downsampled Signal') xlabel('Radians/Sample') ylabel('Magnitude') hold off
Вы видите в предыдущей фигуре, что форма спектра сохраняется. Спектр прореженного сигнала является расширенной и масштабированной версией спектра исходного сигнала, но нет никакого искажения.
Смотрите также
downsample | fir2 | freqz