ギヤ比

ソースコード
https://github.com/index333/gear

ギヤ比は一般には単に前後の比率で表します。
例えばフロント３０T　ｘ　リヤ２０Tの場合
３０➗２０で１．５０

欧米では伝統的にダルマ型自転車の前輪径に換算します。
３０ｘ２０の場合
２７ｘ１．５＝４１インチギヤ
のように表現します。
欧米式の利点は小径車にも統一的に使えることです。
上の場合
２０ｘ１．５＝３０インチギヤとなります。
逆にロードレーサーの４１インチギヤと同じ速度を小径車で得るには
３０ｘ１５あるいは３９ｘ２０という前後の組み合わせが必要となります。

---

ケーデンスを固定した時の時速で表したほうが直感的と思います。

まずタイヤ周長を選択します。

データはキャットアイのページからコピーしました。

$ cat tires
18-622 700x18C 2070 207
19-622 700x19C 2080 208
20-622 700x20C 2086 209
23-622 700x23C 2096 210
25-622 700x25C 2105 211
28-622 700x28C 2136 214
30-622 700x30C 2146 215
32-622 700x32C 2155 216
xx-xxx  Tubular 2130 213
35-622 700x35C 2168 217
38-622 700x38C 2180 218
40-622 700x40C 2200 220
42-622 700x42C 2224 222
44-622 700x44C 2235 224
45-622 700x45C 2242 224
47-622 700x47C 2268 227

＄　cat tires | runghc selectItem.hs > tmp
$ cat tmp
23-622 700x23C 2096 210

---

Round モジュール、MySpinBox　モジュールが必要です。
$ cat tmp | runghc gear.hs

ｐｃｄ，歯底距離

チェーンリングのPCD（ピッチ・サークル・ダイアメーター）を計算します。
チェーンリングのPCDとはチェーン用ギヤにチェーンを巻きつけたときチェーンピンが形作る仮想円の直径のことです。

これは１辺の長さが１２．７ｍｍの正N角形に外接する円の直径の計算に還元できます。

(
クランクのPCDと区別するため、あちらはBCD（ボルト・サークル・ダイアメーター）と表現します

)

次に歯底距離を計算します。

歯底距離を測ればチェーンリングの歯数がわかります。

$ ghci
Prelude> :l Ch.hs
[1 of 1] Compiling Ch               ( Ch.hs, interpreted )
Ok, modules loaded: Ch.
*Ch> rd0 43
(43,166)
*Ch> rd0 44
(44,170)

歯底距離が１６６ｍｍであれば歯数は４３T
歯底距離が１７０ｍｍであれば歯数は４４T

チェーリングのPCDとリヤスプロケットのPCDがわかれば、シングル・バイクのチェーン長が計算できます。

芯間距離（チェーンステイ長）４０５ｍｍ、前スプロケット歯数４８Ｔそして後スプロケット歯数１６Ｔの場合

*Ch> chlen 48 16 405
1226.491724535387
*Ch> it / p
96.57415153821945

９７リンク。（半コマ使用時）

９８リンク。（半コマ不使用時）

---

チェーン長計算は下図を参考にしてください。

SpinBoxで入力を楽に

初期データは「tmp]に

$ cat tmp
[48.0,16.0,405.0]




スプロケを交換したときに
リヤアクスルが何ミリ移動するか計算します。

ランダム法

「Haskell98標準のrandomパッケージは、1.長くない周期、2.Haskell98由来の古いAPIの制約で内部状態を毎回コピーする、3.そのため遅い、4.CPU時間と現在時間で初期化(ポケモンかな？？？)など、過去に標準だったため多くの環境で使用出来ると考えられる、という点以外に余りお勧めできる所が有りません。」
（Haskellの乱数事情より）
https://qiita.com/philopon/items/8f647fc8dafe66b7381b

だそうですがここでは便利さをとり
 System.Random
モジュールを使います。

 System.Randomモジュールはcabalコマンドでパッケージを導入しなければなりません。
まず

sudo apt install cabal-install

でcabalをインストール。

$ cabal　update　-ｖ
$ cabal install random　-ｖ

でImportできるようになります。
結構時間がかかるので　-v オプションをつけて経過を見ながら。

なお表題の「ランダム法」はランダムを使っているという意味であり「モンテカルロ法」とは関係ありません。

基本のコードは

---

これで、未知数が一つ、解の範囲があらかじめ推測できるという条件の方程式を解くことができます。

上のコードに

root n x = evalStateT (try (pf n x)) (1,x) >>= return
pf n x g = g^n < x
main = root 2 2 >>= print

を付け加えて実行しますと

$ runghc -XFlexibleContexts  randomMethod.hs
(1.4142135623728276,1.4142135623731693)

---

3角形の３辺の長さがわかっていれば、その形と大きさが決まることから面積を計算できるはずです。
ヘロンの公式です。
ちなみにヘロンの公式はHaskellでは極めて簡潔にこう書けます。
heron a b c = sqrt (s*(s-a)*(s-b)*(s-c)) where s = (a+b+c) / 2

余談ですがどんな複雑な形の土地でも３角形に分割することで巻き尺１本で面積計算ができます。
（ただし平面地）

ヘロンの公式が思い出せなくてもランダム法で直感的に解くことができます。

三角形の３辺の長さが　a,b,c で

ｂを複素平面上の実数軸においたとき、

三角形の頂点は極座標で

（a,θ）

で表せます。

このθをランダム法で探ればよいのです。

---

$ runghc -XFlexibleContexts  heron.hs
6.0

円に内接する正６角形と、同じ円に外接する正４角形を考えると、
円周率は３と４の間に存在するとわかります。
３と４の間の実数のどれかはπの近似値になるはずです。
では、それをランダム法で求められるでしょうか。
残念ながらできません。円周率を求める方程式が存在しないからです。

複素数を使えば間接的にπを求めることはできます。
複素数　ｃ＝-１＋0i　を極座標で表しますと　（1、π）
位相を未知数θとして　θを　３と　4　の間で、さまざまに変化させ　ｃ　と近似する値をみつければよい。

でもこれはチートなのでは？