Top/EBUG勉強会/20120218_Audacityで信号処理
EBUG 第41回会合
2012年2月18日 @(株)銀座堂
川俣吉広、kaw@on.rim.or.jp

発端 - 最初はまともな動機だったのに.....

なんかもう意地(雀を撃つのに大砲を持ち出す)

で、実際どんなもんなのか? ... 船舶気象通報で耳試し

船舶気象通報

・海上保安庁が運用

・周波数、空中線電力: 1669kHz, 50W

・変調形式: H3E (USB) - 上側波帯全搬送波振幅変調

・各局数分ずつ放送 - 一時間で一巡

XmitSite.jpg
Location.jpg
(送信所の写真はこちらから借用させて貰いました)

というわけで、受信能力を試すには結構よさげな対象。

ラジオを聴くのにAudacityを持ち出す

セッティング

Setups.jpg

Audacity - フリーでクロスプラットフォームな音声エディタ

通常は音楽などの編集・加工に用いられることが多いと思われるツールですが、今回はこのAudacityを用いて信号の変換や測定/解析などを行ってみることにします。

#img(): File not found:
Diagram.gif
系統図

調整

大原則 - レベルオーバーさせない。

調整方法

振幅変調についておさらい

振幅変調とは、正弦波との乗算(周波数変換も)。

AMsignal.gif
時間領域

三角関数の加法定理から

sin((ωc+ωb)*t) = sin(ωc*t)*cos(ωb*t)+cos(ωc*t)*sin(ωb*t)
sin((ωc-ωb)*t) = sin(ωc*t)*cos(ωb*t)-cos(ωc*t)*sin(ωb*t)
sin(ωc*t)*cos(ωb*t) = (sin((ωc+ωb)*t)+sin((ωc-ωb)*t))/2
sin(ωc*t)*(C+cos(ωb*t)) = (C*sin(ωc*t)+sin((ωc+ωb)*t)+sin((ωc-ωb)*t))/2

被変調波に直流成分が含まれていると、変調波では搬送波が現れる。

AMsignal_Spect.gif
周波数領域

要するに正弦波を掛け算するとその周波数だけ、元の信号のスペクトラムが平行移動する。
あと、直交成分が加わった場合の位相表現を複素平面上で云々とかあるけど、とりあえず今回そういうのは置いておく方向で...

試行錯誤の道筋

記録してみよう

受信機側の作業

PC側の作業

考えられる幾つかの方法;

  1. Audacity立ち上げて録音ボタンをポチッ
    1. リアルタイムでグラフィカルに波形とか表示しているので、結構重い。
    2. 非圧縮で記録してるのでディスク容量を食う。
    3. 一定の間隔でファイルを生成する、とかが面倒。
    4. タイムコードとかも付いてないので、あとで実時間との対応をとるのも面倒。

  2. dd if=/dev/audio of=outputfile.raw とか
    1. やっぱり、非圧縮で記録してるのでディスク容量を結構食う。

  3. lameを用いて、MP3にエンコードしながら記録
    1. 不可逆圧縮なので、測定精度に影響がでるかも....

lameを用いる方法でやってみた。2つのシェルコードを走らせる。

while :; do
    lame -r -s 48 -m m -a --preset 32 /dev/audio `date +monitor-%Y%m%d%a-%H%M-1996kHz.mp3`
done
while :; do
    sleep 3600
    pkill lame
done

そして、例えば次の正時が13:00だとすると、

echo 'pkill -f "sleep 3600"' | at 1300

とする。これで毎正時に出力するファイルが更新される。
lameを上記のオプションで実行すると;

Autoconverting from stereo to mono. Setting encoding to mono mode.
Resampling:  input 48 kHz  output 22.05 kHz
Using polyphase lowpass filter, transition band:  8269 Hz -  8535 Hz
Encoding /dev/audio to /dev/null
Encoding as 22.05 kHz single-ch MPEG-2 Layer III (11x) average 32 kbps qval=3
    Frame          |  CPU time/estim | REAL time/estim | play/CPU |    ETA 
   400/400   (100%)|    0:00/    0:00|    0:10/    0:10|   17.415x|    0:00 
  8 [  0] 
 16 [  0] 
 24 [ 75] *****************
 32 [320] *********************************************************************
 40 [  5] **
 48 [  0] 
 56 [  0] 
 64 [  0] 
 80 [  0] 
 96 [  0] 
112 [  0] 
128 [  0] 
144 [  0] 
160 [  0] 
-------------------------------------------------------------------------------
   kbps       mono %     long  %
   30.6      100.0       100.0

解析してみよう  

まずは、受信信号の強度を測定してみる。

受信信号の強度を測定するには、ノイズや側帯波を除外する処理を行う ... フィルタリング

加工してみよう  

Modulation.gif
  1. Toremolo - 4000Hz
  2. HPF - 4600Hz
  3. Toremolo - 4600Hz
  4. LPF - 3500Hz

まとめ  


Top/EBUG勉強会/20120218_Audacityで信号処理
recent(50)
2024-03-19 2024-02-17 2024-01-31 2023-12-24 2023-12-20 2023-11-05 2023-10-04 2023-08-28 2023-06-21 2023-06-19 2023-06-03 2023-05-30 2023-05-20 2023-05-01 2023-04-30 2023-04-23 2023-04-13 2023-04-12 2023-03-17 2023-02-24 2023-02-23 2023-02-22 2022-12-23 2022-12-13 2022-11-17 2022-11-15 2022-11-14 2022-11-13 2022-11-05 2022-09-28 2022-08-08 2022-04-11

Attach file: filecarrier-1hr-mag.png 1836 download [Information] filecarrier-daily.png 1898 download [Information] filecarrier-hourly.png 1664 download [Information] filecarrier-1hr.png 1770 download [Information] fileModulation.gif 1663 download [Information] fileAMsignal_Spect.gif 1866 download [Information] fileAMsignal.gif 1864 download [Information] fileLocation.jpg 1812 download [Information] fileSSBtune.3g2 1558 download [Information] fileDiagram.gif 1874 download [Information] fileSetups.jpg 1797 download [Information] fileXmitSite.jpg 1704 download [Information] fileAntenna.gif 1866 download [Information] fileBoom.jpg 1767 download [Information] filebalun.jpg 1774 download [Information] fileAntenna.jpg 1868 download [Information] fileft757.jpg 1818 download [Information] fileNoisyEnv.3g2 1529 download [Information] fileradio.jpg 1949 download [Information]

Front page   Edit Freeze Diff History Attach Copy Rename Reload   New Page list Search Recent changes   Help   RSS of recent changes
Last-modified: 2012-03-03 (Sat) 22:22:06