Exponenta Banner

Повышающая дискретизация — обработка изображений артефактов

Скрипт Open Live Script

В этом примере показано, как сверхдискретизировать сигнал и как повышающая дискретизация может привести к изображениям. Повышающая дискретизация сигнала сокращает спектр. Например, сверхдискретизировав сигнал 2 результатами в сокращении спектра фактором 2. Поскольку спектр сигнала дискретного времени 2π- периодический, сокращение может заставить копии спектра обычно за пределами основной полосы появляться в интервале [-π,π].

Создайте сигнал дискретного времени, основополосная спектральная поддержка которого [-π,π]. Постройте спектр величины.

F = [0 0.250 0.500 0.7500 1];
A = [1.0000 0.5000 0 0 0];
Order = 511;
B = fir2(Order,F,A);
[Hx,W] = freqz(B,1,8192,'whole');
Hx = [Hx(4098:end) ; Hx(1:4097)];
omega = -pi+(2*pi/8192):(2*pi)/8192:pi;

plot(omega,abs(Hx))

Сверхдискретизируйте сигнал 2. Постройте спектр сверхдискретизированного сигнала.

y = upsample(B,2);
[Hy,W] = freqz(y,1,8192,'whole');
Hy = [Hy(4098:end) ; Hy(1:4097)];

hold on
plot(omega,abs(Hy),'r','linewidth',2)
xlim([-pi pi])
legend('Original Signal','Upsampled Signal')
xlabel('Radians/Sample')
ylabel('Magnitude')
text(-2,0.5,'\leftarrow Imaging','HorizontalAlignment','center')
text(2,0.5,'Imaging \rightarrow','HorizontalAlignment','center')
hold off

Вы видите в предыдущей фигуре, что сокращение спектра вовлекло последующие периоды спектра на интервал [-π,π].

Смотрите также

| |

Документация Signal Processing Toolbox
Поддержка