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])/2Out[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]=
