書き出しのしかた javax.sound.sampled.AudioSystem のwriteを使えば書けるようです。 書き出す前に長さを決めてデータを流し込むイメージで ちょっと調べてみた感じだとstopまで書き出すみたいな使い方は出来なそうです。 内部的には指定した長さまでfor分…

先日WAVが読み込めるようになったので いろいろやってみるためにもまずはビューアの作成を進めます。 特殊なことはなにもないのですが グラフの表示とか拡大縮小とかマウスの処理がただただめんどくさいし 設計のようなことをしないで作っているのでバグだら…

久しぶりに時間が出来たので 音声ファイルを触ってみます。 フォーマットの確認 とりあえずwindows7のログオンの時の音を読み込んでみる AudioInputStream input = AudioSystem.getAudioInputStream(file);AudioFormat format = input.getFormat(); でフォー…

少子高齢化といわれていますが あまり実感がわかないのでちょっと確認してみます データはこちらからもらいました。 統計局ホームページ/日本の統計 書けって書いてあるのでこれも表記しておきましょう 出典 総務省統計局「日本の統計 2014」 ぐぐれば表はい…

お金の話ばっかり書いてますが これを海外のデータと比べて見たいのでまずは日本のデータを羅列。 海外やる前に飽きそうな不安。 例によって情報ソースはこちら 民間給与実態統計調査結果｜長期時系列データ｜国税庁 大企業は高収入なのか 結果は明らかなの…

平均収入の30年の推移 収入分布の30年の推移 を確認したので 今度は年齢別に見てみます。 データは前回同様 民間給与実態統計調査結果｜長期時系列データ｜国税庁 ここから持ってきたデータをグラフにしただけです。 2012年時点での平均値は 全体平均408万円…

前回平均収入の推移を見てみたので 今度は分布を見てみます。 今回も1世代前、父親世代との比較ということで 30年を目安とします。 データの出展は国税庁。 民間給与実態統計調査結果｜長期時系列データ｜国税庁 収入の分布データを割合に直して見ました。 …

研修で調べ物をしてるついでにいろいろ気になったので 気になったものを調べて見ることに。 情報ソース 割と国の機関で情報を公開してるのですが 一応ここが窓口になっている模様。 政府統計の総合窓口 平均年収はここ30年でどう変わったか まずは国内の平均…

飽きたのと仕事が忙しくなったので しばらく放置でしたが 研修でいろいろ調べてていたところ 割といろいろなデータが公開されていることを知りました。 面白いのでちょっとこのあたりを調べながら 書いておこうと思います。 最近はパズドラもツムツムも飽き…

ちょっと前のことですが 友達がパズドラのサポートをするようなアプリを作っていたので 僕も楽しげなところだけ作らせてもらいました。 公開する気はないようなのですが GUIまでわりとちゃんとつくっていました。 よくある大学の課題みたいなものなので 思っ…

中国に行ってきました。 場所は合肥上海から内陸に行ったあたりです。 僕の印象は三国志で魏と呉が戦った場所で 合肥の戦いで有名な場所です。 合肥の戦い - Wikipedia PM2.5 やっぱり多少の影響はあるようで 微妙に曇っています。 微妙に喉というか胸が痛い…

自分が非常に困ったので誰かの参考になれば。 学生の方向けに。 私は英語がかなり苦手です 就職活動で初めて受けたTOEICが４５０くらい。 大学受験くらいまで 英語はどっちかっていうと文系科目だろーくらいに思ってました。 実際は全くそんなことなくて 研…

内容は以前書いたものと同じですが 前回LPFだけだったのでちょっとまとめます。 条件 入力はFs=48kHz f=1kHzの矩形波 これにFFTをかけて特定の周波数を除去したあとに 逆フーリエ変換をかけて波形の確認をします。 f0を振って波形を確認。 LPF f0=1000のとき…

どうもエッジを検出するには フーリエ変換して周波数軸で低周波を削除したあとに 逆変換して隣のピクセルと正負が反転しているところが エッジにみなせるように思えてきた。 理論を調べるのはつまらないのであとでやるとして といあえず実験してみます。 元…

二次元FFTは画像を画素の行列にして 縦にFFTしたあとに横にFFTすればいい。 戻す場合は横に逆FFTしてから縦に逆FFT これで元にもどりそうだということは前回できたので 周波数軸で操作して遊んでみます。 元画像 例によってこれ。 LPF 中心付近が高周波なの…

いまさらながら画像を載せたいときは 編集画面に画像をドロップするんじゃないということに気づきました。 重い理由はこれだったんですね。 解決。

ゲームタイムが終わったのでつづき。 先ほどの画像を ↓方向でFFTかけたあとに→方向でFFTしてみます。 ほぼ真っ白。 よくわからない。 全ピクセルを調べて最高値が255になるようにスケーリングしてみます。 真っ黒。 でもよく見ると一番左上だけ白い。 1000倍…

目的はとくにないのですが 途中からこれがやってみたくていろいろ遊んでいました。 とりあえずいつもの元画像。 これが640x480なのでFFTをかけるには4^nの数にする必要があります。 周期関数という前提があるので周りの色は同じ色がよさそう。 ということで…

FFTが時間軸の信号を周波数軸に変換するのであれば FFTしたあとの信号から特定の周波数を取り除くと フィルターとして動作するんじゃないだろうか？？ たしか教科書とか読むとするって書いてあるんだけど 読んでも理解できないのでやってみることにします。 …

絵に描いたような三日坊主の理由は 買ったゲームが届いたのでそっちに集中してたからです。 さて 先日FFTのモジュールができたので マイク入力の音声信号をFFTしてみることにします。 条件 サンプリング周波数はDVD同等の48kHz ビット深度は16bit 24bitでや…

めげずにもうちょっと調整 このあたりが限界か！？ というかはてなブログすごーく重い。 画像乗せられないんだろうか。

前回ハイパスフィルタを両方向からかけると エッジ検出が出来そうなことがわかったので ちょっとその方向で実験しました。 方針 １．横方向に左からHPFをかける ２．横方向に右からHPFをかける ３．隣のピクセルと比べて正負が反転してたらエッジとみなす ４…

前回画像にフィルタをかけていたら →のあと←にかけると同じになったように見えたので 一応確認してみることしました。 フィルタ 例によってこいつを こいつにかけてみます。 結果 まず順方向と逆方向のグラフ 赤い方はフィルターのレジスタをリセットしたあ…

2次の伝達関数を持つローパスフィルタを作った場合 位相特性はf0（＝カットオフとか遮断周波数）のときをだいたい90度にして 180度まで回ります。 と大学の教科書に書いてありました。 位相の実測は難しそうなのでやりません。 180度というのは振幅の半分ず…

作ったフィルタを画像にかけてみます。 元画像 これを素材にします。 よく見るといろいろ写っていますがあまり気にしない。 1次元データなので赤青緑の成分に分けて その後3枚の絵を重ねます。 赤 青 緑 一見すると単色化して色を付けただけのように見えます…

逆変換はFFT後の値に下記処理を行う (とWikiに書いてある) ・各値の複素共役をとる ・FFTをかける ・出力の複素共役を取ってn(要素数)で割る。 のでやってみる。 前回のFFT結果 逆変換後の波形と元の波形の比較 完全に一致する。 虚数成分も微妙に値を持って…

フーリエ変換の時間軸を離散化したものが 離散フーリエ変換(discrete fourier transform)で その対象性を利用して早く計算できるようにしたものが 工学的には非常に良く見る高速フーリエ変換(Fast Fourier Transform)。 理屈はいまいち良くわかっていません…

LPFの過渡応答を確認。 教科書どおりなので動作は正しそうです。 LPF 2次のローパスフィルタを作ってQの値を変えてみます。 Qが1/√２を超えるとゲインを持ちますが 矩形波を入力するとQ=0.7でもオーバーシュートを持つようです。 HPF 続いてHPF 矩形波を入れ…

ゲインと位相を出したいのですが 理論値の計算がちょっとむずかしそうなので 入力したサイン波と出力したサイン波のピーク値を比べて ゲインを実測することにしました。 http://www.musicdsp.org/files/Audio-EQ-Cookbook.txt をそのまま実装。 まずはLPF そ…

何はともあれグラフが見れないと話にならないので作ります。 こいつらを吐けるようにしました。 見えないと何がなんだかわからない。 まずはサイン波 ちょっとうねうねしてるけどあまり細かいところは気にしない。 次は方形波 動作はよさげ。 波形表示そのも…