Main Content

最新のリリースでは、このページがまだ翻訳されていません。 このページの最新版は英語でご覧になれます。

ダウンサンプリング — エイリアシング

この例では、信号のダウンサンプリングに際してエイリアシングを回避する方法を示します。離散時間信号の基底域スペクトルのサポートが幅 2π/M ラジアンの区間に限定されない場合、M によるダウンサンプリングはエイリアシングをもたらします。エイリアシングは、オーバーラップしている信号のスペクトルの複製が加算されて一緒になることで生じる歪みです。信号の基底域スペクトルのサポートが 2π/M ラジアンを超えるに従い、エイリアシングはより著しくなります。2 によるダウンサンプリングを行った信号のエイリアシングを示します。信号の基底域スペクトルのサポートは、その幅において π ラジアンを超過します。

基底域スペクトルのサポートが幅 3π/2 ラジアンの信号を作成します。fir2 を使用して信号を設計します。信号のスペクトルをプロットします。信号の基底域スペクトルのサポートは、[-π/2,π/2] を超過します。

f = [0 0.2500 0.5000 0.7500 1.0000];
a = [1.00 0.6667 0.3333 0 0];

nf = 512;
b1 = fir2(nf-1,f,a);
Hx = fftshift(freqz(b1,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 の係数で信号をダウンサンプリングして、ダウンサンプリングした信号のスペクトルを元の信号のスペクトルとともにプロットします。スペクトルの振幅のスケーリングに加えて、オーバーラップしているスペクトルの複製の重ね合わせによって、元のスペクトルの歪みが |ω|>π/2 において生じています。

y = downsample(b1,2,0);
Hy = fftshift(freqz(y,1,nf,'whole'));

hold on
plot(omega/pi,abs(Hy))
hold off
legend('Original','Downsampled')
text(2.5/pi*[-1 1],0.35*[1 1],{'\downarrow Aliasing','Aliasing \downarrow'}, ...
    'HorizontalAlignment','center')

信号の基底域スペクトルのサポートを [-7π/8,7π/8] に増やして、信号に対して 2 によるダウンサンプリングを行います。元のスペクトルをダウンサンプリングした信号のスペクトルと共にプロットします。スペクトルの幅を増やした結果、ダウンサンプリングした信号のスペクトルにおいてエイリアシングがより顕著になります。より多くの信号エネルギーが [-π/2,π/2] の外側にあるためです。

f = [0 0.2500 0.5000 0.7500 7/8 1.0000];
a = [1.00 0.7143 0.4286 0.1429 0 0];

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

plot(omega/pi,abs(Hx))
grid
xlabel('\times\pi rad/sample')
ylabel('Magnitude')

y = downsample(b2,2,0);
Hy = fftshift(freqz(y,1,nf,'whole'));

hold on
plot(omega/pi,abs(Hy))
hold off
legend('Original','Downsampled')

最後に、基底域スペクトルのサポートを [-π/2,π/2] に制限した信号を作成します。2 の係数で信号をダウンサンプリングして、元の信号とダウンサンプリングした信号のスペクトルをプロットします。ダウンサンプリングした信号はフル バンドです。ダウンサンプリングした信号のスペクトルは、元のスペクトルを引き伸ばして縮尺を変更したバージョンですが、形状は維持されます。スペクトルの複製がオーバーラップしていないためです。エイリアシングはありません。

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

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

plot(omega/pi,abs(Hx))
grid
xlabel('\times\pi rad/sample')
ylabel('Magnitude')

y = downsample(b3,2,0);
Hy = fftshift(freqz(y,1,nf,'whole'));

hold on
plot(omega/pi,abs(Hy))
hold off
legend('Original','Downsampled')

参考

| |