言葉の通り信号を処理することだが、もっぱらディジタル信号処理の意味で使われる→DSP
どもです。 およそ1年半前に公開したFIR_x2ですが、処理の非同期設計を極力なくす作業を実行中です。 github.com audio-diy.hatenablog.com audio-diy.hatenablog.com その中で「公開したのは良いけど、これ自分の設計したFIRフィルタの組み込み方法書いてないな」ということに気づいたので、FIR_x2に組み込む形式のFIRフィルタ係数データを生成するサンプルプログラムを追加しました。 このブログではその使用方法について簡単に記したいと思います。 生成AIすげぇな ここに書く内容の英語版は、GitHub上にまとめてあります。 github.c…
私が作ったコンテンツの中でダントツで一番人気なものが、Savitzky-Golayフィルタの記事です (検索でのアクセスの20%近く)。しかも、それはこのブログだけじゃなくて、私のYouTubeチャンネル (@kiyonolab) でも、世界的に同じ現象が見られます。 正直、「えっ、みんなそんなにSavitzky-Golayフィルタに興味があるの⁉」とびっくりしています (Savitzky-Golayフィルタ以外では、興味を持たれない私のコンテンツもどうかと思いますが...)。それだけ興味を持ってくれる人が多いのであれば、みなさんの理解の手助けになるコンテンツを作る必要があると感じました。とい…
パワースペクトルがあれば、線形確率過程や弱定常過程の動的特性をバッチリ (完全に)特徴づけることができます。つまり、線形確率過程では、 パワースペクトルは、時系列データの「DNA」 みたいなもので、パワースペクトルは時系列を生成するルールのすべてを記述しています (線形確率過程の場合だけですよ)。 とはいえ、現実世界の時系列データが、そんな理想的な線形確率過程のモデルにピタッと当てはまるかというと、そうでもありません。現実は非線形です。でもまあ、最初の一歩としてパワースペクトルを推定してみるのは、悪くないアイデアです。 ゆらぎの霧に包まれた真のパワースペクトル 現実のデータには限りがあります。…
時系列のパワースペクトルを理解したり，フィルタの特性を理解したり，デルタ関数の性質を理解したりする基礎知識として，矩形関数とsinc関数の関係を理解しておく必要があります．ということで，今回は，矩形関数とsinc関数についてのお話です． 下の絵のように，矩形関数は，時系列の切り出し (信号と矩形関数の積)や，移動平均フィルタ (矩形関数との畳み込み)などに応用されています．そのような信号処理の理解には，sinc関数 (しんくかんすう)の知識が不可欠です．加えて，「掛け算」と「畳み込み」の関係を理解する必要もあります． 時系列の一部を切り出すことは矩形関数を掛けることと同じ (左)，時系列の移動…
私が，Savitzky-Golayフィルタと出会ってから，もうすぐ１０年になります．私は，Savitzky-Golayフィルタをとても気に入っていましたが，残念ながら，さよならするときが来ました． 最近，Savitzky-Golayフィルタは理想のフィルタではないことを認識しましたので，Savitzky-Golayフィルタに変わるFIRフィルタを探していきたいと思います． Savitzky-Golayフィルタの欠点については，前の記事 chaosmemo.com を参考にしてください． 今回は２次Savitzky-Golayフィルタの欠点を改善した，「お手上げポーズ」フィルタを紹介します．世界…
下の図のように，時系列解析では，パワースペクトルとか，コヒーレンスとか，DFAなどのスケーリング解析とかで，いくつかの結果のグラフをまとめて，平均的なグラフを求めたいことがあります．今回は，この処理をするRスクリプトを紹介します． x軸の各点で複数のグラフの平均を計算したい． データフレームにすべての結果を並べて平均 グラフの平均を計算するＲスクリプトの例を掲載しておきます．この例では，５本のホワイトノイズ時系列のパワースペクトルの平均を計算しています． # 何回分析を繰り返すか n.samp <- 5 ### for(i in 1:n.samp){ # 時系列の例：白色ノイズ y <- rn…
これまで、当方は、信号処理において、その計算は、三角関数の実関数の加算と乗算で構成してきました。 しかしながら、公知例の中には、信号処理を複素数で計算し、その実部だけの値を、DA変換器に出力しているものがあるのを発見しました。 これは、ちょっとしたカルチャーショックです。 高校の時に、「虚数は存在しない。」と数1で習ったことが、実は旧文部省、文部科学省の教育指導要項に基づく学問上の論理ミスであることが判明るまで、ずいぶんと年数がかかり、大きなカルチャーショックを受けましたが、それ以来のかなりのショックです。 見直し内容： How to generate SSB by DSP/CPU
どもです。 昨日は"FIR_x2"について実装の工夫点をいくつかお話しましたが、今回はシミュレーションと実際の使用方法についてまとめていきたいと思います。 audio-diy.hatenablog.com 2025/7/5追記 係数生成方法もまとめました。 audio-diy.hatenablog.com 追記ここまで ディレクトリ構成 ディレクトリ構成は以下の通りとなります。 github.com FIR_x2のGitHubディレクトリ構成 01_DPRAM_CONTシンプルデュアルポートRAMをリングバッファとして機能させるためのコントローラ 02_DATA_BUFFERSPRAM_CON…
どもです。 8月の夏休みから4ヶ月ほどかけてやっと動作確認までできたので公開します。 サンプリング周波数選択もできたぁぁぁぁぁぁぁぁぁデータも問題無し！ pic.twitter.com/Gek3xQQBYP — AUDIY (@AUDIY14) November 29, 2023 PCMデータを2倍オーバーサンプリングするFIRデジタルフィルタIPです。 github.com きっかけ AUDIYがFPGAに興味を持ったのは以下のビデオがきっかけです。 CHORDというオーディオメーカーがFPGAとフリップフロップを使用してDAコンバータを実現している、ということを知り、自分でもDAコンバータ…
※この記事は同日投稿の"One Equals One Plus One Minus One"を翻訳したものです。 NTTドコモ R&D Advent Calendar 2022 の15日目の記事です。 1=1 私は DOCOMO Innovations, Inc. の Senior Research Engineer としてモバイル通信ネットワークの最適化に携わっています。ネットワーク最適化を考える上で，複雑な計算に頭を悩ませることがよくあります。ネットワークは相互に接続された多くのシステムで構成される複雑なシステムですが，単純なものから理解することで複雑なものの理解につながります。例えば次の…
www.ic.is.tohoku.ac.jp
aidiary.hatenablog.com
laputa.cs.shinshu-u.ac.jp
www.ic.is.tohoku.ac.jp
www.ic.is.tohoku.ac.jp
engineer.fabcross.jp
aogiri-m2d.hatenablog.com
av.watch.impress.co.jp
fabcross.jp
