遊び tokidoki 仕事

数学と音楽と教育と遊び

| おしごと - きょういく - がくせい - ゼミ - イベント | すうがく - おんがく - 数理音楽 - DTM - かがく - scratch
| Art - photo - おきにー - Tips - ものもう - あれこれ | About - Top

バランスをとる数(1)―入り口は1+2=3


いまや当大学における絶滅危惧レッドリストの筆頭に挙げられる数学同人Sigmaでのある一幕.
https://aue-sigma.jimdo.com/
www.aichi-edu.ac.jp

久しぶりに教育実習から返ってきたメンバーが「実習中にこんなこと見つけた」と鼻息荒く話していた話題なんだが,それは例えば
\begin{align}
1+2&=3\\
4+5+6&=7+8\\
9+10+11+12&=13+14+15\\
16+17+18+19+20&=21+22+23+24
\end{align}
といった連続整数の和についての等式だった.つまり一般に

\displaystyle
\displaystyle\sum_{k=0}^a (a^2+k)=\sum_{k=a+1}^{(a+1)^2-1}(a^2+k)
ということだ.あるいは
\begin{equation}
\displaystyle2\sum_{k=0}^a (a^2+k)=\sum_{k=0}^{(a+1)^2-1}(a^2+k)
\end{equation}
と言い換えても良い.そしてこれは任意の自然数aについて確かに正しいことが計算すれば分かる.


さて,本当の遊びはここからだ.では一般に
\begin{equation}
\displaystyle2\sum_{k=a}^x k=\sum_{k=a}^b k
\end{equation}
となるような自然数a\le x\le bの組はどんなものだろうか?計算すれば
\begin{equation}
2(x-a+1)(x+a)=(b-a+1)(b+a)
\end{equation}
すなわち
\begin{equation}
2(2x+1)^2=(2a-1)^2+(2b+1)^2
\end{equation}
となる.X=2x+1,A=2a-1,B=2b+1と置けば奇数A,X,Bの組で
\begin{equation}
A^2+B^2=2X^2
\end{equation}
を満たすものを探せということになる.あるいはs=\frac{A}{X},t=\frac{B}{X}と置けば,
\begin{equation}
s^2+t^2=2
\end{equation}
を満たす有理点(s,t)で,分母分子共に奇数になっているものを探せ,ということになる.


そこでとりあえず有理点(s,t)を探すわけだが,例えば(s,t)=(\pm1,\pm1)などがすぐ見つかる.
図形的に見やすい(-1,-1)を足がかりに,有理数m>0を傾きとする直線 t=m(s+1)-1s^2+t^2=2 の交点を考えよう.
すでに(-1,-1)が交点,すなわちこの二式の連立の解なので,2次方程式
\begin{equation}
s^2+(m(s+1)-1)^2=2
\end{equation}
(s+1)で括られ,残りの解が
\begin{equation}
\displaystyle s=\frac{1+2m-m^2}{m^2+1},t=\frac{m^2+2m-1}{m^2+1}
\end{equation}
と得られる.あとはm=\frac{q}{p}と既約分数で表示したとき,s,t共に\frac{\text{奇数}}{\text{奇数}}となっているかを確認すればよい.計算すると
\begin{equation}
s=\displaystyle\frac{p^2+2pq-q^2}{p^2+q^2},t=\frac{q^2+2pq-p^2}{p^2+q^2}
\end{equation}
となり,p,q いずれか一方のみが奇数なら分母分子共に奇数,p,qいずれも奇数ならば
\begin{equation}
p^2+q^2\equiv p^2+2pq-q^2\equiv q^2+2pq-p^2\equiv 2\pmod{4}
\end{equation}
となるため,約分すればやはり分母分子共に奇数と分かる.


こうして任意の有理数mに対する組(s,t)が求めるもの全てとなり,更にp,qいずれか一方のみ奇数ならば
\begin{equation}
A=p^2+2pq-q^2,B=q^2+2pq-p^2,X=p^2+q^2,
\end{equation}
p,q共に奇数ならば
\begin{equation}
A=\displaystyle \frac{p^2-q^2}{2}+pq,B=\frac{q^2-p^2}{2}+pq,X=\frac{p^2+q^2}{2}
\end{equation}
と得られる.


もっとも,元の問題に立ち戻れば自然数和を考えていたので,
\begin{equation}
A=p^2+2pq-q^2=(p+q)^2-2q^2\ge 0
\end{equation}
\begin{equation}
B=q^2+2pq-p^2=(p+q)^2-2p^2\ge 0
\end{equation}
という条件,すなわち
\begin{equation}
(\sqrt{2}-1)p\le q\le (\sqrt{2}+1)p
\end{equation}
p,qに課されている.例えばこの条件下にあるq=p+1ならば,
\begin{equation}
A=2p^2-1,B=2(p+1)^2-1,X=2p^2+2p+1,
\end{equation}
すなわち
\begin{equation}
a=p^2,b=(p+1)^2-1,x=p^2+p
\end{equation}
となって,冒頭の学生が発見した等式に戻る.

また(p,q)=(3,5)とすると(A,B,X)=(7,23,17)すなわち(a,b,x)=(4,11,8)となり,新たな等式
\begin{equation}
4+5+6+7+8=9+10+11
\end{equation}
が見つかる.あるいは(p,q)=(4,7)とすると(A,B,X)=(23,89,65),すなわち(a,b,x)=(12,44,32)となって,別の新たな等式
\begin{equation}
12+13+\cdots+32=33+34+\cdots+44
\end{equation}
が得られる.q=2pも条件下にあり,このときは
\begin{equation}
A=p^2,B=7p^2,X=5p^2
\end{equation}
と綺麗な形になり,奇数pに対して新たな系列
\begin{equation}
a=\frac{p^2+1}{2},b=\frac{7p^2-1}{2},x=\frac{5p^2-1}{2}
\end{equation}
が得られる.



使っているのは高校数学程度なんだが,結構面白い.
おやおや,なんだかミイラ取りがミイラになりかかっているぞ...

今年もガッサガッサしてみた

この時期は毎年,担当講義である統計とコンピュータにてガッサガッサする実験をする.
tokidoki.hatenablog.jp

今年で3年目となるこの実験の元ネタは,啓林館のページにあった授業実践だった.
www.shinko-keirin.co.jp

ただし,気をつけないと隣のクラスからうるさいと言われかねない.
(実際,今年は隣で講義をされていた先生が何事かとそっと教室を見に来た.)
f:id:okiraku894:20150601112752j:plain

3年目ともなると結構な数の実験になる.一人50回×3年間実験参加者数272=13600回だ.
そうして得られた結果が以下.
f:id:okiraku894:20170626121052p:plain
ここまでくると コイン>ピン>キャップ のこの順が入れ替わるなんてことは
とんでもなく低い確率になろう.
実際,10×標準偏差分の幅を取ってやっとコインとピンが入れ替わる可能性が出てくる.
そしてその確率は7.7\times10^{-23}程度だ.
そんな評価ができるのも中心極限定理のおかげなんだが,
やっぱり学生にちっとも伝わらないのがこの現象.
もちろん今年も正方形ダーツボードの実験はしてもらうし,
なかなか軽い動作でシミュレーションできるこの回のEXCEL資料は
我ながら気に入っているのでまた掲載するよ↓
f:id:okiraku894:20150607135853g:plain
f:id:okiraku894:20150607145228p:plain

今年はこれに加えて「上か下か」のランダムウォークもEXCEL上で見せようと思った.
f:id:okiraku894:20170626124705p:plain
5本程度,1000ステップ分を描くと↑のようになるが,この程度ではふーんで終わる.
もっと本数を描いてみよう.例えば100本↓
f:id:okiraku894:20170626124706p:plain
こうなってくると形,いわゆる放物線がうっすら見えてくるのではなかろうか.
そしてこれらの軌跡は中心0に近いほど密度が高くなっている.
いるのだけどこれが平面では見せられない.
ということで今度は OpenProcessing にて立体的にこれを見せられるようにしてみた↓

上部は密度を立体的に表示し,その足元の影がランダムウォークの軌跡を表している.
因みに,パスの数が5000以下のときのみランダムウォークを表示するようにした.
それ以上の本数になるとヒストグラムと軌跡両方の表示はどうやら重すぎるので.


ところで全く関係ないが,とある事情によりイベントが急遽行われた.
この場を借りて,あざっす!
f:id:okiraku894:20170623122352j:plain

私の身の回りのべき乗則(2)―名字編

今クールNHKが放映している名字番組
「人名探究バラエティー 日本人のおなまえっ!」ってのがある.
www4.nhk.or.jp
帰宅がてら見るのにちょうど良い番組なのだけど,
レア名字編ではハンコタワー1つ5000種の印鑑で人口の80%がカバーできるとか,
10万種の印鑑の在庫を持つハンコ屋さんがあってそれでも人口の98%だとか,
って話を聞いてると,「おお!冪乗則!」ってなって,早速検索.
直ぐに名字ランキングサイトが見つかる.
以下のグラフはこのページから参照.
myoji-yurai.net
幸いにもその名字におよそ何人いるか,のデータも合わせて書かれてあったので,
冪乗則が調べられる.さてやってみたのが↓
f:id:okiraku894:20170609153627p:plain
EXCELが吐き出す近似曲線がy\fallingdotseq4\cdot 10^7x^{-1.14}で決定係数R^2が0.98を超える.
やはりなかなか良い冪乗則になっている様子.
\log-\log plot が↓
f:id:okiraku894:20170609153628p:plain

それでも何だか二次関数っぽくもあったので,\log_{10}-\log_{10} plotでの二次近似を行うと
X=\log_{10} x,Y=\log_{10} yとして
 Y= -0.189X^2 - 0.0647X + 6.1599,決定係数 R^2 = 0.9997 となった.
因みに,人数が一人,すなわち Y=0 となるのが X\fallingdotseq 5.54 で,
すなわちランキングとしては x=10^X\fallingdotseq 347014 となって,
名字が約29万種あると言われているところからすると,ちょっと多めの評価になる.
もっとも果たして「一人」の状態までカウントして良いのやら,というところでもある.
f:id:okiraku894:20170609173723p:plain


日本の名字は分家制で名字の派生が起こったり同じ読みのまま漢字を変えたり,と
名字の生成/消滅の過程は集団遺伝の力学系のようでもあって,
名字の力学モデルなんてのも考えても面白いのかもしれない.

私の身の回りのべき乗則

それは気のせいなだけなのかもしれないが,
何だかパソコンの動作がモッサリしてきたなぁ,
と思ったときにClean upした後デフラグをかけて,
モッサリがなくすという作業を定期的にしている.

で,そのとき分析結果として各ファイルの断片化数が報告されるわけだが,
それを分断数が多い順に並べ替えると
いわゆる「べき乗則」に従っているように以前から感じていた.
f:id:okiraku894:20170530130948p:plain

で,とある機会に実際にグラフにしてみたところ,なかなか見事な結果になった↓
f:id:okiraku894:20170526180652p:plain
ただし,デフラグのログが取り出せなかったので手打ちでデータを書き写したのだが,
分断数6以下のものは大量すぎて作業から外した.
EXCELによって与えられる近似曲線はy=5557.4x^{-0.95}で,
世にいうZipf則が成り立っているように見える.

それでも図に騙されている気がしたので両側対数軸 \log y/\log x で描いたのが↓
f:id:okiraku894:20170526180653p:plain
なかなか良い感じだ.

実はこの記事を書こうと思ったのは他でもない,
予算が削減されつつある大学の電子ジャーナルを契約更新するか否かの判断材料として
ジャーナルごとのアクセス数のデータをもらったことに始まる.
早速アクセス数の多いジャーナル順にソートして描いてみたら見事なべき乗則↓
f:id:okiraku894:20170526180655p:plain
f:id:okiraku894:20170526180654p:plain
おやおや.

この結果に気を良くして,では世間で言われているweb siteのアクセス数についても
べき乗則が見えるのかと,例えばこのブログの1日あたりアクセス数を
多い順にソートして並べてEXCELで近似曲線まで描かせると
どうもべき乗則ではなく指数分布のほうが決定係数が大きくなる.
(まぁ,決定係数が8種類もあってEXCELがどれを採用しているのか,
 とかいった厄介な話はこの際考えない.)
f:id:okiraku894:20170529121850p:plain

あれれ,何が起こっているのか?と思い両対数グラフを描いたら
\log x=2.5 あたりから急に落ち込んでいた↓
f:id:okiraku894:20170529121851p:plain

ならば\log x=2.5に該当する20アクセス未満の日を切り落とせば良いのか,と描いてみたのが↓
y=253.8x^{-0.408}と今度はべき乗則のほうが決定係数が大きくなった.
f:id:okiraku894:20170529121852p:plain
まぁ,不都合な部分を切り落としたのだから,そりゃべき乗則のほうが良くなるでしょ,
ともいえるのだけど.
そういえば最初のデフラグデータも分断数が小さいものは調べなかったから,
順位の大きい先のほうは切り落としたことになっているわけで,
世にいう「べき乗則」は順位の比較的小さい部分で起こっている話,と理解するものなのかなぁ.

数理音楽の風景(1)-不協和度曲線

2012年度卒業生から始めた当ゼミでの数理音楽なる分野,
その風景を数回シリーズでまとめてみよう.

準周期系の研究において連分数にまつわる力学系を扱うことが多かったことと,
はるか以前から音楽の仕組みに数理的仕組みが見え隠れしていると感じていたので,
ずっと取り扱いたいとは思っていたことだった.
そもそも音階の作り自体,繰り込み的な要素が多分にあるし,
その辺の感覚は連分数力学系そのものに由来する.
愛知教育大学学術情報リポジトリ

さて,その数理音楽のゼミでの取り扱い開始は2012年度7代目の卒業研究から.
ドレミファそーする?―音律の数理と音階の構成―
不協和度曲線から始まって,連分数展開による音階構成の仕組みまで触れた卒論.
とあるページでもこの卒論が紹介されているのをしばらく前に発見した.
sites.google.com
この卒論,元になった藤沢+クックの論文*1
和音性の計算法と曲線の絵描き方 ―不協和度・緊張度・モダリティ―
の真似事から始まって,
不飽和度曲線で起こっている解析的な現象をちょっと数学的にきちんと書いてみた,
というものだった.
ただ,不協和度曲線の協和点でグラフが尖ってるといった特長があればいいものを,
そうもなっていなかったために,何らかの方法で協和点を解析的に特徴付けねば,
とでっち上げたのが第二種不協和度曲線だった.
tokidoki.hatenablog.jp

不協和度曲線の数理モデルのアイディアはこうだ.
まず,周波数の異なるsin波を同時に2つ鳴らし,
心理的な快不快の度合いを複数の被験者に採点してもらう.
その結果丁度半音階程度ずれたときが最も不快で,
そこを離れるに従い不快感は減っていったので,
この採点結果を適当な関数でモデル化する.それは例えば
f:id:okiraku894:20170129104409p:plain
であり,純音のドに対する別の純音の不協和度が以下.
f:id:okiraku894:20170128144954p:plain

しかし実際の楽器による楽音は複数の倍音が鳴っているため,
複合音の不協和度は純音の不協和度曲線をシフトして重ね合わせたものになるだろう,
つまり協和・不協和感覚に線形性が成り立っているだろう,という仮定を課すのである.
f:id:okiraku894:20170129105131p:plain

その結果,以下のような複合音に対する不協和度曲線が得られる.
例えば楽音のドに対する別の楽音の不協和度は↓
f:id:okiraku894:20170128144955p:plain
そして経験的に人々(正確には西洋音楽圏の人々*2)が協和すると感じる場所で
不協和度が下がっている,という尤もらしい説明ができる,というわけだ.
なお当たり前のことなのだが,この協和点の集合はFarey列にlog2したものに他ならない.

特に完全五度である3/2倍音は常にオクターブに次いで最も協和度が高い.
そのことによってピタゴラスは現在ピタゴラス音律と呼ばれる,
音階づくりのルールを定めたのだった.
そしてそこから連分数力学系が絡む素敵な世界が広がるのだが,その話は次回以降へ.

音律と音階の科学―ドレミ…はどのようにして生まれたか (ブルーバックス)

音律と音階の科学―ドレミ…はどのようにして生まれたか (ブルーバックス)

*1:更にその大元にはPlomp と Leveltによる音楽心理実験があるのだ. R. Plomp and J. M. Levelt "Tonal Consonance and Critical Bandwidth" Journal of the Acoustical Society of America, 1965.

*2:実際,ハーモニーという概念を持たない民族圏の人々は例えば半音差で二音鳴らした場合と完全五度で鳴らした場合とで快不快感に特段の違いは起こらなかった,という実験結果がある.単旋律文化圏であった他ならぬ明治以前の日本人もそうだったのではなかろうか.