私が10年以上前に、仕 - 数学@ふたば保管庫

数学@ふたば保管庫 [戻る]




私が10年以上前に、仕事で取り組んだ設備性能の問題について、下記の数式を解いた事がありまして、自分なりに色々な知見が得られました。
独学なので、「そんなの学校で習うよ」というレベルだと思いますが、何か新しい事が分かるかもなので宜しければ、以下の問題についてお付き合い下さい。
-------------------------------------------------------
条件
1.ある設備の番号を、工程順にY=1、Y=2・・・とする。
2.あるパラメーターとして、W=b/(b+c)、V=c/(b+c)とする。
3.その設備の作業回数をXとする。
4.Yと比べて、Xは圧倒的に大きい(例えば500倍など)
-------------------------------------------------------
以上とした場合、下記の式を求めよ。
@Y=1の時、作業回数X回目の設備では
V(1+W+W^2+W^3+・・・+W^(X-1))の計算値を求める必要が出てきた。
 しかしXが大き過ぎるため手計算が困難である。
 因数分解等して、計算の容易な式に変換せよ。
 (ただし、近似式は不可とする)

  AY=2の時、作業回数X回目の設備では
(V^2)(1+2W+3W^2+4W^3+・・・+XW^(X-1))の計算値を求める必要が出てきた。
 以下、前述の@と同じ。

BY=3の時、作業回数X回目の設備では
(V^3)(1+3W+6W^2+10W^3+・・・+(X(X+1)/2)W^(X-1))の計算値を求める必要が出てきた。
 以下、前述の@と同じ。

CY=4の時、作業回数X回目の設備では
(V^4)(1+4W+10W^2+20W^3+・・・+(X(X+1)(X+2)/6)W^(X-1))の計算値を求める必要が出てきた。
 以下、前述の@と同じ。

  なお、今の段階では条件2.のb、cはまだ関係ありません。

  なお、今の段階では誰からもまだレスはありません。

  なんかワラタ

  @ 1 - (b / (b + c))^X
 
xが偶数ならn^x=(n^(x/2))^2
xが奇数ならn^x=n*n^(x-1)
とすればべき乗での乗算の回数は Log2N〜2*Log2N になる。
 
でどうか。

  S[Y,X]=(V^Y)Σ[k=1,X]a[Y,k]W^(k-1)
a[Y,k]=Σ[t=1,k]a[Y-1,t]
a[1,k]=1
だろうか

(1-W)S[Y,X]=(V^Y){(Σ[k=1,X]a[Y-1,k]W^(k-1))-a(Y,X)W^X}でうまくいくのかなと思ったが
{(1-W)^2}S[Y,X]を計算してみて投げた

  >>No.90121 さん
一気に計算しようとせず
@を計算してから、@をもとにAを求め、AをもとにBを求め・・・とすると求め易いです。

ちなみに@については、1-W^X となります。
>>No.90119さんの1行目と一致します。
1-W^X まで得られれば、今回のスレ的にはOKです。

つまり、
@V(1+W+W^2+W^3+・・・+W^(X-1))=1-W^X となります。

  b,cにもよるけど、Xが大きいとなると手計算は大変だと思うんだが、実際はどうやったの?

  すみません。ここで言う、手計算が大変というのは
@の場合、Xが500とすると
V(1+W+W^2+W^3+W^4+W^5+W^6・・・+W^499)
と、Wの乗数を499回足していくのが大変。という意味です。
W^5とかW^499とかは、単純に電卓で求めました。

上式は因数分解すると1-W^Xとなり、1-W^500を計算するだけで済むのでWの乗数は1回求めるだけなので、500倍ほど計算が楽になったという意味です。誤解させてすみません。

  ですので、
@V(1+W+W^2+W^3+・・・+W^(X-1))=1-W^X
A(V^2)(1+2W+3W^2+4W^3+・・・+XW^(X-1))=○○○
B(V^3)(1+3W+6W^2+10W^3+・・・+(X(X+1)/2)W^(X-1))=△△△
C(V^4)(1+4W+10W^2+20W^3+・・・+(X(X+1)(X+2)/6)W^(X-1))=□□□

以上の○○○、△△△、□□□を、簡易な数式になる様に因数分解等して求めましょう。という問題です。

  http://www.wolframalpha.com/input/?i=sum+%7Bk%28k%2B1%29%28k%2B2%29%2F6%7Dr%5Ek+from+k%3D1+to+n

http://www.wolframalpha.com/input/?i=sum+%7Bk%28k%2B1%29%28k%2B2%29%28k%2B3%29%2F24%7Dr%5Ek+from+k%3D1+to+n++

http://www.wolframalpha.com/input/?i=sum+%7Bk%28k%2B1%29%28k%2B2%29%28k%2B3%29%28k%2B4%29%28k%2B5%29%28k%2B6%29%2F5040%7Dr%5Ek+from+k%3D1+to+n++

のようにここではYに相当するものが増えるに従って
r^(n+1),r^(n+2),r^(n+3),……と
どんどんごちゃごちゃした項が増えてくようだからねえ

  そんな公式があったんですね。ありがとう勉強になります。
ただ、意外と答えはすっきりした形に収まります。

Aについては、(V^2)(1+2W+3W^2+4W^3+・・・+XW^(X-1))
=1-(W^X)(1+VX)
となります。
Wに0.8とか0.5など、1以外の数字を入れると検算し易いかと思います。

  21834 B
どこをどうやったらそんな公式になるのかと
だいたい全体にかかっていたV^2はどこいった一目おかしいとわかるが計算サイトに突っ込んどいたhttp://www.wolframalpha.com/input/?i=%28V%5E2%29%28sum+sw%5E%28s-1%29+from+s%3D1+to+X%29

  1+2W+3W^2+4W^3+・・・+XW^(X-1)が、実は
=(1-(W^X)(1+VX))/(V^2)となりまして、ということは

(V^2)(1+2W+3W^2+4W^3+・・・+XW^(X-1))
=1-(W^X)(1+VX)
となります。

  本当にすまん、V+W=1を忘れていた

  いえ、こちらこそ説明不足でした。
今後の事も考えて、条件の2番と4番に追記します。

条件
1.ある設備の番号を、工程順にY=1、Y=2・・・とする。
2.あるパラメーターとして、W=b/(b+c)、V=c/(b+c)とする。
 (W+V=1となります)
3.その設備の作業回数をXとする。
4.Yと比べて、Xは圧倒的に大きい(例えば500倍など)。
 また、bもcと比べて圧倒的に大きい。

  @をどうやって求めたか、書いてみます。

V(1+W+W^2+W^3+・・・+W^(X-1))の先頭のVは後で考えるとして、まずは
1+W+W^2+W^3+・・・+W^(X-1)から先に考えます。
---------------------------------------------------------
(1+W)(1-W)=1-W^2 であるから、(1+W)=(1-W^2)/(1-W)となる。

次に、1+W+W^2を計算する場合、(1+W)+W^2=(1-W^2)/(1-W)+W^2
を計算すれば良く、これを展開して計算すると、(1-W^3)/(1-W)となる。
従って、(1+W+W^2)=(1-W^3)/(1-W)となる。

次に、1+W+W^2+W^3を計算する場合、(1+W+W^2)+W^3=(1-W^3)/(1-W)+W^3
を計算すれば良く、これを展開して計算すると、(1-W^4)/(1-W)となる。
従って、(1+W+W^2+W^3)=(1-W^4)/(1-W)となる。
---------------------------------------------------------

  以上を繰り返すことによって、最終的に
1+W+W^2+W^3+・・・+W^(X-1)=(1-W^X)/(1-W)となる。
ここで、条件2より、1-W=Vであることから
1+W+W^2+W^3+・・・+W^(X-1)=(1-W^X)/Vとなる。

ここで、元の@の式では先頭にVが付いていたので計算
すると、
V(1+W+W^2+W^3+・・・+W^(X-1))=V(1-W^X)/V となり、
V(1+W+W^2+W^3+・・・+W^(X-1))=1-W^X となる。

以上です。

  (1+W)(1-W)=1-W^2 になる事を利用すれば良いと気付いた時は、本当に嬉しかった思い出があります。
また、Aは@を応用して求めたのですが、このAをどうやって求めるか気付いた時も嬉しく、その日は眠れませんでした。
ほんと数学が好きで良かった。と思った、過去の良い思い出です(ただの因数分解ですがw)

さて、Aの解法に行く前に少し話題を変えて、Wの左側に付いている係数を書いてみたいと思います。

  Wの左隣の系数について
@Y=1の時:1, 1, 1, 1, 1, 1・・・
AY=2の時:1, 2, 3, 4, 5, 6・・・
BY=3の時:1, 3, 6, 10, 15, 21・・・
CY=4の時:1, 4, 10, 20, 35, 56・・・
DY=5の時:1, 5, 15, 35, 70,126・・・
EY=6の時:1, 6, 21, 56,125,252・・・



となります。
1.係数を縦方向に見ていくと、横方向も同じ並びがある。
2.『係数』と『右上の係数』を足すと『右隣の係数』になる。
3.係数を右下方向に斜めに見ていくと、『パスカルの三角形』の数列になる。
などの事が分かります。面白いですよね。

  DEについては勢い余って、書いてしまいました。すみません。

  3.について、訂正します><すみません。

3.係数を『右上方向』に斜めに見ていくと、『パスカルの三角形』の数列になる。
です。

  ただの高二でやる級数じゃん

  母関数

  やはり高校レベルの内容でしたか。すみませんでした。

では、途中経過をはしょって結果だけ下記にまとめます。
A(V^2)(1+2W+3W^2+4W^3+・・・+XW^(X-1))=1-(W^X)(1+VX)
B(V^3)(1+3W+6W^2+10W^3+・・・+(X(X+1)/2)W^(X-1))=1-(W^X)(1+VX+(V^2)X(X+1)/2)
C(V^4)(1+4W+10W^2+20W^3+・・・+(X(X+1)(X+2)/6)W^(X-1))
=1-(W^X)(1+VX+(V^2)X(X+1)/2+(V^3)X(X+1)(X+2)/6)

となります。引き続き考えると、
D(V^5)(1+5W+15W^2+35W^3+・・・+(X(X+1)(X+2)(X+3)/4!)W^(X-1))
=1-(W^X)(1+VX+(V^2)X(X+1)/2+(V^3)X(X+1)(X+2)/3!+(V^4)X(X+1)(X+2)(X+3)/4!)



などとなりますね。

  E以降については長くなり過ぎるので、階乗を使って記載します。
例えば、X(X+1)(X+2)については階乗を使って記載すると、(X+2)!/(X-1)!となります。

E(V^6)(1+6W+21W^2+56W^3+・・・+(X+4)!/(5!(X-1)!))W^(X-1))
=1-(W^X)(1+VX+(V^2)(X+1)!/(2!(X-1)!)+(V^3)(X+2)!/(3!(X-1)!)+(V^4)(X+3)!/(4!(X-1)!)+(V^5)(X+4)!/(5!(X-1)!))
となります。

以上から、F,G,H(Y=7,8,9・・・)もどうなるか推測できますね。

  以上を元に、Yを使って記載すると


(V^Y)(1+WY+(W^2)(Y+1)!/(2!(Y-1)!)+(W^3)(Y+2)!/(3!(Y-1)!)+・・・+(W^(X-1))(Y+X-2)!/((X-1)!(Y-1)!))

=1-(W^X)(1+VX+(V^2)(X+1)!/(2!(X-1)!)+(V^3)(X+2)!/(3!(X-1)!)+・・・+(V^(Y-1))(X+Y-2)!/((Y-1)!(X-1)!))
となります。

  2534 B
Σを使ってまとめると、添付画像の様になります。
http://www.wolframalpha.com/input/?i=%28%28V%5EY%29+sum+%28%28W%5En%29%28Y-1%2Bn%29%21%2F%28%28Y-1%29%21n%21%29%29+from+n%3D0+to+%28X-1%29%29%3D1-%28%28W%5EX%29+sum+%28%28V%5Em%29%28X-1%2Bm%29%21%2F%28%28X-1%29%21m%21%29%29+from+m%3D0+to+%28Y-1%29%29

  2163 B
組み合わせの、『C』を使うと添付画像の様にすっきりまとまります。
『C』は『重複を持たない組合せ』なのですが、『重複組合せ』で、『H』というのがあります。リンク先のサイトは『H』に対応していない様なので、添付画像ではカッコの中が(Y-1+n,n)となっていますが、『H』が使える場合は(Y,n)となり、もっとすっきりまとまります。http://www.wolframalpha.com/input/?i=%28%28V%5EY%29+sum+%28%28W%5En%29c%28Y-1%2Bn%2Cn%29%29+from+n%3D0+to+%28X-1%29%29%3D1-%28%28W%5EX%29+sum+%28%28V%5Em%29c%28X-1%2Bm%2Cm%29%29+from+m%3D0+to+%28Y-1%29%29http://ja.wikipedia.org/wiki/%E7%B5%84%E5%90%88%E3%81%9B_(%E6%95%B0%E5%AD%A6)

  以降についてですが、

A.この式の近似式を求めてみる。
B.高校で習うレベルとの事で釈迦に説法だと思い、>>No.90151 は途中経過をはしょったのですが、解き方を書いてみる。
C.そもそも、どこから(V^3)(1+3W+6W^2+10W^3+・・・+(X(X+1)/2)W^(X-1))等が出て来たのか、解説する。

以上の順でやってみようと思いますが、先に書いて欲しい順などあれば教えて下さい。