Wolfram言語

計算音声

AMラジオ

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

In[1]:=
Click for copyable input
a = BandpassFilter[ AudioNormalize@ AudioResample[Import["ExampleData/rule30.wav"], 96000], {Quantity[ 10, "Hertz"], Quantity[5000, "Hertz"]}, 201]
Out[1]=
In[2]:=
Click for copyable input
Periodogram[a, 2000, ImageSize -> Small, PlotRange -> {{0, 10000}, All}]
Out[2]=

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

In[3]:=
Click for copyable input
amsignal = ((.5 + a) AudioGenerator[{"Sin", Quantity[22050, "Hertz"]}, Duration@a])/2
Out[3]=
In[4]:=
Click for copyable input
Periodogram[amsignal, 2200, ImageSize -> Small]
Out[4]=

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

In[5]:=
Click for copyable input
result = BandpassFilter[ amsignal*AudioGenerator[{"Sin", Quantity[22050, "Hertz"]}, Duration@amsignal], {Quantity[100, "Hertz"], Quantity[5000, "Hertz"]}, 201]
Out[5]=

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

In[6]:=
Click for copyable input
BandpassFilter[# - Mean@# &@ MaxFilter[amsignal, 4], {Quantity[100, "Hertz"], Quantity[5000, "Hertz"]}, 201]
Out[6]=

関連する例

de en es fr ko pt-br ru zh