コグノスケ

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

過去から未来へ表示

もっと前

2014年11月20日 >>> 2014年11月7日

もっと後

2014年11月20日

permalink

編集する

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)

コメント一覧

コメントはありません。

この記事にコメントする

2014年11月15日

permalink

編集する

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 変換を調べていたら、どうも、図のような補完でリンギングが起きているわけではなさそうです。信号処理は難しいな…。

この記事にコメントする

もっと前

2014年11月20日 >>> 2014年11月7日

もっと後

管理用メニュー

記事を新規作成

RSSを更新

<	2014					>
<<	<	11			>	>>
日	月	火	水	木	金	土
-	-	-	-	-	-	1
2	3	4	5	6	7	8
9	10	11	12	13	14	15
16	17	18	19	20	21	22
23	24	25	26	27	28	29
30	-	-	-	-	-	-

本日の日記へ

最近のコメント5件

26年7月14日
hdkさん (07/22 08:37)
「冷却水漏れなんて、修理するまでは補充して...」
26年1月23日
すずきさん (01/29 09:48)
「おおー、そんな昔からなんですね。歴史感じ...」
26年1月23日
hdkさん (01/27 19:53)
「#! はUNIX v8からだったってWi...」
24年12月9日
すずきさん (01/18 15:45)
「Thank you for your i...」
24年12月9日
Up2Uさん (01/15 12:57)
「Hi I also find the p...」

もっとみる

最近の記事20件

23年5月15日
すずき (07/22 03:04)
「[車 - まとめリンク] 目次: 車三菱 FTO GPX '95の話。群馬県へのドライブ1群馬県へのドライブ2将来車を買い替え...」
26年7月22日
すずき (07/22 03:04)
「[ジャガーさんの冷却水漏れの修理見積もり] 目次: 車先週、冷却水がほぼなくなってしまい修理工場にレッカーしたジャガーさんです...」
26年7月14日
すずき (07/22 02:37)
「[ジャガーさんの冷却水漏れ] 目次: 車朝、出勤しようと思ってエンジンをかけたところ冷却水不足と表示されました。わぁ、ついに出...」
26年7月19日
すずき (07/20 04:04)
「[ゲームを買ったら遊びましょう4] 目次: ゲーム前回の振り返り（2025年12月22日の日記参照）から半年経ちました。Swi...」
21年12月28日
すずき (07/20 04:03)
「[ゲーム - まとめリンク] 目次: ゲームNintendo DSを買ったパネルでポンDS最近の朝はパネポンDS聖剣伝説DSチ...」
26年7月13日
すずき (07/19 23:33)
「[Flightradar24のcontributorになる方法] 目次: 自宅サーバーADS-Bを受信できる環境があればFli...」
23年6月1日
すずき (07/19 23:31)
「[自宅サーバー - まとめリンク] 目次: 自宅サーバーこの日記システム、Wikiの話。カウンターをPerlからPHPに移植日...」
23年9月11日
すずき (07/15 12:18)
「[Windows - まとめリンク] 目次: WindowsWindows XPのブリッジ機能colinuxとWindowsの...」
25年1月31日
すずき (07/15 12:18)
「[Windowsの日本語フォントはみな長生き] 目次: WindowsGNOME48からフォントが変わるニュースを見ていて、フ...」
26年7月11日
すずき (07/12 01:26)
「[Timberborn、トロフィーコンプ] 目次: ゲームTimberbornのトロフィーをコンプリートしました。アーリーアク...」
24年2月7日
すずき (07/11 13:05)
「[複数の音声ファイルのラウドネスを統一したい] 目次: PythonPCやデジタル音楽プレーヤーで音楽を聞いていると、曲によっ...」
26年6月27日
すずき (07/08 01:45)
「[ジャガーさんのバンパー修理] 目次: 車以前（2025年11月21日の日記参照）、ジャガーさんの左前をぶつけてバンパーをガリ...」
26年7月3日
すずき (07/08 00:42)
「[ADS-Bを受信して飛行機の位置を見よう] 目次: 自宅サーバー空を飛ぶ飛行機は離陸直後や着陸寸前でもない限り飛んでいる姿は...」
26年7月1日
すずき (07/04 02:37)
「[GPSは世界一正確な時計、その3 - ntpsecとgpsd] 目次: 自宅サーバー以前（2015年5月8日の日記参照）、G...」
26年6月26日
すずき (07/03 05:23)
「[ジャガーさんの不思議なシフト] 目次: 車ジャガーXE（前期型）はシフトがダイヤル型になっていて、P, R, N, D, S...」
20年2月22日
すずき (06/26 02:03)
「[Zephyr - まとめリンク] 目次: Zephyr導入、ブート周りHello! Zephyr OS!!Hello! Ze...」
26年6月18日
すずき (06/26 02:02)
「[ZephyrのOut-of-treeアプリケーションその6 - シリアル出力] 目次: Zephyr前回はシリアルに文字を出...」
26年6月23日
すずき (06/26 01:27)
「[ANA国内線予約サイトが悲惨なことに] 4月くらいにANAの国内線予約サイトが国際線と共通のクラウドシステム（Amadeus...」
26年6月11日
すずき (06/26 01:10)
「[ZephyrのOut-of-treeアプリケーションその5 - OpenOCDとGDBで実行] 目次: Zephyr前回Ze...」
26年4月29日
すずき (06/24 01:10)
「[ぽこあポケモンをクリア] 目次: ゲーム1日だけやって放置していたぽこあポケモンをクリア（＝スタッフロールが流れるイベントを...」