コグノスケ

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

過去から未来へ表示

もっと前

2014年11月23日 >>> 2014年11月10日

もっと後

2014年11月20日

permalink

編集する

DA変換

前回の日記（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も計算可能ですが、長くなってきたので、また後日にでも載せます。

編集者:すずき(2014/11/22 17:13)

コメント一覧

コメントはありません。

この記事にコメントする

2014年11月15日

permalink

編集する

USB-DACと矩形波

前回は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）

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

編集者:すずき(2019/09/08 12:45)

コメント一覧

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

この記事にコメントする

もっと前

2014年11月23日 >>> 2014年11月10日

もっと後

管理用メニュー

記事を新規作成

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件

24年5月16日
すずきさん (05/21 11:41)
「あー、確かにdpkg-reconfigu...」
24年5月16日
hdkさん (05/21 08:55)
「システム全体のlocale設定はDebi...」
24年5月17日
すずきさん (05/20 13:16)
「そうですねえ、普通はStandardなの...」
24年5月17日
hdkさん (05/19 07:45)
「なるほど、そういうことなんですね。Exc...」
24年5月17日
すずきさん (05/19 03:41)
「Standardだと下記の設定になってい...」

もっとみる

最近の記事20件

24年5月19日
すずき (06/04 00:44)
「[Yocto - まとめリンク] 目次: YoctoHello YoctoYoctoのセットアップスクリプトとビルドディレクト...」
24年5月31日
すずき (06/04 00:43)
「[Bitbakeのクラス] 目次: YoctoYocto Scarthgap（5.0.1）のメモです。前回同様、コードを読んで...」
24年5月30日
すずき (06/04 00:35)
「[Bitbakeのレイヤー] 目次: YoctoYocto Scarthgap（5.0.1）のメモです。Yoctoの使い方は以...」
24年5月29日
すずき (06/04 00:09)
「[Yoctoのビルドディレクトリとテンプレートディレクトリ] 目次: YoctoYocto Scarthgap（5.0.1）の...」
24年5月28日
すずき (06/04 00:02)
「[Yoctoのセットアップスクリプトとビルドディレクトリ] 目次: YoctoYocto Scarthgap（5.0.1）のメ...」
20年7月11日
すずき (06/03 23:59)
「[STATIONflow実績コンプリート] 目次: ゲームSTATIONflowの実績をコンプリートしました。「ラッキーセブン...」
20年7月1日
すずき (06/03 23:59)
「[STATIONflowランク100] 目次: ゲームSTATIONflowのランクが100になりました。何か実績と紐づいてい...」
20年6月30日
すずき (06/03 23:59)
「[STATIONflowの駅の評価] 目次: ゲーム以前（2020年5月28日の日記参照）STATIONflowで速度3にする...」
20年6月28日
すずき (06/03 23:59)
「[STATIONflowまさかの実績解除方法] 目次: ゲームSTATIONflowのしょうもない小技その2です。実績の解除...」
20年6月27日
すずき (06/03 23:59)
「[STATIONflowプレイ日記] 目次: ゲームSTATIONflowの基本は理解したつもりなので、実績解除に挑んでますが...」
20年6月26日
すずき (06/03 23:58)
「[STATIONflow小技] 目次: ゲームSTATIONflowのしょうもない小技。将来的に、どこに駅の入り口と電車の乗り...」
20年5月28日
すずき (06/03 23:58)
「[STATIONflowのバグ] 目次: ゲームSTATIONflowバグってますね……。ゲーム進...」
20年5月27日
すずき (06/03 23:58)
「[STATIONflowランク20] 目次: ゲームゲーム進行的に1つの区切りと思われる、ランク20を超えました。駅が広くなる...」
20年5月19日
すずき (06/03 23:58)
「[STATIONflow始めました、超えろ、新宿駅] 目次: ゲームSteamで新たなゲームを買いました。2020年4月ローン...」
21年12月28日
すずき (06/03 23:58)
「[ゲーム - まとめリンク] 目次: ゲーム一覧が欲しくなったので作りました。Wizardry（囚われし亡霊の街）敵が強すぎる...」
20年7月13日
すずき (06/03 23:57)
「[STATIONflow - まとめリンク] 目次: ゲームに統合。...」
23年9月18日
すずき (06/03 23:49)
「[一覧の一覧 - まとめリンク] 一覧の一覧、まとめのまとめが欲しくなったので作りました。OS、アーキテクチャ系。目次: An...」
24年5月27日
すずき (06/03 23:46)
「[Hello Yocto] 目次: YoctoYoctoのメモです。なおYoctoのバージョンはScarthgap（5.0.1...」
24年5月21日
すずき (05/23 23:19)
「[Linux 6.1からLinux 6.6に手抜き更新したらハマった] 目次: 自宅サーバーLinux Kernelのlong...」
23年6月1日
すずき (05/23 22:50)
「[自宅サーバー - まとめリンク] 目次: 自宅サーバーこの日記システム、Wikiの話。カウンターをPerlからPHPに移植日...」