「微分」 現場数学(その1)

更新日

投稿日

♦ 数学の教科書に書いてある通りにプログラムしてもだめ!では計算できない!

 サイン・コサイン何になる…とか嘯(うそぶ)いていた学生時代。しかし「現場に出て数学の重要性をひしひしと感じているのに、確実な理解が出来ていない…」と不安を抱いている技術者の皆さん。そこで、今回から18回にわたり「本当に意味のある“現場数学”とは何か」について私ともう一度考えてみませんか?

 

1.“現場数学”では、hは有限でなければならない

 高校で微分とか積分を習った時、先生が今まで見たことのない如何(いか)にも新しい式であると黒板に書く様に感激した経験を思い出します。しかし、これでは本当は計算できないのです。「現場数学」では、hは有限でなければなりません。となると…小さい、と思われる数、例えば0.00001とかにすれば良いということでしょうか?いえいえ、全くそうではないのです。この辺りから疑問を感じていた人は多いのです。でも、高校の先生に聞くと、ムニャムニャとした返事。それで数学が嫌いになった人がいるほどです。どうなっているのでしょうか?

 実はhは小さいどころか無限大でも良いのです!これは極端な話です。でも、高校数学で重要なことを教えていないことは確かなのです。また、大学に入ると微積分は途端にε-δの話になり、余計に誤魔化されてしまいます。実は、最初に考えなければならないのは対称性の問題なのです。通常の微分係数の記述は、

    

 です。しかし、これでは、xの位置から正方向のみの変化を見ることになり、その位置で増加している関数なら当然、x点での値より大きくなります。それにも関らず なのだから、正しい値に近づくはずと習うのです。これは一種の詭弁(きべん)です。二次関数y=x2にこの微分係数の式を当てはめてみましょう。

    

 となります。hだけの誤差があり、 で正しい2xになるという論法です。実は、この公式は前進公式と呼ばれます。xから正方向だけをみているからです。x点での微分係数を求めたいのですから、-h方向もみないといけません。それは後退公式と呼ばれる次式です。

    

 この式の誤差は-hです。ということは(2)式と(3)式の平均をとれば、hの大きさに依(よ)らず微分係数は2xと厳密に求まります。

    

 つまり微分係数はhに依らず、極端にいえばhは無限大でも良い式が作れたのです。この例では、2次関数の微分は一次関数ですから、hに依らず誤差が消えてしまいました。3次関数でも、何次関数でも同様に誤差を消してしまう公式が作れます。ぜひ、自分の手を動かして試して下さい。

2. より本質理解し、数式を上手に扱う

 この例題は、教科書に書いてある通りの式をそのままプログラムして、数値計算してはいけないことを如実に示しています。「現場数学」では、数値をなるべく正確にしかも効率良く計算しなければなりません。それを身につけるためには、このような高校の教科書に載っていなかった方策を勉強する必要があります。これは、頭の良...

♦ 数学の教科書に書いてある通りにプログラムしてもだめ!では計算できない!

 サイン・コサイン何になる…とか嘯(うそぶ)いていた学生時代。しかし「現場に出て数学の重要性をひしひしと感じているのに、確実な理解が出来ていない…」と不安を抱いている技術者の皆さん。そこで、今回から18回にわたり「本当に意味のある“現場数学”とは何か」について私ともう一度考えてみませんか?

 

1.“現場数学”では、hは有限でなければならない

 高校で微分とか積分を習った時、先生が今まで見たことのない如何(いか)にも新しい式であると黒板に書く様に感激した経験を思い出します。しかし、これでは本当は計算できないのです。「現場数学」では、hは有限でなければなりません。となると…小さい、と思われる数、例えば0.00001とかにすれば良いということでしょうか?いえいえ、全くそうではないのです。この辺りから疑問を感じていた人は多いのです。でも、高校の先生に聞くと、ムニャムニャとした返事。それで数学が嫌いになった人がいるほどです。どうなっているのでしょうか?

 実はhは小さいどころか無限大でも良いのです!これは極端な話です。でも、高校数学で重要なことを教えていないことは確かなのです。また、大学に入ると微積分は途端にε-δの話になり、余計に誤魔化されてしまいます。実は、最初に考えなければならないのは対称性の問題なのです。通常の微分係数の記述は、

    

 です。しかし、これでは、xの位置から正方向のみの変化を見ることになり、その位置で増加している関数なら当然、x点での値より大きくなります。それにも関らず なのだから、正しい値に近づくはずと習うのです。これは一種の詭弁(きべん)です。二次関数y=x2にこの微分係数の式を当てはめてみましょう。

    

 となります。hだけの誤差があり、 で正しい2xになるという論法です。実は、この公式は前進公式と呼ばれます。xから正方向だけをみているからです。x点での微分係数を求めたいのですから、-h方向もみないといけません。それは後退公式と呼ばれる次式です。

    

 この式の誤差は-hです。ということは(2)式と(3)式の平均をとれば、hの大きさに依(よ)らず微分係数は2xと厳密に求まります。

    

 つまり微分係数はhに依らず、極端にいえばhは無限大でも良い式が作れたのです。この例では、2次関数の微分は一次関数ですから、hに依らず誤差が消えてしまいました。3次関数でも、何次関数でも同様に誤差を消してしまう公式が作れます。ぜひ、自分の手を動かして試して下さい。

2. より本質理解し、数式を上手に扱う

 この例題は、教科書に書いてある通りの式をそのままプログラムして、数値計算してはいけないことを如実に示しています。「現場数学」では、数値をなるべく正確にしかも効率良く計算しなければなりません。それを身につけるためには、このような高校の教科書に載っていなかった方策を勉強する必要があります。これは、頭の良くなる本などで、特別な例を覚えておいて計算を早くする、という類のものではありません。より本質を理解して、数式を上手に扱う技法なのです。計算をするのに、そろばんを使おうとスーパーコンピューターを使おうと、どれも桁数は有限です。 は何倍精度を持ってしても厳密には成り立たないのです。そうではなく「現場数学」では、そもそも誤差を発生させなくする定式化、なるべく小さくする定式化をすることが肝要なのです。

 

 次回に続きます。

   続きを読むには・・・


この記事の著者

川添 良幸

市販の材料設計シミュレーションプログラムでは満足できない御社技術者に、本当に意味のある物理の基本に基づいた設計法を伝授します。我々の計算と実験結果が合わない場合は実験の方に問題があると言えるレベルを達成しています。

市販の材料設計シミュレーションプログラムでは満足できない御社技術者に、本当に意味のある物理の基本に基づいた設計法を伝授します。我々の計算と実験結果が合わな...


「SQC一般」の他のキーワード解説記事

もっと見る
~ 無次元化 現場数学(その14)

  ♦ 限りある時間で必要な桁数の解を求める策を採用 1.単位系と自然  力学で使われる量は質量、長さと時間が基本で、1k...

  ♦ 限りある時間で必要な桁数の解を求める策を採用 1.単位系と自然  力学で使われる量は質量、長さと時間が基本で、1k...


~公式とは? 現場数学(その4)

   力学、流体力学、電磁気学、光学、熱力学、統計力学、量子力学…、さらには遺伝子工学まで、多くの数学公式が現れ、その数値解...

   力学、流体力学、電磁気学、光学、熱力学、統計力学、量子力学…、さらには遺伝子工学まで、多くの数学公式が現れ、その数値解...


~マルチスケール模型での計算 現場数学(その6)

  1.  計算の限界  如何(いか)にスーパーコンピューターを活用しても、計算出来る範囲はやはり限られています。例えば、タン...

  1.  計算の限界  如何(いか)にスーパーコンピューターを活用しても、計算出来る範囲はやはり限られています。例えば、タン...


「SQC一般」の活用事例

もっと見る
第1種の誤りと第2種の誤り

 「あわて者の誤り」と「ぼんやり者の誤り」をご存知ですか、あわて者の誤りは正式には『第一種の誤り』(Type Ⅰ error)と呼称し、ぼんやり者の誤りは...

 「あわて者の誤り」と「ぼんやり者の誤り」をご存知ですか、あわて者の誤りは正式には『第一種の誤り』(Type Ⅰ error)と呼称し、ぼんやり者の誤りは...


飛行機事故とセンサー・フィードバックの分布

   ボーイング 737 Max 8機が短い期間に2回の墜落死亡事故を起こし、多くの航空会社が同機種の運航を停止する事態となっています。多くの...

   ボーイング 737 Max 8機が短い期間に2回の墜落死亡事故を起こし、多くの航空会社が同機種の運航を停止する事態となっています。多くの...


統計手法による 2018年LPGA公式記録データの解析事例(その1)

【統計手法による 2018年LPGA公式記録データの解析 連載目次】 統計手法による 2018年LPGA公式記録データの解析事例(その1) 統...

【統計手法による 2018年LPGA公式記録データの解析 連載目次】 統計手法による 2018年LPGA公式記録データの解析事例(その1) 統...