投稿

1月, 2020の投稿を表示しています

万有引力の一次元問題

万有引力の働く二物体の運動を調べる。 衝突するまでに掛かる時間や、途中の運動の様子を運動方程式から解析的に求める。 質量m,Mの質点がt=0で距離R0だけ離れて静止していたとする。この二物体には互いに万有引力のみが働き、その他の力は働かないものとする。また、万有引力定数をGとする。 この二物体は直線上で運動するから、位置をそれぞれx,Xとし、x(0)=R0,X(0)=0,x˙(0)=X˙(0)=0とすれば運動方程式は、 {mx¨=GMm(xX)2MX¨=GMm(xX)2 となる。(但し、ドットは時間微分を表す。) 二式を重心運動と相対運動に分けると、 ddt(mx˙+MX˙)=0d2dt2(xX)=G(M+m)(xX)2 となる。 第一式を初期条件に注意して0からtで定積分すると、 m(x˙(t)x˙(0))+M(X˙(t)X˙(0))=0mx˙(t)+MX˙(t)=0 であり、更に同じ区間で定積分すると、 m(x(t)x(0))+M(X(t)X(0))=0mx(t)+MX(t)=mR0 を得る。 第二式でR:=xXと定義すると、R(0)=R0,R˙(0)=0であり、 \begin{align} \ddot{R}=-\dfrac{G(M+m)}{...

ポッホハマー記号と多項式

無限級数を超幾何級数で表すには、nの多項式をポッホハマー記号を用いて書き直す手続きが必要である。 今回は、ガンマ関数を経由して多項式をポッホハマー記号で表す方法を紹介する。 ポッホハマー(Pochhammer)記号 aC,nZに対し、 (a)n:=Γ(a+n)Γ(a)で定める。 n+aをポッホハマー記号で表す。 ガンマ関数の性質 zΓ(z)=Γ(z+1) を用いて多項式をうまくポッホハマー記号だけで表したい。そこで分母分子にΓ(n+a)を掛けてみると、 n+a=(n+a)Γ(n+a)Γ(n+a)=Γ(n+a+1)Γ(n+a)=Γ(n+a+1)Γ(a+1)Γ(a)Γ(n+a)Γ(a+1)Γ(a)=(a+1)n(a)nΓ(a+1)Γ(a)=a(a+1)n(a)n を得る。 これを用いると、無限級数を超幾何級数で表せて、そこから一般化出来たり、様々な変換公式を利用することが出来ることがある。 Tweet

log²(1-x)の冪級数展開

表題の展開係数は調和数を用いて表せることを紹介する。 東京大学2005年前期大問1と本質的には同じ問題である。 f(x):=log2(1x)と定め、原点を中心とした冪級数を求める。収束性などを無視すれば形式的に、 f(x)=n=0f(n)(0)n!xn であるから、f(n)(0)が求まればよい。 まず、nN2f(n)(x)=2(n1)!(1x)n(Hn1log(1x))()となることを数学的帰納法で示す。 但し、Hnは調和数で、 Hn:=k=1n1kである。 f(x)を順に微分していくと、 f(1)(x)=2log(1x)1x f(2)(x)=2(1x)2(1log(1x)) であり、n=2で成立。(f(1)(0)=0も分かる) Nnで成立しているとする。 すると、 \begin{align} &f^{(n+1)}(x)\ &=\dfrac{d}{dx}\left\{\dfrac{2(n-1)!}{(1-x)^n}\left(H_{n-1}-\log(1-x)\right)\right\}\ &=\dfrac{2(n-1)!}{(1-x)^{2n}}\left\{\dfrac{1}{1-x}(1-x)^n+n(1-x)^{n-1}(H_{n-1}-\log(1-x)\right\}\ &=\dfrac{2(n-1)!}{(1-x)^{...

ジョルダン標準形(1) : ジョルダン細胞と冪零行列

このシリーズでは、ジョルダン標準形について扱う。 今回はジョルダン細胞を定義し、その累乗を求めるための冪零行列について紹介する。 ジョルダン細胞 λCを用いてn正方行列のジョルダン細胞を以下で定義する。 Jn(λ):=[λδi,j+δi+1,j]=λIn+Jn(0) 但し、Inn次単位行列、δi,jはクロネッカーのデルタである。 定理1.1 任意のmNに対して、Jnm(0)=[δi+m,j]である。() 証明 m=1のとき、定義そのものからJn1(0)=Jn(0)=[δi+1,j]であるから成立。 Nm()が成立しているとすると、 Jnm+1(0)=Jn1(0)Jnm(0)=[δi+1,j][δi+m,j]=[l=1nδi+1,lδl+m,j] ここで、和に0でない項があるとしたら、i+1=l, l+m=jでなければならないので、i+m+1=jの項のみ残る。 したがってこの和はクロネッカーのデルタを使ってδi+m+1,jと書けるので、m+1でも成立し、数学的帰納法により()は示された。 冪零行列 n次正方行列Nに、mNがあって、 Nm=O が成立するとき、Nは冪...