平均が0でなく相関がある正規分布の比(商)Z=X/Yの分布 Fieller-Hinkleyの式（平均が0ならコーシー分布）をカシオの高精度計算サイトkeisan.casio.jpにあげていたが式にミスがあったのでやり直し。Pythonでモンテカルロして全然違う値になったのでもっと新しい論文探した。
ずっと前に
というのをカシオの高精度計算サイトkeisan.casio.jpに自作式としてUPしていた。
たまたま正規分布の比が必要となって、Pythonで相関があって分母の平均が0でないような正規分布を乱数で作って
（np.random.default_rng().multivariate_normalを使う）
その自作式と比べたら…形は似ているけど絶対値が全然違う！なんか式が古いので誤記とかあるのかな…と思ってもっと新しいのを探したらあった。
こちら：
Confidence Interval for the Ratio of Two Normal Variables(an Application to Value of Time)
Hisa MORISUGI, Jane ROMERO, and Takayuki MORIGUCHI
Interdisciplinary Information Sciences Vol. 15, No. 1 (2009) 37–43
https://www.jstage.jst.go.jp/article/iis/15/1/15_1_37/_pdf
この式を使った。
でPythonで
Ratio distribution
にあった事例
をやってみることにした。結果がこちら。オレンジが乱数で作ったモンテカルロシミュレーションのヒストグラム、
青線が上に書いた式の計算結果だが、見事に重なっている。早速keisan.casio.jpの自作式も修正しておいた。
検証に使ったPythonのプログラムはこちら。
« 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(5) 中二病のような名前のブルースカイカタストロフィ(Blue-sky catastrophe)を生じるGavrilov Shilnikov modelをドルマン・プリンス法（ode45)で計算する。テキストベース（アスキーアート）で！ | トップページ | 原著で持っているダンズモアさんのHandbook of Microwave Component Measurements With Advanced VNA Techniques(2nd ed.) が邦訳されてマイクロ波部品測定ハンドブック(第2版)として出版されてる！Pozarさんのマイクロ波工学に続き、最近名著の翻訳が増えてるのはいい傾向！ »
「パソコン・インターネット」カテゴリの記事
- Excelで1つのセルに複数の値を入れられるリストやセル内配列、ネストした配列という新機能が使えるようになっていた。例にあったデータ数が行ごとに違うときに使うのがよさそうなのでネイピア数の級数での計算をやってみた。(2026.10.07)
- JevのAPI互換のtev1(Ollama経由で利用）を使って、バンドパスフィルタの測定Sパラメータを読み込みスペックと比較して、単にG／NG判定だけでなく重大度や容認度をスコアとして出してくれるアプリrf_jev_judge.pyを作った。(2026.10.05)
- RF Weekly Digest 2026/9/28-10/4 (Codex(GPT-6 Astra）によるRF情報週刊まとめ)(2026.10.04)
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2026年10月9日 IEEE Microwave Magazineにリコンフィギュラブルフロントエンド記事、Keysightがレイヤー1オートメーションスイッチ発表、GSMAがインドで3つのレポート発行、MinewのSOSボタンにNordicのSoC採用など(2026.10.09)
- 高周波・RFニュース 2026年10月8日 SkyworksがQorvoとの合併を完了、iFixitがApple Watch Series12とUltra4を分解・60GHzアンテナあり、TDKが2GHz対応車載インダクタ発表、KeysightがADS等にAIエージェント搭載、Mini-Circuitsのバトラーマトリクス解説(2026.10.08)
- 1980～2026年のノーベル物理学賞を素粒子・宇宙と物性・基礎に分けてプロット。昔は2つの分野は交互に獲る、という感じだったが去年まで5年連続物性・基礎だった。今年はアイスキューブニュートリノ観測所での成果ということで久しぶりに素粒子・宇宙。(2026.10.06)
- 高周波・RFニュース 2026年10月6日 QorvoがBAWベースのCバンドフロントエンドとXバンドフロントエンド発表、太陽誘電が低直流抵抗のメタル系パワーインダクタ発表、アンリツとキャンドックスシステムズが60GHz評価システム発表、QualcommとHuaweiが特許ライセンス契約(2026.10.06)
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2026年10月9日 IEEE Microwave Magazineにリコンフィギュラブルフロントエンド記事、Keysightがレイヤー1オートメーションスイッチ発表、GSMAがインドで3つのレポート発行、MinewのSOSボタンにNordicのSoC採用など(2026.10.09)
- 高周波・RFニュース 2026年10月8日 SkyworksがQorvoとの合併を完了、iFixitがApple Watch Series12とUltra4を分解・60GHzアンテナあり、TDKが2GHz対応車載インダクタ発表、KeysightがADS等にAIエージェント搭載、Mini-Circuitsのバトラーマトリクス解説(2026.10.08)
- 1980～2026年のノーベル物理学賞を素粒子・宇宙と物性・基礎に分けてプロット。昔は2つの分野は交互に獲る、という感じだったが去年まで5年連続物性・基礎だった。今年はアイスキューブニュートリノ観測所での成果ということで久しぶりに素粒子・宇宙。(2026.10.06)
- 高周波・RFニュース 2026年10月6日 QorvoがBAWベースのCバンドフロントエンドとXバンドフロントエンド発表、太陽誘電が低直流抵抗のメタル系パワーインダクタ発表、アンリツとキャンドックスシステムズが60GHz評価システム発表、QualcommとHuaweiが特許ライセンス契約(2026.10.06)
« 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(5) 中二病のような名前のブルースカイカタストロフィ(Blue-sky catastrophe)を生じるGavrilov Shilnikov modelをドルマン・プリンス法（ode45)で計算する。テキストベース（アスキーアート）で！ | トップページ | 原著で持っているダンズモアさんのHandbook of Microwave Component Measurements With Advanced VNA Techniques(2nd ed.) が邦訳されてマイクロ波部品測定ハンドブック(第2版)として出版されてる！Pozarさんのマイクロ波工学に続き、最近名著の翻訳が増えてるのはいい傾向！ »
コメント