XScale

2012/02/21 20:00:13みたいな書式は人間にとっては良いけど、パソコンで扱う（ソートするとかプロットするとか）には不便なのでタイムスタンプにするのが良い。2012/02/21 20:00:13をタイムスタンプに変換するには次のようにする t1 = '2012/02/21 20:00:13'; formatIn = 'yyyy/mm/dd HH:MM:SS'; serial_date_number = datenum(t1, formatIn); プロットしたいときはtimeseriesとplotを使う。ただし、timeseriesの時間データ(time)の単位は（デフォルトでは）…

マルチバンド処理用のフィルタバンク実装でバターワースフィルタによるハイパスフィルタを使おうとすると，7次くらいから周波数応答が不安定になりました． ハイパスフィルタ 周波数応答をプロット してみると，scipy.signal.butterworthのハイパスフィルタは7次くらいから低域の位相が不安定になっています． おそらくフィルタ係数の量子化誤差などが関係していると考えられます。 カットオフ周波数を低くしても，それでも9次くらいから位相が歪んでいます． ローパスフィルタ 一方、ローパスフィルタのほうが比較的高次でも振幅・位相に歪みは目視では確認できません． 補完フィルタ（7次，8バンド） 高…

以前の記事では，木構造のような信号処理フローのマルチバンド処理を紹介しました． この記事では、実装・利用しやすさに着目して信号処理フローを再設計・実装します。 設計 以前の記事では木構造の処理フローでしたが，今回は低域のカットオフ周波数から徐々に算出するような構造とします。 オールパスフィルタの伝達関数をまとめると以下のようになり，各帯域の信号はそれぞれの最後のオールパスフィルタで位相が同期されることで、回路全体が元信号の振幅成分を変化させないオールパスフィルタ回路となります。 なお、各オールパスフィルタのフィルタ係数はこちらの記事にあるようにローパス・ハイパスフィルタのフィルタ係数から簡単に…

ggsurvplot()を使ってKaplan-Meier生存曲線を描く 信頼区間をつける 色パレットを変更する 線の種類を変更する 打ち切りマークを変更する 中央値を示す補助線をつける Risk tableをつける 軸について色々 軸の目盛り間隔を変更する 軸の単位を変更する 軸の範囲を指定する 凡例を変更する 日本語の文字化けを直す おわりに 参考資料 基本の作図関数plot()ではなくsurviminerパッケージを使って生存曲線を描く方法のまとめ。 まずは必要なパッケージを準備する。 library(survival) library(survminer) library(tidyver…

棒グラフ 棒グラフ コード全体 function display(count1, count2) { // グラフのwidth, height const width = window.innerWidth / 3; const height = 400; const axisHeight = 20; const axisWidth = 25; const svg = d3 .select('#count') .append('svg') .attr('width', width) .attr('height', height); const xScale = d3 .scaleBand() …

こんにちは。takapy（@takapy0210）です。 本エントリは言語処理100本ノック2020の6章を解いてみたので、それの備忘です。 途中まで解いて放置していました（）が、続きをやる機会を得たので簡単な解説をつけながら紹介していきます。 nlp100.github.io 例によってコードはGithubに置いてあります。 github.com 第6章: 機械学習 50. データの入手・整形 51. 特徴量抽出 52. 学習 53. 予測 54. 正解率の計測 55. 混同行列の作成 56. 適合率，再現率，F1スコアの計測 57. 特徴量の重みの確認 58. 正則化パラメータの変更 59…

最終更新： 2021-06-06 GDB (The GNU Project Debugger) にはターゲット CPU 用の命令セット シミュレータが内蔵可能です。これを使用すると、ホスト上でクロスコンパイルしたバイナリ オブジェクトを GDB 単体で実行可能となります。本記事では、WSL の Ubuntu-20.04 上で GDB 10.2 を ARM CPU 用の内蔵シミュレータ付きでソースコードからビルドする手順を紹介します。 尚、本デバッガの使用方法は【こちらの別記事】で紹介しています。 目次： 本デバッガで出来ること： GDB Sim プロジェクトの Wiki ページ： ソースコード…

ある音に対して、人間が感じるうるささと音圧レベルの大きさは必ずしも一致していません。 以下の等ラウドネス曲線のように、音圧レベルと周波数によって人間の聴覚の感度は異なっています。 https://ja.wikipedia.org/wiki/音の大きさ そのため、音信号の分析、特に騒音の分析では、A特性、B特性、C特性といった周波数重み付けがよく利用されます。 すなわち、人間の聴覚の感度に合うように、対象となる音信号の音圧レベルを補正して分析を行います。 Wikipedia の A-weighting の記事では、以下の4つの特性とその算出方法が記載されています。 A特性: 等ラウドネス曲線*1…

黒体放射の式Planckの法則(Planck's law)の計算と描写 今回は波長と周波数の関数を使用。 Planck's law コード 定数の定義 波長の関数 関数の定義 描写 周波数の関数 関数の定義 描写 参考文献 Planck's law プランクの法則とは、 物質と熱平衡状態の放射が出すエネルギー密度（放射輝度）分布を表す法則 である(引用:天文学辞典)。 発見者はドイツの物理学者マックス・プランクである。 プランクの法則は、黒体の輝度(radiance) を表す式であり、温度と波長or波数or周波数の関数となっている。 式は以下のように表される(放射輝度:radiance[W m…

簡単な多重重点サンプリングの検証の続き。 ma38su.hatenablog.com 早く綺麗な絵をレイトレで描こうとした時に、ハマって確認した内容を残しておく。 直接レイトレは扱わないが、暗いシーンに限られた光源が置かれているようなシーンにおいては、収束を早くするため、ライト方向の拡散や反射を重点サンプリングすることがある。 https://raytracing.github.io/books/RayTracingTheRestOfYourLife.html#samplinglightsdirectly 参考にGLSLに書き下ろしたもの。 https://www.shadertoy.com/…

Matplotlibには、簡単にヒストグラムを描く便利な機能plt.hist()が備わっています。 とはいっても、plt.hist()は初心者の方には使い方が少し難しいかもしれません。 そこで、今回は… Matplotlibでヒストグラムを作るにはどうしたらいいの？ データ区間の変更や見た目の設定方法がわからない…。 という方のために、 Matplotlibでヒストグラムを作成する機能plt.hist()の使い方 plt.hist()を使ってデータ区間の設定やヒストグラムの見た目を設定する方法 について紹介します。 plt.hist()をマスターすれば、エクセルよりも簡単にヒストグラムが作れる…