コグノスケ


link 未来から過去へ表示  link 過去から未来へ表示(*)

link もっと前
2014年11月11日 >>> 2014年11月20日
link もっと後

2014年11月15日

USB-DACと矩形波

目次: PC

前回はSin波を見ました(2014年10月18日の日記参照)が、今回は矩形波を見てみようと思います。

サンプリング周波数の1/2の周波数(以降、1/2 fsと書きます)の矩形波のLPCMデータを下記のように表すとします。


周波数1/2 fsの矩形波のLPCMデータ

矩形波と書きましたが、このデータをアナログに変換しても矩形波にはなりません。下記のような22kHzのSin波になります。


周波数1/2 fsの矩形波をアナログに変換

私もそれほど詳しくはありませんが、標本化定理と言って、ある波形をサンプリング周波数44kHzでサンプリングした際に完全に元の波形に戻せるのは1/2 fs未満、つまり22kHzまでの周波数を持った波形です。それ以上は情報が失われ、元に戻せません。

見た目は周波数22kHzの波形に見える矩形波でも、波形の立上がり、立下りの部分に高い周波数を含んでいて、この部分が失われてしまうため、元に戻らないのです。

同じ理屈で1/4 fsの矩形波、つまり44kHzのサンプリング周波数で、11kHzの矩形波を表そうとしても、立上がり、立下りの部分が失われ、11kHzのSin波になります。LPCMデータの最大値は変わりませんが、22kHzの時より振幅が大きくなるのが面白いですね。


周波数1/4 fsの矩形波をアナログに変換

見た目が矩形波らしくなるのは1/6 fsか1/8 fsくらいでしょうか。つまり44kHzのサンプリング周波数で、7.4kHzか5.5kHzの矩形波を表そうとしたときです。


周波数1/6 fsの矩形波をアナログに変換


周波数1/8 fsの矩形波をアナログに変換

実際やってみましょう

では、下記のLPCMデータを使ってUSB-DACで再生してみます。

各種揃えましたが、訳あって、この中からサンプリング周波数44kHzのデータを使うことにします。まずは1/2 fsから。


ONKYO SE-U33GXV2の22kHz矩形波出力(fs = 44.1kHz)

綺麗なSin波です。次は1/4 fsです。


ONKYO SE-U33GXV2の11kHz矩形波出力(fs = 44.1kHz)

これも綺麗なSin波ですね。振幅が先ほどより大きいのがわかると思います。次は1/6 fsと1/8 fsです。


ONKYO SE-U33GXV2の7.4kHz矩形波出力(fs = 44.1kHz)


ONKYO SE-U33GXV2の5.5kHz矩形波出力(fs = 44.1kHz)

角の部分はまだガタガタしていますが、矩形波らしくなってきたのがわかります。

編集者:すずき(2024/10/25 02:03)

コメント一覧

  • すずきさん(2014/11/19 02:30)
    DA 変換を調べていたら、どうも、図のような補完でリンギングが起きているわけではなさそうです。信号処理は難しいな…。
open/close この記事にコメントする



2014年11月20日

DA変換

目次: PC

前回の日記(2014年11月15日の日記参照)では、DA変換を理解できていなくて、いい加減な図を描いてしまったので、今日はLPCMのサンプル値から連続時間信号を生成した場合の波形を計算して描いてみようと思います。

調べたところによれば、サンプル値x[n] を離散時間フーリエ変換し、周波数領域に変換して、得られた周波数領域のうち角周波数 -π/Ts < ω < π/Tsのみを残すローパスフィルタを掛け(※)、残った周波数領域を逆フーリエ変換すると連続時間信号になるそうです。

(※)Tsはサンプリング周期、つまり何秒に1回サンプリングするかを表す記号です。

計算をすっ飛ばして結果だけ言うと、時間をTsずつずらしながら各サンプル値x[n]とsinc関数(理想的なローパスフィルタを逆フーリエ変換するとsinc関数になる)の積を取り、全て足し合わせれば連続時間信号になります。何故なのか説明できるほど理解できていませんが、ひとまずそういうモンとして計算します。


周波数1/2 fsの矩形波のLPCMサンプルデータと連続時間信号

サンプル値は-400から400で、サンプリング周期Tsは20 [s]、サンプリング周波数fsは1/20 = 0.05 [Hz] とし、0〜960 [s] までを計算しています。横軸tの単位は [s] です。

これだけだと「ふーん……?」で終わるので、わかりやすいように1個目のサンプルと、2個目のサンプルのみで、連続時間信号を計算してみます。


各LPCMサンプルデータと対応するsinc関数

赤色(系列2)が1個目のサンプル、緑色(系列3)が2個目のサンプル、青色(系列1)が総和です。1/2 fsかつ、ピーク値がLPCMのサンプル値とほぼイコールのSin波が計算されるのがわかると思います。

じゃあ1/4 fsは?こうなります。


周波数1/4 fsの矩形波のLPCMサンプルデータと連続時間信号

前回の適当な図とは違い1/4 fsのSin波が出力され、ピーク値はLPCMのサンプル値を超えてしまう様子もよくわかると思います。ちなみに連続時間信号のピーク値はt=409(真ん中の山の頂点)で563.7です。

もちろん1/6 fsや1/8 fsも計算可能ですが、長くなってきたのでまた後日にでも載せます。

編集者:すずき(2024/10/25 02:02)

コメント一覧

  • コメントはありません。
open/close この記事にコメントする



link もっと前
2014年11月11日 >>> 2014年11月20日
link もっと後

管理用メニュー

link 記事を新規作成

<2014>
<<<11>>>
------1
2345678
9101112131415
16171819202122
23242526272829
30------

最近のコメント5件

  • link 24年10月1日
    すずきさん (10/06 03:41)
    「xrdpで十分動作しているので、Wayl...」
  • link 24年10月1日
    hdkさん (10/03 19:05)
    「GNOMEをお使いでしたら今はWayla...」
  • link 24年10月1日
    すずきさん (10/03 10:12)
    「私は逆にVNCサーバーに繋ぐ使い方をした...」
  • link 24年10月1日
    hdkさん (10/03 08:30)
    「おー、面白いですね。xrdpはすでに立ち...」
  • link 14年6月13日
    2048player...さん (09/26 01:04)
    「最後に、この式を出すのに紙4枚(A4)も...」

最近の記事3件

  • link 24年10月28日
    すずき (10/30 23:49)
    「[Linuxからリモートデスクトップ] 目次: Linux開発用のLinuxマシンの画面を見るにはいろいろな手段がありますが、...」
  • link 23年4月10日
    すずき (10/30 23:46)
    「[Linux - まとめリンク] 目次: Linux関係の深いまとめリンク。目次: RISC-V目次: ROCK64/ROCK...」
  • link 24年10月24日
    すずき (10/25 02:35)
    「[ONKYOからM-AUDIOのUSB DACへ] 目次: PCかれこれ10年以上(2013年3月16日の日記参照)活躍してく...」
link もっとみる

こんてんつ

open/close wiki
open/close Linux JM
open/close Java API

過去の日記

open/close 2002年
open/close 2003年
open/close 2004年
open/close 2005年
open/close 2006年
open/close 2007年
open/close 2008年
open/close 2009年
open/close 2010年
open/close 2011年
open/close 2012年
open/close 2013年
open/close 2014年
open/close 2015年
open/close 2016年
open/close 2017年
open/close 2018年
open/close 2019年
open/close 2020年
open/close 2021年
open/close 2022年
open/close 2023年
open/close 2024年
open/close 過去日記について

その他の情報

open/close アクセス統計
open/close サーバ一覧
open/close サイトの情報

合計:  counter total
本日:  counter today

link About www2.katsuster.net
RDFファイル RSS 1.0

最終更新: 10/30 23:49