【数列】漸化式から求める近似式:微分近似/積分近似

更新日時:2020/10/29

数学数列数列Y番外編近似式

GOODされた数:126 回

前提など

プロローグ

どうも、安田です。
私は現在数列$Y$と名付けた難解な漸化式を持つ数列について研究をしているのですが、解法を探している過程で数列の近似方法を発見しました。
本記事ではその近似方法、微分近似と積分近似について書き留めます。

数列$Y$で微分近似と積分近似を使った記事

定義しておくこと

本記事では以下を数列$a$と定義しておきます。

\begin{equation} \begin{cases} a_{n+1} = a_n + f(a_n) \\ a_1 = x \end{cases} \end{equation}

積分近似

積分近似の概要

積分近似では、その名から察せる通り、数列を面積的な比較により近似をします。
ここで、$y = f(x)$ のグラフを考えます。$x=a_k$ から次の項に進むとき、次の項は$ x + f(x)$ で表されますので、これをグラフに表すと以下のようになります。

画像が読み込めませんでした
この時、$y=0, y=f(x), x=a_k, x=a_{k+1}$で囲まれる図形と$y=0, x=a_k, x=a_{k+1}$また$(a_k, f(a_k))$と$(a_{k+1}, f(a_{k+1}))$を結んだ線分で囲まれる図形の面積を比べる。
図から双方が近い値を取る事が分かるので、これらを比べると \begin{eqnarray*} \int_{a_k}^{a_{k+1}} f(x)dx &\simeq& \frac{1}{2} \left( a_{k+1} - a_k \right)\left( f(k+1) + f(k) \right) \\ \int_{a_k}^{a_{k+1}} f(x)dx &\simeq& \frac{1}{2} f(a_k)\left( f(k+1) + f(k) \right) \\ \int_{a_s}^{a_{t}} f(x)dx &\simeq& \frac{1}{2} \sum_{k=s}^{t-1} f(a_k)\left( f(k+1) + f(k) \right) \end{eqnarray*} である。

このように求める事が出来ました。また、積分近似では$f''(x)$の絶対値が小さいほど、傾きの変化が少なくなることから近似が一致していきます

解決済み数列での積分近似の利用

積分近似で近似をすると漸化式(の近似式)が複雑になりますので、積分近似を漸化式の近似解解法として利用するのは難しいです。
そこで、あらかじめ解決された数列において積分近似を使います。

$a_n = \sqrt{n}$ の時を考える。
\[ f(a_n) = a_{n+1} - a_n = \sqrt{n+1} - \sqrt{n} \] である。また、 \begin{eqnarray*} a_{n+1}^2 &=& a_n^2 + 1 \\ a_{n+1} &=& \sqrt{a_n^2+1} \end{eqnarray*} より \[ f(a_n) = a_{n+1} - a_n = \sqrt{a_n^2 + 1} - a_n \] なので、 \[ f(x) = \sqrt{x^2 + 1} - x \] である。この数列で積分近似を用いると \begin{eqnarray*} \int_{\sqrt{s}}^{\sqrt{t}} \left( \sqrt{x^2 + 1} - x \right)dx &\simeq& \frac{1}{2} \sum_{k=s}^{t - 1} \left(\sqrt{k+1} - \sqrt{k} \right)\left( \sqrt{k+2} - \sqrt{k} \right) \\ 2\int_{\sqrt{s}}^{\sqrt{t}} \sqrt{x^2 + 1}dx - \left[x^2\right]_{\sqrt{s}}^{\sqrt{t}} &\simeq& \sum_{k=s}^{t-1} \left(\sqrt{k+1} - \sqrt{k} \right)\left( \sqrt{k+2} - \sqrt{k} \right) \\ 2\int_{\sqrt{s}}^{\sqrt{t}} \sqrt{x^2 + 1}dx - t + s &\simeq& \sum_{k=s}^{t-1} \left(\sqrt{k+1} - \sqrt{k} \right)\left( \sqrt{k+2} - \sqrt{k} \right) \end{eqnarray*} ここで、 \begin{equation} \int \sqrt{x^2+1} dx = \frac{1}{2} \left( x\sqrt{x^2+1} + \log\left(x+\sqrt{x^2+1}\right) \right) + C \end{equation} なので、 \begin{eqnarray*} \left[ x\sqrt{x^2+1} + \log\left(x+\sqrt{x^2+1}\right) \right]_{\sqrt{s}}^{\sqrt{t}} - t + s &\simeq& \sum_{k=s}^{t-1} \left(\sqrt{k+1} - \sqrt{k} \right)\left( \sqrt{k+2} - \sqrt{k} \right) \\ \sqrt{t^2 + t} - t - \sqrt{s^2 + s} + s + \log\left(\frac{\sqrt{t}+\sqrt{t+1}}{\sqrt{s}+\sqrt{s+1}}\right)&\simeq& \sum_{k=s}^{t-1} \left(\sqrt{k+1} - \sqrt{k} \right)\left( \sqrt{k+2} - \sqrt{k} \right) \end{eqnarray*} である。

このように、不等式を算出する事ができました。左辺右辺を入れ替えると、シグマの計算の近似となっている事が分かります。

\begin{eqnarray*} &&\sum_{k=s}^{t-1} \left(\sqrt{k+1} - \sqrt{k} \right)\left( \sqrt{k+2} - \sqrt{k} \right) \\ &\simeq& \sqrt{t^2 + t} - t - \sqrt{s^2 + s} + s + \log\left(\frac{\sqrt{t}+\sqrt{t+1}}{\sqrt{s}+\sqrt{s+1}}\right) \end{eqnarray*}

この近似は$s$が大きくなるほど精度が上がりますが、$s=1, t=10000$の時ですら、左辺または右辺な値が$4.5$程度なのに対して誤差は$0.005$に満たないほど精度が高いです

また、この数列において、$n$を無限に飛ばす事を考えます。

\[ f'(x) = \frac{2x}{2\sqrt{x^2+1}} - 1 = \frac{x}{\sqrt{x^2+1}} - 1 \] より \[ \lim_{x \to \infty} f'(x) = 1 - 1 = 0 \] である。よって、$x$を無限に飛ばすと傾きが収束するので、面積の近似は完全に一致する。ゆえに \begin{equation} \lim_{s \to \infty}\lim_{t \to \infty} \left\{ \sum_{k=s}^{t-1} \left(\sqrt{k+1} - \sqrt{k} \right)\left( \sqrt{k+2} - \sqrt{k} \right) - \left( \sqrt{t^2 + t} - t - \sqrt{s^2 + s} + s + \log\left(\frac{\sqrt{t}+\sqrt{t+1}}{\sqrt{s}+\sqrt{s+1}}\right) \right) \right\} = 0 \end{equation} である。

このように積分近似では、シグマ計算の極限を求める事も可能となります。
以上が積分近似でした。

微分近似

微分近似の概要

微分近似では、積分近似のようにグラフの形を見るのではなく、式に着目します。

\[ g(xh) = a_x \]

となる関数$g$について考えます($h$は定数)。

数列$a$の漸化式から、 \begin{eqnarray*} g((x+1)h) &=& g(xh) + f(g(xh)) \\ g(xh + h) - g(xh) &=& f(g(xh)) \\ \frac{g(xh + h) - g(xh)}{h} &=& \frac{1}{h} f(g(xh)) \end{eqnarray*} ここで、$t=xh$と置くと、 \begin{eqnarray*} \frac{g(t + h) - g(t)}{h} &=& \frac{1}{h} f(g(t)) \\ \frac{d}{dt} g(t) &\simeq& \frac{1}{h} f(g(t)) \\ \frac{1}{f(g(t))} dg(t) &\simeq& \frac{dt}{h} \end{eqnarray*}

このように関数$g$が満たす微分方程式の近似式を得る事ができます。これを逆に考えると、以下の事が言えます。

\[ \frac{1}{f(g(t))} dg(t) = \frac{dt}{h} \] を満たす関数$g$を用いて \[ a_n \simeq g(nh) \] と近似出来る。

この時、$h \to 0$ で完全一致するように見えますが、実はそうはなりません

$g$が満たす微分方程式において $t = xh$ である事を考慮すると、右辺は$x$の関数となり、左辺は$g(xh)$の関数となるため、$g(xh)$は$h$に依存しない事が分かります。 すなわち、$g(nh)$は$h$に依存しないので、$h=1$の時と$h \to 0$の時の近似も変わらない事が分かります。よって一致するという事はあり得ません。

また、これを活かすと以下のように言い換えられます。

\[ \frac{1}{f(g(x))} dg(x) = dx \] を満たす関数$g$を用いて \[ a_n \simeq g(x) \] と近似できる。

また、この微分方程式を見ると分かりますが、この近似方法は$\frac{1}{f(u)}$の積分が可能であれば$g(x)$を求める事が出来ます
微分近似を用いて漸化式の解にアプローチした例ば以下の記事にあります。

【数列】数列Yとの闘い p4:近似式【シリーズ記事】

微分近似の特徴

微分近似は視覚的な近似方法ではなく式からの近似方法ですので、どういった場合に近似の精度が上がるのか分かりずらいです。
そこで、微分近似の式と漸化式を、差分記号(定義はここを参照)を用いて表して視覚的に比較できるようにします。

微分近似を言い換えると \begin{equation} \begin{cases} \Delta a_n = f(a_n) \\ \frac{d}{dx} g(x) = f(g(x)) \end{cases} \Rightarrow a_n \simeq g(n) \end{equation}

ここから、微分近似は$g'(x)$の変化が少ない方が近似の精度が向上するという事が分かります。$g'(x)$は$f(g(x))$と等しいですので、

\[ \frac{d}{dx} f(g(x)) = f'(g(x))g'(x) \]
の絶対値が小さい場合に精度が向上する事が分かります。 また、$g'(x)$は数列$a$の変化と少しは対応していますので、すなわち$|f(a)|$の大きさと対応します。
以上の事から、$|f(x)|$と$|f'(x)|$が小さい時に精度が向上するという事が分かります。
以上、微分近似でした。

微分近似や積分近似は、離散的なものを連続的なものと比べて近似するというものでした。是非皆さんもこのような近似を探してみてると面白いと思います。
それでは、お疲れさまでした。
この記事いいね!


コメント送信フォームまで飛ぶ

この記事には103件のコメントがあります。

【数列】漸化式から求める近似式:微分近似/積分近似|Mathlize
gckyissklg http://www.g1j05xd0j039mhrpj9i26y694sn281nds.org/
http://www.g1j05xd0j039mhrpj9i26y694sn281nds.org/">agckyissklg
[url=http://www.g1j05xd0j039mhrpj9i26y694sn281nds.org/]ugckyissklg[/url]
【数列】漸化式から求める近似式:微分近似/積分近似|Mathlize
http://www.g71k75266m02i7d6tpu41zw3tj7pj2tys.org/">allciqoohsk
[url=http://www.g71k75266m02i7d6tpu41zw3tj7pj2tys.org/]ullciqoohsk[/url]
llciqoohsk http://www.g71k75266m02i7d6tpu41zw3tj7pj2tys.org/
2.4 Amp Charger https://www.lashuacashotel.com/blog/best-usb-c-pd-chargers-for-2022-cnet/
ルイヴィトン財布コピー激安代引き https://travisxiax95575.educationalimpactblog.com/37027016/ルイヴィトン財布コピー激安代引き
Infared Sauna Room https://www.yanasibathroom.com/infared-sauna-room.html
ルイヴィトン財布スーパーコピー https://dominicka4ym8.blogdemls.com/17054501/ルイ-ヴィトン財布コピー-for-dummies
Fully Adjustable Silicone Vibrating Penis Ring https://www.cnt-works.com/fully-adjustable-silicone-vibrating-penis-ring.html
Pp Webbing Transportation https://www.syycupmachine.com/blog/glass-onion-a-knives-out-mystery-review-hilarious-sequel/
Ceramic Honeycomb Filter https://www.surf-duonuo.de/blog/the-pitchs-holiday-gift-guide-2022/
ルイヴィトン財布スーパーコピー https://travist3cx4.laowaiblog.com/17136343/little-known-facts-about-ルイ-ヴィトン財布コピー激安代引き
EV Charging Cable https://www.mustang-tec.com/
ルイヴィトン財布コピー激安代引き https://tysonb6oh7.blognody.com/17290627/ルイ-ヴィトンコピー財布-things-to-know-before-you-buy
50mm Glass Marbles https://www.notprops.it/blog/12-dry-shampoos-that-are-safe-to-use-amid-massive-recalls/
Automobile Clutch Release Bearing 0K20116510 https://www.xybbearings.com/automobile-clutch-release-bearing-0k20116510.html
イヴィトン財布コピー激安代引き https://jaredviry97531.dailyblogzz.com/19460567/ルイヴィトン財布コピー激安代引き
Android Code Scanner Impinj Multireader https://www.legourmetdelivery.it/blog/rfid-news-roundup/
M12 Connector https://www.hyfdwaterproofconnector.com/m12-connector
ルイヴィトン財布コピー激安代引き https://griffinjcpx86420.look4blog.com/54477221/ルイヴィトン財布コピー激安代引き
755nm 808nm 1064nm Diode Laser https://www.winhere.es/blog/MoXl-cureus-gingival-pigmentation-removal-with-a-high-power-diode-laser/
Visual Ear Cleaner Otoscope https://www.ypincheng.com/visual-ear-cleaner-otoscope
ルイヴィトン財布コピー激安代引き https://donovankvrh20753.dbblog.net/46064374/ルイヴィトン財布コピー激安代引き
Auto Spray Machine https://www.greatglasswasher.com/blog/4pxe-mercado-de-equipamentos-de-secagem-por-pulveriza-o-no/
Flexible Silicone Electronic Products Vacuum Casting https://www.shburycnc.com/flexible-silicone-electronic-products-vacuum-casting.html
ルイヴィトン財布コピー激安代引き https://dhjmjp.exblog.jp/
Adult Luxury Freestanding Bathtub https://www.yanasibathroom.com/adult-luxury-freestanding-bathtub.html
Whisky Bottle http://cmdalayout.com/blog/HYwR-irish-coffee-recipe/
ルイヴィトン財布コピー http://b17.chip.jp/yui37560
Big Spot Size Hybrid Diode Laser https://www.wdmedicshop.de/blog/rpcj-sensus-healthcare-introduces-silk™-laser-hair-removal/
617 Series C x C x FI Tee https://www.ecfitfitting.com/617-series-c-x-c-x-fi-tee.html
ルイヴィトンコピー財布 http://b17.chip.jp/yui37560
Educational Building Blocks Toys https://www.baxinpneumatic.com/blog/psfa-robot-cars-motorcycles-hoped-for-this-christmas-|-imperial-republican/
Corrugated Graphite Tapes https://www.seal-china.com/corrugated-graphite-tapes.html
30 Hp Soft Starter https://www.eoesupplyla.es/blog/n5vx-stackpath/
ルイヴィトン財布コピー激安代引き http://gifu.areablog.jp/ghthtg
Security PCB https://www.pcb-yf.com/security-pcb.html
イヴィトン財布コピー激安代引き http://htyr446.jugem.jp/
Rotavirus Adenovirus (Feces) Combo Rapid Test Kit https://www.bailimedical.com/rotavirus-adenovirus-feces-combo-rapid-test-kit.html
Aftermarket Exhaust Manifold https://www.e-gatesupply.es/blog/nmcj-the-first-cut-is-the-deepest-bimmerlife/
Floral Print Bowknot Satin Ribbon Long Satin Ribbon Headband Hair Rope https://www.nbunionsource.com/floral-print-bowknot-satin-ribbon-long-satin-ribbon-headband-hair-rope.html
ルイヴィトン財布スーパーコピー https://plaza.rakuten.co.jp/adrjkl1/
Polytetrafluoroethylene And Copper Powder Guide Ring https://www.prohuakuang.de/blog/cy9o-oil-and-gas-seal-market-to-rebound-according-to/
Metal 5 - Person Seating Group https://www.ymoutdoor.com/metal-5-person-seating-group.html
ルイヴィトン財布コピー激安代引き http://b17.chip.jp/yui37560
Pvc Foam Sheet Machine https://www.prohuakuang.de/blog/3ejz-learn-the-7-types-of-additive-manufacturing-|-techtarget/
Furniture Screw https://www.neville-taylor.com/blog/types-of-screws-[uses-facts-advantages]-engineerine/
G-Spot Vibrators https://www.zmamil.com/g-spot-vibrators
グッチ指輪コピー https://www.zhu555.jp/Gucci_Rings/

変換


※上の変換器は、TeXが正しいかどうかの確認に使ってください。
※TeXを入力する場合は、コメント本文に、$\$$ または、$\verb|\|$[, $\verb|\|$]で囲った中にTeX表示をそのまま挿入してください。
※URLは、自動的にハイパーリンクに変形されます。
※他のコメントに返信する場合は、「#コメント番号」を挿入してください。
※日本語を含まないコメントはスパムと認定されます。
※Comment including no japanese will be regarded as spam.