雑感等

音楽,数学,語学,その他に関するメモを記す.

扇型に内接する楕円

扇型に内接する楕円を求めた ――かったが、実質的には二等辺三角形に内接する楕円の式が求まった。 扇型は、半径を1、中心角をとする。 楕円は、一方(扇型の径方向)の半径を、もう一方の半径をとする。 各図形は図の座標に配置し、とを既知の値とすると、は…

丸かっこで囲まれた部分以外にマッチする正規表現

丸かっこで囲まれた部分以外にマッチする正規表現 [^\(\)]+(?=\([^\(\)]+\))|[^\(\)]+$|^[^\(\)]+ テキストエディタMeryで動作確認済み。 説明 この正規表現は、orを意味するパイプ記号"|"と、3つの部分(下記ABC)で構成されている。 A | B | C ↓ [^\(\)]+(?=…

物理モデルで自励振動させようとした

いろいろ調整してみたら自励振動させられてないっぽい。管体部分でループ(ハウリングみたいな状態)していただけだった。 物理モデリング音源の特許(だいたい権利失効してる)の実施例にある、励振部をpythonで実装した。 一応音は出るが、パラメータの調整…

Rustで配列(Vec)の要素数を変数で指定:vec!マクロ

出典 https://users.rust-lang.org/t/attempt-to-use-a-non-constant-value-in-a-constant/32112#:~:text=Sep%20%2719-,Use%20vec!%5BDefault%3A%3Adefault()%3B%20n%5D%20or%20Vec%3A%3Awith_capacity(n)%20and%20push%20to%20it.,-Solution # 配列の要素数…

複数フォルダをそれぞれzipに圧縮(バッチファイル)

↓このあたりのソースを改変して、圧縮後のzipファイルの保存先を指定できるようにした。 rezv.net oshiete-suitman.com batプログラムソースは以下。 【保存先フォルダ(フルパス)】のところに、保存先フォルダのパスを書く。 ZIP_PATHには7zipの実行ファイ…

二次元の波動方程式を陰解法+ガウス・ザイデル法で解く(Rust)

バグを直した&メッシュ数に応じて出力画像の比率を変えるようにした。 ■出力画像 ■Rustソース use std::time::{Duration, Instant}; use plotters::prelude::*; const PLT_RATIO: usize = 20; const T_MAX: usize = 1000; const X_MAX: usize = 10; const Y…

二次元の波動方程式を陰解法+ガウス・ザイデル法で解く(Rust) ※バグあり

※バグあり メッシュ数を縦≠横にするとおかしくなる とりあえずソースと結果を貼っておく。 陰解法と書いたが、実際は陰解法と陽解法を平均した差分方程式を使っている。(クランクニコルソン法っぽくした) ↓の記事で導出した式と、陽解法を組み合わせている。…

波動方程式の有限差分法を半自動で導出(sympy)

sympyで差分化した2次元の波動方程式を定義して、式の整理を自動で実行する。 jupyterで実行時、最終的な実行結果は下記。 uの添え字は、1番目が時間、2番目と3番目がxとy。 左辺がu_{1,0,0}なので、次の時間ステップにおける波の振幅を表す。 右辺は、u{0,?…

波形のメモ

aは指数的に変化させる 縦軸y 横軸t

lilypondコンパイルエラーの対処 error: syntax error, unexpected SCM_TOKEN, expecting '='

コンパイルエラーの修正 lilypondファイルをコンパイルすると、以下のエラーメッセージが出力。 ././01-byrd-a3-0-score.ly:36:62: error: syntax error, unexpected SCM_TOKEN, expecting '=' \override StaffGrouper #'staff-staff-spacing #'padding = #4…

波動方程式を有限差分法の陰解法で解く(Julia)

波動方程式を↓の文献にある式で数値計算する。 https://www.uni-muenster.de/imperia/md/content/physik_tp/lectures/ws2016-2017/num_methods_i/wave.pdf https://www.uni-muenster.de/Physik.TP/archive/fileadmin/lehre/part2_hypebolic/node20.html 記号…

Karpus-Strong AlgorithmをJuliaで実装

Karpus-Strongアルゴリズム(K-Sアルゴリズム)を↓ページのブロック線図に基づきJuliaで実装した. https://ccrma.stanford.edu/~jos/pasp/Karplus_Strong_Algorithm.html 実行すると生成波形のグラフが表示され音が再生される. ↓生成波形例 K-Sアルゴリズ…

Python+PyQt5で4次元図形の描画

HSP(Hot Soup Processor)で4次元図形を回転させて描画 - 雑感等 これの二番煎じだけどクラスの設計はMVCを意識したつもり. 4次元図形の回転は上の記事と同様 https://ch.nicovideo.jp/4dimensions/blomaga/ar969330 https://www.researchgate.net/publicati…

pythonでmidiメッセージ送信

pythonのライブラリmidoを使用してリアルタイムにmidiメッセージを送信する. 以下のプログラムでは,ランダムな音高・音価・サステインペダルの踏み度合いを送信する. pianoteqのstandaloneとloopMIDIも合わせて使っている. import mido import time impo…

変分法のオイラー方程式の変形

Rudan, Massimo. Physics of Semiconductor Devices. 2015. Springer. にあった変形 オイラー方程式 (1.4) から(1.5)の導出が分からなかった. # の微分の順番を入れ替えて,(無条件で入れ替えられるかは不明) (1.4)' ここでに多変数関数の偏微分の連鎖律…

プログラムが終了したらメールを送信するpythonプログラム

windows 7や10で動作するはず. 使い方 "監視する実行ファイル名.exe"がすでに実行されているときに, 以下のプログラムを開始する. time.sleep(10)のとおり,10秒おきにプログラムが終了したか確認する. プログラムが起動していなければ(終了していれば…

格子ボルツマン法の局所平衡分布関数を自動で展開するプログラム

格子ボルツマン法(Lattice Boltzmann Method, LBM)の局所平衡分布関数(local equilibrium distribution)を自動で展開するプログラム. 数値計算のループを実行する際に,展開された(ベクトル,行列を使わない)式を使うことで計算の高速化が期待できるかも…

JupyterLabでコード整形(フォーマッタはautopep8)を導入

下記のサイトのコードだと"black"とかいうフォーマッタを使った例になってるらしい(一般名詞で大変紛らわしい) https://jupyterlab-code-formatter.readthedocs.io/en/latest/how-to-use.html { "shortcuts": [ { "command": "jupyterlab_code_formatter:b…

scipyのsolve_ivpでローレンツ方程式を解く

以下のページを参考にした. https://org-technology.com/posts/ordinary-differential-equations.html https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.solve_ivp.html from scipy.integrate import solve_ivp import matplotlib.py…

フィボナッチ数列の兎の個体数増加をグラフで

フィボナッチ数列的な個体数の変化を,下記リンクの図のように表したかった. https://commons.wikimedia.org/wiki/File:FibonacciRabbit.svg 生成される図は下記のようになる. 上記の図は以下のプログラムで生成した. import matplotlib.pyplot as plt fr…

CNNでのBatch Normalizationの平均・分散

畳み込みニューラルネットを使ってるときにバッチ正規化したら, BN層のパラメータが(ch数×4)になった. これはchごとに正規化していて,各chに含まれる各要素(各画素や各時刻の電圧)をそれぞれ正規化しているわけではない. BNの解説では,ミニバッチに含…

secret pure love(2号.)のコード進行

https://youtu.be/BEiyIYEMqIM?t=179 | G- C- F- Bb- | CΔ ~ m7で五度圏下って変終止

ブル調のII-V

IIm7→V7にII Dorian使うのってIV Lydian使うのと一緒か

対数確率関数

の係数が,の係数をで割ると.

ベイズ推論

モチベーション 点推定よりも多くのことがわかる 分布全体がわかる 推定量の自信のなさがわかる ベイズ推論と「ベイズの定理」の関係 ベイズの定理は確率分布と確率分布の組み合わせ方を示すだけ 事前分布(常識とか前知識)と尤度関数(観測したこと)とを…

ベイズ的な式の表記

ベイズでありがちな確率とか積分とか総和とかが入り混じった式を陽的に表現する方法はないのか

neural ODEを理解したい

https://arxiv.org/abs/1806.07366 neural ODEがわからない.この理由として考えるのは, のfがどうなってるのかわからない tがどのようにfの中でパラメータ化されてるかがわからない neural ODEがResNetやニューラルネットの拡張であることを考えれば, fに…

pythonで並列(joblib)使った際のパーミッションエラー

https://github.com/joblib/joblib/issues/806 albertcthomas commented on 3 Dec 2018 @robintwhite for now, if this is possible in your case, you can disable memmaping by setting max_nbytes=None in Parallel Parallel(n_jobs=-1, max_nbytes=None)…

PRMLの確率的主成分分析の式とPCA/ZCA白色化

https://www.microsoft.com/en-us/research/people/cmbishop/prml-book/ https://cbrnr.github.io/2018/12/17/whitening-pca-zca/ PRMLのp.574の確率的主成分分析(PPCA)に関する式(12.45) において, とした(確率的ではない普通のPCAとした)上で, とすれ…

集合による数の表現

自然数 0 = {} k+1 = { k , { k } } 名前 展開 実体 0 {} 1 { 0, {0} } { {} , { {} } } 2 { 1, {1} } { { {} , {{}} } , { { {} , {{}} } } } 3 { 2, {2} } { { { {} , {{}} } , { { {} , {{}} } } } , { { { {} , {{}} } , { { {} , {{}} } } } } } 整数 n…