2017年12月2日

映画 gifted/ギフテッド

Tom Flynn(脚本) 2017年 米国 

このブログの100作品目です!

ガウス積分

亡き姉の子で数学の才能を持つメアリーを育てているフランクは,母から彼女に英才教育を受けさせるよう勧められますが,彼女を普通に育てようとして奮闘する話です.

大学の講義室の黒板に次の内容が書かれてあり,メアリーが間違いを指摘します.
Problem
Show that $\displaystyle\int_{-\infty}^{\infty}e^{x^2/2\sigma^2}dx=\sqrt{2\pi}|\sigma|$
Hint: First show that
$\displaystyle\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{(x^2+y^2)/2\sigma^2}dxdy=2\pi\sigma^2$ 
メアリー「You forgot the negative sign on the exponent.(指数に-をつけるの忘れてる)」
シーモア「Mary, why didn't you say anything?(メアリー,どうして言わなかったの?)」
メアリー「Frank says I'm not supposed to correct older people. Nobody likes a smart-ass.(フランクが年上の間違いを正すなって.嫌われるから)」
最後の"smart-ass”は字幕にありませんでした.直訳すると「賢い尻」ですが,「知ったかぶり」というような意味です.

メアリーは最初の式の $\displaystyle\int_{-\infty}^{\infty}e^{x^2/2\sigma^2 }dx$ に-を付け足して $\displaystyle\int_{-\infty}^{\infty}e^{-x^2/2\sigma^2}dx$ に訂正していましたが,その下の式の中の$(x^2+y^2)$の前にも-が必要です.

ガウス積分
この証明すべき等式で$\sigma=\frac{1}{\sqrt{2}}$のとき,すなわち$$\displaystyle\int_{-\infty}^{\infty}e^{-x^2}dx=\sqrt{\pi}$$をガウス積分といい,関数$y=e^{-x^2}$とx軸とで挟まれる部分の面積が$\sqrt{\pi}$=1.77245であることを示しています.この証明は他のサイトでも見つかりますが,この映画に出てきた$\sigma$のついたままの等式の導出を確認してみましょう.

z=exp(-(x^2+y^2))
求める積分の値を $I=\displaystyle\int_{-\infty}^{\infty}e^{-x^2/2\sigma^2}dx$ とおき,$I=\sqrt{2\pi}|\sigma|$を示します.
\begin{align}
I^2&=\displaystyle\int_{-\infty}^{\infty}e^{-x^2/2\sigma^2}dx\int_{-\infty}^{\infty}e^{-x^2/2\sigma^2}dx\\
&=\displaystyle\int_{-\infty}^{\infty}e^{-x^2/2\sigma^2}dx\int_{-\infty}^{\infty}e^{-y^2/2\sigma^2}dy\\
&=\displaystyle\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{-(x^2+y^2)/2\sigma^2}dxdy
\end{align} このとき,$I^2$は右上図の盛り上がった曲面とxy平面で挟まれた部分の体積を表しています.ここで $x=r\cos\theta$,$y=r\sin\theta$ と置換すると,\begin{align}
I^2&=\displaystyle\int_{0}^{\infty}\int_{0}^{2\pi}e^{-r^2/2\sigma^2}rd\theta dr\\
&=2\pi\displaystyle\int_{0}^{\infty}e^{-r^2/2\sigma^2}rdr\\
&=2\pi \left[ -\sigma^2 e^{-r^2/2\sigma^2} \right] _0^{\infty}\\
&=2\pi\sigma^2
\end{align}よって,この等式を示すことができました.$$I=\displaystyle\int_{-\infty}^{\infty}e^{-x^2/2\sigma^2}dx$=\sqrt{2\pi}|\sigma|$$ところでこの式は,平均0で標準偏差$\sigma$の正規分布(ガウス分布)の$\sqrt{2\pi}|\sigma|$倍を表す式になっています.これを,平均μ,標準偏差σの正規分布を表す式に変形してみましょう.

平均0で標準偏差1の標準正規分布
まず,確率密度関数の-∞から∞までの積分は1にならなければいけないので,両辺を$\sqrt{2\pi}|\sigma|$で割ります.$$\frac{1}{\sqrt{2\pi}|\sigma|}\displaystyle\int_{-\infty}^{\infty}e^{-x^2/2\sigma^2}dx=1$$これを平均μになるように平行移動すると,$$\frac{1}{\sqrt{2\pi}|\sigma|}\displaystyle\int_{-\infty}^{\infty}e^{-(x-\mu)^2/2\sigma^2}dx=1$$となり,平均μで標準偏差σの正規分布を表す式に変形できました.

2017年11月15日

ドラマ 相棒 season16 第4話

2017年11月8日放送 テレビ朝日

Φ オイラーの関数 zero slash 漸化式 二重積分 確率


数学の才能を持つコンビニ店員の若い男(健次郎)が,大学入試問題漏洩に関与させられたあげく殺されてしまった事件の謎を,杉下右京が解き明かすという話です.

漸化式の解法

健次郎がニセ学生となって講義を聴き,終了後,教授に質問をする場面です.
「この証明,別の方法もあると思うんです.僕ちょっと考えたんですけど,まずですね,両辺を$(-1)^{n+1}$で割るんです.そうするとこれ,ただの等差数列になるんです.」 
この場面では数学に関する台詞はこれだけでした.別解を考えついて正しいかどうかを確認しに来たものと思われます.これは漸化式の解法のひとつです.注意深く録画したビデオを見てみると,板書には,漸化式から一般項を推測し,それを数学的帰納法で証明する解答が書かれてありました.その板書の一部は次の通りです.
$a_1=1, \space \space a_{n+1}=-a_n+(-1)^n$ 
$a_1=1, \space \space a_2=-2, \space \space a_3=3, \space \space \cdot \cdot \cdot \cdot, \space \space a_n=n \cdot (-1)^{n+1}$
実際にこの漸化式の両辺を$(-1)^{n+1}$で割ると,$\frac{a_{n+1}}{(-1)^{n+1}}=\frac{a_n}{(-1)^{n}}-1$となります.すると数列$\lbrace{\frac{a_n}{(-1)^{n}}}\rbrace$は,初項${\frac{a_1}{(-1)}}=-a_1=-1$,公差$-1$の等差数列となり,$\frac{a_n}{(-1)^{n}}=-1+(n-1)\cdot(-1)=-n$なので,一般項は$a_n=-n\cdot(-1)^n=n\cdot(-1)^{n+1}$と分かります.

受験生にはよく知られた解法ですが,仮に未習だったとしても,教授が「驚きました.彼は,我々には考えつかないユニークな発想を持っていたんです.」と言っていたのは少しオーバーな気がしました.

積分順序の交換

健次郎が二重積分の問題を解くシーンがありました.
$I=\displaystyle \int_0^{2a} dx \displaystyle \int_{\sqrt{2ax-x^2}}^{\sqrt{4ax}}f(x,y)dy$ 
の積分順序を交換せよ
これは,$I=\displaystyle \int\int_D f(x,y)dydx$(Dは右図の放物線と円と直線x=2aで囲まれた領域)とも書けます.先にyで積分してx軸に垂直な切り口を求め,それをxで積分して,領域Dと曲面z=f(x,y)に挟まれた部分の体積を求める式ですが,積分順序の交換だけなので$f(x,y)$は何でもOKです.

この手順を交換するのですから,先にxで積分してy軸に垂直な切り口を求め,それをyで積分する式をつくるとこのようになります.
$\displaystyle \int_0^a\int_\frac{y^2}{4a}^{a-\sqrt{a^2-y^2}}f(x,y) dx dy+\displaystyle \int_0^a\int_{a+\sqrt{a^2-y^2}}^{2a} f(x,y) dx dy+\displaystyle \int_a^{2\sqrt{2}a}\int_\frac{y^2}{4a}^{2a}f(x,y) dx dy$

道順の確率

「選択肢に正解がなかった」という台詞が何度もありました.どんな問題だったのでしょう.録画したビデオからその問題が分かりました.
[問題4] xy平面上の点Aが,原点(0, 0)から点(n, n)(nは3以上の自然数)まで以下のルールで動く.
1. 点(k, l)にあり,k<n,l<nならば,点(k+1, l)か点(k, l+1)に1/2の確率で動く.
2. 点(n, l)(l<n)なら点(n, l+1)へ,点(k, n)(k<n)なら点(k+1, n)に確率1で動く.
このとき点Aが(n, 2)を通過する確率P(n, 2)は次のうちどれか.
①$(n+1)(\frac{1}{2})^{n+1}$ ②$(n^2+1)(\frac{1}{2})^{n+2}$ ③$(n^2+5n)(\frac{1}{2})^{n+3}$ ④$\frac{1}{(n^2-4n+1)^n}$ ⑤$\frac{n}{2^n+1}$
調べてみたらこの出典は,東工大入試問題の1987年第5問の改題と分かりました.元の問題には選択肢がないので難しいですが,この改題なら,n=3の時の解,すなわちP(3, 2)を求め,選択肢にn=3を代入してこの値になるものを選べば良いわけです.(3, 3)へ行くには(3, 2)または(2, 3)を通るしかないので,対称性よりP(3, 2)=P(2, 3)=1/2となります.ところがn=3を選択肢①~⑤のどれに代入しても1/2になりません.なので「選択肢に正解がなかった」というわけです.

[Reference]
Math-Station 東工大入試研究
http://kubojie.net/titech.html

2017年11月5日

NEWS 太陽系外からの彗星発見か

2017年10月26日 日本経済新聞他 
双曲線 楕円
[ワシントン=共同] 国際天文学連合小惑星センターは25日、太陽系外から飛んできた可能性がある彗星(すいせい)を発見したと発表した。確認されれば、恒星間の軌道を飛行する初の 「恒星間彗星」となる。
 彗星は、米ハワイ大の望遠鏡が発見した「C/2017U1」。現在は地球の軌道と火星の軌道の間を飛んでいるとみられる。世界各地の天文台が30回以上観測した結果、太陽系の外からやってきた可能性があることが分かった。
 同センターは「大きな双曲線軌道を描いているようだ」としており、太陽に近づくのは1度きりで戻ってこないとみられる。オーストラリアのメディアは「太陽から25光年と比較的近くにある恒星、こと座のベガのある方向から来たようにみえる」と報じた。
 彗星は太陽の周りを回る楕円軌道を描いたり、惑星からの力を受けて太陽系外にはじき飛ばされたりするものがこれまで発見されている。
彗星や惑星は楕円軌道で,太陽はその楕円の持つ2つの焦点のひとつだと習ったので,彗星に双曲線軌道のものが存在するのは意外でした.楕円は「2点(焦点)からの距離の和が一定である点の軌跡」,双曲線は「2点(焦点)からの距離の差が一定である点の軌跡」です.いちばん上の図では軌道が急カーブしている少し内側に太陽があるので,そこが焦点ということになります.

まとめて二次曲線と呼ばれる円,楕円,放物線,双曲線は,円錐を切断する方向を変えるとこれらの形の切り口が得られるので,円錐曲線とも呼ばれています.同じ仲間といえますが,かなり形は違いますね.円とどれだけ近い形か,かけ離れた形かを示す値eを離心率といい,円の離心率は0,楕円は0<e<1,放物線はe=1,双曲線は円からかけ離れた形をしているので1<eとなります.地球の軌道離心率は0.0167なので,ほとんど円に近いということが分かります.一方,あの周期75年といわれるハレー彗星の軌道離心率は0.967なので,かなり放物線に近い,細長い楕円ということができます.

双曲線軌道なら,右上の図でいえばこの彗星は遥か遠くからやってきて,太陽(焦点)の近くの頂点 (a, 0) で急カーブして方向を変え,また遥か彼方に去って行くというイメージです.もし仮に双曲線の式(標準形は$\frac{x^2}{a^2}-\frac{y^2}{b^2}=1$)が $x^2-y^2=1$ならa=b=1なので,頂点は(1, 0),焦点$(\sqrt {a^2+b^2}, 0)$は$(\sqrt {2}, 0)$になります.双曲線の最も易しい例は,中1で学習する反比例のグラフ$y=\frac {k}{x}$ですね.$y=\frac {1}{x}$のとき,頂点は(1, 1)ですから,この軌道なら太陽は$(\sqrt {2}, \sqrt {2})$という位置にあるということになります.

<余談1> 回転軸に平行でない直線が回転すると,一葉双曲面ができます.兵庫県神戸市のポートタワーがそのようなつくりになっています.一度近くまで行って,確かめてみてください.

<余談2> 私が小中高と育った兵庫県伊丹市の市章は双曲線に似ています.一度検索して見つけてみてください.

[Refference]
For the first time, astronomers are tracking a distant visitor streaking through our solar system
http://www.sciencemag.org/news/2017/10/first-time-astronomers-are-tracking-distant-visitor-streaking-through-our-solar-system
軌道離心率
https://ja.wikipedia.org/wiki/%E8%BB%8C%E9%81%93%E9%9B%A2%E5%BF%83%E7%8E%87

2017年10月18日

小説 ラメルノエリキサ

渡辺 優 2016年 集英社

アシンメトリー

自分が傷つけられたら必ず復讐をするという女子高校生小峰りなの話です.
懐かしい。 駅からバスで15分ほど先にある、絶妙に不便な立地の図書館を前にして、私はまずそう感じた。小学生の時、遠足で訪れたことがある。それから中学の時にも一度、何かの行事で来ていたはずだ。質感の違う素材がモザイクのように組み合わされた、濃いグレイを基調とした建物。表のプレートによると、築52年にもなるそうだけれど、それにしては古めかしい感じはしない。正面ゲートからまっすぐ見据えると、ビルの建ち並ぶ街の中心から外れ、賛沢な土地の使い方をした3階建ての本館が、アシンメトリーな曲線を描いているのがわかる。
アシンメトリー(asymmetry)は非対称という意味なので,「アシンメトリーな曲線」とは対称でない曲線という意味です.アシンメトリーな部分のある建物の図書館を探してみたら,近いイメージなのがこの写真でした.この小説とは全く関係ありません(笑).発音記号はは æsímətri or eisímətri の2通りあります.ネイティブの発音ではアシリ(またはエイシリ)と聞こえます.

ここで数学の話が出てきたわけではないのですが,アシンメトリーは数学用語のひとつでもあるので少し見てみましょう.

[確率分布のグラフにおける対称性と非対称性]

左右対称なグラフといえば正規分布ですが,その峰が左右にずれた場合,アシンメトリーな曲線になります.分布の非対称性を示す指標を歪度(わいど skewness)といい,峰が右にずれて左が低い場合(J型)の歪度は負,左右対称であれば歪度は0,峰が左にずれて右が低い場合(L型)の歪度は正になります.

因みに,正規分布の場合,平均(mean),中央値(median),最頻値(mode)は同じ値になりますが,他の場合はそれぞれの値がずれていきます.正規分布から少しずれた非対称な分布では,平均と最頻値との差は平均と中央値との差の約3倍になることが知られていて,カール・ピアソンの経験則と呼ばれています.
    Mean–Mode≒3(Mean–Median)
    または 3Median≒2Mean+Mode

[二項関係における対称関係と非対称関係]

二項関係の最も簡単な例は2つの数の関係です.aからbへRという関係があるとき,aRbと表します.例えば3>2などの関係です.

「aからbへの関係が成り立つならば,bからaへの関係も成り立つ」という場合,対称関係といいます.すなわち,$$\ aRb\;\Rightarrow bRa$$例えばRが=という関係のとき,a=bならばb=aも成立するので対称になります.

「aからbへの関係が成り立つならば,bからaへの関係が成り立たない」という場合,非対称関係といいます.すなわち,$$\ aRb\;\Rightarrow \lnot (bRa)$$例えばRが>という関係のとき,a>bならばb>aは成立しないので非対称になります.

小説のタイトルの意味を言いたくてしょうがないのですが,ネタバレになるのでここでは控えておきます(笑).

[Reference]
Asymmetric relation
https://en.wikiedia.org/wiki/Asymmetric_relation
Sample Mean
http://mathworld.wolfram.com/SampleMean.html
二項関係
https://ja.wikipedia.org/wiki/%E4%BA%8C%E9%A0%85%E9%96%A2%E4%BF%82

2017年8月16日

小説 偏差値70の野球部 レベル3 守備理論編

松尾清貴 2012年 小学館文庫

中学時代にピッチャーで全国大会準優勝した新(あらた)真之介(しんのすけ)が,なぜか東大合格者数全国1位の超進学校に入学し,甲子園を目指すという話です.
 
「角運動量保存則を利用した等速円運動打法」を考案し,「強いのは野球のセオリーを作っているチームだ」という才女のヒカルさんが,打球の落下地点を知るための方法として,重力,摩擦抗力,揚力,粘性率,レイノルズ数などを考慮した流体力学の理論を話しました.
ヒカルさん「変数は,打球の水平面に対しての角度と速度と角速度.一応,抗力係数は0.6に設定してる.レイノルズ数がはっきりすればいいんだけど,信頼のおける実験結果が得られなかったから.計算上は105程度になるから,どっちみちナビエ-ストークスの方程式は使えない.(後略)」
ヒカルさんが画面を覗き込むようにしながらマウスをクリックすると,新たに浮き出たドットが,さらに複数の相対的にきれいな二次関数を描き出した.頂点を挟んで右側,つまり下り曲線が上り曲線よりも急角度だった.
この数式を期待したのですが,このままでは複雑すぎるということで,その後,以下のように外的な力を無視し,単なる二次関数でボールの到達距離と時間を予測していました.話の流れからすると残念です.しかもその中の数式に1か所ミスがありました(読んだのは第3刷なので,その後訂正されているかも知れません).
ヒカルさん「初速度v0,水平面との角度θの放物運動する物体は,水平方向x鉛直方向yの二次元座標において,x=v0cosθ・t,y=v0sinθ・t-gtで表されるよね.ボールが到達する地点では y=0 だから,まずy式にそれを代入する.到達時間t2はもちろん0ではないから,式 t2=2v0sinθ/g が得られるでしょ.水平到達距離をx軸上のD点として,先に求めたt解をx式に代入する.だから,具体的には,D=v0cosθ・tを,v0cosθ・2v0sinθ/g に変換すると,これは v02sin2θ/g になるね.」
式中の$g$は重力加速度9.8m/秒2,$t$は時間(秒)を表します.はじめの式 $$ y=v_0 \sinθ \cdot t-gt^2 $$ は,正しくは次式になりますね.$$y=v_0\sin\theta \cdot t-\frac{1}{2}gt^2\tag{1}$$この場合,速度は $$(\frac{dx}{dt}, \frac{dy}{dt})=(v_0 \cosθ, v_0 \sinθ-gt)$$ なので,これを$t$で積分すると,位置は次式になります(整式の積分なので高校数学Ⅱの知識でOK).$$(x, y)=(v_0 \cosθ・t, v_0 \sinθ・t-\frac{1}{2}gt^2)$$ この$(x, y)$から$t$を消去すると次式になります.$$y=\tan\theta \cdot x-\frac{g}{2v_0^2\cos^2\theta}\cdot x^2\tag{2}$$例え空気抵抗などを無視したとしても,相手がボールを打ってフライを上げた瞬間に,この計算をして外野手に知らせるなんて,超人的な能力ですね(笑).

初速度や打出角度が変わればどうなるかをGeoGebraで作ってみました.青のグラフは横軸が時間$t$で式(1)を,紫のグラフは横軸が水平飛行距離$x$で式(2)を表しています.

なお,空気抵抗を加味した放物運動は,映画「ST赤と白の捜査ファイル」ですでに紹介しています.

2017年8月11日

随筆 墨汁一滴

正岡子規(1867-1902) 岩波文庫 岩波書店
1927年12月第1刷 1998年1月第35刷  

正岡子規は「柿くへば鐘が鳴るなり法隆寺」という俳句で有名です.36才で亡くなる前の年の1901年1月から7月まで、正岡子規が記者として勤めていた当時の新聞『日本』に連載されていた随筆の一説ですが,この文章は東京大学予備門(のちの第一高等学校、戦後は東京大学教養学部)に入学した18才の頃を回顧しています.
しかし余の最も困つたのは英語の科でなくて数学の科であつた。この時数学の先生は隈本(有尚)先生であつて数学の時間には英語より外の語は使はれぬといふ制規であつた。数学の説明を英語でやる位の事は格別むつかしい事でもないのであるが余にはそれが非常にむつかしい。つまり数学と英語と二つの敵を一時に引き受けたからたまらない、とうとう学年試験の結果幾何学の点が足らないで落第した。(六月十四日) 
余が落第したのは幾何学に落第したといふよりもむしろ英語に落第したといふ方が適当であらう。それは幾何学の初にあるコンヴアース、オツポジトなどといふ事を英語で言ふのが余には出来なんだのでそのほか二行三行のセンテンスは暗記する事も容易でなかつた位に英語が分らなかつた。落第してからは二度目の復習であるから初のやうにない、よほど分りやすい。コンヴアースやオツポジトを英語でしやべる位は無造作に出来るやうになつたが、惜しい事にはこの時の先生はもう隈本先生ではなく、日本語づくめの平凡な先生であつた。しかしこの落第のために幾何学の初歩が心に会得せられ、従つてこの幾何学の初歩に非常に趣味を感ずるやうになり、それにつづいては、数学は非常に下手でかつ無知識であるけれど試験さへなくば理論を聞くのも面白いであらうといふ考を今に持つて居る。これは隈本先生の御蔭かも知れない。(六月十五日)
幾何学でコンヴアース(converse),オツポジト(opposite)が登場するものを推測してみましょう.Converseは,条件文P⇒Q(PならばQ)に対してその逆,Q⇒P(QならばP)を意味します.Oppositeは「向かい側の」とか「反対側の」という意味でよく使われます.

■平行四辺形の性質「平行四辺形の2組の対角(opposite angles)は等しい」の逆(converse)は「2組の対角(opposite angles)が等しい四角形は平行四辺形である」

■平行四辺形の性質「平行四辺形の2組の対辺(opposite sides)は等しい」の逆(converse)は「2組の対辺(opposite sides)が等しい四角形は平行四辺形である」

■三平方の定理「直角三角形の斜辺の2乗は他の2辺の2乗の和に等しい」の逆(converse)は「三角形で,ひとつの辺の2乗と他の辺の2乗の和が等しいとき,最初の辺の対角(opposite angle)は直角になる」

■二等辺三角形の定理「二等辺三角形の底角は等しい」の逆(converse)は「三角形の2角が等しいとき,それらの対辺(opposite sides)も等しい」

■三角形の辺と角の大小「三角形の2辺が異なるとき,長い方の辺の対角(opposite angle)は短い方の辺の対角より大きい」の逆(converse)は「三角形の2角が異なるとき,大きい方の角の対辺(opposite side)は小さい方の角の対辺より長い」

以上はいずれも今の日本の中学程度の内容ですが,もともとユークリッドの「原論」に載っていたものです.正岡子規が存在していた明治時代初期の頃,幾何学はユークリッドの「原論」が主に学習されていたようなので,その中の初歩といえば,このような内容だったと思われます.

[Reference]
明治前期の日本において教えられ,学ばれた幾何(数学史の研究)

2017年7月28日

映画 Hidden Figures (邦題「ドリーム」)本編

2016年 米国 20th Century Fox

Euler's Method

日本では2017年9月公開予定の映画です.主役のキャサリン・ジョンソンは、宇宙船を軌道から戻すための計算にオイラーの方法(Euler's Method)を使います.これは常微分方程式 ordinary differential equations (ODEs)の厳密解(解析解)が求められないときに近似解(数値解)を求める方法の一つです.


解こうとする微分方程式が $y'=f(t, y(t))$,初期値が$y(t_0)=y_0$のとき,まずtを幅hで区切って$t_0=0, t_1=h, t_2=2h, t_3=3h,$ …とし,最初はt=0における接線で近似,その後は,$y_{n+1}=y_n+h・f(t_n,y_n)$として,2点$(t_n, y_n)$, $(t_{n+1},y_{n+1})$を結ぶ直線で近似していくという方法です.

簡単な例として,厳密解(解析解)が$y=e^t$になる微分方程式$y'=y$を見てみましょう.初期値は$y_0=1$になります.この例では$f(t_n,y_n)=y_n$なので,$y_{n+1}=y_n+h・y_n$として$(t_n,y_n)$を次々に求めて行き,折れ線グラフを作っていくというイメージです.

$y_1=y_0+h\cdot y_0=1+h\cdot 1=1+h$
$y_2=y_1+h\cdot y_1=(1+h)+h\cdot(1+h)=(1+h)^2$
$y_3+h\cdot y_3=(1+h)^2+h\cdot(1+h)^2=(1+h)^3$

例えば,h=0.5なら,$y_1=1.5$,$y_2=2.25$,$y_3=3.375$,$y_4=5.0625$,…となって,右図のようになります.
(オレンジ色がEuler's methodによる近似解,青色が$y=e^x$)

hの値(分割の幅)が小さいほどより良い近似になりますが,その分,計算は煩雑になります.hの値を変えたらどうなるかわかるものをGeoGebraで作ってみましたので,試してみてください.

因みにEuler's Methodは,世界中最も多くの国で普及している高校カリキュラム「国際バカロレアディプロマプログラム(IB Diploma Program)」のMathematics Higher Level(主に理系向き)のOption科目"Calculus"に登場します。

余談ですが,日本語の方のWikipediaで,1か所間違いを見つけたので訂正しておきました.これでWikipediaの記事を訂正したのは4~5回目ぐらいになのですが,以前に訂正した内容を記録してなくて思い出せないので,今後履歴が残るように今回初めてWikipediaにアカウント登録してみました.

<reference>
Euler method
https://en.wikipedia.org/wiki/Euler_method
オイラー法
https://ja.wikipedia.org/wiki/%E3%82%AA%E3%82%A4%E3%83%A9%E3%83%BC%E6%B3%95

2017年7月22日

漫画 Q.E.D. 証明終了

 加藤元浩 1997-2014年 講談社

以前から数学の話題が頻繁に登場する漫画であることを知っていたのに,なぜかここに書かなかったという作品のひとつです.久しぶりに第1巻を含むいくつかを読んで,ようやく書く気になりました.

理工系大学の世界最高峰である米国のMIT(マサチューセッツ工科大学)を15歳で卒業したが,普通の高校生活を体験したいと,日本の高校に入学してきた燈馬想が,同級生の体育会系水原可奈の協力のもと,次々と難事件を解決していくという漫画です.Q.E.D.はラテン語でQuod Erat Demonstrandum.タイトル通り「証明終了」という意味です.

第1巻/第2話「銀の瞳」

数学の苦手な水原可奈に燈馬想が数学を教える場面です.
燈馬想「だから交点が1になる場合を考えて判別式が0になる値を求めるんです.」
水原可奈「じぇんじぇんわかりましぇん」
問題が掲載されていなかったので推理してみましょう.このころはまだ2人とも高校1年生なので,高校数学Ⅰの問題なら,放物線と直線の共有点が1つ,すなわち接するときの式中の定数を求める問題でしょう.そうだとすれば,「交点が1になる場合を考えて」よりも「共有点が1つになる場合だから」という表現のほうが適切ですね.

[簡単な例題] 放物線$y=x^2$と直線$y=ax-3$が接するとき,定数$a$の値を求めよ.(正解は±2√3)

第7巻/第1話「Serial John Due」
ロキ「じゃア,カメーネフの死は"π"に符合してるとして…,劉の方は?」
燈馬想「自然数eだよ」
ロキ「e? あのオイラーの数か?」
水原可奈「eって?」
燈馬想は「自然数」,ロキは「オイラーの数」と言っていますが,eは「自然対数の底」,または「ネイピア数」と呼ばれることが多いです.「自然数」だと正の整数(または負でない整数)の意味を持つnatural numberと誤解しやすくなります.また,eという文字で表すのはオイラーが最初だったので「オイラーの数」でもいいのですが,他に「オイラーの定数(Euler's constant)」というのがあるので,混同しないようにしたいですね.
燈馬想「円周率πは円周を円の直径で割った数.自然数eはネイピアの考案した対数という概念から生まれた数.そして虚数iは2乗して-1になる想像上の数として生まれました.この3つの記号はそれぞれ人類の歴史の中でなんの関連もなく生み出されたものです.ところが大数学者オイラーはこの3つの数に一つの関係を見つけ出した.
$e^{πi}=-1$
これがオイラーの公式,人類の数学史上最も美しい式と呼ばれるものです.」
ここでもeを「自然数」と呼んでいますが,さらにオイラーの公式(Euler's formula)$$e^{θi}=\cosθ+i\sinθ$$に$θ=π$を代入して得られるオイラーの等式(Euler's identity)$$e^{πi}=-1$$を「オイラーの公式」と呼んでいるのが気になりますね.

第33巻/第2話「推理小説家殺人事件」
燈馬想「『グラフで囲まれた部分の面積を求めろ』ってことなので,まず2つの交点のxの値を出して範囲を決めます.yが同値になるxってことだから…」
燈馬想が水原可奈に数学を教えている場面ですが,台詞はこれだけで図もありません.これは積分の問題でしょう.中には公式を使って交点のxの値を出さずに解く方法があります.以前,映画「容疑者xの献身」で述べました.

第38巻/第2話「十七」

xのn乗根がn=1から4まで紹介されていましたが,n=4のときの解(1の4乗根)は,$x^4-1=0$を解いて$x=±1$,$±i$になります.間違ってn=4のところにn=5のときの解(1の5乗根)が書かれてあります.と思ったらよく見ると値も違っていました.正しくは次のようになります.他のサイトでいくつも見つかりますから,確認しみてください.$$\frac{(-1+\sqrt{5})±i\sqrt{10+2\sqrt{5}}}{4}, \quad \frac{(-1-\sqrt{5})±i\sqrt{10-2\sqrt{5}}}{4}$$

第44巻/第2話「Question!」

トップのタイトルページの数式です.右上の図より,円$x^2+y^2=1$と,(-1,0)を通る直線$y=m(x+1)$の,(-1,0)でない方の交点の座標を求める問題のようです.横書きの数式なのに,板書が右から始まるのが変ですね.話は右上の図から始まって,右下の連立させた方程式を変形する過程まではいいのですが,そのあとの左の計算式(燈馬想に隠れて見えにくいですが)が少しおかしいですね.右下の式を整理すると次式になります.
$(1+m^2)x^2+2m^2x+(m^2-1)=0$
この左辺をを因数分解するのに,解のひとつがx=-1とわかっているので,x+1がひとつの因数になります.その割り算は右図のようになるはずなので,$m^2>1$や,$x√$というのは写し間違いだと思われます.

他の巻にも以下のように数学の話題が頻出します.数学の話ではなくても十分楽しめる知的な漫画なので,どなたにも強くお勧めします.
  9巻 ケーニヒスベルグの橋
13巻 クラインの壺
15巻 デデキントの切断
20巻 カントールの無限集合
23巻 リーマン予想
29巻 ポアンカレ予想
44巻 フェルマー予想 ゴールドバッハ予想 ABC予想
47巻 P≠NP問題

2017年6月11日

テレビ番組 とくダネ! 中国の過酷すぎる受験対策

2017年6月9日放送 約940万人が挑む一発勝負
中国の全国大学統一入試は「高考(ガオカオ)」と呼ばれ、約940万人が受験に臨み,2800校中超一流と呼ばれる80校を目指します.15校が志望でき、ランクごとに5校ずつ選べます.学力のみの一発勝負で,試験の成績によって入学校が決まり,成績トップ者は新聞一面で紹介されます.大学ごとの入試は行なわれません.東京大学の倍率約3倍に対して、北京大学の倍率は150倍にもなります.(番組情報より)
今年広東省で出題された数学の問題の前半に出てきた比較的易しい問題がひとつ紹介されていました.少し分かりにくい問題文でしたが,要するに太極と呼ばれる図形の黒い部分の面積の,外接する正方形の面積に対する割合を求める問題です.正解の選択肢は
A 1/4  B π/8  C 1/2  D π/4

東大の建築学科卒で数学に関する著作もある司会の女性タレントがこの問題の正解をすぐに答えていました.「どう考えても1/4じゃないから…,Bですよね」 つまり計算はせずに,半分よりは小さく,1/4よりは大きいからその間のπ/8を正解にしたようです.

一応確認してみましょう.円の半径をrとすると,黒い部分は白い部分と全く同じ形なので円の面積の半分で$\pi r^2/2$,正方形の面積は$(2r)^2=4r^2$,前者を後者で割ると$\pi /8$になります.

テレビの情報番組やクイズ番組などで数学の問題が出るときは,このようにすぐに正解が得られる問題しか出てきませんね.特にクイズ番組では,ただ覚えている知識を答えるだけの問題が多いので,物足りなく感じます.

他の問題も気になったので探してみました.2017広東省数学問題解答です.上の問題の原文も見つかりました.「理科数学」は「理系数学」という意味です.せっかく見つけたので少し時間のかかる最後の記述式問題も紹介しておきます.
[理科数学 2] (上の問題)如图,正方形ABCD内的图形来自中国古代的太极图正方形內切圓中的黑色部分和白色部分关干正方形的中心成中心対称.在正方形内随机取一点,剣此点取自黒色部分的概率是.
[理科数学 23] 己知函数$f(x)=-x^2+ax+4$, $g(x)=|x+1|+|x-1|$.
(1) 当$a=1$时,求不等式$f(x)\ge g(x)$的解集.
(2) 若不等式$f(x)\ge g(x)$的解集包含$[-1, 1]$,求$a$的取値范围.
[理系数学 23](和訳)関数$f(x)=-x^2+ax+4$,$g(x)=|x+1|+|x-1|$がある.
(1) $a=1$のとき,不等式$f(x)\ge g(x)$を解け.
(2) 不等式$f(x)\ge g(x)$の解が$[-1, 1]$に含まれるとき,$a$の取り得る値の範囲を求めよ.

正解は  (1) -1≦x≦(-1+√17)/2   (2) -1≦a≦1

[Reference]
2017年广东高考真题及答案解析
http://www.gaokao.com/e/20170607/59375d1f5f355.shtml

2017年6月4日

NEWS 大学新テスト記述式問題例 

2017年5月16日に発表された独立行政法人大学入試センター「大学入学共通テスト(仮称)」記述式問題のモデル問題例が,翌日の新聞各紙に掲載されました.「どんな問題だろう」と解いてみた人も多いのではないでしょうか.その中の4問目(数学の2問目),銅像が最もよく見える位置を考察する,すなわち銅像の頭Aと足元Bまでを見込む角が最大になる視点Pの位置を考えるという問題が少し気になりました.問(1)は具体的な数字を与えられて∠APBを求めるもので特に問題を感じませんでしたが,この後の展開が気になりました.
問(2) 銅像を見込む角が最大となるときの,見る人の足元の位置を「ベストスポット」と呼ぶことにする.この「ベストスポット」について,太郎さんは次のように考えた.
[太郎さんの考え] 3点A,B,Pを通る円の半径をRとすると,ABの長さは常に一定であることから,∠APBが鋭角ならば,∠APBが最大となるのは,Rが最小のときである.
ということは,Rがどんなときに最小になるかをこれから考えていくんだなと思ったら,次はこんな問いでした.
(i) ∠APBが鋭角であることを確かめる方法を,△APBの3辺の長さAB,AP,BPについての式を用いて説明せよ.
ここでわざわざ「∠APBが鋭角である」ことを確かめる必要があるのでしょうか.銅像の足元を視点より低くして,よっぽど近寄らないと見込む角は鈍角にはなりません.銅像の足元Bは視点Pより高いので,∠APBが鋭角であることは明らかです.この状況を図に描いたら,100人中100人は∠APBが鋭角になるでしょう.
(ii) [太郎さんの考え] が正しいことは,sin∠APB,AB,R を用いたある関係式と,「∠APBが鋭角のとき,∠APBが大きくなるほどsin∠APB の値は大きくなる」ことからわかる.その関係式を答えよ.
これだけ誘導されていたら正弦定理とすぐに分かりますね.sin∠APB=AB/(2R)なので,「Rが最小」⇔ 「sin∠APBが最大」⇔ 「∠APBが最大」といえます.
(iii) 二人は [太郎さんの考え] について先生に相談したところ,Rが最小になるのは,3点A,B,Pを含む平面上において,3点A,B,Pを通る円と点Pを通り直線ABに垂直な直線が接するときであることを教えてもらった.この考え方に基づいて,目の高さが1.5mの花子さんが,高さ6.5mの台座の上に乗せた高さ4mの銅像を見る場合の最小R,最大∠APB,ベストスポットの位置を求めよ.
ようやくRが最小になるときを考えるのかなと思ったら,いきなり先生に相談です.自分たちで考える問題なのに「先生に教えてもらったこと」を前提にして話を次に進めています.この「先生に教えてもらったこと」を理解せずに次を考えるのは気持ち悪くないのでしょうか(この解説はこちらのサイトにあります).

この解説の別解で書かれているように, 「先生に教えてもらったこと」は,数学Ⅰまでの知識なら円周角を考えれば説明できます.他の方法で説明しようとすれば,正接の加法定理と相加平均・相乗平均の関係か方べきの定理,または微分を利用してもできますが,数学Ⅰの範囲を超えてしまいます.

日本の試験ではまだ使えませんが,グラフ電卓やそれに類似するソフト・アプリを使えば,数学Ⅰまでの知識でもこのことを確認することはできます.

まずA(0,9), B(0,5), P(x,0)とし,3点A,B,Pを通る円の中心をC(c,7)として,cをxの関数で表します.AC=BC=CPなので,
c^2+2^2=(x-c)^2+7^2
整理すると
c=(x^2+45)/(2x)
あとはこの関数を未習であっても,グラフ電卓等でグラフを描かせて最小値を表示させれば,ベストスポットの位置の近似解6.7が得られます.

国際バカロレア,米国のAPやSAT,英国のA Levelなどでは,グラフ電卓を使える試験と使えない試験の併用が当たり前のように実施されています.

視点をいろいろ変えたらどうなるかをGeogebraで作ってみました.点Pを動かしてみてください.

2017年3月5日

映画 君に届け


椎名軽穂(原作) 2010年 東宝

北海道の高校が舞台です.数学の授業の最後に問題を出しておいてノートに解答を書くよう指示して先に退室した先生から,黒沼爽子がノートを集めて持ってくるよう指示された…のではないかと思われます(こんな先生,あまりいないと思いますが…).
黒沼爽子「数学のノートを集めますので,持って来てください」 
かなり板書が読みづらかったのですが,こんな内容でした.
三角形の性質
<問題>
∠C>∠Bである△ABCで,∠Aの2等分線と辺BCとの交点をPとし,点Aから辺BCに下した垂線をAQとすると,2∠PAQ=∠C-∠Bとなることを証明せよ.
(以下の解答部分は自分で解いた後にドラッグして見てください)

<解答>
APは∠Aの2等分線より,∠BAP=∠CAP
∠APB=∠CAP+∠C … ①
∠APC=∠BAP+∠B … ②
①②より,∠APB-∠APC=∠C-∠B
仮定より,∠APB=180°-∠APC
     ∠APC=90°-∠PAQ
よって
∠C-∠B=∠APB-∠APC
    =180°-2∠APC
    =180°-2(90°-∠PAQ)
    =2∠PAQ
以上は板書の左半分に書かれてあった問題と解答ですが,右半分にあった生徒に出題されたであろう問題は読み取れませんでした.

上の問題はあまり見たことがなかったので,検索してみたらこちらにひとつ見つかりました.

この結果がどんな三角形でも成り立つことを確認できるものをGeogebraで作ってみました.頂点を動かしてみてください.

2016年8月23日

ドラマ 受験のシンデレラ

和田秀樹(原作) 2016年 NHK

分数 三角関数 積分 実数 直円錐

人気塾講師五十嵐透が貧しい高校生遠藤真紀を東大に合格させようとする話です.1/2+1/3=2/5と答えた真紀に対し,お好み焼きを6等分して1/2+1/3=5/6を説明していました.

第3話
オープニングクレジット(Opening credits)にいくつかの問題が一瞬だけ表示されます.

①sinθ-cosθ=14となるとき,sinθcosθの値を求めよ.
この問題は,第4話から「sinθ-cosθ=1/4となるとき」に訂正されていました.高校数学Ⅱの三角関数の合成を習うと,sinθ-cosθ=√2sin(θ-π/4)と変形でき,この値の絶対値は√2より小さいことが分かりますから,sinθ-cosθ=14はあり得ません.たぶん,視聴者からの指摘があって訂正したのではないかと思われます.訂正された問題の正解は15/32になります.解いてみてください.
第4話

他にもかすかに問題が読み取れました.

②曲線y=x^2+x+2-2xと直線y=2mx+2で囲まれる部分の面積が最小値となる定数mを求めよ.
この曲線の式は同類項xと-2xがまとめられていないのが既におかしいですね.まとめると曲線y=x^2-x+2となりますが,このまま解くと面積の最小値が0(このときm=-1/2)になってしまって,間違いではないですが少しおかしな問題になります.

③実数t,x,y,zが下の条件を満たす.x,y,zのうち,2つが等しい値となるときのtの値を求めよ.
x+y+z=t
x^2+y^2+z^2=t^2-t+10
x^3+y^3+z^3=t^3-3t^2+4t+12
いろいろやってみましたが,すっきりした値が出ないので,zをyに置き換えてWolframAlphaに
solve x+2y=t,x^2+2y^2=t^2-t+10,x^3+2y^3=t^3-3t^2+4t+12
と入力してみたら,有効数字6桁のtの値が5つ出てきました.問題のどこかが間違っているような気がします.

④直円錐の頂点をP,底面の直径の両端をA,B,線分APの中点をQとする.このとき,点BからQに至る最短距離を求めよ.
これは具体的な値が与えれていないので,文字だけで計算することになります.母線APの長さをR,底面の半径をrとして,展開図の側面の扇形を考え,線分BQの長さを求めればいいですね.側面の扇形の中心角をθとすれば,θ=2πr/Rなので,△PQBの∠P=θ/2=πr/Rとなり,PB=R, PQ=R/2なので,余弦定理より
$BQ^2=R^2+\left(\frac{R}{2}\right)^2-2\cdot R\cdot\frac{R}{2}\cos\frac{\pi r}{R}$
よって,BQの長さは次式になります.
$BQ=R\sqrt{\frac{5}{4}-\cos\frac{\pi r}{R}}=\frac{R}{2}\sqrt{5-4\cos\frac{\pi r}{R}}$

このドラマの原作者は東大卒の精神科医で,大学受験対策本や啓蒙書も多数執筆されていて,中には「なぜ数学が得意な人がエグゼクティブになるのか」とか「数学は暗記だ!」など,ちょっと気になる本も出されています.

<2017年1月3日追記>
最終回(2016年8月28日放送)の中で2017年1月15日センター試験「数学Ⅰ・数学A」の問題が少しだけ映っていました.もちろんフィクションなので本物ではありません.実際は,第1問から第3問が必答で,第4問~第6問から2問選択ですが,ドラマでは必答にすべき2次関数の問題が選択になっていました.読み取った問題を掲載しますので解いてみてください.

(ドラマの中の問題)
数学Ⅰ・数学A
第5問(選択問題)(配点 20)
2次関数
      y=-x^2+2x+1  …①
のグラフの頂点の座標は(ア,イ)である.また
      y=f(x)
はxの2次関数で,そのグラフは,①のグラフをx軸方向にp,y軸方向にqだけ平行移動したものであるとする.
(1) 下の(ウ)(オ)には次の⓪~④のうちから当てはまるものを一つずつ選べ.ただし,同じものを繰り返し選んでもよい.
      ⓪> ①< ②≧ ③≦ ④≠
2≦x≦4におけるf(x)の最大値がf(2)になるようなpの値の範囲は
      p(ウ)(エ)
であり,最小値がf(2)になるようなpの値の範囲は
      p(オ)(カ)
である.
(2) 2次不等式f(x)>0の解が-1≦x≦5になるのは
      p=(キ) q=(ク)
のときである.
(3) 2次不等式f(x)>1の解が-1≦x≦5になるのは
      p=(ケ) q=(コ)
のときである.

二次関数の最大値/最小値を求めたり,二次不等式を解くだけなら易しいのですが,このように,解が先に分かっていて元の関数の式を決定する問題は少し難しくなります.ドラマの中での主人公の解答は(カ)と(コ)が間違っていましたが,全科目900点満点中728得点だったので良かったですね.

(正解)
(1) (ア,イ)(1, 2)   (ウエ)≦1   (オカ)≧2
(2) (キ)1   (ク)7   (ケ)1   (コ)8

GeoGebraで正解を確認できるものを作ってみました.pとqの値を変えてみてください.
https://www.geogebra.org/m/CT4gKJvV

2016年8月22日

映画 Hidden Figures (予告編)

2016年 20th Century Fox

等辺 台形 二等辺 正四面体 4次方程式

まだ人種差別法が存在していた1960年代、NASA(アメリカ航空宇宙局)でマーキュリー計画に数学で貢献した実在のアフリカ系アメリカ人女性Katherine Johnsonたちを描いた映画です.2016年8月現在,まだ公開されていないので,予告編を見ただけでこれを書いています.
Katherine "equilateral, trapezoid, isosceles, tetrahedlon,....."
Teacher "I've never seen a mind like the one your daughter has."
最初の4語は,等辺,台形,二等辺,正四面体です.equilateral triangleなら正三角形,equilateral quadrilateralなら等辺四角形すなわちrhombus(菱形),isosceles triangleなら二等辺三角形になります.その次が,小学生の時の先生が両親にKatherineの才能について話した台詞です.飛び級していたのでしょうか,Katherineが高校生のクラスで,2次式×2次式=0という4次方程式を解く場面がありました.

$(x^2+6x-7)(2x^2-5x-3)=0$
$(x+7)(x-1)(2x+1)(x-3)=0$
$x=-7,\quad  1,\quad  -1/2,\quad  3$

日本語のタイトルは何だろうと思って探しましたが見つからないので機械翻訳してみると,「隠された数」「隠された図」などが出てきたのですが,figureには人物の意味もあり,栄光の陰で支えた人たちを描いているので,「陰で活躍した人たち」がいいのではないかと思います(検索したら中国語のサイトには「隐藏人物」と出てきました).

2016年8月9日

アニメ 山田君と7人の魔女 第1話

吉河美希(原作) 2014年 講談社

ベクトル 微分積分

身体を入れ替ることができる不良生徒山田竜と優等生白石うらら,他に数名が超常現象研究部をつくって,同じ高校にいる7人の魔女探しをするという話です.冒頭でいきなり数学Ⅱの小テストの答案が現れ,できなかった山田竜を先生が叱っている場面から始まります.
先生「君は学校を何だと思っているのかね.この私立の名門,朱雀高校開校以来初めてだよ.毎日のように遅刻・早退を繰り返し,あちこちで喧嘩三昧,揚句,授業中に居眠りだと!? 山田竜!」
この小テスト,変なところががいくつかありました.
・小テストなのにやたら問題数が多い.
・「解答」と書くべきところがすべて「回答」となっている.
・どの解答欄も極端に狭い.
・なぜか横に緑茶の入った湯呑がある.
・数学Ⅱなのに数学Bや数学Ⅲの範囲の問題が出ている.
・しかも小テストにしてはけっこう計算の面倒な問題がある.例えば, 「次の曲線や直線で囲まれた部分の面積を求めよ」
(1) y=(x-e)logx, y=0
(2) (3) 略
(4) y=sinx, y=sin2x (0≦x≦π)
解いてみたらこうなりました.
(1) -1/4e^2+e-1/4(≒0.621)
(4) 5/2(=2.5)
(数学Ⅲを既習の人は確かめてみてください)
それにしても数学Ⅱの小テストなのですから,生徒の立場からすれば,出題範囲を越境しないでほしいですね.


2016年8月8日

小説 ラプラスの魔女

東野圭吾 2015年 角川書店

素数 ナビエ–ストークス方程式 ラプラスの悪魔
桐宮玲「もし,この世に存在するすべての原子の現在位置と運動量を把握する知性が存在するならば,その存在は,物理学を用いることでこれらの原子の時間的変化を計算できるだろうから,未来の状態を完全に予知できる――」
「ラプラスは,このような仮説を立てました.その存在のことは後年,ラプラスの悪魔と呼ばれるようになります.」
「ラプラスの悪魔」と呼ばれる知性を持つという主人公が女の子なので,このタイトルになったようです.数学では,ある物体の「現在位置と運動量を把握する」ということは「任意の点で導関数が満たす微分方程式を作る」ということを意味し,この微分方程式を解くことで,ある関数が求められ,それによってその物体の過去や未来の位置が計算できるようになります.

ナビエ–ストークス方程式は,流体(液体や気体など)の運動を表す非常に複雑な2階非線型偏微分方程式です. 「2階」は第2次導関数まで式の中にあるという意味で,「非線型」は「線型」でない,すなわち次の条件を満たさないという意味です.

        加法性: 任意のx,yに対してf(x+y)=f(x)+f(y)
        斉次性: 任意のx,αに対してf(αx)=αf(x)       

「偏微分方程式」は2変数以上の関数において,ある1つの変数について微分して得られる導関数を含んだ方程式という意味です.例えば,y=f(x)のとき,$\frac{dy}{dx}$はyをxで微分したもの,z=f(x,y)のとき,$\frac{∂z}{∂x}$はzを(x以外は定数と見做して)xだけで微分したものという意味です(∂は「ラウンド」とか「デー」などいくつか読み方があります).ナビエ–ストークス方程式は,難しすぎて未だに解決されず,クレイ研究所が「解けたら100万ドル」を贈るという懸賞付きの問題のひとつです.

ところで,$\frac{dy}{dx}$をもう一度微分すると,$\frac{d^2y}{dx^2}$になります.なぜ$\frac{dy^2}{dx^2}$ではなく,分子は$d^2y$で分母が$dx^2$と書くのかというと,$\frac{dy}{dx}$は$\frac{d}{dx}y$とも書けるので,$\frac{dy}{dx}$をもう一度微分すると,$\frac{d}{dx}\left ( \frac{dy}{dx}\right )$となり,分子はdが2個でyが1個,分母はdxが2個になるからです.(注)$(dx)^2$は$dx^2$とも書きます.$(AB)^2$を$AB^2$と書くのと同様です.

2016年7月31日

小説 未来方程式 -fate equation-

神世 希 2011年 講談社

英数字 三次元空間座標 緯度 経度

予知能力があるとされる少女初音未来を中心に起こる数々の事件が描かれた話です.
皇(すめらぎ)警部 「違う……そもそも8桁のパスワードが適当に打って当たるハズがないだろう.全部で218兆3401億0558万4896通りあるんだぞ」
波多野巡査 「そんなにあるんすか!?」
英数字(alphanumeric)8桁ということは,英字(alphabetic)26文字と数字(numeric)10文字で合計36文字が8桁なので,重複順列 36Π8=36^8=2兆8211億0990万7456 (通り) のはずです.しかし値が違うので,逆に218兆3401億0558万4896の8乗根を計算してみると62になりました.62=26×2+10なので,ここでは大文字と小文字を区別しているということが分かりました.それにしてもこんな数字を覚えているはずがないのにすらすらと言えるのが不思議ですね.
所長「20XX年に打ち上げられた人工衛星"おまわり参号(仮称)"は,明日―――世界時において6月6日正午には,地球を中心とする3次元空間上の座標,(x,y,z)=(171,1789,6841).緯度経度では,北緯75.3度,東経146.1度の,高度およそ7000メートルに存在することになる.これも立派な未来情報の一つだ」
この座標から緯度,経度,高度を計算してみましょう.原点をO,点(171,0,0)をA,点(171,1789,0)をB,点(171,1789,6841)をPとします.

[緯度] まず北緯のほうを確認してみます.$$OB=\sqrt{171^2+1789^2}≒1797$$なので,北緯は次の値になり,上の台詞の値と一致しました.$$\arctan\frac{BP}{OB}=\arctan\frac{6841}{1797}≒75.3°$$

[経度] OBx軸とのなす角は,$$\arctan\frac{AB}{OA}=\arctan\frac{1789}{171}≒84.5°$$になります.東経146.1度なので,146.1-84.5=61.6° が子午線と交わる赤道円の半径とx軸とのなす角になります.すなわち,x軸が東経61.6度の経線と交わっているとすれば,この座標の位置は東経146.1度になります.

[高度] まずOPの長さ(地球の中心から人工衛星までの距離)を求めます.$$OP=\sqrt{171^2+1789^2+6841^2}≒7073$$ここから地球の赤道半径約6378kmを引くと,7073-6378=695 (km) なので,高度は上の台詞の「およそ7000メートル」ではなく,およそ700kmということになります.飛行機の平均高度が約10000m (10km) なので,人工衛星は飛行機より約70倍高いところにあることがわかります.

因みに人工衛星の高度,速度,周期を調べてみたところ,重力加速度をg,地球の赤道半径をRとし,高度h(km)で回っている円軌道の場合,その速度v(km/sec)は次の式で計算できることがわかりました.
$$v=\sqrt{\frac{gR^2}{R+h}}$$
従って人工衛星の速度は、高度だけで決まり,高いところほど遅いということになります.軌道の円周2π(R+h)をこの速度vで割れば周期 t(秒)が計算できます.
$$t=\frac{2\pi(R+h)}{v}$$

高度        速度        周期
200km    7.8km/sec   1時間28分
500km     7.6km/sec    1時間34分
1,000km    7.4km/sec    1時間45分
36,000km    3.1km/sec    23時間56分(静止軌道)
380,000km    約1.0km/sec    約27日間

問題
上に出てきた高度700kmの人工衛星の速度,周期を求めよ.

(答)約7.5km/sec, 約1時間39分

参考
「人工衛星の高度と速度」宇宙航空研究開発機構(JAXA)

2016年7月23日

小説 数学的帰納の殺人

草上 仁 2009年 ハヤカワ・ミステリワールド

ピタゴラス 完全数 フィボナッチ数列 嘘つき村 ガロア カルダン・グリル ステヴィン 懸垂線(カテナリー) ワイルズ メルセンヌ素数 フェルマーの最終定理 谷山・志村予想 楕円関数 円周率 モンティ・ホール問題 虚点

フォト・ジャーナリストの女性と大学教授の男性が,カルト教団による殺人連鎖の謎を追う話です.この殺人連鎖を数学的帰納法に例えてこのタイトルにしたようです.

2段組みで400頁を越える長編の随所に数学の話題が多数登場します.中には話の流れに乗らない余分な感じのするところもありましたが,謎解きは面白く,最後も少し驚かされたりして,けっこう楽しめました.

作品中に少し気になった点やもっと説明のほしいところがありました.
「三角形の下にぶら下がる球はロゴでは半円形に見えるが,実際には懸垂線つまりカテナリー曲線を描く.式にすると$y = a \cosh(x/a)$だね.この双曲線はバランス的に中立だから…(以下略)」
双曲線は反比例などの2つの同じ形をした曲線をいいます.懸垂線は双曲線関数(hyperbolic function)のひとつですが,双曲線関数は双曲線の媒介変数表示に使うからそう呼ぶのであって,形は双曲線ではありません.円関数(三角関数)は円の媒介変数表示に使いますが,形は円ではないというのと同様です.従ってここは「この双曲線は」ではなく,「この曲線は」でいいと思います.

懸垂線がなぜこの式になるかこちらにまとめてみました。
「二つの円は,常に四点で交わる.体験的にわかるのは二点までだ.しかし体系を普遍化していくと,そちらのほうが特殊な例であることがわかる.より一般的には,見えない二点――虚点を想定したほうが,多くの場合に成り立つのだ.」
2つの図形を表す連立方程式が実数解を持たない場合は$xy$平面上では交わりませんが,複素平面上で虚数解を座標とする点(虚点)で交わると見なすことができます.例えば,$y=x^2$と$y=2x-5$は実数解を持たないので$xy$平面上では交わりませんが,複素平面上の虚点$(1±2i, -3±4i)$で交わると考えることができます.

2つの円の場合は$x^2$と$y^2$の係数が等しいので,連立させると$x^2$と$y^2$の項が消えてしまい,交点が2つしか出てきません.従って複素平面上で考慮しても4つの交点はないのですが,さらに無限遠点を含む射影平面上で考えれば,平行な直線でも1点で交わり,どんな2つの円も4点で交わると見なすことができます.これをベズーの定理(Bézout's theorem)といいます.

2016年5月3日

映画 落ちこぼれの天使たち Stand and Deliver


1987年 米国

部分積分

荒れた高校の生徒たちと,そこへ赴任してきた教師がだんだん信頼し合っていくという話です.1987年の映画ですが,部分積分が登場することを最近知りました.生徒が黒板に書いた数式はこれです.
$\int{x^2 \sin x}dx$
$u=\sin x$           $dv=x^2dx$
$du=\cos xdx$     $v=\frac{1}{3}x^3$
$\frac{1}{3}x^3\sin x-\int{\frac{1}{3}x^3\cos x}dx$

 生徒は$du=\sin xdx$,$v=x^2$にすれば良かったのですが,逆にしてしまいました.教師が書いた解き方は,日本ではあまり見ない方法です.


これは表を使う方法(tabuler method)で解いています.この方法を使うと,繰り返しの必要な部分積分の計算が速くできます.

表をじっくり見てみましょう.左の列は$x^2$を次々に微分したもの,中央の列は$\sin x$を次々に積分したもの,右の列は符号です.折れ線でつながるものを与式の右辺に書き出せば正解が得られます.$$\int{x^2 \sin x}dx=-x^2 \cos x+2x \sin x+2 \cos x+C$$


2016年1月9日

漫画 和算に恋した少女 第2巻

天元術 算額

▼ 「はああ…… 天元術かあ……」
 「ごらん律,これが算額だ.自分で考えた問題を額にして神社に奉納するんだよ.父さんの算額はこれだ.答えが書いてないのは『遺題』といってね,『誰か解ける人はいるかい?』って訪ねているんだよ」
 「じゃあ,あたしが解く!」
 「うん,そうしたらおまえの答えを算額にしてこの隣に掲げなさい.」

江戸時代の数学を和算といいます.天元術とは一言でいうと方程式で問題を解く方法です.これまでに解いた算額の問題は三平方の定理を用いて円の半径などを求める問題が多かったので,この「父さんの算額」の問題も同じように解けると思ったらこれはかなりの難問でした.

他サイトで類題があり,①デカルトの円定理を使った解答と②反転法を使った解答が見つかったのですが,ここではピタゴラスの定理を使って解く方法を考えてみました.

計算を簡単にするため,図のように大円$O$の中心を原点,半径を$6$として,円$C$の中心を$(3, 0)$,円$C_0$の中心$(a_0, b_0)=(-3, 0)$とすれば,円$C_0$の半径 $r_0=3$となります.それ以外の円${C_n}$の中心を$(a_n, b_n)$,半径を$r_n$としましょう.
$C_{n-1}C_n$を斜辺とする直角三角形で三平方の定理より$$(a_{n-1}-a_n)^2+(b_{n-1}-b_n)^2=(r_{n-1}+r_n)^2\tag{1-n}$$ $OC_n$を斜辺とする直角三角形で三平方の定理より$$a_n^2+b_n^2=(6-r_n)^2\tag{2}$$ $CC_n$を斜辺とする直角三角形で三平方の定理より$$(a_n-3)^2+b_n^2=(3+r_n)^2\tag{3}$$ 以上の式を使って$r_n$をいくつか求めます

■まず$r_1$を求めましょう.この場だけ円$C_1$の中心$(a_1, b_1)=(a,b)$,半径$r_1=r$として計算します.式(1-n)でn=1とすると$$(-3-a)^2+(0-b)^2=(3+r)^2\tag{1-1}$$ 式(2)より$$a^2+b^2=(6-r)^2\tag{2}$$ 式(3)より$$(a-3)^2+b^2=(3+r)^2\tag{3}$$ まず(2)(3)から$a$と$b$を$r$で表すと$$a=6-3r,\quad b^2=24r-8r^2\tag{4}$$ (4)を式(1-1)に代入して$r$についての方程式を解くと$r=2$.これを(4)に代入して$(a, b)=(0, 4)$.すなわち,円$C_1$の半径と中心は次式になります.$$r_1=2,\quad (a_1, b_1)=(0, 4)$$ ■次に$r_2$を求めましょう.計算を簡単にするためにまた同じ文字を使います.この場だけ円$C_2$の中心$(a_2, b_2)=(a,b)$,半径$r_2=r$として計算します.式(1-n)でn=2とすると$$(0-a)^2+(4-b)^2=(2+r)^2\tag{1-2}$$ 同様に(4)を(1-2)に代入して$r$についての方程式を解くと$r=1, r=3$となりますが,$r$は$2$より小さいので$r=1$.式(1-2)が$r_2$の前後の$r$を求められる式になっているので,もうひとつの解である$3$は$r_0$の値になっています.$r=1$を(4)に代入して$(a, b)=(3, 4)$.すなわち,円$C_2$の半径と中心は次式になります.$$r_2=1,\quad (a_2, b_2)=(3, 4)$$ ■次に$r_3$を求めましょう.計算を簡単にするためにまた同じ文字を使います.この場だけ円$C_3$の中心$(a_3, b_3)=(a,b)$,半径$r_3=r$として計算します.式(1-n)でn=3とすると$$(3-a)^2+(4-b)^2=(1+r)^2\tag{1-3}$$ 同様に(4)を(1-3)に代入して$r$についての方程式を解きましょう.
$(3r-3)^2+(4-\sqrt{24r-8r^2})^2=(r+2)^2$
$9r^2-18r+9+16-8\sqrt{24r-8r^2}+24r-8r^2=r^2+2r+1$
$4r+24=8\sqrt{24r-8r^2}$
$r+6=2\sqrt{24r-8r^2}$
$r^2+12r+36=4(24r-8r^2)$
$33r^2-84r+36=0$
$11r^2-28r+12=0$
$(11r-6)(r-2)=0$
$r=\frac{6}{11}, r=2$
となりますが,$r$は$1$より小さいので$r=\frac{6}{11}$.$r=\frac{6}{11}$を(4)に代入して$(a, b)=\left(\frac{48}{11}, \frac{36}{11}\right)$.すなわち,円$C_3$の半径と中心は次式になります.$$r_3=\frac{6}{11},\quad (a_3, b_3)=\left(\frac{48}{11}, \frac{36}{11}\right)$$ ■次に$r_4$を求めましょう.計算を簡単にするためにまた同じ文字を使います.この場だけ円$C_4$の中心$(a_4, b_4)=(a,b)$,半径$r_4=r$として計算します.式(1-n)でn=4とすると$$\left(\frac{48}{11}-a\right)^2+\left(\frac{36}{11}-b\right)^2=\left(\frac{6}{11}+r\right)^2\tag{1-4}$$ 同様に(4)を(1-4)に代入して$r$についての方程式を解くと(計算略)$r=\frac{1}{3}, r=1$となりますが,$r$は$\frac{6}{11}$より小さいので$r=\frac{1}{3}$.すなわち,円$C_4$の半径は $r_4=\frac{1}{3}$ になります.

ここで求めた$r_n$を並べてみると,$2, 1, \frac{6}{11}, \frac{1}{3}, \cdot\cdot\cdot$
分子を6にすると,$\frac{6}{3}, \frac{6}{6}, \frac{6}{11}, \frac{6}{18}, \cdot\cdot\cdot$
分母は,$3, 6, 11, 18, \cdot\cdot\cdot$
この階差数列は,$3, 5, 7, \cdot\cdot\cdot$となり,その第k項は2k+1と推測できるので,分母の第n項は$$3+\sum_{k=1}^{n-1}(2k+1)=n^2+2$$ 従って円$C_n$の半径は次式になると推測できます.$$r_n=\frac{6}{n^2+2}\tag{5}$$ あとは数学的帰納法で$$r_{n+1}=\frac{6}{(n+1)^2+2}\tag{6}$$が成り立つことを示せばOKです.式(4)の添え字をnにして式(1-n)に代入して整理すると次式を得ます.$$12(r_{n+1}-r_n)^2+3r_n^2r_{n+1}^2-4r_nr_{n+1}(r_n+r_{n+1})=0$$ これに(5)を代入して整理すると(6)を得られて終了なのですが,この計算は大変なので,WolframAlphaにしてもらいました.($x=r_n, y=r_{n+1}$として計算しています)

これで式(5)が正しいことが分かったので,これを(2)と(3)に代入すると円$C_n$の中心の座標は次式になります.$$(a_n, b_n)=\left(\frac{6(n^2-1)}{n^2+2},\quad \frac{12n}{n^2+2}\right)\tag{7}$$
さらに式(5)も式(7)もnに0と負の整数を代入しても成り立ちますから,円Cの周りのすべての円の中心と半径が求められたことになります.確かに$n\rightarrow\pm\infty$のとき,中心の座標は$(6,0)$に,半径は$0$に近づいていくことが分かります.

<2016年12月4日追記>

上図の半径3の円CとC0と半径6の大円Oとで囲まれる部分をArbelos(アルベロス)と呼ぶことが分かりました.古代ギリシアで使われていたといわれる「靴屋のナイフ(Arbelos)」がこの形と似ているのでこういう名前になったそうです.