パスワードを忘れた? アカウント作成

tetsuyaさんのトモダチの日記。 あなたのタレコミがスラドを支えます。さぁ、タレこめ!

13583848 journal
日記

akiraaniの日記: 本日のねためも

日記 by akiraani

「漫画村」のロゴ入り漫画、Amazon Kindleストアで無断販売される 集英社「至急対応する」(ITmedia)
 ……まあ、そういう底の浅いことをするやつも出てくるよな。問題は、やらかしたやつをきっちり捕まえてとっちめられる買って事なんだけど、Amazonはその辺ちゃんと保証できるだけの仕事はしてるのか? とりあえず、続報に注目かな。

NTTグループ4社、海賊サイトブロッキング実施を正式発表 「法制度が整備されるまでの短期的な緊急措置として」(ITmedia)
 うーん、ブロッキングするまでもなく対象のサイトが壊滅しているというのが茶番をステキに演出してるな。この反応は、ここまで既定路線で、水面下でいろいろ動いてたんだろうかね。

中学校の性教育で大論争、東京都議vs教育現場それぞれの言い分(ASCII.jp)

筆者は古賀都議に質問をし、回答を得た。内容を要約すれば以下の3点となる。

(1)「性教育は、各々の子どもの成長段階に応じて行うべきもので、中学生の授業で、生徒に対し一律に『性交』だけでなく、『避妊』『人工妊娠中絶』を教えるのは適切とは言えない。

 いやー、なに寝言言ってるんですかねこの人……。避妊教えるな、はいくらなんでもないだろう。

13583604 journal
日記

taggaの日記: 意識が朦朧としてるから、拡散方程式から正規分布を出す

日記 by tagga

;; 書類が……、というより、その前の調整でガリガリと……。

∂u/∂t = D ∂^2u/∂x^2, u(0, x)=0 when x!=0, ∫(-∞,∞)u(0,x)dx = 1 っていうのを解きたい。

要するにフーリエ変換するのだけど、 u≈Σ(k=-∞,∞)exp(i 2πk/T x)φ(t, 2πk/T) * 2π/T * 1/2π と書いてみる。

x で微分すると、 ∂u/∂x≈Σ(k=-∞,∞)(i 2πk/T)exp(i 2πk/T x)φ(t, 2πk/T) * 2π/T * 1/2π.

もう1度 x で微分すると、 ∂^2u/∂x^2≈Σ(k=-∞,∞)-(2πk/T)^2exp(i 2πk/T x)φ(t, 2πk/T) * 2π/T * 1/2π.

t で微分すると、 ∂u/∂t≈Σ(k=-∞,∞)(i 2πk/T)exp(i 2πk/T x)∂φ(t, 2πk/T)/∂t * 2π/T * 1/2π.

微分方程式から、∂φ(t, 2πk/T)/∂t ≈ - D(2πk/T)^2 * φ(t, 2πk/T).

2πk/T=ωとして、T→∞にすると∂φ(t, ω)/∂t = - Dω^2 * φ(t, ω).

ところで、のっけの式を t=0のところで、-T/2からT/2まで∫すると、 ∫(-T/2, T/2) u(0, x)dx ≈ φ(0, 0).

ということで、初期条件から、φ(0, 0)=1 なので、φ(t, ω)=exp(- Dω^2 t).

のっけの式に入れて、T→∞にすると、 u= 1/(2π) ∫(-∞,∞) exp(i ω x)exp(- Dω^2 t) dω.

あとは指数部分を平方完成して、地道に計算すると u = 1/(√(2π)√(2Dt)) * exp(- x^2 / (2 * 2Dt)).

ということで、分散 2Dt の正規分布.

とりあえず、正規分布の出し方でメジャーなやつは説明したはず。 次からストレスで何な状態になったら、何を書くのだろう……。

13582939 journal
Firefox

alpの日記: 4/24 通勤音楽

日記 by alp
  1. 4/24 通勤音楽
    • Ravel, Violinsonate #1 & #2, Tzigane, Berceuse Sur Le Nom De Gabriel Fauré, Ibragimova Tiberghien
    • 同, Klavierkoncert für linke hand D-Dur, Klavierkoncert G-Dur, Pierre-Laurent Aimard Pierre Boulez Cleveland Orchestra
  2. 先週に Pale Moon 27.9.0 が来ましたが、動きおかしいな。かなり重いのと、描画が乱れるサイト (Google だ) がある。ハングアップも既に数回起きています。Security fix はない模様なのでまだ更新していない人は様子見た方がいいかも。
13582915 journal
日記

taggaの日記: 気分がはれないので、雑な中心極限定理

日記 by tagga

中心極限定理で正規分布を出すには、 キュムラント母関数を使うのが簡単。 だけど、2種類あるキュムラント母関数を出すには、 特性関数か積率母関数を出す必要がある。 という具合で、先が長い。

ということで、特性関数。

密度関数を波の合成と思えば、
f(x)≈Σ[k=-∞, ∞] 1/(2π) * exp(-ix* 2πk/T) * φ(2πk/T) * 2π/T.

逆向きにφ(2πk/T) ≈ ∫(-T/2, T/2) exp(ix * 2πk/T) f(x) dx.

そこでフーリエ変換とその逆変換で、φ(t) = ∫(-∞, ∞) exp(itx) f(x) dx, f(x) = 1/(2π) * ∫(-∞, ∞) exp(-itx) φ(t) dt.

前者が特性関数だけど、それを期待値の式とみると、φ(t) = E[exp(itX)].

マクローリン展開をすると、 φ(t)=E[1]/0!*(it)^0 + E[X]/1!*(it)^1+E[X^2]/2!*(it)^2+...

何度も微分して t=0 を代入すると、 φ^(k)(0) = i^k E[X^k] と原点まわりのモーメントが出せる。

独立な確率変数の和については、 φX+Y(t)= E[exp(it(X+Y))] = E[exp(itX) * exp(itY)] = E[exp(itX)] * E[exp(itY)] = φX(t) * φY(t).

特性関数の対数が(第2)キュムラント母関数で、HX(t)=logφX(t) = log E[exp(itX)].

独立な確率変数の和については、 HX+Y(t)=HX(t)+HY(t) と和なのがうれしい。

マクローリン展開して、 HX(t) = 〈X〉c/1! * (it) + 〈X^2〉c / 2! * (it)^2 + 〈X^3〉c / 3! * (it)^3 + ....

〈X^k〉c のところが k 次のキュムラント。 何度も微分して 0 を代入すると出てくる。 しかも、〈X〉c=μ, 〈X^2〉c=σ^2.

確率変数に定数を足したとき、 HX+b(t) = log E[exp(it(X+b))] = HX(t) + itb
となり、〈X+b〉c=〈X〉c+b, 〈(X+b)^k〉c=〈X^k〉c (k>1).

確率変数を定数倍したとき HaX(t)=E[exp(itaX)] = HX(at) なので、 〈(aX)^k〉c=a^k〈X^k〉c.

なお、正規分布のばあい、 logφ(t)=log(exp(iμt-(σ^2)/2*t^2)=iμt-(σ^2)/2*t^2.

ようやく到着。

X1, ..., Xnが独立で同一分布、 それぞれ平均μ、分散σ^2のとき、
X=((X1+...+Xn)/n - μ)/(σ/√n)の分布が、 n→∞で、 平均 0、分散 1 に近付く。

Xを簡単にすると X = 1/(σ√n){X1+...+Xn}-μ√n/σ.

平均こと第1次のキュムラントを計算すると、 1/(σ√n) * (n * μ) - μ√n/σ = 0.

分散こと第2次のキュラントを計算すると、 1/(σ^2 n) * (n * σ^2) = 1.

3次より上のキュムラントは、1/(σ√n)^k * (n * 〈X^k〉c)→0.

ということで、キュムラント母関数が正規分布のそれに近付くので、 特性関数も正規分布のに近付き、 密度関数も近付く。

独立だったら、平均や分散が違っていても、 モーメントが高次でも定義できてれば、同じ感じ。

13581945 journal
日記

taggaの日記: 入居して10年かあ

日記 by tagga

引越して10年らしい。 管理会社は、系列にリフォーム会社があるので、 リフォームを勧めるメールが来た。

そんなに経ったのか、 元相方は元気なんだろうか。 僕はぼろぼろだけど。

13581620 journal
日記

akiraaniの日記: 御殿場高原への交通アクセスメモ

日記 by akiraani

あなげ部合宿の会場までのルートについて個人的なメモ

参考:御殿場高原ビール 交通アクセス

御殿場高原へ公共交通機関を利用していく場合、シャトルバス(PDF注意)を利用することになる。
ルートは主に3つ。

・JR三島駅からのシャトルバス
 14:05発-14:40着
 三島駅は新幹線の駅だけあって、駅前にショッピングモールがあって、飲食店がいくつかある。新幹線で来る場合、社内で駅弁食べれば良いので利用することはまれ。

・JR御殿場駅からのシャトルバス
 14:25発-14:50着
 →東から来る場合、13:56御殿場着の沼津行きからの連絡になる。御殿場駅前にはある程度食べ物屋があるので、もう1~2本早めの電車で来て、御殿場駅前で昼食をとるのもあり。ただし、さわやか御殿場IC店は激混みで2時間以上平気で待たされるので無理。
 →小田急の新松田からJR御殿場線の松田に乗り換えの場合、松田発は13:24。なんでもいいならここで昼食をとるのもあり。ただし、松田駅前で何か食べようにも選択肢はあまりない。
 →小田急ロマンスカーふじさん号を使う場合、新宿発10:40-12:24御殿場着となる。ただし、シャトルバスとの連絡が最悪で、御殿場駅前で2時間待ち、昼食をとってなお時間があまる。

・JR岩波駅からのシャトルバス(もしくは徒歩)
 14:35発-14:40着
 岩波駅から御殿場高原ビールまでは約2.5kmなので、徒歩でいけないこともない。三島駅からのシャトルバスが岩波駅にも寄ります。まあ、これ利用するくらいなら三島か御殿場からのシャトルバス使うと思いますが。

typodupeerror

UNIXはただ死んだだけでなく、本当にひどい臭いを放ち始めている -- あるソフトウェアエンジニア

読み込み中...