AMラジオ

帯域制限信号から開始する．

In[1]:=

a = BandpassFilter[ AudioNormalize@ AudioResample[Import["ExampleData/rule30.wav"], 96000], {Quantity[ 10, "Hertz"], Quantity[5000, "Hertz"]}, 201]

Out[1]=

In[2]:=

Periodogram[a, 2000, ImageSize -> Small, PlotRange -> {{0, 10000}, All}]

Out[2]=

22,050Hzの正弦曲線で振幅を変調し，信号のパワースペクトルを高い位置にシフトする．結果は，ほとんどの人には，すでに聞こえない．

In[3]:=

amsignal = ((.5 + a) AudioGenerator[{"Sin", Quantity[22050, "Hertz"]}, Duration@a])/2

Out[3]=

In[4]:=

Periodogram[amsignal, 2200, ImageSize -> Small]

Out[4]=

結果に，位相が等しい別の22,050Hzの正弦曲線を掛けて，AM信号を復調する．

In[5]:=

result = BandpassFilter[ amsignal*AudioGenerator[{"Sin", Quantity[22050, "Hertz"]}, Duration@amsignal], {Quantity[100, "Hertz"], Quantity[5000, "Hertz"]}, 201]

Out[5]=

非線形フィルタを使って，AM信号を復調する．

In[6]:=

BandpassFilter[# - Mean@# &@ MaxFilter[amsignal, 4], {Quantity[100, "Hertz"], Quantity[5000, "Hertz"]}, 201]

Out[6]=