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

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

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

f = [0 0.250 0.500 0.7500 1];
a = [1.0000 0.5000 0 0 0];

nf = 512;
b = fir2(nf-1,f,a);
Hx = fftshift(freqz(b,1,nf,'whole'));

omega = -pi:2*pi/nf:pi-2*pi/nf;
plot(omega/pi,abs(Hx))
grid
xlabel('\times\pi rad/sample')
ylabel('Magnitude')

Figure contains an axes. The axes contains an object of type line.

Сверхдискретизируйте сигнал 2. Постройте спектр сверхдискретизированного сигнала. Сокращение спектра вовлекло последующие периоды спектра на интервал [-π,π].

y = upsample(b,2);
Hy = fftshift(freqz(y,1,nf,'whole'));

hold on
plot(omega/pi,abs(Hy))
hold off
legend('Original','Upsampled')
text(0.65*[-1 1],0.45*[1 1],["\leftarrow Imaging" "Imaging \rightarrow"], ...
    'HorizontalAlignment','center')

Figure contains an axes. The axes contains 4 objects of type line, text. These objects represent Original, Upsampled.

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

| |