Этот пример показывает, как улучшить сигнал и как увеличение дискретизации может привести к изображениям. Увеличение дискретизации сигнала сокращает спектр. Для примера увеличение дискретизации сигнала на 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')
Увеличьте значение сигнала на 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')